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.

TACTGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCC  >  NC_000913/2229028‑2229323
                                                                                                                                                      |                                                                                                                                                   
tACTGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCa                                                                                                                                                     >  2:679910/1‑150 (MQ=255)
  cTGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGACTATTTCAGGCTCTCATCGTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCTCTGGCGCATGGCGCGCGCggg                                                                                                                                                  >  1:795615/1‑150 (MQ=255)
  cTGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCCgg                                                                                                                                                  >  1:1678253/1‑150 (MQ=255)
  cTGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAgg                                                                                                                                                  <  2:2224477/150‑1 (MQ=255)
  cTGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAgg                                                                                                                                                  >  1:1834195/1‑150 (MQ=255)
   tggtggTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAggg                                                                                                                                                 <  1:3197651/150‑1 (MQ=255)
   tggtggTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCAGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTTCATTGCTTCATATAATCACTGGCGCAATGAGCGCGCAggg                                                                                                                                                 <  1:3008341/150‑1 (MQ=255)
   tggtggTGACTGGAGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAggg                                                                                                                                                 <  1:3008314/150‑1 (MQ=255)
    ggtggtGACTGTCGTATTTAAGATTATATTTTAATATTAATTAATGTCATCAGGTCCGAAAATAAAGAGAATATTTCTGTCTCTCATCCTTTTGCGCTCCTGTCATGTGCATTGCTTAATATAATCACTGGCGCAAGGAGCGCGCAgggg                                                                                                                                                <  1:2079931/150‑1 (MQ=255)
     gtggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAggggg                                                                                                                                               >  2:1281629/1‑150 (MQ=255)
      tggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGCgggggg                                                                                                                                              >  1:36455/1‑149 (MQ=255)
      tggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCCgggggg                                                                                                                                              >  2:2257573/1‑149 (MQ=255)
      tggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCAGCCTGTTGCGCTCCTGTCATGTGCCTTGCTTCATATAATCACTGGCGCAAGGAGCGCGGCgggggg                                                                                                                                              >  1:313216/1‑149 (MQ=255)
         tgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGCGGGgggggc                                                                                                                                           >  2:907760/1‑150 (MQ=255)
         tgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCCGGGggcggc                                                                                                                                           >  2:1091436/1‑150 (MQ=255)
           aCTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCCGGGGGGGGccc                                                                                                                                         >  2:708930/1‑149 (MQ=255)
             tGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGGTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAAt                                                                                                                                       >  2:2929369/1‑150 (MQ=255)
             tGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCCGGGGGCGGcccaa                                                                                                                                       >  2:2886203/1‑147 (MQ=255)
             tGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCaaa                                                                                                                                       >  2:187205/1‑149 (MQ=255)
             tGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCACGGGGGGGCCaaa                                                                                                                                       >  2:2886219/1‑149 (MQ=255)
             tGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGAGCGCCGGGGGCGGcccaa                                                                                                                                       >  2:2886211/1‑147 (MQ=255)
               gCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGGGGCCAAtag                                                                                                                                     >  1:510479/1‑148 (MQ=255)
                cGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGcccaaccc                                                                                                                                    >  2:2097242/1‑144 (MQ=255)
                cGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAAAcgc                                                                                                                                    >  2:2938550/1‑150 (MQ=255)
                 gTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGGCGCCAAAcgcc                                                                                                                                   >  2:810651/1‑150 (MQ=255)
                       aaGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGGGGCCAATCgccgccgcc                                                                                                                             >  2:2382552/1‑150 (MQ=255)
                       aaGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCgccgccgcc                                                                                                                             <  1:2617440/150‑1 (MQ=255)
                         gatgatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGAGGGGGCGGGCAAACGCCGCCGcccc                                                                                                                           >  2:2890375/1‑150 (MQ=255)
                         gatgatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGcccc                                                                                                                           <  2:1899816/150‑1 (MQ=255)
                           tgatATGTTCAAATTCATTCATGTCATCAGGTCCGAAAATAACGAGAAGAGGTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCt                                                                                                                         <  2:297601/150‑1 (MQ=37)
                            gatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGAGGGGGGGGGCCAAAccccccccccccag                                                                                                                         >  2:578980/1‑148 (MQ=255)
                            gatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCCGCCAAAAGCccccccccccag                                                                                                                        >  2:460948/1‑148 (MQ=39)
                                ttttAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCCGGGGGGGGGCCAACCCCCccgccccccgcccc                                                                                                                    >  1:1833785/1‑150 (MQ=255)
                                   tAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAATGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCAccccc                                                                                                                 <  2:2945228/150‑1 (MQ=255)
                                     aaaTTAATTATTGTCATCTGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCGGTCATGTGATTTGCTTCATATAATCACTGGCGCATGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCgg                                                                                                               <  2:655750/150‑1 (MQ=37)
                                     aaaTTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGCGGGGGGGGGCCAACGCCCCCCCCCCCCCCAccccccc                                                                                                               >  1:1404991/1‑148 (MQ=38)
                                     aaaTTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGGGGGGGGGGGCAAAACCCCCCCCCCCCAGCcccccccc                                                                                                               >  1:2828995/1‑148 (MQ=255)
                                         taattaatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGAGGGGGCGGGCCAAAGccccccccccccccccccccccgcgc                                                                                                           >  1:1953334/1‑144 (MQ=255)
                                         taattaatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCCGGGGGGGGGCAAAAGCCCCCCCCCCCCCCAcccccccgcgc                                                                                                           >  1:1084954/1‑144 (MQ=37)
                                          aattaatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGctct                                                                                                          <  2:2749360/150‑1 (MQ=255)
                                                tGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGaa                                                                                                    <  2:1795902/150‑1 (MQ=255)
                                                tGTCATCAGGTCCGAAAATAACGAGAATATTTAAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGAATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGaa                                                                                                    <  2:111927/150‑1 (MQ=255)
                                                ggTCAGCAGGTCCGACACTCACGAGAATCTTTCAGGCGCTCATCCTGGTGCGCTCCTGTCAGGTGCATTGCTTCATATACTCACGGGCGCAAGGCGCGCGCAGGGGGCGGCCCATCGCCGCCGCCCCCTGCACCCCCGGGCTCGGGCGaa                                                                                                    <  2:6104/149‑1 (MQ=21)
                                                  tcatcaGGTCCGAAAATAACGAGAATATTTCAGTGTCTGATCCTGTTGCGGTCCTGTCATGTGCATTGCTTCATATATTCACTGGCGCAAGGAGCGCGGGGGGGGGGGGGAATTGGCGGCGGCCCCGGGACCCCCGGGGTGTgggggacg                                                                                                  >  2:2901954/1‑144 (MQ=255)
                                                       aGGTCCGAAAATAACGAGAATATTTCAGTCGCTCAGCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCCGGGGGCGGCCCAACCCCCCCCCCCCCCGCCCCCCCCCGCCCCGGCGCACcaaaac                                                                                             >  2:2594189/1‑148 (MQ=255)
                                                         gTCCTAAAATAACGTTAATTTTTCTGTCTCTCATCTTGTTGCGATCCTGTCTTGTGCTTTGATTCATATAATCAGTGGCGCACGGCGCGCGCTGGGGGCGGCAAAGCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATcgc                                                                                           <  1:28477/150‑1 (MQ=18)
                                                         gTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCCGGGGGCCGCCCAACCCCGCCCCCCCCCCCCCCCCCCCGCCCGGGCGCaaaaaaaagc                                                                                           >  2:2325276/1‑146 (MQ=25)
                                                         gTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCATTCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGATCAAAATcgc                                                                                           <  2:2917496/150‑1 (MQ=255)
                                                         gTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATcgc                                                                                           <  2:729505/150‑1 (MQ=255)
                                                          tCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATcgcc                                                                                          <  2:2292781/150‑1 (MQ=255)
                                                                  tAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGcgcg                                                                                  <  2:2119779/150‑1 (MQ=255)
                                                                       agaCTATTTCAGTCTCTCAGCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCACGGGGCGGCCCATCGCCCCCCCCCCCGGCACCCCCCGGCTCAGGCGCCCCAAAACGCCGCCGCGCGGTGcc                                                                             >  2:535021/1‑150 (MQ=21)
                                                                            ttttCCGTCTCTAATCCTGTGGCGCTCCTGTCTGGTGATTTGCTTCTTATAATCACGGGCGCACGGCGCGCGCGGGGGGGGGCCAACCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCgg                                                                        <  1:3226726/149‑1 (MQ=17)
                                                                            ttttCCGTCTCAAATCCTGTGGCGCCCCGGTCTTGTGAATTGCTTAATATAATCACTGGCGCATGGCGCGCGCTGGGGGCGGCCATGCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCgg                                                                        <  1:3226755/149‑1 (MQ=17)
                                                                              ttCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAAAACCCCCCCCCCCCACCACCCCCCCGCACAGGCGAAAAAAAAACACCCAGCGCGGTGCCCACGGca                                                                      >  2:488180/1‑149 (MQ=16)
                                                                                 aGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCACGGGGCGGCCAAACCCCCCCCCCCCCACCACCCCCCGGCACAGGCGAACAAAAAAGCCCCGGCCCGGTGCCCTCGGCGTAt                                                                   >  2:1842516/1‑150 (MQ=18)
                                                                                   tctctcATCCTGTTGGGCTCCTGTCTTGTGATTTGCTTCATATAATCCCTGGCGCAAGGTGCGCGCTGGGGGCGGCCACTCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATcc                                                                 <  1:161041/150‑1 (MQ=17)
                                                                                      ctcATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCtt                                                              <  2:2037512/150‑1 (MQ=255)
                                                                                       tcATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTa                                                             <  2:811431/150‑1 (MQ=255)
                                                                                         aTCCTGTTGCGCTCTTGTCTTGTGATTTGCTTAATATAATCACGGGCGCAGGGCGCGCGCTGGGGGGGCCAAATCCCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACg                                                           <  1:1138599/150‑1 (MQ=17)
                                                                                             tGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCCAACGCCGCCCCCCCCAGCACCCCCCGGCGCAGGCGAACAAAAAAGCCCCGGCGCGGTGCCCACCGCGTAACCCGTAAGGcaa                                                       >  2:2409722/1‑148 (MQ=17)
                                                                                                 gggcTCGTGTCTTGTGCATTGCTTCATTTAATCACTGGCGCATGGAGCGCGCGGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCgg                                                   <  1:1528780/148‑1 (MQ=17)
                                                                                                 gcgcTCCTGTCTTGTGCTTTGATTCATATAATCACTGGCGCATGGCGCGCGCTGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCgg                                                   <  1:2025497/150‑1 (MQ=17)
                                                                                                 gcgcTCCTGTCATGTGCATTGCTTCATATAATCCCTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCgg                                                   <  2:2575553/150‑1 (MQ=255)
                                                                                                 gcgcTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCgg                                                   <  2:2575381/150‑1 (MQ=255)
                                                                                                   gcTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAAACGTCGCCGCCCCCAGCCCCCCCGGGCGCTGGCCAACAAAATCGCCGCGGCGCGGTGCCCGCGGCATATCCCCTACCGCTACCCGGt                                                 >  2:2853397/1‑150 (MQ=17)
                                                                                                    cTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTc                                                <  2:1914242/150‑1 (MQ=255)
                                                                                                       cTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGGGGCCAAAAGCCGCCGCCCCCAGCACCCCCCGGCGCAGGCGAACAAAAACCCCCCGGCGCGGTGCCCACCGCTTATCCCATAACGCAACCCCGTCggg                                             >  2:3197422/1‑150 (MQ=17)
                                                                                                       cTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCggg                                             <  2:1461583/150‑1 (MQ=255)
                                                                                                        tGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATAGCCGCCGCCCCCTGCTCCCCCGGGCTCTGGCGAACAAAATCGCCGGTGCGCGGTGCGCTCGGCTTATCCCTTTCGGCTACCGGGTCgggg                                            >  2:2421618/1‑149 (MQ=17)
                                                                                                           cATGTGCATTGCTTCATATAATCACTGGCGCACGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGcgcg                                         <  2:1458061/150‑1 (MQ=255)
                                                                                                            aTGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGGGGCCAATAGCCGGCGCCCCCAGGACCCCCGGGGTCTTTCGAACAAAATCGCCGCTTCGCGGTTCCCTCGGCTTATACCTTAAGGCTACCGGGTCggggggga                                        >  2:1873881/1‑145 (MQ=17)
                                                                                                             tgtgCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAg                                       <  2:1546650/150‑1 (MQ=255)
                                                                                                             tgtgCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAg                                       <  2:315740/150‑1 (MQ=255)
                                                                                                               tgctttgcttCATATAATCACTGGCGCAAGGAGCGCGCCGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGt                                     <  2:1414303/150‑1 (MQ=255)
                                                                                                                gCATTGCTTCATATAATCAGTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTa                                    <  2:2105750/150‑1 (MQ=255)
                                                                                                                gCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTa                                    <  2:2034768/150‑1 (MQ=255)
                                                                                                                 cATTGCTTCATATAACCCCTGGCGCAAGGAGCGCGCCGGGGGCGGCCAACCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTaa                                   <  1:1200639/150‑1 (MQ=255)
                                                                                                                    tGCTTCATATAATCACTGGCGCAAGGAGCGCGCTGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACAt                                <  1:708930/150‑1 (MQ=255)
                                                                                                                      ctacATATAATCACTGGCGCATGGAGCGCGCAGGGGGCGGCCACTCGCCGCCGCGCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTATCATcc                              <  1:1932036/147‑1 (MQ=11)
                                                                                                                      ctacATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAAACGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATcc                              <  1:2901954/147‑1 (MQ=255)
                                                                                                                      cTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATcc                              <  2:1625754/150‑1 (MQ=40)
                                                                                                                      cTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATcc                              <  2:1272712/150‑1 (MQ=40)
                                                                                                                        tCATATACTCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCACTCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCt                            <  2:2828875/150‑1 (MQ=255)
                                                                                                                        tCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCt                            <  1:506418/150‑1 (MQ=39)
                                                                                                                        tCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCt                            <  1:2169783/150‑1 (MQ=39)
                                                                                                                        tCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCt                            <  1:1389020/150‑1 (MQ=39)
                                                                                                                        tCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCt                            <  2:2801809/150‑1 (MQ=39)
                                                                                                                              aaTCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTaaaa                      <  1:2284784/150‑1 (MQ=34)
                                                                                                                               aTCACTGGCGCACGGAGCGCGCCGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAc                     <  2:1456369/150‑1 (MQ=18)
                                                                                                                               aTCACTGGCGCAAGGAGCGCGCTGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATTCCTGTAAAAc                     <  1:1091436/150‑1 (MQ=18)
                                                                                                                               aTCACTGGCGCAAGGAGCGCGCTGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATTCCTGTAAAAc                     <  1:907760/150‑1 (MQ=18)
                                                                                                                                 cACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAcgc                   <  1:454812/150‑1 (MQ=32)
                                                                                                                                   cTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACgcgc                 <  2:647041/150‑1 (MQ=32)
                                                                                                                                          aaGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGc          <  1:126968/150‑1 (MQ=32)
                                                                                                                                          aaGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGc          <  2:1794874/150‑1 (MQ=32)
                                                                                                                                            ggAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGccc        <  2:434177/150‑1 (MQ=32)
                                                                                                                                             gAGCGCGCTGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCca       <  1:3045943/150‑1 (MQ=17)
                                                                                                                                             gAGCGCGCGGGGTGCGGCCTATCGCCGCCGCCCCCTGCACCCCCGGGCTATGGCGAACTAAATAGCCGCTGCGCGGTGCCGTCGGCTTATCCCTTACGGCTACGGGGTCGGGCGCGAGGTAACTTCCCTGTAATACGAGCGCTCAGccga       <  1:2473212/150‑3 (MQ=11)
                                                                                                                                                  cgcAGGGGGCGGCAAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACAtgg  <  2:1404991/150‑3 (MQ=17)
                                                                                                                                                      |                                                                                                                                                   
TACTGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGCAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCC  >  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: