Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,743,517:1 +TTG 100% coding (61/255 nt) rpsQ → 30S ribosomal subunit protein S17

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,743,5171.T81.3% 226.2 / 48.9 75coding (61/255 nt)rpsQ30S ribosomal subunit protein S17
Reads supporting (aligned to +/- strand):  ref base . (0/14);  new base T (61/0);  total (61/14)
Fisher's exact test for biased strand distribution p-value = 1.78e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
*minE2,743,5172.T81.3% 221.8 / 48.9 75coding (61/255 nt)rpsQ30S ribosomal subunit protein S17
Reads supporting (aligned to +/- strand):  ref base . (0/14);  new base T (61/0);  total (61/14)
Fisher's exact test for biased strand distribution p-value = 1.78e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.07e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
*minE2,743,5173.G81.3% 234.3 / 43.9 75coding (61/255 nt)rpsQ30S ribosomal subunit protein S17
Reads supporting (aligned to +/- strand):  ref base . (0/14);  new base G (61/0);  total (61/14)
Fisher's exact test for biased strand distribution p-value = 1.78e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

ATAAAATCCGTACTCTGCAAGGTCGCGTTGTTAGCGACAAAATGGAGAAATCCA‑‑‑TTGTTGTTGCTATCGAACGTTTTGTGAAACACCCGATCTACGGTAAATTCATCAAGCGTACGAC  >  minE/2743464‑2743581
                                                      |||                                                                
aTAAAATCCGTACTCTGCAAGGTCGCGTTGTTAGCGACAAAATGGAGAAATCCA‑‑‑TTGTTGTTGCTAt                                                     <  1:1041644/67‑1 (MQ=255)
aTAAAATCCGTACTCTGCAAGGTCGCGTTGTTAGCGACAAAATGGAGAAATCCA‑‑‑TTGTTGTTGCTAt                                                     <  1:1732514/67‑1 (MQ=255)
aTAAAATCCGTACTCTGCAAGGTCGCGTTGTTAGCGACAAAATGGAGAAATCCA‑‑‑TTGTTGTTGCTAt                                                     <  1:900719/67‑1 (MQ=255)
aTAAAATCCGTACTCTGCAAGGTCGCGTTGTTAGCGACAAAATGGAGAAATCCA‑‑‑TTGTTGTTGCTAt                                                     <  1:1370233/67‑1 (MQ=255)
aTAAAATCCGTACTCTGCAAGGTCGCGTTGTTAGCGACAAAATGGAGAAATCCA‑‑‑TTGTTGTTGCTAt                                                     <  1:2040781/67‑1 (MQ=255)
aTAAAATCCGTACTCTGCAAGGTCGCGTTGTTAGCGACAAAATGGAGAAATCCA‑‑‑TTGTTGTTGCTAt                                                     <  1:2123027/67‑1 (MQ=255)
                   aagTCGCGTTGTTAGCGACAAAATGGAGAAATCCATTGTTGTTGTTGCTATCGAACGTTTTGTGaaa                                     >  1:2100938/3‑67 (MQ=255)
                   aGGTCGCGTTGTTAGCGACGAAATGGAGAAATCCATTGTTGTTGTTGCTATCGAACGTTTTGTGaaa                                     >  1:608992/1‑67 (MQ=255)
                   aGGTCGCGTTGTTAGCGACAAAATTGAGAAATCCATTGTTGTTGTTGCTATCGAACGTTTTGTGcaa                                     >  1:1137160/1‑67 (MQ=255)
                   aGGTCGCGTTGTTAGCGACAAAATGGAGAAATCCATTGTTGTTGTTGCTATCGATCGTTTTGTGaaa                                     >  1:216235/1‑67 (MQ=255)
                   aGGTCGCGTTGTTAGCGACAAAATGGAGAAATCCATTGTTGTTGTTGCTATCGAACGTTTTGTGaaa                                     >  1:1936198/1‑67 (MQ=255)
                   aGGTCGCGTTGTTAGCGACAAAATGGAGAAATCCATTGTTGTTGTTGCTATCGAACGTTTTGTGaaa                                     >  1:23914/1‑67 (MQ=255)
                   aGGTCGCGTTGTTAGCGACAAAATGGAGAAATCCATTGTTGTTGTTGCTATCGAACGTTTTGTGaaa                                     >  1:220132/1‑67 (MQ=255)
                   aGGTCGCGTTGTTAGCGACAAAATGGAGAAATCCATTGTTGTTGTTGCTATCGAACGTTTTGTGaaa                                     >  1:2151458/1‑67 (MQ=255)
                   aGGTCGCGTTGTTAGCGACAAAATGGAGAAATCCATTGTTGTTGTTGCTATCGAACGTTTTGTGaaa                                     >  1:2105167/1‑67 (MQ=255)
                   aGGTCGCGTTGTTAGCGACAAAATGGAGAAATCCATTGTTGTTGTTGCTATCGAACGTTTTGTGaaa                                     >  1:518008/1‑67 (MQ=255)
                   aGGTCGCGTTGTTAGCGACAAAATGGAGAAATCCATTGTTGTTGTTGCTATCGAACGTTTTGTGaaa                                     >  1:1928881/1‑67 (MQ=255)
                   aGGTCGCGTTGTTAGCGACAAAATGGAGAAATCCATTGTTGTTGTTGCTATCGAACGTTTTGTGaaa                                     >  1:191553/1‑67 (MQ=255)
                   aGGTCGCGTTGTTAGCGACAAAATGGAGAAATCCATTGTTGTTGTTGCTATCGAACGTTTTGTGaaa                                     >  1:1911306/1‑67 (MQ=255)
                   aGGTCGCGTTGTTAGCGACAAAATGGAGAAATCCATTGTTGTTGTTGCTATCGAACGTTTTGTGaaa                                     >  1:1881638/1‑67 (MQ=255)
                   aGGTCGCGTTGTTAGCGACAAAATGGAGAAATCCATTGTTGTTGTTGCTATCGAACGTTTTGTGaaa                                     >  1:1852108/1‑67 (MQ=255)
                   aGGTCGCGTTGTTAGCGACAAAATGGAGAAATCCATTGTTGTTGTTGCTATCGAACGTTTTGTGaaa                                     >  1:492584/1‑67 (MQ=255)
                   aGGTCGCGTTGTTAGCGACAAAATGGAGAAATCCATTGTTGTTGTTGCTATCGAACGTTTTGTGaaa                                     >  1:503390/1‑67 (MQ=255)
                   aGGTCGCGTTGTTAGCGACAAAATGGAGAAATCCATTGTTGTTGTTGCTATCGAACGTTTTGTGaaa                                     >  1:185829/1‑67 (MQ=255)
                   aGGTCGCGTTGTTAGCGACAAAATGGAGAAATCCATTGTTGTTGTTGCTATCGAACGTTTTGTGaaa                                     >  1:546934/1‑67 (MQ=255)
                   aGGTCGCGTTGTTAGCGACAAAATGGAGAAATCCATTGTTGTTGTTGCTATCGAACGTTTTGTGaaa                                     >  1:653919/1‑67 (MQ=255)
                   aGGTCGCGTTGTTAGCGACAAAATGGAGAAATCCATTGTTGTTGTTGCTATCGAACGTTTTGTGaaa                                     >  1:662641/1‑67 (MQ=255)
                   aGGTCGCGTTGTTAGCGACAAAATGGAGAAATCCATTGTTGTTGTTGCTATCGAACGTTTTGTGaaa                                     >  1:796673/1‑67 (MQ=255)
                   aGGTCGCGTTGTTAGCGACAAAATGGAGAAATCCATTGTTGTTGTTGCTATCGAACGTTTTGTGaaa                                     >  1:809940/1‑67 (MQ=255)
                   aGGTCGCGTTGTTAGCGACAAAATGGAGAAATCCATTGTTGTTGTTGCTATCGAACGTTTTGTGaaa                                     >  1:831957/1‑67 (MQ=255)
                   aGGTCGCGTTGTTAGCGACAAAATGGAGAAATCCATTGTTGTTGTTGCTATCGAACGTTTTGTGaaa                                     >  1:851946/1‑67 (MQ=255)
                   aGGTCGCGTTGTTAGCGACAAAATGGAGAAATCCATTGTTGTTGTTGCTATCGAACGTTTTGTGaaa                                     >  1:96981/1‑67 (MQ=255)
                   aGGTCGCGTTGTTAGCGACAAAATGGAGAAATCCATTGTTGTTGTTGCTATCGAACGTTTTGTGaaa                                     >  1:970569/1‑67 (MQ=255)
                   aGGTCGCGTTGTTAGCGACAAAATGGAGAAATCCATTGTTGTTGTTGCTATCGAACGTTTTGTGaaa                                     >  1:971147/1‑67 (MQ=255)
                   aGGTCGCGTTGTTAGCGACAAAATGGAGAAATCCATTGTTGTTGTTGCTATCGAACGTTTTGTGaaa                                     >  1:972126/1‑67 (MQ=255)
                   aGGTCGCGTTGTTAGCGACAAAATGGAGAAATCCATTGTTGTTGTTGCTATCGAACGTTTTGTGaaa                                     >  1:975007/1‑67 (MQ=255)
                   aGGTCGCGTTGTTAGCGACAAAATGGAGAAATCCATTGTTGTTGTTGCTATCGAACGTTTTGTGaaa                                     >  1:1502058/1‑67 (MQ=255)
                   aGGTCGCGTTGTTAGCGACAAAATGGAGAAATCCATTGTTGTTGTTGCTATCGAACGTTTTGTGaaa                                     >  1:1120339/1‑67 (MQ=255)
                   aGGTCGCGTTGTTAGCGACAAAATGGAGAAATCCATTGTTGTTGTTGCTATCGAACGTTTTGTGaaa                                     >  1:1130198/1‑67 (MQ=255)
                   aGGTCGCGTTGTTAGCGACAAAATGGAGAAATCCATTGTTGTTGTTGCTATCGAACGTTTTGTGaaa                                     >  1:1182326/1‑67 (MQ=255)
                   aGGTCGCGTTGTTAGCGACAAAATGGAGAAATCCATTGTTGTTGTTGCTATCGAACGTTTTGTGaaa                                     >  1:1201528/1‑67 (MQ=255)
                   aGGTCGCGTTGTTAGCGACAAAATGGAGAAATCCATTGTTGTTGTTGCTATCGAACGTTTTGTGaaa                                     >  1:1277199/1‑67 (MQ=255)
                   aGGTCGCGTTGTTAGCGACAAAATGGAGAAATCCATTGTTGTTGTTGCTATCGAACGTTTTGTGaaa                                     >  1:10973/1‑67 (MQ=255)
                   aGGTCGCGTTGTTAGCGACAAAATGGAGAAATCCATTGTTGTTGTTGCTATCGAACGTTTTGTGaaa                                     >  1:1312412/1‑67 (MQ=255)
                   aGGTCGCGTTGTTAGCGACAAAATGGAGAAATCCATTGTTGTTGTTGCTATCGAACGTTTTGTGaaa                                     >  1:1337145/1‑67 (MQ=255)
                   aGGTCGCGTTGTTAGCGACAAAATGGAGAAATCCATTGTTGTTGTTGCTATCGAACGTTTTGTGaaa                                     >  1:1350593/1‑67 (MQ=255)
                   aGGTCGCGTTGTTAGCGACAAAATGGAGAAATCCATTGTTGTTGTTGCTATCGAACGTTTTGTGaaa                                     >  1:1374279/1‑67 (MQ=255)
                   aGGTCGCGTTGTTAGCGACAAAATGGAGAAATCCATTGTTGTTGTTGCTATCGAACGTTTTGTGaaa                                     >  1:109045/1‑67 (MQ=255)
                   aGGTCGCGTTGTTAGCGACAAAATGGAGAAATCCATTGTTGTTGTTGCTATCGAACGTTTTGTGaaa                                     >  1:1384908/1‑67 (MQ=255)
                   aGGTCGCGTTGTTAGCGACAAAATGGAGAAATCCATTGTTGTTGTTGCTATCGAACGTTTTGTGaaa                                     >  1:1395111/1‑67 (MQ=255)
                   aGGTCGCGTTGTTAGCGACAAAATGGAGAAATCCATTGTTGTTGTTGCTATCGAACGTTTTGTGaaa                                     >  1:1481199/1‑67 (MQ=255)
                   aGGTCGCGTTGTTAGCGACAAAATGGAGAAATCCATTGTTGTTGTTGCTATCGAACGTTTTGTGaaa                                     >  1:1059672/1‑67 (MQ=255)
                   aGGTCGCGTTGTTAGCGACAAAATGGAGAAATCCATTGTTGTTGTTGCTATCGAACGTTTTGTGaaa                                     >  1:180812/1‑67 (MQ=255)
                   aGGTCGCGTTGTTAGCGACAAAATGGAGAAATCCATTGTTGTTGTTGCTATCGAACGTTTTGTGaaa                                     >  1:1730545/1‑67 (MQ=255)
                   aGGTCGCGTTGTTAGCGACAAAATGGAGAAATCCATTGTTGTTGTTGCTATCGAACGTTTTGTGaaa                                     >  1:1707221/1‑67 (MQ=255)
                   aGGTCGCGTTGTTAGCGACAAAATGGAGAAATCCATTGTTGTTGTTGCTATCGAACGTTTTGTGaaa                                     >  1:1706683/1‑67 (MQ=255)
                   aGGTCGCGTTGTTAGCGACAAAATGGAGAAATCCATTGTTGTTGTTGCTATCGAACGTTTTGTGaaa                                     >  1:1647995/1‑67 (MQ=255)
                   aGGTCGCGTTGTTAGCGACAAAATGGAGAAATCCATTGTTGTTGTTGCTATCGAACGTTTTGTGaaa                                     >  1:1639830/1‑67 (MQ=255)
                   aGGTCGCGTTGTTAGCGACAAAATGGAGAAATCCATTGTTGTTGTTGCTATCGAACGTTTTGTGaaa                                     >  1:1560526/1‑67 (MQ=255)
                   aGGTCGCGTTGTTAGCGACAAAATGGAGAAATCCATTGTTGTTGTTGCTATCGAACGTTTTGTGaaa                                     >  1:1608113/1‑67 (MQ=255)
                   aGGTCGCGTTGTTAGCGACAAAATGGAGAAATCCATTGTTGTTGTTGCTATCGAACGTTTTGTGaaa                                     >  1:1603455/1‑67 (MQ=255)
                   aGGTCGCGTTGTTAGCGACAAAATGGAGAAATCCATTGTTGTTGTTGCTATCGAACGTTTTGTGaaa                                     >  1:1602085/1‑67 (MQ=255)
                   aGGTCGCGTTGTTAGCGACAAAATGGAGAAATCCATTGTTGTTGTTGCTATCGAACGTTTTGTGaaa                                     >  1:1594214/1‑67 (MQ=255)
                   aGGTCGCGTTGTTAGCGACAAAATGGAGAAATCCATTGTTGTTGTTGCTATCGAACGTTTTGTGaaa                                     >  1:1585350/1‑67 (MQ=255)
                   aGGTCGCGTTGTTAGCGACAAAATGGAGAAATCCATTGTTGTTGTTGCTATCGAACGTTTTGTGaaa                                     >  1:1569830/1‑67 (MQ=255)
                   aGGTCGCGTTGTTAGCGACAAAATGGAGAAATCCATTGTTGTTGTTGCTATCGAACGTTTTGTGaa                                      >  1:732458/1‑66 (MQ=255)
                   aGGTCGCGTTGTTAGCGACAAAATGGAGAAATCCATTGTTGTTGTTGCTATCGAACGTTTTGTGaa                                      >  1:828545/1‑66 (MQ=255)
                                                   ccA‑‑‑TTGTTGTTGCTATCGAACGTTTTGTGAAACACCCGATCTACGGTAAATTCATCAAGCGTACGAc  <  1:1082645/67‑1 (MQ=255)
                                                   ccA‑‑‑TTGTTGTTGCTATCGAACGTTTTGTGAAACACCCGATCTACGGTAAATTCATCAAGCGTACGAc  <  1:101772/67‑1 (MQ=255)
                                                   ccA‑‑‑TTGTTGTTGCTATCGAACGTTTTGTGAAACACCCGATCTACGGTAAATTCATCAAGCGTACGAc  <  1:92148/67‑1 (MQ=255)
                                                   ccA‑‑‑TTGTTGTTGCTATCGAACGTTTTGTGAAACACCCGATCTACGGTAAATTCATCAAGCGTACGAc  <  1:746645/67‑1 (MQ=255)
                                                   ccA‑‑‑TTGTTGTTGCTATCGAACGTTTTGTGAAACACCCGATCTACGGTAAATTCATCAAGCGTACGAc  <  1:543138/67‑1 (MQ=255)
                                                   ccA‑‑‑TTGTTGTTGCTATCGAACGTTTTGTGAAACACCCGATCTACGGTAAATTCATCAAGCGTACGAc  <  1:2039938/67‑1 (MQ=255)
                                                   ccA‑‑‑TTGTTGTTGCTATCGAACGTTTTGTGAAACACCCGATCTACGGTAAATTCATCAAGCGTACGAc  <  1:2014831/67‑1 (MQ=255)
                                                   ccA‑‑‑TTGTTGTTGCTATCGAACGTTTTGTGAAACACCCGATCTACGGTAAATTCATCAAGCGTACGAc  <  1:1001083/67‑1 (MQ=255)
                                                      |||                                                                
ATAAAATCCGTACTCTGCAAGGTCGCGTTGTTAGCGACAAAATGGAGAAATCCA‑‑‑TTGTTGTTGCTATCGAACGTTTTGTGAAACACCCGATCTACGGTAAATTCATCAAGCGTACGAC  >  minE/2743464‑2743581

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