Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,229,177 A→C 18.9% intergenic (+112/+261) yohP → / ← dusC uncharacterized protein/tRNA‑dihydrouridine synthase C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,229,1770AC18.9% 132.8 / 7.6 79intergenic (+112/+261)yohP/dusCuncharacterized protein/tRNA‑dihydrouridine synthase C
Reads supporting (aligned to +/- strand):  ref base A (14/36);  new base C (10/2);  total (32/47)
Fisher's exact test for biased strand distribution p-value = 6.83e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TACTGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATG  >  NC_000913/2229028‑2229326
                                                                                                                                                     |                                                                                                                                                      
tACTGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCa                                                                                                                                                        <  1:2776663/150‑1 (MQ=255)
  cTGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAgg                                                                                                                                                      <  2:2084597/150‑1 (MQ=255)
  cTGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAgg                                                                                                                                                      <  2:2084642/150‑1 (MQ=255)
   tggtggTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCCggg                                                                                                                                                     >  2:1610024/1‑150 (MQ=255)
     gtggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCgggggggg                                                                                                                                                   >  2:792868/1‑150 (MQ=255)
     gtggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAggggg                                                                                                                                                   <  1:1611545/150‑1 (MQ=255)
       ggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGggc‑g                                                                                                                                                <  2:2133046/150‑1 (MQ=255)
       ggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGggc‑g                                                                                                                                                <  2:360960/150‑1 (MQ=255)
         tgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCCCGGggc‑ggc                                                                                                                                              >  2:1500748/1‑150 (MQ=255)
           aCTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGC‑GGCCa                                                                                                                                            >  2:1729733/1‑150 (MQ=255)
             tGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGCCGGGGC‑GGcccaa                                                                                                                                          >  2:589238/1‑147 (MQ=255)
             tGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGC‑GGCCaaa                                                                                                                                          >  2:2848831/1‑149 (MQ=255)
             tGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGC‑GGCCaaa                                                                                                                                          >  2:708951/1‑149 (MQ=255)
             tGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGC‑GGCCCAt                                                                                                                                          >  2:1572518/1‑150 (MQ=255)
             tGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGC‑GGCCAAt                                                                                                                                          <  1:1720203/150‑1 (MQ=255)
             tGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGC‑GGCCAAt                                                                                                                                          <  1:1788695/150‑1 (MQ=255)
             tGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGC‑GGCCAAt                                                                                                                                          <  1:1788483/150‑1 (MQ=255)
              ggCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGC‑GGCCAATc                                                                                                                                         <  1:2554349/150‑1 (MQ=255)
               gCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGC‑GGCCAATcg                                                                                                                                        >  1:2541946/1‑150 (MQ=255)
                cGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGC‑GGCCAATcgc                                                                                                                                       <  1:1939167/150‑1 (MQ=255)
                 gTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCCCCGGGC‑GGccccacccc                                                                                                                                      >  2:2059885/1‑143 (MQ=255)
                   aTTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGC‑GGCCAATcgccgc                                                                                                                                    <  2:2179308/150‑1 (MQ=255)
                     ttAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGC‑GGCCAATcgccgccg                                                                                                                                  <  1:1437093/150‑1 (MQ=255)
                     ttAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGC‑GGCCAATcgccgccg                                                                                                                                  <  1:883219/150‑1 (MQ=255)
                           tgatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGC‑GGCCAATCGCCGCCGCCCCCt                                                                                                                            <  1:681353/150‑1 (MQ=255)
                           tgatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGC‑GGCCAATCGCCGCCGCCCCCt                                                                                                                            <  1:677362/150‑1 (MQ=255)
                            gatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGCCGGGGC‑GGCCCAACcccccccccccccc                                                                                                                           >  2:861260/1‑148 (MQ=255)
                                ttttAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCAcc                                                                                                                       <  1:424467/150‑1 (MQ=255)
                                 tttAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGGGGGGGC‑GGCCCCACGCCCCCCCCCCCCGCcccc                                                                                                                      >  1:2760723/1‑150 (MQ=38)
                                 tttAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCGGGGGGC‑GGGCCCGCGCCCCCCCCCCCCCCAccc                                                                                                                      >  1:604896/1‑150 (MQ=38)
                                 tttAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGC‑GGCCCAACGCCGCCGCCCCCGGCAccc                                                                                                                      >  2:2187754/1‑150 (MQ=255)
                                      aattaattAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCGGGGGGG‑GGCCCAACGCCGCCGCCCCCCCCAcccccccg                                                                                                                 >  1:51431/1‑147 (MQ=38)
                                       attaattaATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGGGGGGGG‑GGGCCCACGCCGCccccccccccccccccccgc                                                                                                                >  2:63412/1‑146 (MQ=255)
                                       attaattaATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCCGGGGGCGGGCCAAAC‑CCCCCCCCCCCCGCACCCCCCGGc                                                                                                                >  2:2058694/1‑150 (MQ=255)
                                         taattaatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGTGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGGGGGGGG‑GGGCCAGCGGCGCCGCCCCCCCCACCCCCcgggcc                                                                                                              >  2:2750800/1‑147 (MQ=255)
                                         taattaatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGGGGGGGG‑GGGCCCGCGCCCCCCCCCCCCCAAccccccccccc                                                                                                              >  1:1300759/1‑144 (MQ=255)
                                         taattaatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGGGGGGGCGGGCCAAAC‑CCCCCCCCCCCCCCAccccccccccc                                                                                                              >  1:545681/1‑144 (MQ=255)
                                         taattaatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCCCGGGGC‑GGCCCAACGCCGCCCCCCCCAGCCCCCCCCGGctc                                                                                                              >  2:918009/1‑150 (MQ=37)
                                         taattaatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGC‑GGCCCAACGCCCCCCCCCCCCCCAACCCCCGGcgc                                                                                                              >  2:917960/1‑148 (MQ=37)
                                             taatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAACGCGCAGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGc                                                                                                          <  2:2599072/150‑1 (MQ=255)
                                                    atcaGGTCCTAAAATAACTATAATATTTCTGTCTCTCATCCTGTTGCGCTCCTGTCTTGTGCTTTGCTTCATATAATCCCGGGCGAATGGAGGGGGCTGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACaaa                                                                                                   <  1:1686154/150‑1 (MQ=21)
                                                    atcaGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCCGGGGGG‑GGGCCATCGCCCCCCCCCCCCGCACCCCCCGGCGCAGGCGGACcaa                                                                                                   >  1:2899495/1‑150 (MQ=255)
                                                     tcaGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCGCGGGGC‑GGCCAAACGCCCCCCCCCCCCGCCCCCCCCGGCGGTGGCGGACacaa                                                                                                  >  2:2521178/1‑150 (MQ=25)
                                                         gTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGG‑GGCCAATCGCCCCCGCCCCCACCACCCCCCGGCGCTGGCGAACaaaaaaga                                                                                              >  1:1205470/1‑146 (MQ=25)
                                                                 aTAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTgcgc                                                                                      <  2:2609070/150‑1 (MQ=255)
                                                                    actttaaTTTTTCTGTCTCTCTTGCTGGTGCGCTCGTGTCTTGTGCTTTGCTTAATATAATCACGGGCGAATGGCGCGCGCTGGGGGC‑GGCAAACGGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGt                                                                                   <  1:2403056/145‑1 (MQ=17)
                                                                         aaTATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCACAGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCt                                                                              >  1:2749699/1‑150 (MQ=255)
                                                                            atgtctgTCTCTCATCCTGTTGCGCTCCGGTCATGGGCTTTGATTCATATAAGCCCTGGCGCAAGGAGGGCGCCGGGGGG‑GGCCAATCGCGGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCgg                                                                           <  2:1110585/144‑1 (MQ=17)
                                                                             tttCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGCGCGCGCAGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGc                                                                          <  1:1937702/150‑1 (MQ=255)
                                                                                cAGTCTCTCAGCTTGTTGCGCTCCTGTCTTGTGATTTGCTTATTATAACCAGGGGCGCATGGCGCGCGCGGGGGGC‑GGCCACCCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAAAAAAATCGCCGCTGCGCGGTGCCCTCGGCTTa                                                                       <  1:1292113/150‑1 (MQ=17)
                                                                                   tctctcATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATcc                                                                    <  2:1625162/150‑1 (MQ=255)
                                                                                   tctctcAGCCCGTTGCGCTCCGGTCATGCGCAGTGCTACATATACTCACTGGCGCAAGGAGCGCGCAGGGGGC‑GGCCACTCGCCGCCGCCCCCTGCACCCCCGGGCTCGGGCGCACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATcc                                                                    <  2:2227104/150‑1 (MQ=17)
                                                                                   cctctcAGCCTGGAGCGCTCCTGTCATGGGCAGTGCTGCATATACCCACTGGCGCAAGGAGCGCGCAGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATcc                                                                    <  2:2227088/149‑1 (MQ=14)
                                                                                       tcttcttGTTGCGCTCCGGTCTTGTGCATTGCTTAATATAATCAGGGGCGCATGGAGCGCGCAGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTa                                                                <  1:2556349/144‑1 (MQ=17)
                                                                                       tcATCCTGTTGCGCTCCTGTCATGTGCTTTGCTTCATTTAATCACGGGCGCATGGAGCGCGCCGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTa                                                                <  1:454285/150‑1 (MQ=17)
                                                                                         aTCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACg                                                              <  2:1894459/150‑1 (MQ=255)
                                                                                         aTCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACg                                                              <  2:2541946/150‑1 (MQ=255)
                                                                                                 gcgcTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCgg                                                      <  2:2301472/150‑1 (MQ=255)
                                                                                                    cTCCTGTCTTGTGCTTTGATTCATATAATCACTGGCGCAAGGAGCGCGCGGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTc                                                   <  1:1336287/150‑1 (MQ=255)
                                                                                                      ccTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCgg                                                 <  2:1435437/150‑1 (MQ=255)
                                                                                                      ccTGTCATGTGCATTGCTTCATATAATCACTGGCGAATGGAGCGCGCTGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCgg                                                 <  1:63412/150‑1 (MQ=255)
                                                                                                           cttgtgCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGcgcg                                            <  1:162342/148‑1 (MQ=255)
                                                                                                              gtgCATTGCTTCATATAATCACTGGCGCATGGAGCGCGCTGGGGGG‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAgg                                         <  1:302524/150‑1 (MQ=255)
                                                                                                                   ttGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGC‑GGCCAATAGCCCCCCCCCCCAGCACCCCCCGGCTCTGGCGCAAAAAAACCCCGCGGCGCGGGGCCCGCGGCGTATACCATACGGCTACCCGGGCGGGCGCCCGGGAACa                                    >  2:1872621/1‑150 (MQ=9)
                                                                                                                       ttCATATAATCACTGGCGCAAGGAGCGCGCTGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATccc                                <  1:589238/150‑1 (MQ=255)
                                                                                                                       ttCATATAATCACTGGCGCAAGGAGCGCGCCGGGGGC‑GGCCCAAACCCCCCCCCCCCCGCAACCCCCGGCTCAGGCGAACAAAAAAGCCCCGGCCCGGTGCCCGCCGCGTAACCCATACGGCTACCCGGGCGGGCGCGCGGTAAAAAccc                                >  2:2622312/1‑150 (MQ=255)
                                                                                                                       ttCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGC‑GGCCCAACGCCGCCCCCCCCAGCCCCCCCCGGCTGGGGCGGACAAAAACGCCCCAGCGCCGTGCCCCACGCGTATCCCATAACGCTACCGGGTCGGGCGCGCGGGAACATccc                                >  2:2403272/1‑150 (MQ=9)
                                                                                                                       ttCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGC‑GGCCAAAACCCCCCCCCCCCAGCCCCCCCCGGCGCAGGCGAACAAAAAACCCCCGGCGCGGGGCCCGCGGCTTAAACCCTACCGCTACCCGGTCGGGCGCGCGGTAACaaacc                                >  2:2403198/1‑146 (MQ=9)
                                                                                                                        tCATATACGCCCTGGCGCAAGGTGCGCGCTGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCt                               <  1:2061327/150‑1 (MQ=11)
                                                                                                                               aTCACTGGCGCAAGGAGCGCGCAGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAc                        <  1:2059885/150‑1 (MQ=33)
                                                                                                                                  aCTGGCGCAAGGAGCGCGCCGGGGGG‑GGCCCAACGCCGCCGCCCCCGGCACCCCCCGGCGCGGGCGCACAAAAACGCCCCCGCGCGGGGCCCGCCGCGGATCCCAGAAGGCTACCGGGGCGGGCGCGGGGGGACAAACCAGTAAAAcccg                     >  2:1859066/1‑150 (MQ=2)
                                                                                                                                   cTGGCGCAAGGAGCGCGCAGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACgcgc                    <  2:1452226/150‑1 (MQ=32)
                                                                                                                                    tGGCGCAAGGAGCGCGCTGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGcc                   <  1:2032327/150‑1 (MQ=17)
                                                                                                                                    tGGCGCAAGGAGCGCGCAGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGcc                   <  1:1610024/150‑1 (MQ=32)
                                                                                                                                     ggCGCAAGGAGCGCGCAGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGccc                  <  2:694279/150‑1 (MQ=32)
                                                                                                                                         cacggAGCGCGCAGGGGGC‑GGCCACTCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAg              <  1:2622237/147‑1 (MQ=17)
                                                                                                                                          aaGGAGCGCGCAGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGc             <  1:2797178/150‑1 (MQ=32)
                                                                                                                                                 gcgcAGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATc      <  1:1492770/150‑1 (MQ=32)
                                                                                                                                                  cgcAGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATcc     >  2:2836984/1‑150 (MQ=32)
                                                                                                                                                  cgcAGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATcc     >  2:1658614/1‑150 (MQ=32)
                                                                                                                                                   gcAGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCa    <  2:2281034/150‑1 (MQ=32)
                                                                                                                                                     aGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAtg  >  2:2911081/1‑150 (MQ=32)
                                                                                                                                                     |                                                                                                                                                      
TACTGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATG  >  NC_000913/2229028‑2229326

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: ‑