Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 3,908,439 | Δ1 bp | intergenic (‑73/+110) | pstB ← / ← pstA | phosphate ABC transporter ATPase/phosphate ABC transporter permease |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,908,439 | 0 | G | . | 97.2% | 152.2 / ‑1.4 | 36 | intergenic (‑73/+110) | pstB/pstA | phosphate ABC transporter ATPase/phosphate ABC transporter permease |
Reads supporting (aligned to +/- strand): ref base G (0/0); major base . (18/17); minor base C (0/1); total (18/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.25e-01 |
GCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTACTCGGGGCAGTTTCAACCATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCG > NC_000913/3908304‑3908568 | gCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTACTCGGGGCAGTTTCAACCATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCAC‑CCGCATCCGGCa > 7:26811/1‑147 (MQ=255) gtagAAGTTCAAATTACGAACCTGAATTTTACTCGGGGCAGTTTCAACCATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCAC‑CCGCATCCGGCAAAGTGtt > 8:92449/1‑148 (MQ=255) gtagAAGTTCAAATTACGAACCTGAATTTTACTCGGGGCAGTTTCAACCATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCAC‑CCGCATCCGGCAAAGTGTTa < 8:26811/149‑1 (MQ=255) tagAAGTTCAAATTACGAACCTGAATTTTACTCGGGGCAGTTTCAACCATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCAC‑CCGCATc < 3:94584/135‑1 (MQ=255) tagAAGTTCAAATTACGAACCTGAATTTTACTCGGGGCAGTTTCAACCATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCAC‑CCGCATc > 4:94584/1‑135 (MQ=255) tagAAGTTCAAATTACGAACCTCAATTATACTCGGCGCAGTTTCAACCATACTCATTGCAAACTAATCCACATCGTGCCTGATCCACTACGCTAATCAGGCTTACCGTAGCCCGGATCACGCGATCAC‑CCGCATCCGGCAAAGAGTTaa > 1:133876/1‑149 (MQ=255) cAAATTACGAACCTGAATTTTACTCGGGGCAGTTTCAACCATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCAC‑CCGCATCCGGCAAAGTGTTAACCGTGTa < 3:117559/148‑1 (MQ=255) cAAATTACGAACCTGAATTTTACTCGGGGCAGTTTCAACCATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCAC‑CCGCATCCGGCAAAGTGTTAACCGTGTa < 4:143841/148‑1 (MQ=255) ataTTACGATCCTGAATTTTACTCGGGGCAGTTTCAACCATACTCATTGCAATCTTTTCCTCTTCGTGGCTGATGCACTTCGCTAATGTGGCTTACCGTAGGCCGGATCTGGCGTTCAC‑CCGCATCCGGGAATGTGTTAACCGTGTACa < 2:133876/147‑1 (MQ=255) aTTACGAACCTGAATTTTACTCGGGGCAGTTTCAACCATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCAC‑CCGCATCCGGCAAAGTGTTAACCGTGTACa > 8:153522/1‑147 (MQ=255) gAACCTGAATTTTACTCGGGGCAGTTTCAACCATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCAC‑CCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCa > 6:189688/1‑148 (MQ=255) cTGAATTTTACTCGGGGCAGTTTCAACCATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCAC‑CCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAg > 5:10425/1‑149 (MQ=255) ggggCAGTTTCAACCATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCAC‑CCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTcggtcggt > 8:157493/1‑149 (MQ=255) gggCAGTTTCAACCATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCAC‑CCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCggtcggtc > 4:121358/1‑149 (MQ=255) ggCAGTTTCAACCATACTCATTGCAATCTTTTCCTCATCGTGCCTGCTGCACTTCGCTAATCATGATTACAGTAGTCCGGATCAGTCGTTCAC‑CCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTcc > 2:129850/1‑149 (MQ=255) ggCAGTTTCAACCATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCAC‑CCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTcc > 2:133217/1‑149 (MQ=255) tCAACCATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCAC‑CCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGc < 5:189688/148‑1 (MQ=255) aaCCATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCAC‑CCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGcccc > 7:92686/1‑149 (MQ=255) aCCATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCAC‑CCGCATCCGGCAAAGTGtt < 4:163078/102‑1 (MQ=255) aCCATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCAC‑CCGCATCCGGCAAAGTGtt > 3:163078/1‑102 (MQ=255) cTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCAC‑CCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTgggg < 8:92686/149‑1 (MQ=255) gCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCAC‑CCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAgg < 7:157493/148‑1 (MQ=255) cAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCAC‑CCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTggg < 1:157430/141‑1 (MQ=255) cAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCAC‑CCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTggg > 2:157430/1‑141 (MQ=255) tttCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCAC‑CCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCaaaaa < 4:54601/104‑1 (MQ=255) tttCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCAC‑CCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCaaaaa > 3:54601/1‑104 (MQ=255) ttCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCAC‑CCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAgg < 8:53363/146‑1 (MQ=255) ttCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTC‑CCCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAgg < 7:92449/146‑1 (MQ=255) aTCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCAC‑CCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAggg < 7:153522/147‑1 (MQ=255) aTCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCAC‑CCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAggg < 6:10425/147‑1 (MQ=255) ccTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCAC‑CCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCtcggtcg < 6:50365/103‑1 (MQ=255) aCTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCAC‑CCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGc > 8:77372/1‑148 (MQ=255) ttCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCAC‑CCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCaa > 7:114140/1‑148 (MQ=255) ttACCGTAGGCCGGATCAGGCGTTCAC‑CCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTa > 5:67606/1‑149 (MQ=255) tAGGCCGGATCAGGCGTTCAC‑CCGCATCCGGCAAAGTGTTAACCTTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTAttcttc < 3:121358/149‑1 (MQ=255) ggCCGGATCAGGCGTTCAC‑CCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCg < 7:77372/148‑1 (MQ=255) | GCCGTAGTAGAAGTTCAAATTACGAACCTGAATTTTACTCGGGGCAGTTTCAACCATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCG > NC_000913/3908304‑3908568 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |