Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 16,749 T→A intergenic (+192/+2) insL1 → / ← hokC IS186/IS421 family transposase/protein HokC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091316,7490TA88.2% 35.0 / ‑0.1 18intergenic (+192/+2)insL1/hokCIS186/IS421 family transposase/protein HokC
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base A (12/3);  total (14/4)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.66e-01

GATTTCCCCCCCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGT‑CATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCC‑TTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAAT‑CATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACT  >  NC_000913/16516‑16994
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                      
gATTTCCCCCCCAGAAGGGCCGGAACCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGT‑CATCTGAAGCTTACAGGTGGCCATTCGTGGGACAGTAGCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGACACCGAGGCGTAGCCTTCTTCAGCTGGTTCAACAC‑TAAAGTCAGCGCTTATGGGATCACTACCAGCCGTTGCAATTACTCTGATTCGTa                                                                                                                                                                                                                                        >  2:82568/1‑250 (MQ=11)
      cccccaaaaaggaccGGATCCAAAAAAAAGATCAAACTCGTTGTGGAGATAAACAAAAAGGGT‑CATCGGGAGCTTACGGGTGGCCATTCGGGGGACAGTCCCCCTGACAGCCTACAAAACTCAATTGAAAAACGCGAGGCATCGTCTTAACGAGGCACCGTGGCGTCGTATCCTGGAGTTGGTCCACCCC‑TAAAGTGAGTTTTTAGGGGACGTCTCCCCGCCGTTGCACTTACTCTGATTCGTAAGCCGTg                                                                                                                                                                                                                                 <  1:129756/238‑1 (MQ=2)
          caaGAAGTGCCGGATCAGAAAAGAAGAACTAACTCGTTGGGGAGAAAAAAAAAAATGGT‑CATCGGGAGCTTACAGGTGGCCATTCGGGGGACAGTACCCCTGACAGCCTAAAAAACGCAATGGAAGACCGCGAGGCATCGTTTTAACGAGGCCCCGGGGCGTCGCATTCTTCAGATGGTTAACCCC‑TTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCACTTACTCGGATTCGTAAGCCGTGaaaa                                                                                                                                                                                                                             <  1:208955/249‑1 (MQ=255)
                         cggaaaaaaaaaaaaactgcGTTGTGGAAAATAA‑AAAAAAGGTACATCTGAAATTTACAGGTGGCTATTCGGGGGACAGTACCCCTGACAGCCAAAAAACCGAAATGGAAAACCGCGAGGCATCGTCTTAACGAGGCCCCGAGGCGACGCATTCTTCA‑AAGGTTTAACCCCTTAAGTTAGCGCTTTTGGGACCACTCCCCGCCGTTGCACTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCt                                                                                                                                                                                                              <  1:40794/232‑1 (MQ=255)
                             aaaGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGT‑CATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCC‑TTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCACTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGcc                                                                                                                                                                                                          >  1:55964/1‑251 (MQ=255)
                                                   gagaATAACAAAAATGGT‑CATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCC‑TTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCACTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTcac                                                                                                                                                                                      >  1:345543/1‑249 (MQ=255)
                                                             aaaaTGGT‑CATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCC‑TTAAGTTAGCGCTTATGGGATCACTCCACGCCGTTGCACTTAATCGGATTCGTAAGCCGAGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACCGAGGTCTTTtc                                                                                                                                                                          >  1:306536/1‑251 (MQ=255)
                                                                                 tACAGGTGGCCATTCGTGGGGCAGTATCCCTGACATCCTACCAAACGCAATTGAAGGACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCC‑TTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGGGAAATCGGCAACCTCCGTCTGGc                                                                                                                                                                                                           >  2:171141/1‑199 (MQ=18)
                                                                                  aCAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCC‑TTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCACTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACt                                                                                                                                                      >  2:310384/1‑251 (MQ=255)
                                                                                    aGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCC‑TTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCACTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTAc                                                                                                                                                    >  1:299884/1‑251 (MQ=255)
                                                                                          ccATTCGTGGAACAGTATCCCTGACAGCCTAAAAAACGCAATTGAAGCACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCC‑TTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCACTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTAcggcggg                                                                                                                                              >  2:234487/1‑250 (MQ=255)
                                                                                           cATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCC‑TTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCACTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGACAGTTCGGATGTGAACCTCCCAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTg                                                                                                                                             >  1:381555/1‑251 (MQ=255)
                                                                                                                                  aTTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCC‑TTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCACTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGGTACAGAAGACGATCAGGGCGACAATACA‑CGCCTTAtt                                                                                                      >  1:68778/1‑249 (MQ=255)
                                                                                                                                      aagaaCGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCC‑TTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGCGATAAAGGTGACGATCAGGGCGACAAT‑CATACCCCTATgctgct                                                                                                 >  2:35870/1‑251 (MQ=255)
                                                                                                                                                          tAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCC‑TTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCACTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGGGATACAGATGACGGTCAGGGCGACAAT‑CACCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTg                                                                             >  1:214481/1‑251 (MQ=255)
                                                                                                                                                                                                                                   gccgTTGCACTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCAGTCTGGCCAGTTCGGCTGTGAACCTCACAGAGGACTTTTCTCGTTACCAGCGCCGCCACTACGGCGGCGAAAAAGAAGTCGATCAGGGGGACCAT‑CCACGGCTTATGCTACTGCATATCACTTTATTCCTTGACCTTACAGCCAGTAAACGGAACCATACATGTGTTTAGCATAAaggggggccccgggtgttagaaaaatat      >  2:215000/1‑225 (MQ=255)
                                                                                                                                                                                                                                    ccgTTGCACTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAAA‑CATACCCTTTTTCCTCTTTTCTGCTCTCTTCTCCTTACCCTTACAGTCAGTAAACGGCCCCCTACATGTGTTCCACATCTAGGGGGCCTCGGGTTGTTAGTAAAATATcc    >  2:175364/1‑250 (MQ=255)
                                                                                                                                                                                                                                       ttGCCCTTATTCGGATTCGTAACCCTTGAAACCAGCACCCCCCGTCTGCCCAGTTCGGATGTGACCCTCACAGAGGTTTTTTCTCTTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAAT‑CATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGTTGATAGTAAAATATCACt  <  1:215000/250‑1 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                      
GATTTCCCCCCCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGT‑CATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCC‑TTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCAACCTCCGTCTGGCCAGTTCGGATGTGAACCTCACAGAGGTCTTTTCTCGTTACCAGCGCCGCCACTACGGCGGTGATACAGATGACGATCAGGGCGACAAT‑CATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACT  >  NC_000913/16516‑16994

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