| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_000913 | 4,037,458 | G→A | 31.6% | intergenic (+123/‑61) | alaT → / → rrlA | tRNA‑Ala/23S ribosomal RNA of rrnA operon |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 4,037,458 | 0 | G | A | 31.6% | 66.9 / 52.1 | 57 | intergenic (+123/‑61) | alaT/rrlA | tRNA‑Ala/23S ribosomal RNA of rrnA operon |
| Reads supporting (aligned to +/- strand): ref base G (20/19); new base A (9/9); total (29/28) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.22e-01 | |||||||||||
TCCACCATCTCTGTAGTGATTAAATAAAAAATACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGC > NC_000913/4037329‑4037594 | tCCACCATCTCTGTAGTGATTAAATAAAAAATACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGgt < 4:1317294‑M2/132‑3 (MQ=255)tCCACCATCTCTGTAGTGATTAAATAAAAAATACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGgt > 3:1317294‑M2/1‑130 (MQ=255) aTCTCTGTAGTGATTAAATAAAAAATACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAacac < 3:215475/149‑1 (MQ=31) aTCTCTGTAGTGATTAAATAAAAAATACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAacac > 4:215475/1‑149 (MQ=31) tctGTAGTGATTAAATAAAAAATACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACgat < 3:19339/149‑1 (MQ=33) gTGATTAAATAAAAAATACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGaa < 2:1140891/149‑1 (MQ=37) aTTAAATAAAAAATACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCg < 3:820461/149‑1 (MQ=35) aaataaaAAATACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCgaaa < 4:1316889/149‑1 (MQ=38) taaaAAATACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCgaaagaa > 4:2278/1‑149 (MQ=37) aaaaaTACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCg < 4:374215/141‑1 (MQ=35) aaaaaTACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCg > 3:374215/1‑141 (MQ=35) aaTACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACgat > 1:589665/1‑129 (MQ=32) aaTACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACgat < 2:589665/129‑1 (MQ=32) ttCAGAGTGTGCCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAACTTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAGTTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCgg > 4:1271021/1‑149 (MQ=255) gagTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTtg > 1:830696/1‑149 (MQ=38) agTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTtgt > 2:851201/1‑149 (MQ=39) tgtACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGa > 1:726707/1‑149 (MQ=37) gtACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAg < 1:851201/149‑1 (MQ=39) gtACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAg < 2:726707/149‑1 (MQ=38) aCCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACga > 1:475097/1‑113 (MQ=33) aCCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACga < 2:475097/113‑1 (MQ=32) gCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGTTAAg > 1:1163473/1‑149 (MQ=38) aaaGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCgaaa > 4:1102507/1‑120 (MQ=40) aaaGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCgaaa < 3:1102507/120‑1 (MQ=40) aCTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCg > 2:494553/1‑149 (MQ=38) aCTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCg < 2:830696/149‑1 (MQ=38) aCTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCg > 1:1198812/1‑149 (MQ=37) tGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAacac > 2:1110134/1‑95 (MQ=32) tGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAacac < 1:1110134/95‑1 (MQ=32) tGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTc < 4:1381780/121‑1 (MQ=37) tGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTc > 3:1381780/1‑121 (MQ=37) gAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGTTaa > 3:1071494/1‑133 (MQ=33) gAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGTTaa < 4:1071494/133‑1 (MQ=33) aaaCACTGAACAACGAAAGTTGTTCGgtgtttcaattttcagcttgatccagatttttaaagagcaaatatctcaaacatcacccgaa > 2:1278882‑M2/1‑26 (MQ=255) aaaCACTGAACAACGAAAGTTGTTCGgtgtttcaattttcagcttgatccagatttttaaagagcaaatatctcaaacatcacccgaa < 1:1278882‑M2/88‑63 (MQ=255) acTGAACAACGAAAGTTGTTCGgtgtttcaattttcagcttgatccagatttttaaagagcaaatatctcaaacatcacccgaagatgagttttgagatattaa > 4:514553‑M2/1‑22 (MQ=255) acTGAACAACGAAAGTTGTTCGgtgtttcaattttcagcttgatccagatttttaaagagcaaatatctcaaacatcacccgaagatgagttttgagatattaa < 3:514553‑M2/104‑83 (MQ=255) acTGAACAACGAAAGTTGTTCGgtgtttcaattttcagcttgatccagatttttaaagagca > 4:814118‑M2/1‑22 (MQ=255) acTGAACAACGAAAGTTGTTCGgtgtttcaattttcagcttgatccagatttttaaagagca < 3:814118‑M2/62‑41 (MQ=255) acTGAACAACGAAAGTTGTTCAgtgtttcaattttcagcttgatccagatttttaaagagcaaatatctcaaacatcacccgaagatgagttttgagatattaaggcaggtgactttcactcacaaaccagcaagtggcgtcccctagg < 2:1355659‑M2/149‑128 (MQ=255) acTGAACAACGAAAGTTGTTCAgtgtttcaattttcagcttgatccagatttttaaagagcaaatatctcaaacatcacccgaagatgagttttgagatattaaggcaggtgactttcactcacaaac > 4:1236019‑M2/1‑22 (MQ=255) acTGAACAACGAAAGTTGTTCAgtgtttcaattttcagcttgatccagatttttaaagagcaaatatctcaaacatcacccgaagatgagttttgagatattaaggcaggtgactttcactcacaaac < 3:1236019‑M2/128‑107 (MQ=255) acTGAACAACGAAAGTTGTTCAgtgtttcaattttcagcttgatccagatttttaaagagcaaatatctcaaacatcacccgaagatgagttttgagatattaag > 3:742096‑M2/1‑22 (MQ=255) acTGAACAACGAAAGTTGTTCAgtgtttcaattttcagcttgatccagatttttaaagagcaaatatctcaaacatcacccgaagatgagttttgagatattaag < 4:742096‑M2/105‑84 (MQ=255) cTGAACAACGAAAGTTGTTCAgtgtttcaattttcagcttgatccagatttttaaagagcaaatatctcaaacatcacccgaagatgagttttgagatattaaggcaggtgactttcactcacaaaccagcaagtggcgtcccctaggg > 1:1355659‑M2/1‑21 (MQ=255) cTGAACAACGAAAGTTGTTCAgtgtttcaattttcagcttgatccagatttttaaagagcaaatatctcaaacatcacccgaagatgagttttgagatattaaggcaggtgactttcactcacaaaccagcaagtggcgtcccctaggg < 2:1324076‑M2/149‑129 (MQ=255) cTGAACAACGAAAGTTGTTCAgtgtttcaattttcagcttgatccagatttttaaagagcaaatatctcaaacatcacccgaa > 1:1025247‑M2/1‑21 (MQ=255) cTGAACAACGAAAGTTGTTCAgtgtttcaattttcagcttgatccagatttttaaagagcaaatatctcaaacatcacccgaa < 2:1025247‑M2/83‑63 (MQ=255) gAACAACGAAAGTTGTTCAgtgtttcaattttcagcttgatccagatttttaaagagcaaatatctcaaacatcacccgaagatgagttttgagatattaaggcaggtgactttcac > 3:1181920‑M2/1‑19 (MQ=255) gAACAACGAAAGTTGTTCAgtgtttcaattttcagcttgatccagatttttaaagagcaaatatctcaaacatcacccgaagatgagttttgagatattaaggcaggtgactttcac < 4:1181920‑M2/117‑99 (MQ=255) gAACAACGAAAGTTGTTCAgtgtttcaattttcagcttgatccagatttttaaagagcaaatatctcaaacatcacccgaagatgagttttgagatattaag > 4:47890‑M2/1‑19 (MQ=255) gAACAACGAAAGTTGTTCAgtgtttcaattttcagcttgatccagatttttaaagagcaaatatctcaaacatcacccgaagatgagttttgagatattaag < 3:47890‑M2/102‑84 (MQ=255) gAACAACGAAAGTTGTTCAgtgtttcaattttcagcttgatccagatttttaaaga < 3:1229240‑M2/56‑38 (MQ=255) gAACAACGAAAGTTGTTCAgtgtttcaattttcagcttgatccagatttttaaaga > 4:1229240‑M2/1‑19 (MQ=255) aacaacGAAAGTTGTTCAgtgtttcaattttcagcttgatccagatttttaaagagcaaatatctcaaacatcacccgaagatgagttttga < 4:1409426‑M2/92‑75 (MQ=255) aacaacGAAAGTTGTTCAgtgtttcaattttcagcttgatccagatttttaaagagcaaatatctcaaacatcacccgaagatgagttttga > 3:1409426‑M2/1‑18 (MQ=255) cGAAAGTTGTTCAgtgtttcaattttcagcttgatccagatttttaaagagcaaatatctcaaacatcacccgaagatgagttttgagatattaaggcaggtgactttcactcacaaaccagcaagtggcgtcccctaggggattcgaa > 1:1324076‑M2/1‑13 (MQ=255) | TCCACCATCTCTGTAGTGATTAAATAAAAAATACTTCAGAGTGTACCTGCAAAGGTTCACTGCGAAGTTTTGCTCTTTAAAAATCTGGATCAAGCTGAAAATTGAAACACTGAACAACGAAAGTTGTTCGTGAGTCTCTCAAATTTTCGCAACACGATGATGAATCGAAAGAAACATCTTCGGGTTGTGAGGTTAAGCGACTAAGCGTACACGGTGGATGCCCTGGCAGTCAGAGGCGATGAAGGACGTGCTAATCTGCGATAAGC > NC_000913/4037329‑4037594 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |