Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,377,770 A→C 5.8% intergenic (‑350/+45) sspA ← / ← rpsI stringent starvation protein A, phage P1 late gene activator, RNAP‑associated acid‑resistance protein, inactive glutathione S‑transferase homolog/30S ribosomal subunit protein S9

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,377,7700AC 5.8% 670.9 / 24.4 189intergenic (‑350/+45)sspA/rpsIstringent starvation protein A, phage P1 late gene activator, RNAP‑associated acid‑resistance protein, inactive glutathione S‑transferase homolog/30S ribosomal subunit protein S9
Reads supporting (aligned to +/- strand):  ref base A (90/88);  new base C (7/4);  total (97/92)
Fisher's exact test for biased strand distribution p-value = 5.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

Only 100 of 195 total aligned reads displayed.

GGAAAAAGAGATGCGCTTTAGTCTGAAATAGTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTC  >  NC_000913/3377621‑3377919
                                                                                                                                                     |                                                                                                                                                     
ggAAAAAGAGATGCGCTTTAGTCTGAAATAGTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGa                                                                                                                                                       >  2:2523867/1‑150 (MQ=255)
ggAAAAAGAGATGCGCTTTAGTCTGAAATAGTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGATATAAAAAAACCCGCCGa                                                                                                                                                       >  1:1176217/1‑150 (MQ=255)
 gAAAAAGAGATGCGCTTTAGTCTGAAATAGTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGaa                                                                                                                                                      <  2:629432/150‑1 (MQ=255)
    aaaGAGATGCGCTTTAGTCTGAAATAGTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAcgcg                                                                                                                                                   >  2:2520259/1‑150 (MQ=255)
    aaaGAGATGCGCTTTAGTCTGAAATAGTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCg                                                                                                                                                   >  2:2520282/1‑150 (MQ=255)
     aaGAGATGCGCTTTAGTCTGAAATAGTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCgg                                                                                                                                                  <  2:1221129/150‑1 (MQ=255)
        agaTGCGCTTTAGTCTGAAATAGTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGtt                                                                                                                                               <  2:1620655/150‑1 (MQ=255)
         gaTGCGCTTTAGTCTGAAATAGTTGACTTAGTCCCATATTGGCGATGTGGTTTTTGTTTTACCGGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCAGCGTCAGAGTTTAGCTATAAAAAAACCCGCCGAAGCGGGttt                                                                                                                                              >  1:1125945/1‑150 (MQ=255)
          aTGCGCTTTAGTCTGAAATAGTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGtttt                                                                                                                                             <  2:1766891/150‑1 (MQ=255)
           tGCGCTTTAGTCTGAAATAGTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGttttt                                                                                                                                            <  1:2698207/150‑1 (MQ=255)
            gcgcTTTAGTCTGAAATAGTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGtttttt                                                                                                                                           <  2:2470123/150‑1 (MQ=255)
            gcgcTTTAGTCTGAAATAGTTGACTTAGTCCCTTAGTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGtttttt                                                                                                                                           <  2:1511030/150‑1 (MQ=255)
               cTTTAGTCTGAAATAGTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGa                                                                                                                                        <  2:2203560/150‑1 (MQ=255)
                tttAGTCTGAAATAGTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCTTCGGCGGGtttttttata                                                                                                                                       <  1:1698699/150‑5 (MQ=255)
                  tAGTCTGAAATAGTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCGGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGaaac                                                                                                                                     <  2:473837/150‑2 (MQ=255)
                        gAAATAGTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATGTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGttt                                                                                                                               >  2:4711/1‑150 (MQ=255)
                         aaaTAGTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGtttt                                                                                                                              <  1:370216/150‑1 (MQ=255)
                             aGTTGACTTAGTCCCTTATTGGCGTTGTGTTTTTTTTTTTACCTGTCTGTCAGGTGGCTGAAAAAAGCAACTTTCCAGTTTTTACGCTTATACAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGc                                                                                                                          <  1:1018484/150‑1 (MQ=255)
                              gTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGcc                                                                                                                         >  2:2890791/1‑150 (MQ=255)
                              gTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGcc                                                                                                                         >  2:1078776/1‑150 (MQ=255)
                              gTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGcc                                                                                                                         >  1:670466/1‑150 (MQ=255)
                              gTTGACATAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGcc                                                                                                                         <  1:818382/150‑1 (MQ=255)
                                 gACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTGAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGTTTCAGTTTTTTGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTATGCCGGa                                                                                                                      >  1:2235030/1‑150 (MQ=255)
                                 gACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGa                                                                                                                      <  2:229231/150‑1 (MQ=255)
                                 gACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGa                                                                                                                      <  2:1128881/150‑1 (MQ=255)
                                    ttAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACATTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCa                                                                                                                   >  1:732063/1‑150 (MQ=255)
                                        tCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAg                                                                                                               >  2:2432310/1‑150 (MQ=255)
                                          ccTTATTGGCGATGTGGTTTTTGTTTTACCAGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAAACCGCCGAAGCGGGTTTTTTCGAAAATTTTTTTATTCCGGAGCAGAAGcc                                                                                                             >  1:1278645/1‑150 (MQ=255)
                                          ccTTATTGGCGATGTGGTTTTTGTTTTAACGGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCGGATTCAGATTTTAGCTAAAAAAAAAACCGCCGAAGCGGGTTTTTTCGAAAATTTTTTTCTGCCGGAGCAGAAGcc                                                                                                             >  1:2632653/1‑150 (MQ=255)
                                           cTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCa                                                                                                            >  2:231730/1‑150 (MQ=255)
                                           cTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTAAGCTGATTCAGATTTTAGCTATAAAAAAAACCCCCGAAGCGGGTTTTTTAGAAAATTGTTTTATGCCGGAGCAGAAGCCa                                                                                                            >  1:2480202/1‑150 (MQ=255)
                                           cTTATTGGCGATGTGGTTTTTGTTTTAACAGTCTGTCAGGTGGCAGCAAAAAGCAACATTCCAGATTTTACGCTGATTCAGATTTTAGATATAAAAAAACCCGCCGAAGCGGGTTTTTTAGAAAATTGTTTTATGCCGGAGCAGAAGCCa                                                                                                            >  1:1814924/1‑150 (MQ=255)
                                            ttATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCaa                                                                                                           <  1:664293/150‑1 (MQ=255)
                                            ttATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCaa                                                                                                           <  2:1804955/150‑1 (MQ=255)
                                              aTTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAAtt                                                                                                         <  2:1858317/150‑1 (MQ=255)
                                                tGGCGATGTGGTTCTTGTATTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTaa                                                                                                       <  2:2532960/150‑1 (MQ=255)
                                                  gCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACg                                                                                                     >  2:1248162/1‑150 (MQ=255)
                                                  gCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACg                                                                                                     >  2:669217/1‑150 (MQ=255)
                                                   cGATGTGGTTTTTTTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTAAGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTATGCCGGAGCAGAAGCCAATTAACGt                                                                                                    >  1:1810058/1‑150 (MQ=255)
                                                   cGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGt                                                                                                    >  2:1895943/1‑150 (MQ=255)
                                                       gtgGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTgg                                                                                                <  2:1445936/150‑1 (MQ=255)
                                                        tgGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGga                                                                                               <  1:581197/150‑1 (MQ=255)
                                                        tgGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGga                                                                                               <  1:581252/150‑1 (MQ=255)
                                                        tgGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGga                                                                                               <  1:1778631/150‑1 (MQ=255)
                                                        tgGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGga                                                                                               <  1:2352030/150‑1 (MQ=255)
                                                           tttttGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGaa                                                                                            >  2:1396015/1‑150 (MQ=255)
                                                               tGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGc                                                                                        <  1:2573645/150‑1 (MQ=255)
                                                                 ttttACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCgg                                                                                      >  2:330310/1‑150 (MQ=255)
                                                                 ttttACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCgg                                                                                      >  2:401961/1‑150 (MQ=255)
                                                                     aCCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTATCGAAAATTGTTTTCTGCCGGAGCTGAAGGCAATTAACGATTGGAGAACTGCGgacga                                                                                  >  2:1739358/1‑150 (MQ=255)
                                                                       ctgtctgtCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCATATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAATATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGgacgacg                                                                                <  1:2221379/150‑1 (MQ=255)
                                                                       ctgtctgtCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGgacgacg                                                                                >  1:2896932/1‑150 (MQ=255)
                                                                       ctgtctgtCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGgacgacg                                                                                >  1:703616/1‑150 (MQ=255)
                                                                       ctgtctgtCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGgacgacg                                                                                >  2:2808460/1‑150 (MQ=255)
                                                                       ctgtctgtCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGgacgacg                                                                                >  1:2882234/1‑150 (MQ=255)
                                                                       ctgtctgtCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTGTGGAGAACTGCGgacgacg                                                                                >  2:1068475/1‑150 (MQ=255)
                                                                       ctgtctgtCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGGTTGGAGAACTGCGgacgacg                                                                                >  1:2894672/1‑150 (MQ=255)
                                                                       ctgtctgtCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTATACGCTGATACAGATTTTAGCTATAAAAAAAAAAGCCGAAGCGGGATATTTCGAAAATAGTTAACTGCAGTAGCAGAAGCCAATAAACGTTTGGAGAAATGCGgaagacc                                                                                >  2:2166707/1‑149 (MQ=255)
                                                                       ctgtctgtCAGGTGGCAGAAAAAAGCAACTTTCCAGTTTTTACGCTGATTCATATTTTAGCTATAAAAAAAGCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGgacgacg                                                                                <  1:2220662/150‑1 (MQ=255)
                                                                        tgtctgtcAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGgacgacga                                                                               >  2:2771380/1‑150 (MQ=255)
                                                                        tgtctgtcAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGgacgacga                                                                               >  2:2248634/1‑150 (MQ=255)
                                                                         gtctgtcAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGgacgacgac                                                                              <  1:700974/150‑1 (MQ=255)
                                                                         gtctgtcAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGgacgacgac                                                                              <  1:1099885/150‑1 (MQ=255)
                                                                         gtctgtcAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATGTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGgacgacgac                                                                              >  2:1484119/1‑150 (MQ=255)
                                                                            tgtcAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTg                                                                           <  2:1598959/150‑1 (MQ=255)
                                                                                  gTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTAc                                                                     <  1:441875/150‑1 (MQ=255)
                                                                                   tGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTGTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACg                                                                    >  2:501678/1‑150 (MQ=255)
                                                                                     gcagcaAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCa                                                                  >  2:2660072/1‑150 (MQ=255)
                                                                                     gcagcaAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCa                                                                  >  1:1564314/1‑150 (MQ=255)
                                                                                       agcaAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGa                                                                >  2:1710647/1‑150 (MQ=255)
                                                                                          aaaaaGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCg                                                             >  1:1022796/1‑150 (MQ=255)
                                                                                          aaaaaGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCg                                                             >  2:855824/1‑150 (MQ=255)
                                                                                           aaaaGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGa                                                            >  1:187529/1‑150 (MQ=255)
                                                                                           aaaaGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGa                                                            <  1:1895943/150‑1 (MQ=255)
                                                                                            aaaGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGAc                                                           >  1:69119/1‑150 (MQ=255)
                                                                                            aaaGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGAc                                                           >  1:2531129/1‑150 (MQ=255)
                                                                                                cAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTc                                                       <  2:2408494/150‑1 (MQ=255)
                                                                                                 aaCTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCt                                                      >  2:2103928/1‑150 (MQ=255)
                                                                                                      tCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGt                                                 >  2:514687/1‑150 (MQ=255)
                                                                                                          gTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCaa                                             >  1:2794875/1‑150 (MQ=255)
                                                                                                           tttttACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAAc                                            <  1:89317/150‑1 (MQ=255)
                                                                                                           tttttACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAAc                                            <  1:2749134/150‑1 (MQ=255)
                                                                                                                       ttCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTc                                <  1:986489/150‑1 (MQ=255)
                                                                                                                       ttCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTc                                <  2:754219/150‑1 (MQ=255)
                                                                                                                        tCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCa                               >  1:1181944/1‑150 (MQ=255)
                                                                                                                         cAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCAc                              >  2:1524803/1‑150 (MQ=255)
                                                                                                                         cAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCAc                              >  1:1431947/1‑150 (MQ=255)
                                                                                                                                 aGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACg                      <  1:2653572/150‑1 (MQ=255)
                                                                                                                                     ataAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGc                  >  1:2683032/1‑150 (MQ=255)
                                                                                                                                        aaaaaaCCCGCTTCGGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAg               >  1:119284/1‑150 (MQ=255)
                                                                                                                                        aaaaaaCCCGCTTCGGCGGGTTTTTTAGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATAAAACTTTGGAGAACAGCGGACGACGACGTGCTTTACGCAGACCGACATTCTAACGTTCAACCTGACGAGCGTCAAGAGTAACGAAGCCAg               >  1:794616/1‑150 (MQ=255)
                                                                                                                                         aaaaaCCCGCTTCGGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGc              >  1:774285/1‑150 (MQ=255)
                                                                                                                                         aaaaaCCCGCTTCGGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGc              >  1:1759066/1‑150 (MQ=255)
                                                                                                                                           aaaCCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGATGAGCGTCACGAGTAACGAAGCCAGCtt            <  1:2429923/150‑1 (MQ=255)
                                                                                                                                             aCCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTa          <  1:1248162/150‑1 (MQ=255)
                                                                                                                                                 gccgAAGCGGGTTTTTTCGAAAATTGTTTTCTGTCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTAACGTTCAACCTGACGAGCGTCACGAGAAACGAAGCCAGCATTACGCa      >  1:2205366/1‑150 (MQ=255)
                                                                                                                                                  ccgAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAg     <  1:2126903/150‑1 (MQ=255)
                                                                                                                                                  ccgAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAg     >  2:2836446/1‑150 (MQ=255)
                                                                                                                                                     aaGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTc  <  2:897834/150‑1 (MQ=255)
                                                                                                                                                     aaGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCAGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTc  >  1:456952/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGAAAAAGAGATGCGCTTTAGTCTGAAATAGTTGACTTAGTCCCTTATTGGCGATGTGGTTTTTGTTTTACCTGTCTGTCAGGTGGCAGCAAAAAGCAACTTTCCAGTTTTTACGCTGATTCAGATTTTAGCTATAAAAAAACCCGCCGAAGCGGGTTTTTTCGAAAATTGTTTTCTGCCGGAGCAGAAGCCAATTAACGTTTGGAGAACTGCGGACGACGACGTGCTTTACGCAGACCGACTTTCTTACGTTCAACCTGACGAGCGTCACGAGTAACGAAGCCAGCTTTACGCAGTTC  >  NC_000913/3377621‑3377919

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