Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,229,177 A→C 19.7% 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,1770AC19.7% 154.1 / 7.8 114intergenic (+112/+261)yohP/dusCuncharacterized protein/tRNA‑dihydrouridine synthase C
Reads supporting (aligned to +/- strand):  ref base A (26/50);  new base C (18/2);  total (54/60)
Fisher's exact test for biased strand distribution p-value = 7.24e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 116 total aligned reads displayed.

ACTGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATG  >  NC_000913/2229029‑2229326
                                                                                                                                                     |                                                                                                                                                      
aCTGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTACTGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAg                                                                                                                                                       <  1:3965047/150‑1 (MQ=255)
  tggtggTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAggg                                                                                                                                                     >  2:2241644/1‑150 (MQ=255)
    gtggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAggggg                                                                                                                                                   <  2:1900584/150‑1 (MQ=255)
     tggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGACTATTTCAGTCTCTCATACAGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGGGGGggc                                                                                                                                                  >  1:1691547/1‑150 (MQ=255)
     tggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGggc                                                                                                                                                  <  2:2669978/150‑1 (MQ=255)
     tggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGggc                                                                                                                                                  <  2:136935/150‑1 (MQ=255)
      ggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGggcg                                                                                                                                                 >  1:2989547/1‑150 (MQ=255)
       gtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGCGCGCGCGGGGggcgg                                                                                                                                                >  1:270951/1‑150 (MQ=255)
       gtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGCGCGCGCCGGGggcgg                                                                                                                                                >  1:1405535/1‑150 (MQ=255)
        tgACTGTCGTATTTTAGATGATATTTTAAAATTAATTAATTTAATGAGGTCCTTATCTTACGAGAATATTTCTGTCTGTCATCCTGTTGCGCTCGTGTCTTGTGCATTGCTTAATATAATCACTGGCTGAATGATCGCGCAggggtattc                                                                                                                                               <  2:4184098/150‑6 (MQ=255)
        tgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGggcggc                                                                                                                                               <  1:1414255/150‑1 (MQ=255)
        tgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGggcggc                                                                                                                                               <  2:4191081/150‑1 (MQ=255)
          aCTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCa                                                                                                                                             <  1:1046504/150‑1 (MQ=255)
           cTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGCggggggggggca                                                                                                                                            >  1:3280854/1‑146 (MQ=255)
           cTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCaa                                                                                                                                            <  1:3027595/150‑1 (MQ=255)
           cTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTAAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGCGGGGGGGGGccc                                                                                                                                             >  1:3280794/1‑149 (MQ=255)
            tGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGGCGGGGGGGcccca                                                                                                                                           >  2:368618/1‑147 (MQ=255)
            tGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGCgggggggggcgca                                                                                                                                           >  1:2747745/1‑145 (MQ=255)
            tGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGGGGCCaaa                                                                                                                                           >  1:3893777/1‑149 (MQ=255)
            tGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGcccag                                                                                                                                           >  2:2777426/1‑147 (MQ=255)
              gCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCCGGGGGCCGCCaaaac                                                                                                                                         >  2:2073024/1‑147 (MQ=255)
              gCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCCCGGGGCGGcccaacc                                                                                                                                         >  2:4021262/1‑145 (MQ=255)
              gCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAAtag                                                                                                                                         >  2:310242/1‑148 (MQ=255)
              gCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCACGGGgccgcccaaac                                                                                                                                         >  2:4021093/1‑145 (MQ=255)
                gTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATcgcc                                                                                                                                       <  2:445492/150‑1 (MQ=255)
                gTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATcgcc                                                                                                                                       <  1:3310360/150‑1 (MQ=255)
                  aTTTAAGATGATATTTTAAAATTAATTAATGTCATCGGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATcgccgc                                                                                                                                     <  1:122586/150‑1 (MQ=255)
                   tttAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCGATTggccgcc                                                                                                                                    <  1:2002086/150‑1 (MQ=255)
                      aaGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGCGGGGGCGGGCCAACgcccccccc                                                                                                                                 >  1:3786887/1‑150 (MQ=255)
                      aaGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCGGGGGGGGGGCAAAAgcccccgcc                                                                                                                                 >  2:570778/1‑150 (MQ=255)
                      aaGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGGCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGGGGGGGGGGGCCAGCgcccccccc                                                                                                                                 >  1:1107705/1‑150 (MQ=255)
                       aGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGccc                                                                                                                                <  1:3693046/150‑1 (MQ=255)
                        tatgatATTTTAAAATTAATTAATTTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCGTGTCATGTGCTTTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGcccc                                                                                                                               <  1:3270625/149‑1 (MQ=255)
                        gatgatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGcccc                                                                                                                               <  2:2701651/150‑1 (MQ=255)
                        gatgatATTTTAAAATTAATTAAGGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCCGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCTGGGGGCGGCCAATCGCCGCCGcccc                                                                                                                               <  1:4165814/150‑1 (MQ=255)
                            atatTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCCGGGGGCGGCCAAAACCCGCCGCCCCCAGc                                                                                                                           >  1:321795/1‑150 (MQ=38)
                            atatTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCCATCGCCGCCGCCCCCGGc                                                                                                                           >  1:321892/1‑150 (MQ=255)
                              tttttAAAATTAATTATGTTCATCAGGTCCTAAAATAACGAGAATATTTCAGTCTCTCTTCCTGTTGCGCTCCTTTCTTGTGCTTTGCTTCATATAATCACTGGCGCAATGAGCGCGCCGGGGGCGGCCAATCGCCGCCGCCCCCTGCAc                                                                                                                         <  1:1172848/149‑1 (MQ=25)
                               ttttAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGCGGGGGGGGGCCCAccccccccccccccccccc                                                                                                                         >  1:1939388/1‑145 (MQ=255)
                                tttAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCCGGGGGCGGCCCAACGCCCCCCCCCCCCCCAccc                                                                                                                       >  2:1563970/1‑150 (MQ=38)
                                  tAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGGGGGCCAAAGCCGCCGCCCCCCGCAccccc                                                                                                                     >  2:2549304/1‑150 (MQ=38)
                                  tAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCCAGACCCCCCCCCCCCCGCAccccc                                                                                                                     >  2:2549331/1‑150 (MQ=38)
                                  tAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCGTCATATACTCACTGGCGCAAGGAGCGCGCCGGGGGCGGCCAAAAGCCCCCCCCCCCCCCAccccc                                                                                                                     >  2:2549329/1‑150 (MQ=255)
                                    aaaTTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCCGGGGGGGGCCCAACGCCGCCGCCCCCCGCCCCCCCgg                                                                                                                   >  2:3794012/1‑150 (MQ=38)
                                    aaaTTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCCGGGGGGGGCCCAAAGCCCCCCCCCCCAGCAccccccc                                                                                                                   >  2:2648347/1‑148 (MQ=38)
                                     aattaattAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCggg                                                                                                                  <  2:4127298/150‑1 (MQ=255)
                                      attaattaATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGCGGGGGGGGCCCAACCCCCCCCCCCCCAccccccccccgc                                                                                                                 >  1:1962185/1‑146 (MQ=255)
                                      attaattaATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCGGGGGGCGGCCCAACGCCGCCCCCCCCCGCACCCCCCGGc                                                                                                                 >  1:2472256/1‑150 (MQ=255)
                                      attaataaATGTCATCAGGTCCGAAAATAACGAGAATAATACAGACTATCATCCAGTTGCGCTCATGACATGTGCATTGCTTCATATAATCACTGTCGCAAGGAGCGCGGAGGGCGCGGGCCAACACACCCCCCCCCCCCACCCCCcggg                                                                                                                 >  1:1258343/1‑149 (MQ=255)
                                           ttaatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGGGGGGGGGGGCCATACCCCCCGCCCCCCGCACCCCCCGGcgcggg                                                                                                            >  1:163922/1‑145 (MQ=255)
                                           ttaatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCCGGGGGGGGGCCAACGCCCCCCCCCCCCGCACCCCCCGGCTCggg                                                                                                            >  1:3838555/1‑150 (MQ=255)
                                            taatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCCGGGGGGGGCCCAACCCCCCCCCCCCCCCCCCCCCCCCGCACTGGc                                                                                                           >  2:2065111/1‑150 (MQ=37)
                                            taatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCCGGGGGCGGCCCAACGCCCCCCCCCCCCGCCCCCCCCGGCGCCGGc                                                                                                           >  2:2065043/1‑150 (MQ=37)
                                            taatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCCCGGGGCGGCCACACCCCCCCCCCCCCCGCCCCCCCCGGCGCCGGc                                                                                                           >  2:2065028/1‑150 (MQ=37)
                                             aatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAAAACCCCCCCCCCCCAGCAACCCCCGGCTCTGGCg                                                                                                          >  2:561410/1‑150 (MQ=37)
                                              atGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCGCGGGGGGGCCCAACGCCGCCCCCCCCACCCCCCCCCCGCCCTGGCGa                                                                                                        >  1:3473816/1‑150 (MQ=255)
                                               tGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATATTCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGaa                                                                                                       <  2:1940216/150‑1 (MQ=255)
                                               tGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGaa                                                                                                       <  2:1164053/150‑1 (MQ=255)
                                                  catcaGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCCAACGCCCCCCCCCCCACCACCCCCCGGCGCTGGCGACCaa                                                                                                    >  2:244501/1‑150 (MQ=37)
                                                      tgttCCTAAAATATCTATAATTTTTCAGTCTCTCATCCTGTTGCGCTCCTGTCTTGTGATTTGCTTAATATAATCCGGGGCGCAAGGAGGGCGCTGGGGGCGGCAAAGCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATc                                                                                                <  1:2211668/147‑1 (MQ=18)
                                                       ggTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATcg                                                                                               <  2:2288497/150‑1 (MQ=255)
                                                         tCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCCAGGGGCGGCCCAAACCCCCCCCCCCCCGCACCCCCCCGCCCAGGCGCAACaaaaaccc                                                                                              >  2:2194902/1‑146 (MQ=25)
                                                         tCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGGGGCCCAAAGCCGCCCCCCCCAGAACCCCCCGGCACGGGCGAACaaaaaagcc                                                                                             >  2:571923/1‑145 (MQ=25)
                                                                    cGATAATATTTAAGTCTCTCATCCTGTTGCGCTCCTGTCTTGTGCATTGCTTCATATAATCACTGGCGCATGGAGCGCGCTGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTg                                                                                  <  2:915749/150‑1 (MQ=18)
                                                                    cGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTg                                                                                  <  2:3036537/150‑1 (MQ=255)
                                                                     gagaATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCACGGGGCGGCCCAAAGCCCCCCCCCCCAGCCCCCCCCGGCAAAGGCCAAAAAAAACCCCCCAGCGCCGTGc                                                                                 >  2:2698271/1‑150 (MQ=21)
                                                                          tatTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCg                                                                            <  2:3538902/150‑1 (MQ=255)
                                                                          tatTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCg                                                                            <  2:2194365/150‑1 (MQ=255)
                                                                           atTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCgg                                                                           <  2:2989547/150‑1 (MQ=255)
                                                                            tttCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCACAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGc                                                                          <  2:1844113/150‑1 (MQ=255)
                                                                                  tctctcATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACGGGCGCACGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATcc                                                                    <  2:3930935/150‑1 (MQ=255)
                                                                                  tctctcATCCTGTTGCGCTCCTGTCATGTGCATTGAGTCATATAATCACTGGCGCCAGGAGCGCGCAGGGGGCGGCCATTCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATcc                                                                    <  2:1084219/150‑1 (MQ=255)
                                                                                      tcATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGGGGCCCAAACCCCCCGCCCCCCGCACCCCCCGGCGCGGGCGGACAAAAACCCCGCGGCGCGGTGGCCGCGGCGTATCCCCTa                                                                >  2:1335426/1‑150 (MQ=18)
                                                                                      tcATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGGCAAAAGCCGCCCCCCCCAGCACCCCCGGGCTCTGGCGAACAAAATAGCCCCGGCCCGGTGGCCGCGGCGTATcccaaa                                                                >  2:1347769/1‑147 (MQ=18)
                                                                                                     ccTGTCATGTGAATTGCTTCATATAAGCACTGGCGCATGGAGCGCGCTGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCgg                                                 <  1:368618/150‑1 (MQ=255)
                                                                                                           aTGTGCTTTGATTCATATAATCACTGGCGCACGGAGCGCGCCGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGa                                           <  1:36376/150‑1 (MQ=255)
                                                                                                              tgCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAAAAGCCGCCGCCCCCAGCACCCCCGGGCACTGGCGAACAAAATCGCCGCTGCGCGGTGCCCAAGGCTTATCCCATACCGCTAACCGGTCGGGCGCGAGGt                                        >  2:3192074/1‑150 (MQ=14)
                                                                                                                      ttCATATAATCACTGGCGCAAGGAGCGCGCTGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATccc                                <  1:244501/150‑1 (MQ=255)
                                                                                                                       tCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCCATCGCCGCCGCCCCCAGCACCCCCCGGCTCTGGCGAACAAAAACCCCGCTGCGCGGGGCCCGCGGCTTATCCCATACGGCACCCCGGTCGGGCGCGCGGTAACAAcccg                               >  2:2469697/1‑149 (MQ=14)
                                                                                                                       tCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCt                               <  1:786136/150‑1 (MQ=39)
                                                                                                                         atatAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGt                             <  1:4081494/150‑1 (MQ=38)
                                                                                                                         atatAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCCCCGCCCCCAGCACCCCCCGGCTCAGGCGAACAAAATCGCCGCGGCGCGGTGCCCACCGCGTATACCATACGGCTAACCGGTCCGGCGCGAGGTTACAACCCAGt                             >  2:4214936/1‑150 (MQ=11)
                                                                                                                            taATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTaaa                          <  1:255652/150‑1 (MQ=35)
                                                                                                                             aaTCACTGGCGCAAGGAGCGCGCTGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTaaaa                         <  1:2538518/150‑1 (MQ=21)
                                                                                                                             aaTCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTaaaa                         <  1:3193163/150‑1 (MQ=34)
                                                                                                                             aaTCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAAAAGCCCCCCCCCCCGGCACCCCCCGGCGCTGGCGAACAAAAACGCCCCGGCGCCGGGGCCACGGCGTATACCATACCGCTACCCGGTCGGGCGCCAGGTAAAAACCCAGTaaaa                         >  2:2311631/1‑150 (MQ=2)
                                                                                                                                 aCTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAcgcg                     <  2:3514735/150‑1 (MQ=32)
                                                                                                                                  cTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACgcgc                    <  2:3845465/150‑1 (MQ=32)
                                                                                                                                    ggCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGccc                  <  2:3085726/150‑1 (MQ=32)
                                                                                                                                        cAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAg              <  1:193942/150‑1 (MQ=32)
                                                                                                                                        cAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAg              <  1:292467/150‑1 (MQ=32)
                                                                                                                                        cAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAg              <  2:3955827/150‑1 (MQ=32)
                                                                                                                                           ggAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGccc           <  2:3074802/150‑1 (MQ=32)
                                                                                                                                                 cgcAGGGGGCGGCCACTCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATcc     <  2:413428/150‑1 (MQ=17)
                                                                                                                                                 cgcAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATcc     <  1:2819192/150‑1 (MQ=32)
                                                                                                                                                 cgcAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATcc     <  1:2073081/150‑1 (MQ=32)
                                                                                                                                                 cgcAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATcc     <  2:882577/150‑1 (MQ=32)
                                                                                                                                                    cAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAt   <  1:3869089/150‑1 (MQ=32)
                                                                                                                                                     aGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAtg  >  2:66527/1‑150 (MQ=32)
                                                                                                                                                     aGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAtg  >  2:290355/1‑150 (MQ=32)
                                                                                                                                                     |                                                                                                                                                      
ACTGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATG  >  NC_000913/2229029‑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: