Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,229,177 A→C 17.3% 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,1770AC17.3% 171.0 / 2.7 104intergenic (+112/+261)yohP/dusCuncharacterized protein/tRNA‑dihydrouridine synthase C
Reads supporting (aligned to +/- strand):  ref base A (22/50);  new base C (14/3);  total (39/65)
Fisher's exact test for biased strand distribution p-value = 1.71e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 106 total aligned reads displayed.

TACTGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCC  >  NC_000913/2229028‑2229323
                                                                                                                                                      |                                                                                                                                                   
tACTGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCa                                                                                                                                                     >  2:679910/1‑150 (MQ=255)
  cTGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGACTATTTCAGGCTCTCATCGTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCTCTGGCGCATGGCGCGC‑GCg‑gg                                                                                                                                                  >  1:795615/1‑150 (MQ=255)
  cTGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCC‑gg                                                                                                                                                  >  1:1678253/1‑150 (MQ=255)
  cTGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑gg                                                                                                                                                  <  2:2224477/150‑1 (MQ=255)
  cTGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑gg                                                                                                                                                  >  1:1834195/1‑150 (MQ=255)
   tggtggTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑ggg                                                                                                                                                 <  1:3197651/150‑1 (MQ=255)
   tggtggTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCAGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTTCATTGCTTCATATAATCACTGGCGCAATGAGCGC‑GCA‑ggg                                                                                                                                                 <  1:3008341/150‑1 (MQ=255)
   tggtggTGACTGGAGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑ggg                                                                                                                                                 <  1:3008314/150‑1 (MQ=255)
    ggtggtGACTGTCGTATTTAAGATTATATTTTAATATTAATTAATGTCATCAGGTCCGAAAATAAAGAGAATATTTCTGTCTCTCATCCTTTTGCGCTCCTGTCATGTGCATTGCTTAATATAATCACTGGCGCAAGGAGCGC‑GCA‑gggg                                                                                                                                                <  1:2079931/150‑1 (MQ=255)
     gtggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑ggggg                                                                                                                                               >  2:1281629/1‑150 (MQ=255)
      tggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GGC‑gggggg                                                                                                                                              >  1:36455/1‑149 (MQ=255)
      tggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCC‑gggggg                                                                                                                                              >  2:2257573/1‑149 (MQ=255)
      tggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCAGCCTGTTGCGCTCCTGTCATGTGCCTTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GGC‑gggggg                                                                                                                                              >  1:313216/1‑149 (MQ=255)
         tgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GGC‑GGGgggggc                                                                                                                                           >  2:907760/1‑150 (MQ=255)
         tgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCC‑GGGggcggc                                                                                                                                           >  2:1091436/1‑150 (MQ=255)
           aCTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCC‑GGGGGGGGccc                                                                                                                                         >  2:708930/1‑149 (MQ=255)
             tGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGGTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGGGGCGGCCAAt                                                                                                                                       >  2:2929369/1‑150 (MQ=255)
             tGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCC‑GGGGGCGGcccaa                                                                                                                                       >  2:2886203/1‑147 (MQ=255)
             tGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGGGGCGGCCaaa                                                                                                                                       >  2:187205/1‑149 (MQ=255)
             tGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑CGGGGGGGCCaaa                                                                                                                                       >  2:2886219/1‑149 (MQ=255)
             tGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGAGC‑GCC‑GGGGGCGGcccaa                                                                                                                                       >  2:2886211/1‑147 (MQ=255)
               gCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGGGGGGGCCAAtag                                                                                                                                     >  1:510479/1‑148 (MQ=255)
                cGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGGGGCGGcccaaccc                                                                                                                                    >  2:2097242/1‑144 (MQ=255)
                cGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGGGGCGGCCAAAcgc                                                                                                                                    >  2:2938550/1‑150 (MQ=255)
                 gTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGGGGGCGCCAAAcgcc                                                                                                                                   >  2:810651/1‑150 (MQ=255)
                       aaGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGGGGGGGCCAATCgccgccgcc                                                                                                                             >  2:2382552/1‑150 (MQ=255)
                       aaGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGGGGCGGCCAATCgccgccgcc                                                                                                                             <  1:2617440/150‑1 (MQ=255)
                         gatgatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GGA‑GGGGGCGGGCAAACGCCGCCGcccc                                                                                                                           >  2:2890375/1‑150 (MQ=255)
                         gatgatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGGGGCGGCCAATCGCCGCCGcccc                                                                                                                           <  2:1899816/150‑1 (MQ=255)
                           tgatATGTTCAAATTCATTCATGTCATCAGGTCCGAAAATAACGAGAAGAGGTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGGGGCGGCCAATCGCCGCCGCCCCCt                                                                                                                         <  2:297601/150‑1 (MQ=37)
                            gatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GGAGGGGGGGGGCCAAAccccccccccccag                                                                                                                         >  2:578980/1‑148 (MQ=255)
                            gatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGGGGCCGCCAAAAGCccccccccccag                                                                                                                        >  2:460948/1‑148 (MQ=39)
                                ttttAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCC‑GGGGGGGGGCCAACCCCCccgccccccgcccc                                                                                                                    >  1:1833785/1‑150 (MQ=255)
                                   tAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAATGAGCGC‑GCA‑GGGGGCGGCCAATCGCCGCCGCCCCCTGCAccccc                                                                                                                 <  2:2945228/150‑1 (MQ=255)
                                     aaaTTAATTATTGTCATCTGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCGGTCATGTGATTTGCTTCATATAATCACTGGCGCATGGAGCGC‑GCA‑GGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCgg                                                                                                               <  2:655750/150‑1 (MQ=37)
                                     aaaTTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGCG‑GGGGGGGGCCAA‑CGCCCCCCCCCCCCCCAccccccc                                                                                                               >  1:1404991/1‑148 (MQ=38)
                                     aaaTTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GGG‑GGGGGGGGGCAAAACCCCCCCCCCCCAGCcccccccc                                                                                                               >  1:2828995/1‑148 (MQ=255)
                                         taattaatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GGA‑GGGGGCGGGCCAAAGccccccccccccccccccccccgcgc                                                                                                           >  1:1953334/1‑144 (MQ=255)
                                         taattaatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCC‑GGGGGGGGGCAAAAGCCCCCCCCCCCCCCAcccccccgcgc                                                                                                           >  1:1084954/1‑144 (MQ=37)
                                          aattaatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGctct                                                                                                          <  2:2749360/150‑1 (MQ=255)
                                                tGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGaa                                                                                                    <  2:1795902/150‑1 (MQ=255)
                                                tGTCATCAGGTCCGAAAATAACGAGAATATTTAAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGAATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGaa                                                                                                    <  2:111927/150‑1 (MQ=255)
                                                ggTCAGCAGGTCCGACACTCACGAGAATCTTTCAGGCGCTCATCCTGGTGCGCTCCTGTCAGGTGCATTGCTTCATATACTCACGGGCGCAAGGCGCGC‑GCA‑GGGGGCGGCCCATCGCCGCCGCCCCCTGCACCCCCGGGCTCGGGCGaa                                                                                                    <  2:6104/149‑1 (MQ=21)
                                                  tcatcaGGTCCGAAAATAACGAGAATATTTCAGTGTCTGATCCTGTTGCGGTCCTGTCATGTGCATTGCTTCATATATTCACTGGCGCAAGGAGCGC‑GGG‑GGGGGGGGGGAATTGGCGGCGGCCCCGGGACCCCCGGGGTGTgggggacg                                                                                                  >  2:2901954/1‑144 (MQ=255)
                                                       aGGTCCGAAAATAACGAGAATATTTCAGTCGCTCAGCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCC‑GGGGGCGGCCCAACCCCCCCCCCCCCCGCCCCCCCCCGCCCCGGCGCACcaaaac                                                                                             >  2:2594189/1‑148 (MQ=255)
                                                         gTCCTAAAATAACGTTAATTTTTCTGTCTCTCATCTTGTTGCGATCCTGTCTTGTGCTTTGATTCATATAATCAGTGGCGCACGGCGCGC‑GCT‑GGGGGCGGCAAAGCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATcgc                                                                                           <  1:28477/150‑1 (MQ=18)
                                                         gTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCC‑GGGGGCCGCCCAACCCCGCCCCCCCCCCCCCCCCCCCGCCCGGGCGCaaaaaaaagc                                                                                           >  2:2325276/1‑146 (MQ=25)
                                                         gTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGGGGCGGCCATTCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGATCAAAATcgc                                                                                           <  2:2917496/150‑1 (MQ=255)
                                                         gTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATcgc                                                                                           <  2:729505/150‑1 (MQ=255)
                                                          tCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATcgcc                                                                                          <  2:2292781/150‑1 (MQ=255)
                                                                  tAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGcgcg                                                                                  <  2:2119779/150‑1 (MQ=255)
                                                                       agaCTATTTCAGTCTCTCAGCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑CGGGGCGGCCCATCGCCCCCCCCCCCGGCACCCCCCGGCTCAGGCGCCCCAAAACGCCGCCGCGCGGTGcc                                                                             >  2:535021/1‑150 (MQ=21)
                                                                            ttttCCGTCTCTAATCCTGTGGCGCTCCTGTCTGGTGATTTGCTTCTTATAATCACGGGCGCACGGCGCGC‑GCG‑GGGGGGGGCCAACCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCgg                                                                        <  1:3226726/149‑1 (MQ=17)
                                                                            ttttCCGTCTCAAATCCTGTGGCGCCCCGGTCTTGTGAATTGCTTAATATAATCACTGGCGCATGGCGCGC‑GCT‑GGGGGCGGCCATGCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCgg                                                                        <  1:3226755/149‑1 (MQ=17)
                                                                              ttCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGGGGCGGCCAAAACCCCCCCCCCCCACCACCCCCCCGCACAGGCGAAAAAAAAACACCCAGCGCGGTGCCCACGGca                                                                      >  2:488180/1‑149 (MQ=16)
                                                                                 aGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑CGGGGCGGCCAAACCCCCCCCCCCCCACCACCCCCCGGCACAGGCGAACAAAAAAGCCCCGGCCCGGTGCCCTCGGCGTAt                                                                   >  2:1842516/1‑150 (MQ=18)
                                                                                   tctctcATCCTGTTGGGCTCCTGTCTTGTGATTTGCTTCATATAATCCCTGGCGCAAGGTGCGC‑GCT‑GGGGGCGGCCACTCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATcc                                                                 <  1:161041/150‑1 (MQ=17)
                                                                                      ctcATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCtt                                                              <  2:2037512/150‑1 (MQ=255)
                                                                                       tcATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTa                                                             <  2:811431/150‑1 (MQ=255)
                                                                                         aTCCTGTTGCGCTCTTGTCTTGTGATTTGCTTAATATAATCACGGGCGCAGGGCGCGC‑GCT‑GGGGGGGCCAAATCCCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACg                                                           <  1:1138599/150‑1 (MQ=17)
                                                                                             tGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGGGGCGGCCCAACGCCGCCCCCCCCAGCACCCCCCGGCGCAGGCGAACAAAAAAGCCCCGGCGCGGTGCCCACCGCGTAACCCGTAAGGcaa                                                       >  2:2409722/1‑148 (MQ=17)
                                                                                                 gggcTCGTGTCTTGTGCATTGCTTCATTTAATCACTGGCGCATGGAGCGC‑GCG‑GGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCgg                                                   <  1:1528780/148‑1 (MQ=17)
                                                                                                 gcgcTCCTGTCTTGTGCTTTGATTCATATAATCACTGGCGCATGGCGCGC‑GCT‑GGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCgg                                                   <  1:2025497/150‑1 (MQ=17)
                                                                                                 gcgcTCCTGTCATGTGCATTGCTTCATATAATCCCTGGCGCAAGGAGCGC‑GCA‑GGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCgg                                                   <  2:2575553/150‑1 (MQ=255)
                                                                                                 gcgcTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCgg                                                   <  2:2575381/150‑1 (MQ=255)
                                                                                                   gcTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGGGGCGGCCAAACGTCGCCGCCCCCAGCCCCCCCGGGCGCTGGCCAACAAAATCGCCGCGGCGCGGTGCCCGCGGCATATCCCCTACCGCTACCCGGt                                                 >  2:2853397/1‑150 (MQ=17)
                                                                                                    cTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTc                                                <  2:1914242/150‑1 (MQ=255)
                                                                                                       cTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGGGGGGGCCAAAAGCCGCCGCCCCCAGCACCCCCCGGCGCAGGCGAACAAAAACCCCCCGGCGCGGTGCCCACCGCTTATCCCATAACGCAACCCCGTCggg                                             >  2:3197422/1‑150 (MQ=17)
                                                                                                       cTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCggg                                             <  2:1461583/150‑1 (MQ=255)
                                                                                                        tGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGGGGCGGCCAATAGCCGCCGCCCCCTGCTCCCCCGGGCTCTGGCGAACAAAATCGCCGGTGCGCGGTGCGCTCGGCTTATCCCTTTCGGCTACCGGGTCgggg                                            >  2:2421618/1‑149 (MQ=17)
                                                                                                           cATGTGCATTGCTTCATATAATCACTGGCGCACGGAGCGC‑GCA‑GGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGcgcg                                         <  2:1458061/150‑1 (MQ=255)
                                                                                                            aTGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGGGGGGGCCAATAGCCGGCGCCCCCAGGACCCCCGGGGTCTTTCGAACAAAATCGCCGCTTCGCGGTTCCCTCGGCTTATACCTTAAGGCTACCGGGTCggggggga                                        >  2:1873881/1‑145 (MQ=17)
                                                                                                             tgtgCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAg                                       <  2:1546650/150‑1 (MQ=255)
                                                                                                             tgtgCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAg                                       <  2:315740/150‑1 (MQ=255)
                                                                                                               tgctttgcttCATATAATCACTGGCGCAAGGAGCGC‑GCC‑GGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGt                                     <  2:1414303/150‑1 (MQ=255)
                                                                                                                gCATTGCTTCATATAATCAGTGGCGCAAGGAGCGC‑GCA‑GGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTa                                    <  2:2105750/150‑1 (MQ=255)
                                                                                                                gCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTa                                    <  2:2034768/150‑1 (MQ=255)
                                                                                                                 cATTGCTTCATATAACCCCTGGCGCAAGGAGCGC‑GCC‑GGGGGCGGCCAACCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTaa                                   <  1:1200639/150‑1 (MQ=255)
                                                                                                                    tGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCT‑GGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACAt                                <  1:708930/150‑1 (MQ=255)
                                                                                                                      ctacATATAATCACTGGCGCATGGAGCGC‑GCA‑GGGGGCGGCCACTCGCCGCCGCGCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTATCATcc                              <  1:1932036/147‑1 (MQ=11)
                                                                                                                      ctacATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGGGGCGGCCAAACGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATcc                              <  1:2901954/147‑1 (MQ=255)
                                                                                                                      cTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATcc                              <  2:1625754/150‑1 (MQ=40)
                                                                                                                      cTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATcc                              <  2:1272712/150‑1 (MQ=40)
                                                                                                                        tCATATACTCACTGGCGCAAGGAGCGC‑GCA‑GGGGGCGGCCACTCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCt                            <  2:2828875/150‑1 (MQ=255)
                                                                                                                        tCATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCt                            <  1:506418/150‑1 (MQ=39)
                                                                                                                        tCATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCt                            <  1:2169783/150‑1 (MQ=39)
                                                                                                                        tCATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCt                            <  1:1389020/150‑1 (MQ=39)
                                                                                                                        tCATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCt                            <  2:2801809/150‑1 (MQ=39)
                                                                                                                              aaTCACTGGCGCAAGGAGCGC‑GCA‑GGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTaaaa                      <  1:2284784/150‑1 (MQ=34)
                                                                                                                               aTCACTGGCGCACGGAGCGC‑GCC‑GGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAc                     <  2:1456369/150‑1 (MQ=18)
                                                                                                                               aTCACTGGCGCAAGGAGCGC‑GCT‑GGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATTCCTGTAAAAc                     <  1:1091436/150‑1 (MQ=18)
                                                                                                                               aTCACTGGCGCAAGGAGCGC‑GCT‑GGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATTCCTGTAAAAc                     <  1:907760/150‑1 (MQ=18)
                                                                                                                                 cACTGGCGCAAGGAGCGC‑GCA‑GGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAcgc                   <  1:454812/150‑1 (MQ=32)
                                                                                                                                   cTGGCGCAAGGAGCGC‑GCA‑GGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACgcgc                 <  2:647041/150‑1 (MQ=32)
                                                                                                                                          aaGGAGCGC‑GCA‑GGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGc          <  1:126968/150‑1 (MQ=32)
                                                                                                                                          aaGGAGCGC‑GCA‑GGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGc          <  2:1794874/150‑1 (MQ=32)
                                                                                                                                            ggAGCGC‑GCA‑GGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGccc        <  2:434177/150‑1 (MQ=32)
                                                                                                                                             gAGCGC‑GCT‑GGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCca       <  1:3045943/150‑1 (MQ=17)
                                                                                                                                             gAGCGC‑GCG‑GGGTGCGGCCTATCGCCGCCGCCCCCTGCACCCCCGGGCTATGGCGAACTAAATAGCCGCTGCGCGGTGCCGTCGGCTTATCCCTTACGGCTACGGGGTCGGGCGCGAGGTAACTTCCCTGTAATACGAGCGCTCAGccga       <  1:2473212/150‑3 (MQ=11)
                                                                                                                                                  c‑gcA‑GGGGGCGGCAAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACAtgg  <  2:1404991/150‑3 (MQ=17)
                                                                                                                                                      |                                                                                                                                                   
TACTGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCC  >  NC_000913/2229028‑2229323

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