Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009974 3,929,980 (G)7→8 intergenic (‑84/‑145) RPPX_17505 ← / → RPPX_17510 ATPase AAA/type IV secretion protein Rhs

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0099743,929,9731.G88.5% 494.6 / 65.7 156intergenic (‑77/‑152)RPPX_17505/RPPX_17510ATPase AAA/type IV secretion protein Rhs
Reads supporting (aligned to +/- strand):  ref base . (6/12);  new base G (103/35);  total (109/47)
Fisher's exact test for biased strand distribution p-value = 7.37e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 174 total aligned reads displayed.

ACCAGTTTCTCCTGGCCTAGGATCTGGGTAGAAAGAAAAGTGCGCAGCGCGATCAGCGCTTCACGGTGTTCCATCGGCGACGGTTCCTGGTAGGGGGCCTGAGGCGGCGTGAGGCAGCAGCCTTGGCAAGGGGGGATACTTTAATCTATC‑C‑GGGGGGGCGACGACCAATGGCGGCGGGTAAATTTATCGTTGCAGACGGGCCTGCGGGCAAAGAAATGTCCCACGCATGGCTGGCCTTTGCCTGATGGTCCGTGACGGCGTGCTGGCGCCTGCTTGTGCTTTTCGGACAAGCAGGA  >  CP009974/3929823‑3930118
                                                                                                                                                        |                                                                                                                                                 
aCCAGTTTCTCCTGGCCTAGGATCTGGGTAGAAAGAAAAGTGCGCAGCGCGATCAGCGCTTCACGGTGTTCCATCGGCGACGGTTCCTGGTAGGGGGCCTGAGGCGGCGTGAGGCAGCAGCCTTGGCAAGGGGGGATACTTTAATCTATc‑c                                                                                                                                                    <  1:1273115/151‑1 (MQ=255)
  cAGTTTCTCCTGGCCTAGGATCTGGGTAGAAAGAAAAGTGCGCAGCGCGATCAGCGCTTCACGGTGTTCCATCGGCGACGGTTCCTGGTAGGGGGCCTGAGGCGGCGTGAGGCAGCAGCCTTGGCAAGGGGGGATACTTTAATCTATC‑C‑gg                                                                                                                                                 <  2:1631442/151‑1 (MQ=255)
  cAGTTTCTCCTGGCCTAGGATCTGGGTAGAAAGAAAAGTGCGCAGCGCGATCAGCGCTTCACGGTGTTCCATCGGCGACGGTTCCTGGTAGGGGGCCTGAGGCGGCGTGAGGCAGCAGCCTTGGCAAGGGGGGATACTTTAATCTATC‑C‑gg                                                                                                                                                 >  2:591223/1‑151 (MQ=255)
   ggTTTCTCCTGGCCTAGGATCTGGGTAGAAAGAAACGTGCGCAGCGCGATCAGCGCTTCACGGTGTTCCATCGGCGACGGTTCCTGGTAGGGGGCCTGAGGCGGCGTGAGGCAGCAGCCTTGGCAAGGGGGGATACTTTAATCTATC‑C‑ggg                                                                                                                                                <  2:2002501/150‑1 (MQ=255)
     tttCTCCTGGCCTAGGATCTGGGTAGAAAGAAAAGTGCGCAGCGCGATCAGCGCTTCACGGTGTTCCATCGGCGACGGTTCCTGGTAGGGGGCCTGAGGCGGCGTGAGGCAGCAGCCTTGGCAAGGGGGGATACTTTAATCTATC‑C‑ggggg                                                                                                                                              <  1:1246308/151‑1 (MQ=255)
       tctcCTGGCCTAGGATCTGGGTAGAAAGAAAAGTGCGCAGCGCGATCAGCGCTTCACGGTGTTCCATCGGCGACGGTTCCTGGTAGGGGGCCTGAGGCGGCGTGAGGCAGCAGCCTTGGCAAGGGGGGATACTTTAATCTATC‑C‑ggggggg                                                                                                                                            <  1:1847699/151‑1 (MQ=255)
       tccccTGGCCTAGGATCGGGGTAGAAAGAAAAGTGCGCAGCGCGATCAGCGCTTCACGGTGTTCCATCGGCGACGGTTCCTGGTAGGGGGCCTGAGGCGGCGTGAGGCAGCAGCCTTGGCAAGGGGGGATACTTTAATCTATC‑C‑ggggggg                                                                                                                                            <  1:2293705/148‑1 (MQ=255)
        ctcCTGGCCTAGGATCTGGGTAGAAAGAAAAGTGCGCAGCGCGATCAGCGCTTCACGGTGTTCCATCGGCGACGGTTCCTGGTAGGGGGCCTGAGGCGGCGTGAGGCAGCAGCCTTGGCAAGGGGGGATACTTTAATCTATC‑C‑gggggggg                                                                                                                                           >  1:2032116/1‑150 (MQ=255)
           cTGGCCTAGGATCTGGGTAGAAAGAAAAGTGCGCAGCGCGATCAGCGCTTCACGGTGTTCCATCGGCGACGGTTCCTGGTAGGGGGCCTGAGGCGGCGTGAGGCAGCAGCCTTGGCAAGGGGGGATACTTTAATCTATC‑C‑ggggggggcga                                                                                                                                        >  1:419301/1‑147 (MQ=255)
            gggCCGAGGAGCTGGGGAGAAAGGAAAGGGCGCGGCGCGGTCAGCGCTGCACGGTGTTCCATGGGCGACGGTTCCTGGGAGGGGGCCGGAGGCGGCGTGGGGCAGCGGCCTTGGCAAGGGGGGATACTTTAATCTATC‑C‑ggggggggcgac                                                                                                                                       <  2:3592403/150‑6 (MQ=255)
             ggCCTAGGATCTGGGTAGAAAGAAAAGTGCGCAGCGCGATCAGCGCTTCACGGTGTTCCATCGGCGACGGTTCCTGGTAGGGGGCCTGAGGCGGCGTGAGGCAGCAGCCTTGGCAAGGGGGGATACTTTAATCTATC‑CGGGGGGGGcgacg                                                                                                                                       >  2:2101024/1‑151 (MQ=255)
                  aGGATCTGGGTAGAAAGAAAAGTGCGCAGCGCGATCAGCGCTTCACGGTGTTCCATCGGCGACGGTTCCTGGTAGGGGGCCTGAGGCGGCGTGAGGCAGCAGCCTTGGCAAGGGGGGATACTTTAATCTATC‑CGGGGGGGGCGACGACCaa                                                                                                                                  >  2:270420/1‑151 (MQ=255)
                   ggATCTGGGTAGAAAGAAAAGTGCGCAGCGCGATCAGCGCTTCACGGTGTTCCATCGGCGACGGTTCCTGGTAGGGGGCCTGAGGCGGCGTGAGGCAGCAGCCTTGGCAAGGGGGGATACTTTAATCTATC‑CGGGGGGGGCGACGACCAAt                                                                                                                                 >  2:2948733/1‑151 (MQ=255)
                   ggATCTGGGTAGAAAGAAAAGTGCGCAGCGCGATCAGCGCTTCACGGTGTTCCATCGGCGACGGTTCCTGGTAGGGGGCCTGAGGCGGCGTGAGGCAGCAGCCTTGGCAAGGGGGGATACTTTAATCTATC‑CGGGGGGGGCGACGACCAAt                                                                                                                                 >  1:1406856/1‑151 (MQ=255)
                    gATCTGGGTAGAATGTTTAGTGCGCAGCGCGATCAGCGCTTCACGGTGTTCCATCGGCGACGGTTCCTGGTAGGGGGCCTGAGGCGGCGTGAGGCAGCAGCCTTGGCAAGGGGGGATACTTTAATCTATC‑CGGGGGGGGCGACGACCAATg                                                                                                                                >  2:1237383/1‑151 (MQ=255)
                      tCTGGGTAGAAAGAAAAGTGCGCAGCGCGATCAGCGCTTCACGGTGTTCCATCGGCGACGGTTCCTGGTAGGGGGCCTGAGGCGGCGTGAGGCAGCAGCCTTGGCAAGGGGGGATACTTTAATCTATC‑CGGGGGGGGCGACGACCAATggg                                                                                                                              >  1:3444233/1‑150 (MQ=255)
                         ggggAGCAAGAAGTGGGCGGAGCGCGATGGGGGCGTCACGGTGTTCCAGCGGCGACGGGTCGTGGTAGGGGGCCTGAGGCGGCGGGAGGCAGCAGCCTTGGGAAGGGGGGATACTTTAATCTATC‑CGGGGGGGGCGACGACCAATggcggc                                                                                                                           <  2:3064987/151‑1 (MQ=255)
                           ggggaaagaagAGTGCGCAGGGGGATCGGGGCTTCAGGGTGTTCCATCGGGGGGGGTTCGTGGTAGGGGGCCGGGGGCGGCGTGAGGCGGCAGCCTTGGCAGGGGGGGATACTTTAATCTATC‑CGGGGGGGGCGACGACCAATGGcggcgg                                                                                                                         <  2:2217684/148‑1 (MQ=255)
                            tAGAAAGAAAAGTGCGCAGCGCGATCAGCGCTTCACGGTGTTCCATCGGCGACGGTTCCTGGTAGGGGGCCTGAGGCGGCGTGAGGCAGCAGCCTTGGCAAGGGGGGATACTTTATTCTATC‑CGGGGGGGGCGACGACCAATGGCGGCggg                                                                                                                        >  2:722410/1‑151 (MQ=255)
                            gagaaagagAAGGGGGCGGCGCGAGCAGCGCTTCACGGGGTGCCAGCGGGGAGGGTGCCGGGTAGGGGGCCGGAGGCGGCGTGAGGCAGAGGGCTTGGAAAGGGGGGATACTTTAATCTATC‑CGGGGGGGGCGACGACCAATGGCGGCggg                                                                                                                        <  2:4341575/150‑1 (MQ=255)
                             agaaagaaAAGTGCGCAGCGCGATCAGCGCTTCACGGTGTTCCATCGGCGACGGTTCCTGGTAGGGGGCCTGAGGCGGCGTGAGGCAGCAGCCTTGGCAAGGGGGGATACTTTAATCTATC‑CGGGGGGGGCGACGACCAATGGCGGCGGGt                                                                                                                       >  2:4728061/1‑151 (MQ=255)
                                aagaaaAGTGCGCAGCGCGATCAGCGCTTCACGGTGTTCCATCGGCGACGGTTCCTGGTAGGGGGCCTGAGGCGGCGTGAGGCAGCAGCCTTGGCAAGGGGGGATACTTTAATCTATC‑CGGGGGGGGCGACGACCAATGGCGGCGGGTaaa                                                                                                                    >  2:2339491/1‑151 (MQ=255)
                                 agaaaAGTGCGCAGCGCGATCAGCGCTTCACGGTGTTCCATCGGCGACGGTTCCTGGTAGGGGGCCTGAGGCGGCGTGAGGCAGCAGCCTTGGCAAGGGGGGATACTTTAATCTATC‑CGGGGGGGGCGACGACCAATGGCGGCGGGTAAAt                                                                                                                   >  1:1156605/1‑151 (MQ=255)
                                 agaaaAGTGCGCAGCGCGATCAGCGCTTCACGGTGTTCCATCGGCGACGGTTCCTGGTAGGGGGCCTGAGGCGGCGTGAGGCAGCAGCCTCGGCAAGGGGGGATACTTTAATCTATC‑CGGGGGGGGCGACGACCAATGGCGGCGGGTagcc                                                                                                                   >  2:4076673/1‑148 (MQ=255)
                                      ggTGCGCGGGGCGATCAGCGCTTCACGGTGTCCCTCGGGGGCCGGTTCCTGGGAGGGGGCCTGAGGCGGCGTGAGGCAGCAGCCTTGGCAAGGGGGGATACTTTAATCTATC‑CGGGGGGGGCGACGACCAATGGCGGCGGGTAAATTTATc                                                                                                              <  2:4705467/150‑1 (MQ=255)
                                      aGTGCGCAGCGCGATCAGCGCTTCACGGTGTTCCATCGGCGACGGTTCCTGGTAGGGGGCCTGAGGCGGCGTGAGGCAGCAGCCTTGGCAAGGGGGGATACTTTAATCTATC‑CGGGGGGGGCGACGACCAATGGCGGCGGGTAAATTTATc                                                                                                              >  2:1770443/1‑151 (MQ=255)
                                      aGTGCGCAGCGCGATCAGCGCTTCACGGTGTTCCATCGGCGACGGTTCCTGGTAGGGGGCCTGAGGCGGCGTGAGGCAGCAGCCTTGGCAAGGGGGGATACTTTAATCTATC‑CGGGGGGGGCGACGACCAATGGCGGCGGGTAAATTTATc                                                                                                              >  1:4991683/1‑151 (MQ=255)
                                      aGTGCGCAGCGCGATCAGCGCTTCACGGTGTTCCATCGGCGACGGTTCCTGGTAGGGGGCCTGAGGCGGCGTGAGGCAGCAGCCTTGGCAAGGGGGGATACTTTAATCTATC‑CGGGGGGGGCGACGACCAATGGCGGCGGGTAAATTTATc                                                                                                              >  1:151801/1‑151 (MQ=255)
                                       gTGCGCAGCGCGATCAGCGCTTCACGGTGTTCCATCGGCGACGGTTCCTGGTAGGGGGCCTGAGGCGGCGTGAGGCAGCAGCCTTGGCAAGGGGGGATACTTTAATCTATC‑CGGGGGGGGCGACGACCAATGGCGGCGGGTAAATTTATCg                                                                                                             >  1:476597/1‑151 (MQ=255)
                                         gcgcAGCGCGAGCAGCGCGTCAGGGTGTTCCAGCGGGGGGGGTTCCGGGTAGGGGGCGTGAGGCGGCGGGAGGGAGCAGCCTTGGCAAGGGGGGATACTTTAATCTATC‑CGGGGGGGGCGACGACCAATGGCGGCGGGTAAATTTATCGtt                                                                                                           <  1:2971398/151‑1 (MQ=255)
                                          cgcAGCGCGATCAGCGCTTCACGGTGTTCCATCGGCGACGGTTCCTGGTAGGGGGCCTGAGGCGGCGTGAGGCAGCAGCCTTGGCAAGGGGGGATACTTTAATCTATC‑CGGGGGGGGCGACGACCAATGGCGGCGGGTAAATTTATCGTTg                                                                                                          >  1:432207/1‑151 (MQ=255)
                                             aGCGCGATCAGCGCTTCACGGTGTTCCATCGGCGACGGTTCCTGGTAGGGGGCCTGAGGCGGCGTGAGGCAGCAGCCTTGGCAAGGGGGGATACTTTAATCTATC‑CGGGGGGGGCGACGACCAATGGCGGCGGGTAAATTTATCGTTGCAg                                                                                                       >  2:1416790/1‑151 (MQ=255)
                                              gggggATCAGCGCTTCACGGTGGTCCATCGGCGGCGGTTCCGGGTAGGGGGCTGGAGGCGGCGTGAGGCAGCAGCCTTGGCAAGGGGGGATACTTTAATCTATC‑CGGGGGGGGCGACGACCAATGGCGGCGGGTAAATTTATCGTTGCAGa                                                                                                      <  2:5061531/147‑1 (MQ=255)
                                               ggcgATCGGGGCTGCACGGTGTTCCGGCGGCGACGGTTCCTGGGAGGGGGCCTGGGGCGGGGTGGGGCAGCAGCCTTGGCAGGGGGGGATACTTTAATCTATC‑GGGGGGGGGCGACGACCAATGGCGGCGGGTAAATTTATCGTTGCAGAc                                                                                                     <  1:661617/150‑1 (MQ=255)
                                               cgcgATCAGCGCTTCACGGTGTTCCATCGGCGACGGTTCCTGGTAGGGGGCCTGAGGCGGCGTGAGGCAGCAGCCTTGGCAAGGGGGGATACTTTAATCTATC‑CGGGGGGGGCGACGACCAATGGCGGCGGGTAAATTTATGGTTGCAGAc                                                                                                     >  2:4110027/1‑151 (MQ=255)
                                               cccgctCACCGCTCCGCGGTGTTCCATGGGCGACGGTTCCTGGTTGGGGGCCGGTGGCGGCGTGGGGCGGCAGCCTTGGCAAGGGGGGATCCTTTAATCAATC‑CGGGGGGGGCGACGACCAATGGCGGCGGGTAAATTTATCGTTGCAGAc                                                                                                     <  1:1313671/146‑1 (MQ=255)
                                                gcgATCAGCGCTTCTCGGTGTTCCATCGGCGACGGTTCCTGGTAGGGGGCCTGAGGCGGCGTGAGGCAGCAGCCTTGGCAAGGGGGGATACTTTAATCTCTC‑C‑GGGTGGGCGACGACCAATGGCGGCGGGTAAATTTATCGTTGCAGACgg                                                                                                   >  2:3623115/1‑151 (MQ=255)
                                                gcgATCAGCGCTTCACGGTGTTCCATCGGCGACGGTTCCTGGTAGGGGGCCTGAGGCGGCGTGAGGCAGCAGCCTTGGCAAGGGGGGATACTTTAATCTATC‑CGGGGGGGGCGACGACCAATGGCGGCGGGTAAATTTATCGTTGCAGACg                                                                                                    >  2:242387/1‑151 (MQ=255)
                                                   aTCAGCGCTTCACGGTGTTCCATCGGCGACGGTTCCTGGTAGGGGGCCTGAGGCGGCGTGAGGCAGCAGCCTTGGCAAGGGGGGATACTTTAATCTATC‑CGGGGGGGGCGACGACCAATGGCGGCGGGTAAATTTATCGTTGCAGACGGGc                                                                                                 >  1:608749/1‑151 (MQ=255)
                                                     cAGCGCTTCACGGTGTTCCATCGGCGACGGTTCCTGGTAGGGGGCCTGAGGCGGCGTGAGGCAGCAGCCTTGGCAAGGGGGGATACTTTAATCTATC‑CGGGGGGGGCGACGACCAATGGCGGCGGGTAAATTTATCGTTGCAGACGGGCCt                                                                                               >  1:4389382/1‑151 (MQ=255)
                                                      aGCGCTTCACGGTGTTCCATCGGCGACGGTTCCTGGTAGGGGGCCTGAGGCGGCGTGAGGCAGCAGCCTTGGCAAGGGGGGATACTTTAATCTATC‑CGGGGGGGGCGACGACCAATGGCGGCGGGTAAATTTATCGTTGCAGACGGGCCTg                                                                                              >  1:4663454/1‑151 (MQ=255)
                                                      aGCGCTTCACGGTGTTCCATCGGCGACGGTTCCTGGTAGGGGGCCTGAGGCGGCGTGAGGCAGCAGCCTTGGCAAGGGGGGATACTTTAATCTATC‑CGGGGGGGGCGACGACCAATGGCGGCGGGTAAATTTATCGTTGCAGACGGGCCTg                                                                                              >  1:4663305/1‑151 (MQ=255)
                                                      aGCGCTTCACGGTGTTCCATCGGCGACGGTTCCTGGTAGGGGACCTGAGGCGGCGTGAGGCAGCAGCCTTGGCAAGGGGGGATACTTTAATCTATC‑CGGGGGGGGCGACGCCCAATGGCGGCGCGCAAATTTATCGTTGCAGACGGGCCTg                                                                                              >  2:5222815/1‑151 (MQ=255)
                                                         gcTTCACGGTGTTCCATCGGCGACGGTTCCCGGTAGGGGGCCTGAGGCGGCGTGAGGCAGCAGCCTTGGCAAGGGGGGATAGTTTAATCTATC‑CGGGGGGGGCGACGACCAATGGCGGCGGGTACATTTATCGTTGCAGACGGGCCTGCgg                                                                                           >  2:3122228/1‑151 (MQ=255)
                                                          gggcACGGTGTTCCATGGGCGACGGTTCTTGGTAGGGGGCCTGAGGCGGCGGGAGGCAGCAGCCTTGGCAAGGGGGGATACTTTAATCTATC‑CGGGGGGGGCGACGACCAATGGCGGCGGGTAAATTTATCGTTGCAGACGGGCCTGCggg                                                                                          <  2:1840792/148‑1 (MQ=255)
                                                          cTTCACGGTGTTCCATCGGCGACGGTTCCTGGTAGGGGGCCTGAGGCGGCGTGAGGCAGCAGCCTTGGCAAGGGGGGATACTTTAATCTATC‑CGGGGGGGGCGACGACCAATGGCGGCGGGTAAATTTATCGTTGCAGACGGGCCTGCggg                                                                                          >  2:1413059/1‑151 (MQ=255)
                                                               gggTGTTCCAGCGGCGGCGGTTCCTGGGAGGGGGCCTGAGGCGGCGTGAGGCAGCAGCCTTGGCAAGGGGGGATACTTTAATCTATC‑CGGGGGGGGCGACGACCAATGGCGGCGGGTAAATTTATCGTTGCAGACGGGCCTGCGGGCAAAg                                                                                     <  2:1880172/150‑1 (MQ=255)
                                                               cGGTGTTCCATCGGCGACGGTTCCTGGTAGGGGGCCTGAGGTGGTGTGAGGCAGCTGCCTTGGCAAGGGGGGTTACTTTAATCTATC‑CGGGGGGGGTGACGACCAATGGCGGCGGGTAAAGTTATCGTTGCAGACGGGCCTGCGGGCATAg                                                                                     >  2:1316610/1‑151 (MQ=255)
                                                               cGGTGTTCCATCGGCGACGGTTCCTGGTAGGGGGCCTGAGGCGGCGTGAGGCAGCAGCCTTGGCAAGGGGGGATACTTTAATCTATC‑C‑GGGGGGGCGACGACCAATGGCGGCGGGTAAATTTATCGTTGCAGACGGGCCTGCGGGCAAAGa                                                                                    >  2:120967/1‑151 (MQ=255)
                                                                   gtTCCATCGGCGACGGTTCCTGGTAGGGGGCCTGAGGCGGCGTGAGGCAGCAGCCTTGGCAAGGGGGGATACTTTAATCTATC‑CGGGGGGGGCGACGACCAATGGCGGCGGGTAAATTTATCGTTGCAGACGGGCCTGCGGGCAAAGAAAt                                                                                 >  1:4250415/1‑151 (MQ=255)
                                                                      cgaTCGGCGAGGGTTCCTGGAAGGGGGCCGGAGGCGGGGGGGGGCAGCAGCAGTGGAAAGGGGGGATACTTTAATCTATC‑CGGGGGGGGCGACGACCAATGGCGGCGGGTAAATTTATCGTTGCAGACGGGCCTGCGGGCAAAGAAATGTc                                                                              <  1:285333/149‑1 (MQ=255)
                                                                       gggcGGCGACGGTGGCGGGTGGGGGGTCGGAGGGGGCGGGAGGCAGCAGCCTTGGCAGGGGGGGATACTTTGATCTATC‑GGGGGGGGGCGACGACCAATGGCGGCGGGTAAATTTATCGTTGCAGACGGGCCTGCGGGCAAAGAAATGTcc                                                                             <  2:2787689/148‑1 (MQ=255)
                                                                        aTCGGGGAGGGTTGTTGGTAGGGGGCCGGAGGGGGGGGGACGCGGCAGCGTTGGCAGGGGGGGATACTTTAATCTATCGG‑GGGGGGGCGACGACCAATGGCGGCGGGTAAATTTATCGTTGCAGACGGGCCTGCGGGCAAAGAAATGTccc                                                                            <  1:3642437/151‑1 (MQ=255)
                                                                          gggCGACGGTTGCTGGTAGGGGGTCTGAGGCGGGGTGAGGCAGCAGCCTTGGCAAGGGGGGATACTTTAATCTATC‑C‑GGGGGGGCGACGACCAATGGCGGCGGGTAAATTTATCGTTGCAGACGGGCCTGCGGGCAAAGAAATGTCCCACg                                                                         <  2:3244256/150‑1 (MQ=255)
                                                                          cGGCGACGGTTCCTGGTAGGGGGCCTGAGGCGGCGTGAGGCAGCAGCCTTGGCAAGGGGGGATACTTTAATCTATC‑CGGGGGGGGCGCCGACCAATGGCGGCGGGGACATTTATCGTTGCAGACGGGCCTGCGGGCCAAGAAATGTccccc                                                                          >  2:2655556/1‑149 (MQ=255)
                                                                           ggCGACGGTTCCTGGTAGGGGGCCTGAGGCGGCGTGAGGCAGCAGCCTTGGCAAGGGGGGATACTTTAATCTATC‑CGGGGGGGGCGACGACCAATGGCGGCGGGTAAATTTATCGTTGCAGACGGGCCTGCGGGCAAAGAAATGTCCCACg                                                                         >  1:3873342/1‑151 (MQ=255)
                                                                           ggCGACGGTTCCTGGTAGGGGGCCTGAGGCGGCGTGAGGCAGCAGCCTTGGCAAGGGGGGATACTTTAATCTATC‑CGGGGGGGGCGACGACCAATGGCGGCGGGTAAATTTATCGTTGCAGACGGGCCTGCGGGCAAAGAAATGTCCCACg                                                                         >  1:2936640/1‑151 (MQ=255)
                                                                             cGACGGTTCCTGGTAGGGGGCCTGAGGCGGCGTGAGGCAGCAGCCTTGGCAAGGGGGGATACTTTAATCTATC‑CGGGGGGGGCGACGACCAATGGCGGCGGGTAAATTTATCGTTGCAGACGGGCCTGCGGGCAAAGAAATGTCCCACGCa                                                                       >  2:810343/1‑151 (MQ=255)
                                                                                 ggTTCCTGGTAGGGGGCCTGAGGCGGCGTGAGGCAGCAGCCTTGGCAAGGGGGGATACTTTAATCTATC‑CGGGGGGGGCGACGACCAATGGCGGCGGGTAAATTTATCGTTGCAGACGGGCCTGCGGGCAAAGAAATGTCCCACGCAtggc                                                                   >  2:3377211/1‑151 (MQ=255)
                                                                                 ggTTCCTGGTAGGGGGCCTGAGGCGGCGTGAGGCAGCAGCCTTGGCAAGGGGGGATACTTTAATCTATC‑CGGGGGGGGCGACGACCAATGGCGGCGGGTAAATTTATCGTTGCAGACGGGCCTGCGGGCAAAGAAATGTCCCACGCAtggc                                                                   >  2:2387198/1‑151 (MQ=255)
                                                                                 ggTTCCTGGTAGGGGGCCTGAGGCGGCGTGAGGCAGCAGCCTTGGCAAGGGGGGATACTTTAATCTATC‑C‑GGGGGGGCGACGACCAATGGCGGCGGGTAAATTTATCGTTGCAGACGGGCCTGCGGGCAAAGAAATGTCCCACGCAtggct                                                                  >  2:4717373/1‑151 (MQ=255)
                                                                                   taccgggTAGGGGGCGTGAGGCGGCGTGAGGCAGCAGCCATGGCAAGGGGGGATACTTTAATCTATC‑CGGGGGGGGCGACGACCAATGGCGGCGGGTAAATTTATCGTTGCAGACGGGCCTGCGGGCAAAGAAATGTCCCACGCAtggctg                                                                 <  2:409832/146‑1 (MQ=255)
                                                                                     ccTGGTAGGGGGCCTGAGGCGGCGTGAGGCAGCAGCCTTGGCAAGGGGGGATACTTTAATCTATC‑CGGGGGGGGCGACGACCAATGGCGGCGGGTAAACTTATCGTTGCAGACGGGCCTGCGGGCAAAGAAATGTCCCACGCAtggctggc                                                               >  2:405721/1‑151 (MQ=255)
                                                                                     ccTGGTAGGGGGCCTGAGGCGGCGTGAGGCAGCAGCCTTGGCAAGGGGGGATACTTTAATCTATC‑C‑GGGGGGGCGACGACCAATGGCGGCGGGTAAATTTATCGTTGCAGACGGGCCTGCGGGCAAAGAAATGTCCCACGCATGGCTGGcc                                                              >  1:4422807/1‑151 (MQ=255)
                                                                                      gtGGGAGGGGGCCTGGGGCGGCGCGAGGCAGCAGCCTTGGCAAGGGGGGATACTTTAATGTATCGG‑GGGGGGGCGACGACCAATGGCGGCGGGTAAATTTATCGTTGCAGACGGGCCTGCGGGCAAAGAAATGTCCCACGCATGGCTGGcc                                                              <  2:4981611/150‑1 (MQ=255)
                                                                                       tGGTAGGGGGCCTGAGGCGGCGTGAGGCAGCAGCCTTGGCAAGGGGGGATACTTTAATCTATC‑CGGGGGGGGCGACGACCAATGGCGGCGGGTAAATTTATCGTTGCAGACGGGCCTGCGGGCAAAGAAATGTCCCACGCATGGCTGGCCt                                                             >  1:3032891/1‑151 (MQ=255)
                                                                                       ggggaGGGGGCGGGGGGCGGCGTGAGGCAGCAGCCTTGGCAAGGGGGGATACTTTAATCTAGC‑CGGGGGGGGCGACGACCAATGGCGGCGGGTAAATTTATCGTTGCAGACGGGCCTGCGGGCAAAGAAATGTCCCACGCATGGCTGGCCt                                                             <  2:1268778/147‑1 (MQ=255)
                                                                                          tAGGGGGCCTGAGGCGGCGTGAGGCAGCAGCCTTGGCAAGGGGGGATACTTTAATCTATC‑CGGGGGGGGCGACGACCAATGGCGGCGGGTAAATTTATCGTTGCAGACGGGCCTGCGGGCTAAGAAATGTCCCACGCCTGGCTGGCCTTTg                                                          >  2:2329953/1‑151 (MQ=255)
                                                                                          tAGGGGGCCTGAGGCGGCGTGAGGCAGCAGCCTTGGCAAGGGGGGATACCCTAATCTATC‑CGGGGGGGGCGACGACCAATGGCGGCGGGTAAATTTATCCTTGCAGACGGGCCTGCGGGCAAGGGAAAGGTCCACCCCTGGCTGGCCttcg                                                          >  2:963275/1‑149 (MQ=255)
                                                                                           aGGGAGCCTGAGGCGGCGTGAGACAGCATCCTTGACAAGGGGTGATACTTACTCCTACC‑CGGGGGGGGCGACGACCAATGGCGGCGGGTAAATTTATCGTTGCAGACGGGCCTGCGGGCAAAGAAATGTCCCACGCATGGCTGGCCTTTGc                                                         <  2:2768036/151‑1 (MQ=255)
                                                                                             ggggCCTGAGGCGGCGTGAGGCAGCAGCCTTGGCAAGGGGGGATACTTTAATCTATC‑CGGGGGGGGCGACGACCAATGGCGGCGGGTAAATTTATCGTTGCAGACGGGCCTGCGGGCAAAGAAATGTCCCACGCATGGCTGGCCTTTGCCt                                                       >  1:3676804/1‑151 (MQ=255)
                                                                                                  cTGAGGCGACGTGAGGCAGCAGCCTTGGCAAGGGGGGATACTTTAATCTATC‑CGGGGGGGGCGACGACCAATGGCGGCGGGTAAATTTATCGTTGCAGACGGGCCTGCGGGCAAAGACATGTCCCACGCATCGCTGGCCTTTGCCTGACgg                                                  >  2:3119563/1‑151 (MQ=255)
                                                                                                   tGAGGCGGCGTGAGGCAGCAGCCTTGGCAAGGGGGGATACTTTAATCTATC‑CGGGGGGGGCGACGACCAACGGCGGCGGGTAAATCTATCGTTGCAGACGGGCCTGCGGGCAAAGAAATGTCCCACGCATGGCTCGCCTTTGCCTGATggc                                                 >  2:1929539/1‑150 (MQ=255)
                                                                                                    gAGGCGGCGTGAGGCAGCGGCCTTGGCAAGGGGGGATACTTTGAGCGATC‑CGGGGGGGGCGACGACCATTGGCGGCGGGTAAATTTATCGTTGCAGACGGGCCTGCGGGCAAAGAAATGTCCCACGCATGGCTGGCCTTTGCCTGATGGTc                                                <  2:2021017/151‑1 (MQ=255)
                                                                                                     gggcggcGTGGGGCAGCAGCCTGGGCAAGGGGGGATACTTTAATGTATC‑CGGGGGGGGCGACGACCAATGGCGGCGGGTAAATTTATCGTTGCAGACGGGCCTGCGGGCAAAGAAATGTCCCACGCATGGCTGGCCTTTGCCTGATGGTcc                                               <  1:117632/150‑1 (MQ=255)
                                                                                                     agcgggcgggcgGCAGCGGCCGTGGCAAGGGGGGATACTTTAGTCTATC‑CGGGGGGGGCGACGACCAATGGCGGCGGGTAAATTTAGCGTTGCAGACGGGCCTGCGGGCAAAGAAATGTCCCACGCATGGCTGGCCTTTGCCTGATGGTcc                                               <  2:4837836/147‑1 (MQ=255)
                                                                                                         ggcgTGAGGCAGCAGCCTTGGCAAGGGGGGATACTTTAATCTATC‑CGGGGGGGGCGACGACCAATGGCGGCGGGTAAATTTATCGTTGCAGACGGGCCTGCGGGCAAAGAAATGTCCCACGCATGGCTGGCCTTTGCCTGATGGTCCGTGa                                           >  1:3919397/1‑151 (MQ=255)
                                                                                                           cgaggggcggcAGCGTGGGCACGGGGGGCGATGTTAGTCGATCGG‑GGGGGGGCGACGACCAATGGCGGCGGGTAAATTTATCGTTGCAGACGGGCCTGCGGGCAAAGAAATGTCCCACGCATGGCTGGCCTTTGCCTGATGGTCCGTGACg                                         <  1:4885683/146‑1 (MQ=255)
                                                                                                            gTGAGGCAGCAGCCTTGGCAAGGGGGGATACTTTAATCTATC‑CGGGGGGGGCGACGACCAATGGCGGCGGGTAAATTTATCGTTGCAGACGGGCCTGCGGGCAAAGAAATGTGCCACGCATGGCTGGCCTTTGCCTGATGGTCCGTGACgg                                        >  2:105284/1‑151 (MQ=255)
                                                                                                            gTGAGGCAGCAGCCTTGGCAAGGGGGGATACTTTAATCTATC‑CGGGGGGGGCGACGACCAATGGCGGCGGGTAAATTTATCGTTGCAGACGGGCCTGCGGGCAAAGAAATGTCCCACGCATGGCTGGCCTTTGCCTGATGGTCCGTGACgg                                        >  1:4167465/1‑151 (MQ=255)
                                                                                                              gAGGCAGCAGCCTTGGCAAGGGGGGATACTTTAATCTATC‑CGGGGGGGGCGACGACCAATGGCGGCGGGTAAATTTATCGTTGCAGACGGGCCTGCGGGCAAAGAAATGTCCCACGCATGGCTGGCCTTTGCCTGATGGTCCGTGACGGCg                                      >  1:3865381/1‑151 (MQ=255)
                                                                                                               aGGCAGCAGCCTTGGCAAGGGGGGATACTTTAATCTATC‑CGGGGGGGGCGACGACCAATGGCGGCGGGTAAATTTATCGTTGCAGACGGGCCTGCGGGCAAAGAAATGTCCCACGCATGGCTGGCCTTTGCCTGATGGTCCGTGACGGCGt                                     <  2:2900222/151‑1 (MQ=255)
                                                                                                                ggCAGCAGCCTTGGCAAGGGGGGATACTTTAATCTATC‑CGGGGGGGGCGACGACCAATGGCGGCGGGTAAATTTATCGTTGCAGACGGGCCTGCGGGCAAAGAAATGTCCCACGCATGGCTGGCCTTTGCCTGATGGTCCGTGACGGCGTg                                    >  2:3971000/1‑151 (MQ=255)
                                                                                                                   ggcagcCTTGGCAAGGGGGGATACTTTAATCTAGC‑CGGGGGGGGCGACGACCAATGGCGGCGGGTAAATTTATCGTTGCAGACGGGCCTGCGGGCAAAGAAATGTCCCACGCATGGCTGGCCTTTGCCTGATGGTCCGTGACGGCGTGCTg                                 <  1:2686317/150‑1 (MQ=255)
                                                                                                                   agcagcCTGGGAAAGGGGGGAGACTTTAATCTGTCGG‑GGGGGGGCGACGACCAATGGCGGCGGGTAAATTTATCGTTGCAGACGGGCCTGCGGGCAAAGAAATGTCCCACGCATGGCTGGCCTTTGCCTGATGGTCCGTGACGGCGTGCTg                                 <  1:2685333/151‑1 (MQ=255)
                                                                                                                      agcCTTGGCAAGGGGGGATACTTTAATCTATC‑CGGGGGGGGCGACGACCAATGGCGGCGGGTAAATTTATCGTTGCAGACGGGCCTGCGGGCAAAGAAATGTCCCACGCATGGCTGGCCTTTGCCTGATGGTCCGTGACGGCGTGCTGgcc                              >  2:1769055/1‑150 (MQ=255)
                                                                                                                       gcCTTGGCAAGGGGGGATACTTTAATCTATC‑CGGGGGGGGCGACGACCAATGGCGGCGGGTAAATTTATCGTTGCAGACGGGCCTGCGGGCAAAGAAATGTCCCACGCATGGCTGGCCTTTGCCTGATGGTCCGTGACGGCGTGCTGgcgc                             >  2:4209533/1‑151 (MQ=255)
                                                                                                                         tttgctaaGGGGGGATACGTTACGCTCTC‑CGGGGGGGGCGACGACCAATGGCGGCGGGTAAATTTATCGTTGCAGACGGGCCTGCGGGCAAAGAAATGTCCCACGCATGGCTGGCCTTTGCCTGATGGTCCGTGACGGCGTGCTGGCGCCt                           <  1:357206/145‑1 (MQ=255)
                                                                                                                         cTTGGCAAGGGGGGATACTTTAATCTATC‑CGGGGGGGGCGACGACCAATGGCGGCGGGTAAATTTATCGTTGCAGACGGGCCGGCGGGCAAAGAAATGTCCCACGCATGGCTGGCCTTTGCCTGATGGTCCGTGACGGCGTGCTGGCGccc                           >  2:3507843/1‑150 (MQ=255)
                                                                                                                            ggCAAGGGGGGATACTTTAATCTATC‑CGGGGGGGGCGACGACCAATGGCGGCGGGTAAATTTATCGTTGCAGACGGGCCTGCGGGCAAAGAAATGTCCCACGCATGGCTGGCCTTTGCCTGATGGTCCGTGACGGCGTGCTGGCGCCTGCt                        >  1:3193504/1‑151 (MQ=255)
                                                                                                                             gCAAGGGGGGATACTTTAATCTATC‑CGGGGGGGGCGACGACCAATGGCGGCGGGTAAATTTATCGTTGCAGACGGGCCTGCGGGCAAAGAAATGTCCCACGCATGGCTGGCCTTTGCCTGATGGTCCGTGACGGCGTGCTGGCGCCTGCtt                       >  2:4802508/1‑151 (MQ=255)
                                                                                                                                  gggggATACTTTAATCTATC‑CGGGGGGGGCGACGACCAATGGCGGCGGGTAAATTTATCGTTGCAGACGGGCCTGCGGGCAAAGAAATGTCCCACGCATGGCTGGCCTTTGCCTGATGGTCCGTGACGGCGTGCTGGCGCCTGCTTGTGCt                  >  2:1672806/1‑151 (MQ=255)
                                                                                                                                    gggATACTTTAATCTATC‑CGGGGGGGGCGACGACCAATGGCGGCGGGTAAATTTATCGTTGCAGACGGGCCTGCGGGCAAAGAAATGTCCCACGCATGGCTGGCCTTTGCCTGATGGTCCGTGACGGCGTGCTGGCGCCTGCTTGTGCttt                >  1:4395068/1‑151 (MQ=255)
                                                                                                                                     ggATACTTTAATCTATC‑CGGGGGGGGCGACGACCAATGGCGGCGGGTAAATTTATCGTTGCAGACGGGCCTGCGGGCAAAGAAATGTCCCACGCATGGCTGGCCTTTGCCTGATGGTCCGTGACGGCGTGCTGGCGCCTGCTTGTGCtttt               >  1:2652457/1‑151 (MQ=255)
                                                                                                                                      gATACTTTAATCTATC‑CGGGGGGGGCGACGACCAATGGCGGCGGGTAAATTTATCGTTGCAGACGGGCCTGCGGGCAAAGAAATGTCCCACGCATGGCTGGCCTTTGCCTGATGGTCCGTGACGGCGTGCTGGCGCCTGCTTGTGCTTTTc              <  1:619899/151‑1 (MQ=255)
                                                                                                                                         aCTTTAATCTATC‑CGGGGGGGGCGACGACCAATGGCGGCGGGTAAATTTATCGTTGCAGACGGGCCTGCGGGCAAAGAAATGTCCCACGCATGGCTGGCCTTTGCCTGATGGTCCGTGACGGCGTGCTGGCGCCTGCTTGTGCTTTTCggg           >  2:398333/1‑150 (MQ=255)
                                                                                                                                         aCTTTAATCTATC‑CGGGGGGGGCGACGACCAATGGCGGCGGGTAAATTTATCGTTGCAGACGGGCCTGCGGGCAAAGAAATGTCCCACGCATGGCTGGCCTTTGCCTGATGGTCCGTGACGGCGTGCTGGCGCCTGCTTGTGCTTTTCGGa           >  2:4409355/1‑151 (MQ=255)
                                                                                                                                          cTTTAATCTATC‑CGGGGGGGGCGACGACCAATGGCGGCGGGTAAATTTATCGTTGCAGACGGGCCTGCGGGCAAAGAAATGTCCCACGCATGGCTGGCCTTTGCCTGATGGTCCGTGACGGCGTGCTGGCGCCTGCTTGTGATTTTCGGAc          >  2:614221/1‑151 (MQ=255)
                                                                                                                                                tCTATC‑CGGGGGGGGCGACGACCAATGGCGGCGGGTAAATTTATCGTTGCAGACGGGCCTGCGGGCAAAGAAATGTCCCACGCATGGCTGGCCTTTGCCTGATGGTCCGTGACGGCGTGCTGGCGCCTGCTTGTGCTTTTCGGACAAGCag    >  1:3608082/1‑151 (MQ=255)
                                                                                                                                                  tATC‑CGGGGGGGGCGACGACCAATGGCGGCGGGTAAATTTATCGTTGCAGACGGGCCTGCGGGCAAAGAAATGTCCCACGCATGGCTGGCCTTTGCCTGATGGTCCGTGACGGCGTGCTGGCGCCTGCTTGTGCTTTTCGGACAAGCagga  >  2:3052857/1‑151 (MQ=255)
                                                                                                                                                        |                                                                                                                                                 
ACCAGTTTCTCCTGGCCTAGGATCTGGGTAGAAAGAAAAGTGCGCAGCGCGATCAGCGCTTCACGGTGTTCCATCGGCGACGGTTCCTGGTAGGGGGCCTGAGGCGGCGTGAGGCAGCAGCCTTGGCAAGGGGGGATACTTTAATCTATC‑C‑GGGGGGGCGACGACCAATGGCGGCGGGTAAATTTATCGTTGCAGACGGGCCTGCGGGCAAAGAAATGTCCCACGCATGGCTGGCCTTTGCCTGATGGTCCGTGACGGCGTGCTGGCGCCTGCTTGTGCTTTTCGGACAAGCAGGA  >  CP009974/3929823‑3930118

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑