Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,229,177 A→C 23.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,1770AC23.7% 88.0 / 7.3 63intergenic (+112/+261)yohP/dusCuncharacterized protein/tRNA‑dihydrouridine synthase C
Reads supporting (aligned to +/- strand):  ref base A (16/26);  new base C (12/2);  total (31/32)
Fisher's exact test for biased strand distribution p-value = 4.34e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGG‑CGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATC  >  NC_000913/2229035‑2229322
                                                                                                                                               |                                                                                                                                                  
ggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGgg‑cg                                                                                                                                            >  2:402608/1‑150 (MQ=255)
ggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAAGGgg‑cg                                                                                                                                            >  2:402588/1‑150 (MQ=255)
 gtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGgg‑cgg                                                                                                                                           <  2:1069547/150‑1 (MQ=255)
 gtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGgg‑cgg                                                                                                                                           <  2:1469632/150‑1 (MQ=255)
  tgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACCGGCGCAAGGAGCGC‑GGCGGGgg‑gggc                                                                                                                                          >  2:1063006/1‑150 (MQ=255)
   gACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCCGGGGG‑CGGcc                                                                                                                                         >  2:1220665/1‑150 (MQ=255)
   gACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGG‑CGGcc                                                                                                                                         >  2:1953684/1‑150 (MQ=255)
   gACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGG‑CGGcc                                                                                                                                         >  2:1953662/1‑150 (MQ=255)
       ggCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGG‑CGGCCAATc                                                                                                                                     <  2:2221517/150‑1 (MQ=255)
       ggCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGG‑CGGCCAATc                                                                                                                                     <  1:1491334/150‑1 (MQ=255)
        gCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCCGGGGG‑GGGCCaaaac                                                                                                                                    >  2:282144/1‑147 (MQ=255)
        gCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGG‑CGGGCCATcg                                                                                                                                    >  1:1407142/1‑150 (MQ=255)
        gCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGG‑CGGCCAATcg                                                                                                                                    >  1:1133287/1‑150 (MQ=255)
              ttAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGG‑CGGCCAAAcgccgccc                                                                                                                              >  2:392640/1‑149 (MQ=255)
                 aGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGG‑CGGCCAATCGCCGCCGccc                                                                                                                           <  2:1133287/150‑1 (MQ=255)
                     gatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCCGGGGGCCGGCCCAACccccccccccccc                                                                                                                        >  1:822/1‑149 (MQ=255)
                            tcaaaTTAATTAATGTCATCAGGTCCGAAAATAACGATAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGG‑CGGCCAATCGCCGCCGCCCCCTGCAccccc                                                                                                                <  2:612049/148‑1 (MQ=255)
                            tAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGG‑CGGCCAATCGCCGCCGCCCCCTGCAccccc                                                                                                                <  2:612154/150‑1 (MQ=255)
                            tAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGG‑CGGCCAATAGCCGCCGCCCCCTGCcccccc                                                                                                                >  2:1605960/1‑150 (MQ=255)
                               aattaattAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCCGGGGG‑CGGCCCAACGCCGCCGCCCCCAGCACCCCCcgg                                                                                                             >  2:951755/1‑150 (MQ=38)
                               aattaattAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGGCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCCGGGGG‑GGGCCCCACGCCGCCCCCCCCAccccccccccg                                                                                                             >  1:979632/1‑147 (MQ=37)
                                attaattaATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GGGGGGGG‑GGGGGGATTGGCGGCGCCCCCGGCACCCCCcggg                                                                                                            >  1:2125602/1‑149 (MQ=255)
                                attaattaATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GGCGGGGG‑GGGCCCAAAGCCCCCCCCCCCAcccccccccgcc                                                                                                            >  1:2322678/1‑146 (MQ=255)
                                attaattaATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCCGGGGG‑CGGCCCAGcccccccccccccccccccccccccc                                                                                                            >  1:1137861/1‑146 (MQ=255)
                                  taattaatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCCGGGGG‑CGGCCAAACGCCGCCCCCCCCAGCCCCCCCCGGcgc                                                                                                          >  2:1959857/1‑148 (MQ=37)
                                  taattaatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCCGGGGG‑CCGCCAAAACCCCCCCCCCCCCGCccccccccgcgc                                                                                                          >  2:767063/1‑144 (MQ=37)
                                    attaatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGGGGGGGG‑GGGCCAAAcccccccccccccccccccccccccccag                                                                                                         >  2:1013955/1‑143 (MQ=255)
                                    attaatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGG‑CGGCCCAGAGCCCCCCCCCCCCGCCCCCCCCGGcgcag                                                                                                        >  2:319519/1‑146 (MQ=37)
                                    attaatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGG‑CGGCCCAACGCCGCCGCCCCCCGCCCCCCCCGGcgcag                                                                                                        >  2:1524577/1‑146 (MQ=37)
                                    attaatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGG‑CGGCCCAACGCCCCCCCCCCCAGCACCCCCCGGCACTg                                                                                                        >  2:1524527/1‑150 (MQ=37)
                                     ttaatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGG‑CGGCCCAACGCCGCCCCCCCCAGCACCCCCGGGCTCTgg                                                                                                       >  2:2173754/1‑150 (MQ=255)
                                       aatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCGGGGGG‑GGGCCCCACCCCGGCGCCCCCCGCCCCCCCCCGCGCTGGCg                                                                                                     >  2:2293531/1‑150 (MQ=37)
                                       aatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCCGGGGG‑CGGGCCAACGCGGGCCCCCCCCGCCCCCCCCGGCGCTGGCg                                                                                                     >  1:2130304/1‑150 (MQ=37)
                                                aGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGG‑CGGCCAAACGCCCCCCCCCCCCGCACCCCCGGGCTCTGGCGAAAAAAATc                                                                                            >  2:381931/1‑150 (MQ=25)
                                                aGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGAAAGGAGCGC‑GCACGGGG‑CGGCCAAAAGCCGCCGCCCCCAGCACCCCCCGGCTCTGGCGAACaaaaaa                                                                                            >  2:1948796/1‑148 (MQ=255)
                                                        aaaTAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGG‑CGGCCAATCGCCGCCG‑CCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTgcg                                                                                   <  2:1082201/150‑1 (MQ=255)
                                                        aaaTAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGG‑CGGCCAATCGCCGCCG‑CCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTgcg                                                                                   <  2:1082202/150‑1 (MQ=255)
                                                                agaATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCCTATAATCACTGGCGCACGGAGCGC‑GCAGGGGG‑CGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGcc                                                                            <  2:2143436/150‑1 (MQ=255)
                                                                agaATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGG‑CGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGcc                                                                            <  2:1915517/150‑1 (MQ=255)
                                                                   atttttCAGTCTCTCATCCGGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGG‑CGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTc                                                                         <  2:603671/147‑1 (MQ=255)
                                                                     atTTCAGTCTCACATCCGGTTGCGCGCCTGTCTTGTGCTTTGCTACATATAAACACTGGCGCATGGAGCGC‑GCCGGGGG‑CGGCCAATCGCCGCCGCCCCCGGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCgg                                                                       <  1:334886/150‑1 (MQ=17)
                                                                              tctcATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGG‑CGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCt                                                              <  2:1524252/150‑1 (MQ=255)
                                                                                       gTTGCGCTCCTGTCTTGTGCTTTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCTGGGGG‑CGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTAc                                                     <  1:1242138/150‑1 (MQ=255)
                                                                                       gTTGCGCTCCTGTCCTGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCTGGGGG‑CGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTAc                                                     <  1:1242165/150‑1 (MQ=255)
                                                                                       gTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCACGGAGCGC‑GCTGGGGG‑CGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTAc                                                     <  1:1242152/150‑1 (MQ=255)
                                                                                          gcgcTCCTGTATTGTGCATTGCTACATATAATCACTGGCGCACGGAGCGC‑GCAGGGGG‑CGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCgg                                                  <  1:2069710/150‑1 (MQ=255)
                                                                                           cgcTCGTGTCTTGTGCTTTGCTTCATATAATCCGTGGCGCACGGCGCGC‑GCTGGGGG‑CGGCCAAGCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCggg                                                 <  1:223954/150‑1 (MQ=14)
                                                                                           cgcTCCTGTCATGTGCTTTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGG‑CGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCggg                                                 <  1:1365066/150‑1 (MQ=255)
                                                                                           cgcTCCTGTCATGTGCATTGCTTCATATAATCACGGGCGCAAGGAGCGC‑GCAGGGGG‑CGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCggg                                                 <  1:1365107/150‑1 (MQ=255)
                                                                                                cTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCCGGGGG‑CGGGCCAGACCCGCCGCCCCCGCCACCCCCCGGCTCTGGCGCACAAAAAAGCCGCTGCGCGGGGCCCGAGGCGTATCCCATACGGCTACCCGGTCggg                                            >  2:1340688/1‑150 (MQ=255)
                                                                                                cTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGG‑CGGCCCATAGCCGCCCCCCCCTGCAACCCCGGGCTCTGGCGAACAAAAAACCCGCTGCGCGGTGCCCACGGCGTATACCATAAGGCTACCCGGTCggg                                            >  2:1340689/1‑150 (MQ=17)
                                                                                                       gtgCATTGCTTCATATAATCACTGGCGCACGGCGCGC‑GCCGGGGG‑CGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAgg                                     <  1:811612/150‑1 (MQ=255)
                                                                                                                 tCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGG‑CGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCt                           <  2:651503/150‑1 (MQ=39)
                                                                                                                 tCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGG‑CGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCt                           <  1:1192608/150‑1 (MQ=39)
                                                                                                                          cACTGGCGCAAGGAGCGC‑GCAGGGGG‑CGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAcgc                  <  1:1953662/150‑1 (MQ=32)
                                                                                                                          cACTGGCGCAAGGAGCGC‑GCAGGGGG‑CGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAcgc                  <  1:1220665/150‑1 (MQ=32)
                                                                                                                          cACTGGCGCAAGGAGCGC‑GCAGGGGG‑CGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAcgc                  <  1:1953684/150‑1 (MQ=32)
                                                                                                                            cTGGCGCAAGGAGCGC‑GCAGGGGG‑CGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACgcgc                <  2:391796/150‑1 (MQ=32)
                                                                                                                              ggCGCAAGGTGCGC‑GCAGGGGG‑CGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGccc              <  1:879203/150‑1 (MQ=17)
                                                                                                                                    aGGAGCGC‑GCAGGGGG‑CGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGcc        <  2:2082949/150‑1 (MQ=32)
                                                                                                                                      gAGCGC‑GCAGGGGG‑CGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCca      <  1:1959857/150‑1 (MQ=32)
                                                                                                                                      gAGCGC‑GCAGGGGG‑CGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCca      <  1:767063/150‑1 (MQ=32)
                                                                                                                                          gc‑gcAGGGGG‑CGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTCATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATc  >  2:148861/1‑150 (MQ=17)
                                                                                                                                               |                                                                                                                                                  
GGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGG‑CGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATC  >  NC_000913/2229035‑2229322

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