Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_00091316,4900T→G100.0% 0.9 5I349S (ATT→AGT) insL1IS186 transposase
Rejected: E-value exceeds prediction threshold.
Reads supporting (aligned to +/- strand):  new base (0/3):  ref base (2/0):  total (2/3)
Fisher's exact test strand distribution p-value = 1.00e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

GAACCTGAACTCGCGAAAGCGTGGATATTTGCTAATCTACTCGCCGCATTTTTAATTGACGACATAATCCAGCCATCGCTGGATTTCCCCCCCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCA  >  NC_000913/16435‑16781
                                                       |                                                                                                                                                                                                                                                                                                   
gAACCTGAATGCGCGAAAGCGTGGATATTTGCGAATCTAGTCGCCGCATTTTTAAGTGAGTACATAATCCAGCCATAGCTGGAGTTCCCCCCCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATc                                                  <  2:234649/299‑1 (MQ=2)
               aaaGCGTGGATATTTGCTAATCTACTCGCCGCATTTTTAAGTGACGACATAATCCAGCCATCGCTGGATTTCCCCCCCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGctc                                 <  2:112158/301‑1 (MQ=21)
                               aataaagaCTCGCCTCAGTTTGAAGTGACGAAAAAATCAACGCATTCTGGGATTGCCACCCCGGCAGAGACGGATCCGAAAAGAAGAACGATGTCGTGGTGGAGAATAACAAAAATGCTCATCTGTAGCTCACAGGTGGCCATTCGTGGGACAGTATCCCTTACAGCCTACAAAAAGCAATAGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGTGGCGCCGCATTCGCCAGAAGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCAGTTGCTCTTACTCGGATTCGTaa                 <  1:639069/294‑1 (MQ=14)
                                      aCTCGCCGCATTTTTAATTGCCGACATAATCCAGCCATCGCTGGATTTCCCCCACAGAAGTGCCGGATCAGACAAGAAGACCTCACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGCAAGAATCCCTCACCGCCTACAACACGCCATTGAACACCGCGAGGCATCGACTTAACGAGGACCCGAAGCGTCGCCTTCTTCACATGGTTCAACCCTTAAGTTACCGCTTATGTGACCACCACCCGCCGTTCCTCTGCCCCGGCTTCggcccgttga           >  2:639043/1‑291 (MQ=11)
                                              cATTTTTAATTGACGACATAATCCAGCCATCGCTGGATTTCCCCCCCAGAAGTGCCGGATCCGAAAAGACGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGGGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCTTTGCTCTTCCTCGGCTTCTTAATCCTTGAAAACAGCa  >  2:670065/1‑301 (MQ=255)
                                                       |                                                                                                                                                                                                                                                                                                   
GAACCTGAACTCGCGAAAGCGTGGATATTTGCTAATCTACTCGCCGCATTTTTAATTGACGACATAATCCAGCCATCGCTGGATTTCCCCCCCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCA  >  NC_000913/16435‑16781

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.