Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,313,207 G→A 6.1% intergenic (‑29/+135) rbfA ← / ← infB 30s ribosome binding factor/translation initiation factor IF‑2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,313,2070GA 6.1% 420.9 / 9.5 131intergenic (‑29/+135)rbfA/infB30s ribosome binding factor/translation initiation factor IF‑2
Reads supporting (aligned to +/- strand):  ref base G (56/67);  new base A (5/3);  total (61/70)
Fisher's exact test for biased strand distribution p-value = 4.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.85e-01

Only 100 of 135 total aligned reads displayed.

CACCGGAAACGGTGGTCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCTAAAAACCTTAAGCAATGGTACG  >  NC_000913/3313058‑3313356
                                                                                                                                                     |                                                                                                                                                     
cacCGGAAACGGTGGTCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTg                                                                                                                                                       >  2:1670526/1‑150 (MQ=255)
  ccGGAAAAGGGGGTCATGATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGtt                                                                                                                                                     <  2:1040694/150‑1 (MQ=255)
   cGGAAACGGTGGTCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTa                                                                                                                                                    >  1:433548/1‑150 (MQ=255)
    ggAAACGGTGGTCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAg                                                                                                                                                   >  2:670079/1‑150 (MQ=255)
    ggAAACGGTGGTCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATACTCCTGACAAAAAAGGGGCTGTTAg                                                                                                                                                   >  1:255345/1‑150 (MQ=255)
     gAAACGGTGGTCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGc                                                                                                                                                  <  1:886139/150‑1 (MQ=255)
     gAAACGGTGGTCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGc                                                                                                                                                  >  2:124498/1‑150 (MQ=255)
          ggtggtCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCtt                                                                                                                                             >  2:1607184/1‑150 (MQ=255)
          ggtggtCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCtt                                                                                                                                             <  2:647753/150‑1 (MQ=255)
           gtggtCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCttt                                                                                                                                            >  2:115015/1‑150 (MQ=255)
           gtggtCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCttt                                                                                                                                            <  2:742311/150‑1 (MQ=255)
            tggtCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCtttt                                                                                                                                           <  2:1219320/150‑1 (MQ=255)
             ggtCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCttttt                                                                                                                                          <  2:1123406/150‑1 (MQ=255)
             ggtCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCttttt                                                                                                                                          >  1:809236/1‑150 (MQ=255)
             ggtCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCttttt                                                                                                                                          <  1:764910/150‑1 (MQ=255)
             ggtCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTAACAGCCCCttttt                                                                                                                                          <  1:1118547/150‑1 (MQ=255)
                 atcatGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTaaaa                                                                                                                                      <  2:513989/150‑1 (MQ=255)
                  tcatGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAaat                                                                                                                                     >  2:897928/1‑150 (MQ=255)
                  tcatGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAaat                                                                                                                                     >  2:904023/1‑150 (MQ=255)
                   catGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAaatt                                                                                                                                    >  2:1138467/1‑150 (MQ=255)
                       cccAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAaattaatt                                                                                                                                <  1:802793/150‑1 (MQ=255)
                       cccAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAaattaatt                                                                                                                                >  2:122927/1‑150 (MQ=255)
                        ccAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAAttt                                                                                                                               >  1:1489600/1‑150 (MQ=255)
                         cAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTc                                                                                                                              <  2:425221/150‑1 (MQ=255)
                           ggCGAGGATCTTTAATTTCACGCTGGAGGATGAGAGCTATCTCTTTTTGCATTTCCAGCGCTACGCGCTGGGGGCTACCATATTCTTTCGCCATAAAAAATTCTCCCGACAAAAAAGGGGCTGTTATGCCCTTTTTAAAATTAATTTCAg                                                                                                                            <  2:1632330/150‑1 (MQ=255)
                             cGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGt                                                                                                                          <  1:582330/150‑1 (MQ=255)
                             cGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGt                                                                                                                          <  1:634228/150‑1 (MQ=255)
                              gAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTg                                                                                                                         <  2:148510/150‑1 (MQ=255)
                              gAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTg                                                                                                                         >  2:1530793/1‑150 (MQ=255)
                               aGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTgg                                                                                                                        >  2:617925/1‑150 (MQ=255)
                                       tAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACTCGCTGCGGGTGACCAAATTCATTCGACATAATAAATTCTCCTGACAAAAAAGGGGCTGTAAGCCCCATTTTAAAATTAATATCAGGTGGAATGGcag                                                                                                                >  1:141427/1‑148 (MQ=255)
                                         aTTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGtt                                                                                                              >  2:1323090/1‑150 (MQ=255)
                                          tttCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTc                                                                                                             <  1:1627939/150‑1 (MQ=255)
                                             cccGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACg                                                                                                          <  1:1178345/148‑1 (MQ=255)
                                             cACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACg                                                                                                          <  2:826788/150‑1 (MQ=255)
                                              aCGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGt                                                                                                         <  1:1115697/150‑1 (MQ=255)
                                              aCGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGt                                                                                                         >  2:1683061/1‑150 (MQ=255)
                                                gCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTg                                                                                                       >  2:1411875/1‑150 (MQ=255)
                                                 cTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGa                                                                                                      <  2:64230/150‑1 (MQ=255)
                                                    cAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCt                                                                                                   <  2:691407/150‑1 (MQ=255)
                                                    cAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCt                                                                                                   <  1:60298/150‑1 (MQ=255)
                                                    cAGGATGAGAGCGATCTCTTTATGCAATTCCTGCGCTACGCGCGGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAAATTCAGGTGGCAGGGCTGTTCACGTTGACCt                                                                                                   <  2:881585/150‑1 (MQ=255)
                                                      ggATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGACATAATAAATTCTACTGACAAAAAAGGGGCTGTTAGACCCTTTTTAAAATTAATTTCAGGTGGAAGGTCTGTTAACGTTGACCTGa                                                                                                 >  1:120731/1‑150 (MQ=255)
                                                          gagaGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAg                                                                                             >  2:365975/1‑150 (MQ=255)
                                                          gagaGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAg                                                                                             >  1:19918/1‑150 (MQ=255)
                                                            tagCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGAc                                                                                           <  2:1659570/149‑1 (MQ=255)
                                                            gagCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGAc                                                                                           <  1:155912/150‑1 (MQ=255)
                                                              gctaTCTATTTTTGCATTTCCTGCGCTACGCGCTGCGGGAGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGAcgc                                                                                         <  2:1096822/147‑1 (MQ=255)
                                                              gCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGAcgc                                                                                         >  1:1697729/1‑150 (MQ=255)
                                                                gATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACgcgc                                                                                       >  1:127812/1‑150 (MQ=255)
                                                                gATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACgcgc                                                                                       >  2:316845/1‑150 (MQ=255)
                                                                 aTCTCTTTTTGCATTTCCTGCTCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAATAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGcc                                                                                      <  1:1150015/150‑1 (MQ=255)
                                                                    tctTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGc                                                                                   <  1:1256095/150‑1 (MQ=255)
                                                                           gCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACAt                                                                            >  2:30299/1‑150 (MQ=255)
                                                                               ttCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAgg                                                                        >  2:340500/1‑150 (MQ=255)
                                                                                tCCTGCGCTACGCGCGGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGc                                                                       <  2:1505286/150‑1 (MQ=255)
                                                                                     cgcTACGCGCTGCTGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAATTTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATcc                                                                  <  2:346275/150‑1 (MQ=255)
                                                                                       cTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCAt                                                                >  1:687664/1‑150 (MQ=255)
                                                                                              cTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGc                                                                                                                  >  1:267864/1‑93 (MQ=255)
                                                                                              cTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGc                                                                                                                  <  2:267864/93‑1 (MQ=255)
                                                                                               tGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATg                                                        <  1:509418/150‑1 (MQ=255)
                                                                                                gCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGc                                                       <  2:19918/150‑1 (MQ=255)
                                                                                                     cGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGt                                                  <  2:1489600/150‑1 (MQ=255)
                                                                                                     cGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGt                                                  >  1:1106352/1‑150 (MQ=255)
                                                                                                      gACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCTCATCAGGCAATCCATGCCGGATGCAGCGTa                                                 <  2:919074/150‑1 (MQ=255)
                                                                                                      gACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCTATCCATGCCGGATGCAGCGTa                                                 <  1:686218/150‑1 (MQ=255)
                                                                                                         cAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAAc                                              <  2:1028728/150‑1 (MQ=255)
                                                                                                         cAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAAc                                              >  1:357060/1‑150 (MQ=255)
                                                                                                            aTTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGcc                                           <  2:687664/150‑1 (MQ=255)
                                                                                                                  tCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATcc                                     <  1:1152690/150‑1 (MQ=255)
                                                                                                                   cGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTAACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATccc                                    >  1:285415/1‑150 (MQ=255)
                                                                                                                    gCCATAATAACTTCTCCTGACAAAAACGGGGCGGTTAGCCCCTTTTTAAAATTAATTTCAGGTCGAAGGGCTGTTCACGTTTACCTTATAAGACGCGCCATCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCg                                   <  2:667687/150‑1 (MQ=255)
                                                                                                                    gCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCg                                   <  2:809388/150‑1 (MQ=255)
                                                                                                                     ccATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGc                                  >  2:611625/1‑150 (MQ=255)
                                                                                                                      cATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCa                                 >  1:476015/1‑150 (MQ=255)
                                                                                                                       ataataAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCAt                                >  2:692053/1‑150 (MQ=255)
                                                                                                                       ataataAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCAt                                >  2:119359/1‑150 (MQ=255)
                                                                                                                        taataaATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATg                               <  2:433548/150‑1 (MQ=255)
                                                                                                                         aataaATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATgg                              >  1:248079/1‑150 (MQ=255)
                                                                                                                          ataaATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGa                             <  1:1607184/150‑1 (MQ=255)
                                                                                                                            aaaTTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTTGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAAc                           <  2:939321/150‑1 (MQ=255)
                                                                                                                            aaaTTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCTGATCCAGCGTAAACGACTTATCCCGCATGGAAc                           >  2:750983/1‑150 (MQ=255)
                                                                                                                            aaaTTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAAc                           <  2:198685/150‑1 (MQ=255)
                                                                                                                              aTTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAAccc                         <  1:1670526/150‑1 (MQ=255)
                                                                                                                              aTTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGAGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGCATGCAGCGTAAACGCCTTATCCCGCATGGAAccc                         >  2:685598/1‑150 (MQ=255)
                                                                                                                                 ctcCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCTaa                      >  2:514043/1‑150 (MQ=255)
                                                                                                                                   ccTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCTaaaa                    <  2:604915/150‑1 (MQ=255)
                                                                                                                                      ttttaaaaaGGGGCTAACAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCTAAAAAcc                 >  2:323389/5‑150 (MQ=255)
                                                                                                                                      gACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCTAAAAAcc                 <  1:904023/150‑1 (MQ=255)
                                                                                                                                      gACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCTAAAAAcc                 <  1:897928/150‑1 (MQ=255)
                                                                                                                                       tttaaaaaGGGGCTAACAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCTAAAACCCt                >  2:755251/4‑150 (MQ=39)
                                                                                                                                        cAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGAACCCTAAAAACCTTa              <  2:782391/150‑1 (MQ=255)
                                                                                                                                         aaaaaaGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCTAAAAACCTTa              >  2:1707837/1‑150 (MQ=255)
                                                                                                                                           aaaaGTGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATACATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCTAAAAACCTTAAg            <  2:653907/150‑1 (MQ=255)
                                                                                                                                           aaaaGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCTAAAAACCTTAAg            <  2:1495953/150‑1 (MQ=255)
                                                                                                                                           aaaaGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCTAAAAACCTTAAg            >  1:925442/1‑150 (MQ=255)
                                                                                                                                           aaaaGGGGCTAACAGCCCCTTTTTAAATTTAATTTCAGGTGGAAGGGCTGTTCACTTTGACCTGATAAGAAGCGCCAGCGTCACATCAGGCAATCCATTCAGGATGATGCGTAAACGCCTTATCCCGCATGGATCACTATAAACGTTAAg            >  1:929873/1‑150 (MQ=255)
                                                                                                                                            aaaGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCTAAAAACCTTAAGc           >  2:1342389/1‑150 (MQ=255)
                                                                                                                                               ggggCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCTAAAAACCTTAAGCAAt        <  1:1683061/150‑1 (MQ=255)
                                                                                                                                                     gTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCTAAAAACCTTAAGCAATGGTACg  >  1:1470064/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CACCGGAAACGGTGGTCATCATGCCCAGGCGAGGATCTTTAATTTCACGCTGCAGGATGAGAGCGATCTCTTTTTGCATTTCCTGCGCTACGCGCTGCGGGCGACCAAATTCTTTCGCCATAATAAATTCTCCTGACAAAAAAGGGGCTGTTAGCCCCTTTTTAAAATTAATTTCAGGTGGAAGGGCTGTTCACGTTGACCTGATAAGACGCGCCAGCGTCACATCAGGCAATCCATGCCGGATGCAGCGTAAACGCCTTATCCCGCATGGAACCCTAAAAACCTTAAGCAATGGTACG  >  NC_000913/3313058‑3313356

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