Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 933,156 G→A 5.3% intergenic (+67/‑68) lrp → / → ftsK DNA‑binding transcriptional dual regulator Lrp/cell division DNA translocase FtsK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913933,1560GA 5.3% 474.8 / 9.8 152intergenic (+67/‑68)lrp/ftsKDNA‑binding transcriptional dual regulator Lrp/cell division DNA translocase FtsK
Reads supporting (aligned to +/- strand):  ref base G (72/68);  new base A (3/5);  total (77/75)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.24e-01

Only 100 of 158 total aligned reads displayed.

TGCCTGGCGTCAATGACACACGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGCTGA  >  NC_000913/933007‑933305
                                                                                                                                                     |                                                                                                                                                     
  ccTGGCGGCAACGACACACGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGGCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTggg                                                                                                                                                     <  2:296053/150‑1 (MQ=255)
    tGGCGTCAATGACACACGGACATACCTTGTTATGGAAGAAGTCAAGCAGAGTAATAGTCTGGTTATTAAGACGCGCTAACACGGAAAAGGAGCAAAATCGGCGTAATTACATTACACTCATGTTAATCCATACAGCAACAGTACTGGGGt                                                                                                                                                   >  2:2109743/1‑150 (MQ=255)
       cGTCAATGACACACGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAAc                                                                                                                                                >  1:568092/1‑150 (MQ=255)
         tCAATGACACACGGACATACGTTGTTATGGAAGAAGTCAAGGAGATTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCt                                                                                                                                              <  2:1849152/150‑1 (MQ=255)
           actGACAGACGGACATACGGTGTTATGGAAGACGTCAAGCAGAGTAATCGTCTGGGTAGTAAGACGGGCTACCACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTgg                                                                                                                                            <  2:1568731/148‑1 (MQ=255)
           aaTGACACACGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTgg                                                                                                                                            <  2:2376356/150‑1 (MQ=255)
            aTGACACACGGACATACGTTGTTATGGAAGAACTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACATCTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTCGt                                                                                                                                           >  2:1716100/1‑150 (MQ=255)
              gACACACGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTAc                                                                                                                                         <  2:568092/150‑1 (MQ=255)
                 acacGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACCAGGTTACCCCAGTACtgt                                                                                                                                      <  1:415012/150‑1 (MQ=255)
                  cacGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACtgtt                                                                                                                                     >  1:347653/1‑150 (MQ=255)
                   acGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACCAGGTACCCCAGTACtgttgc                                                                                                                                   <  2:1427339/150‑2 (MQ=255)
                      gACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTcc                                                                                                                                 <  1:1784464/150‑1 (MQ=255)
                      gACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACCAGGTTACCCCAGTACtgttgctg                                                                                                                                 <  1:1816967/150‑4 (MQ=255)
                         aTACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGtt                                                                                                                              <  2:64641/150‑1 (MQ=255)
                          tACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGttt                                                                                                                             >  1:128298/1‑150 (MQ=255)
                           aCGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGtttt                                                                                                                            >  2:2082207/1‑150 (MQ=255)
                            cGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTa                                                                                                                           <  1:1937645/150‑1 (MQ=255)
                            cGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTa                                                                                                                           >  1:2145206/1‑150 (MQ=255)
                             gttgttATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAg                                                                                                                          <  1:1804893/150‑1 (MQ=255)
                                gttATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCAt                                                                                                                       >  1:1977472/1‑150 (MQ=255)
                                 ttATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATc                                                                                                                      <  1:421185/150‑1 (MQ=255)
                                    tGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCggg                                                                                                                   >  2:87964/1‑150 (MQ=255)
                                    tGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCggg                                                                                                                   <  1:574756/150‑1 (MQ=255)
                                       aagaagTCAAGCAGATTAATCGTCTGGTTATTAAGACGCGCTAAGACGTAACAGGTTCAAAATCGGCGTATTTTTATTACACTCCTGTTAATCCATACAGCAACAGTACTTGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAg                                                                                                                <  2:2452890/150‑1 (MQ=255)
                                       aagaagTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAg                                                                                                                >  2:2299434/1‑150 (MQ=255)
                                       aagaagGCAAGGAGAGTTATCGTCTGGTTATTAAGACGCGCTAACACGTAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGTGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAg                                                                                                                <  2:1152706/150‑1 (MQ=255)
                                        agaagTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGAAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAgg                                                                                                               >  1:2099760/1‑150 (MQ=255)
                                            gTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGaaaa                                                                                                           <  1:714097/150‑1 (MQ=255)
                                            gTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGaaaa                                                                                                           >  1:1466742/1‑150 (MQ=255)
                                            gTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGaaaa                                                                                                           >  1:19796/1‑150 (MQ=255)
                                              cAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGc                                                                                                         <  1:35668/150‑1 (MQ=255)
                                               aaGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGcc                                                                                                        <  2:1133057/150‑1 (MQ=255)
                                               aaGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGcc                                                                                                        >  1:1157387/1‑150 (MQ=255)
                                               aaGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGcc                                                                                                        <  1:2019423/150‑1 (MQ=255)
                                                aGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCt                                                                                                       >  2:656625/1‑150 (MQ=255)
                                                aGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCt                                                                                                       >  2:711607/1‑150 (MQ=255)
                                                  cAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGCCAGGAAAAGCCTGt                                                                                                     >  2:2408843/1‑150 (MQ=255)
                                                    gagTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTaa                                                                                                   >  2:1740219/1‑150 (MQ=255)
                                                    gagTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTaa                                                                                                   >  2:1028644/1‑150 (MQ=255)
                                                      gTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAAcc                                                                                                 >  1:2143500/1‑150 (MQ=255)
                                                         aTCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTTCACTCGTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTgg                                                                                              <  1:1001975/150‑1 (MQ=255)
                                                           cGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGgag                                                                                            <  1:12903/150‑1 (MQ=255)
                                                            gTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGgaga                                                                                           <  2:337721/150‑1 (MQ=255)
                                                             tCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGagag                                                                                          >  2:993628/1‑150 (MQ=255)
                                                               tGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGcc                                                                                        <  1:521532/150‑1 (MQ=255)
                                                                ggTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCt                                                                                       <  2:347653/150‑1 (MQ=255)
                                                                     ttaAGACGCGCTAACACGGATCAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCtt                                                                                  <  1:2193602/150‑1 (MQ=255)
                                                                      taAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTg                                                                                 >  2:1336034/1‑150 (MQ=255)
                                                                       aaGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGa                                                                                <  2:475349/150‑1 (MQ=255)
                                                                        aGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAg                                                                               >  1:760288/1‑150 (MQ=255)
                                                                        aGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAg                                                                               >  1:577372/1‑150 (MQ=255)
                                                                         gACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCAAACAGCAACAGTACTCGGGTAACCTCGTACTGTTCTCCGTTTTAGCATCGGGCAGGAAAAGCCTCTATCCTCCAGAGCCATTCTTGatc                                                                              >  2:2291145/1‑148 (MQ=255)
                                                                               cTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGaa                                                                        >  2:1621073/1‑150 (MQ=255)
                                                                               cTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTAATGGGGTAACCTGGTACTGTTGTCCGATTTAAAATCGTGCAGGAAAATCCTGTAACATGGAAAACATATCTTGAGCAAGaaa                                                                        >  2:1047709/1‑150 (MQ=255)
                                                                                tAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAAt                                                                       <  1:479848/150‑1 (MQ=255)
                                                                                 aaCACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATa                                                                      <  1:626286/150‑1 (MQ=255)
                                                                                   cacGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACa                                                                    >  1:1231011/1‑150 (MQ=255)
                                                                                         aCAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAg                                                              <  1:1130195/150‑1 (MQ=255)
                                                                                         aCAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAg                                                              <  1:617004/150‑1 (MQ=255)
                                                                                             gTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACaa                                                          >  1:1500582/1‑150 (MQ=255)
                                                                                                cAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAaaga                                                       >  2:553659/1‑150 (MQ=255)
                                                                                                cAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAaaga                                                       >  2:2009745/1‑150 (MQ=255)
                                                                                                  aaaTCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAaagaag                                                     <  1:1587433/150‑1 (MQ=255)
                                                                                                    aTCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTc                                                   <  2:809354/150‑1 (MQ=255)
                                                                                                       ggCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTcaca                                                >  1:2440305/1‑150 (MQ=255)
                                                                                                       ggCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTcaca                                                >  1:2053738/1‑150 (MQ=255)
                                                                                                         cGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACAtt                                              >  2:952028/1‑150 (MQ=255)
                                                                                                          gTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTGGCATCGGGCAGGAAAAGCCTGTAACCTTGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTg                                             <  2:116392/150‑1 (MQ=255)
                                                                                                              tttGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACaa                                         >  1:971977/1‑150 (MQ=255)
                                                                                                                   ttACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTa                                    <  1:1691330/150‑1 (MQ=255)
                                                                                                                   ttACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTa                                    <  2:2167938/150‑1 (MQ=255)
                                                                                                                     acacTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAg                                  <  1:68469/150‑1 (MQ=255)
                                                                                                                      cacTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGt                                 <  2:854675/150‑1 (MQ=255)
                                                                                                                       acTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTa                                <  1:1787766/150‑1 (MQ=255)
                                                                                                                            tGTTAATCCATACAGCAACAGTACTGGGGTATCCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGgc                           <  1:1275259/150‑1 (MQ=255)
                                                                                                                            tGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTCGCGgc                           <  1:336699/150‑1 (MQ=255)
                                                                                                                            tGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTCTAACCTGGAGAGCCTTTCTTGAGCAAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGATAAGTAGCGgc                           >  2:2131706/1‑150 (MQ=255)
                                                                                                                                    cATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCt                   >  1:2403439/1‑150 (MQ=255)
                                                                                                                                      tACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTc                 >  1:1205021/1‑150 (MQ=255)
                                                                                                                                      tACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTc                 >  2:537534/1‑150 (MQ=255)
                                                                                                                                       cggacAACAGTACCAGGTTACCCCAGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCt                >  1:220268/5‑150 (MQ=255)
                                                                                                                                       aCAGCAACAGTACTGGGGTAACCTGTTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCt                <  1:1716100/150‑1 (MQ=255)
                                                                                                                                        ggacAACAGTACCAGGTTACCCCAGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTg               >  1:1662742/4‑150 (MQ=255)
                                                                                                                                        cAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTg               >  1:2335914/1‑150 (MQ=255)
                                                                                                                                        cAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTg               >  2:516510/1‑150 (MQ=255)
                                                                                                                                        cAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTg               >  1:1599795/1‑150 (MQ=255)
                                                                                                                                        cAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTg               >  2:2454187/1‑150 (MQ=255)
                                                                                                                                        cAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTg               >  1:689123/1‑150 (MQ=255)
                                                                                                                                         aGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTgg              <  1:993628/150‑1 (MQ=255)
                                                                                                                                          gCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTgag             >  2:1448676/1‑148 (MQ=255)
                                                                                                                                          gCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGa             >  2:1514777/1‑150 (MQ=255)
                                                                                                                                          gCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGa             <  1:1250477/150‑1 (MQ=255)
                                                                                                                                          gCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGACTGTAAACTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCGCATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGa             >  2:979572/1‑150 (MQ=255)
                                                                                                                                          gCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTATTAGCATCGGGCAGGAAAAGACTGTAACATAGAGAGCATTTATTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGa             >  2:1399570/1‑150 (MQ=255)
                                                                                                                                           cAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAg            >  1:2462622/1‑150 (MQ=255)
                                                                                                                                           cAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAg            <  1:2186896/150‑1 (MQ=255)
                                                                                                                                             aCAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTgg              <  2:785557/146‑1 (MQ=255)
                                                                                                                                             aCAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTgg              >  1:785557/1‑146 (MQ=255)
                                                                                                                                                gTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGtt       <  2:1794560/150‑1 (MQ=255)
                                                                                                                                                     ggggTAACCTGGTACTGTTCTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGCTGa  >  2:1871302/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGCCTGGCGTCAATGACACACGGACATACGTTGTTATGGAAGAAGTCAAGCAGAGTAATCGTCTGGTTATTAAGACGCGCTAACACGGAACAGGTGCAAAATCGGCGTATTTTGATTACACTCCTGTTAATCCATACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGCTGA  >  NC_000913/933007‑933305

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: