Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | U00096 | 3,945,819 | 0 | G | A | 20.0% | 110.3 / 10.1 | 50 | noncoding (2116/2904 nt) | rrlC | 23S ribosomal RNA |
Reads supporting (aligned to +/- strand): ref base G (15/25); new base A (10/0); total (25/25) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.36e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GGAAAGACCCCGTGAACCTTTACTATAGCTTGACACTGAACATTGAGCCTTGATGTGTAGGATAGGTGGGAGGCTTTGAAGTGTGGACGCCAGTCTGCATGGAGCCGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGACCCGTGATCCGGGTTGCGGACAGTGTCTGGTGGGTAGTTTGACTGGGGCGGTCTCCTCC > U00096/3945759‑3945967 | ggAAAGACCCCGTGAACCTTTACTATAGCTTGACACTGAACATTGAGCCTTGATGTGTAGGATAGGTGGGAGGCTTTGAAGTGTGGACGCCAGTCTGCATGGAGCCGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTgg < 2:2892499/149‑1 (MQ=31) ggAAAGACCCCGTGAACCTTTACTATAGCTTGACACTGAACATTGAGCCTTGATGTGTAGGATAGGTGGGAGGCTTTGAAGTGTGGACGCCAGTCTGCATGGAGCCGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTgg < 2:2892498/149‑1 (MQ=31) gAAAGACCCCGTGAACCTTTACTATAGCTTGACACTGAACATTGAGCCTTGATGTGTAGGATAGGTGGGAGGCTTTGAAGTGTGGACGCCAGTCTGCATGGAGCCGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTgg < 1:2892500/148‑1 (MQ=31) gAAAGACCCCGTGAACCTTTACTATAGCTTGACACTGAACATTGAGCCTTGATGTGTAGAAAACGTTGGACGCATTAAAATTAAAAAAAAAATCTGCATGGAGCCGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGa > 2:2892560/1‑149 (MQ=11) aaaGACCCCGTGAACCTTTACTATAGCTTGACACTGAACATTGAGCCTTGATGTGTAGGATAGGTGGGAGGCTTTGAAGTGTGGACGCCAGTCTGCATGGAGCCGCCCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGAc > 1:2892591/1‑149 (MQ=17) aaGACCCCGTGAACCTTTACTATAGCTTGACACTGAACATTGAGCCTTGATGTGTAGAAAACGTTGGACGCATTAAAATTTTAAAAACAGTCTGCATGGAGCCGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGAcc > 2:2892549/1‑149 (MQ=11) aGACCCCGTGAACCTTTACTATAGCTTGACACTGAACATTGAGCCTTGATGTGTAGGATAGGTGGGAGGCTTTGAAGTGTGGACGCCAGTCTGCATGGAGCCGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGAccc < 2:2892501/149‑1 (MQ=32) aGACCCCGTGAACCTTTACTATAGCTTGACACTGAACATTGAGCCTTGATGTGTAGGATAGGTGGGAGGCTTTGAAGTGTGGACGCCAGTCTGCATGGAGCCGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGAccc > 1:2892561/1‑149 (MQ=32) cccccGTGAACCTTTACTATAGCTTGACACTGAACATTGAGCCTTGATGTGTAGAAAAGGTGGGAGGCATTGAAGTGATGACACCAGTCTGCATGGAGCCGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGACCCGt > 2:2892564/2‑149 (MQ=11) ccccGTGAACCTTTACTATAGCTTGACACTGAACATTGAGCCTTGATGTGTAGGATAGGTGGGAGGCTTTGAAGTGTGGACGCCAGTCTGCATGGAGCCGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGACCCGTg < 1:2892502/149‑1 (MQ=32) cGTGAACCTTTACTATAGCTTGACACTGAACATTGAGCCTTGATGTGTAGGATAGGTGGGAGGCTTTGAAGTGTGGACGCCAGTCTGCATGGAGCCGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGACCCGTGATc < 1:2892504/149‑1 (MQ=34) cGTGAACCTTTACTATAGCTTGACACTGAACATTGAGCCTTGATGTGTAGGATAGGTGGGAGGCTTTGAAGTGTGGACGCCAGTCTGCATGGAGCCGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGACCCGTGATc < 2:2892503/149‑1 (MQ=34) gTGAACCTTTACTATAGCTTGACACTGAACATTGAGCCTTGATGTGTAGAAAACGTTGGACGCATTAAAATTAAAAAAAAAGTCTGCATGGAGACGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGACCCGTGATcc > 2:2892580/1‑149 (MQ=11) tGAACCTTTACTATAGCTTGACACTGAACATTGAGCCTTGATGTGTAGAAAACGTTGGACGCATTAACATTAAAAAAAAAGTCTGCATGGAGCCGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGACCCGTGATCCg > 2:2892610/1‑149 (MQ=11) aaCCTTTACTATAGCTTGACACTGAACATTGAGCCTTGATGTGTAGGATAGGTGGGAGGCTTTGAAGTGTGGACGCCAGTCTGCATGGAGCCGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGACCCGTGATcc < 1:2892505/146‑1 (MQ=34) aCCTTTACTATAGCTTGACACTGAACATTGAGCCTTGATGTGTAGGATAGGTGGGAGGCTTTGAAGTGTGGACGCCAGTCTGCATGGAGCCGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGACCCGTGATCCGGGt < 2:2892506/149‑1 (MQ=34) ccTTTACTATAGCTTGACACTGAACATTGAGCCTTGATGTGTAGGATAGGTGGGAGGCTTTGAAGTGTGGACGCCAGTCTGCATGGAGCCGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGACCCGTGATCCGGGtt < 1:2892508/149‑1 (MQ=34) ccTTTACTATAGCTTGACACTGAACATTGAGCCTTGATGTGTAGGATAGGTGGGAGGCTTTGAAGTGTGGACGCCAGTCTGCATGGAGCCGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGACCCGTGATCCGGGtt < 2:2892507/149‑1 (MQ=34) ttACTATAGCTTGACACTGAACATTGAGCCTTGATGTGTAGAAAACGTTGGACGCATTAAAATTAAAAAAAAACTCTGCATGGAGCCGACCTAGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGACCCGTGATCCGGGTTGCg > 2:2892563/1‑149 (MQ=11) aCTATAGCTTGACACTGAACATTGAGCCTTGATGTGTAGGATAGGTGGGAGGCTTTGAAGTGTGGACGCCAGTCTGCATGGAGCCGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGACCCGTGATCCGGGTTGCGGa < 2:2892509/149‑1 (MQ=34) aCTATAGCTTGACACTGAACATTGAGCCTTGATGTGTAGGATAGGTGGGAGGCTTTGAAGTGTGGACGCCAGTCTGCATGGAGCCGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGACCCGTGATCCGGGTTGCGGa < 2:2892510/149‑1 (MQ=34) tataGCTTGCCACTGAACATTGAGCCTTGATGTGTAGGATAGGTGGGAGGCTTTGAAGTGTGGACGCCAGTCTGCATGGAGCCGACCTTGAAATACCACCATTTAATGTTTGATGTTCTAACGTGGACCCGAGATCCGGGTTGCGGACa > 1:2892530/1‑149 (MQ=255) tataGCTTGACACTGAACATTGAGCCTTGATGTGTAGGATAGGTGGGAGGCTTTGAAGTGTGGACGCCAGTCTGCATGGAGCCGCCCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGACCCGTGATCCGGGTTGCGGACa > 1:2892576/1‑149 (MQ=21) tataGCTTGACACTGAACATTGAGCCTTGATGTGTAGGATAGGTGGGAGGCTTTGAAGTGTGGACGCCAGTCTGCATGGAGCCGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGACCCGTGATCCGGGTTGCGGACa > 1:2892567/1‑149 (MQ=34) tataGCTTGACACTGAACATTGAGCCTTGATGTGTAGAAAACGTTGGACGCATTCAATTGTGAACACCAGTCTGCATGGAGCCGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGACCCGTGATCCGGGTTGCGGACa > 2:2892579/1‑149 (MQ=11) tataGCTTGACACTGAACATTGAGCCTTGATGTGTAGAAAACGTTGGACGCATTAAAATTAAAAAACCAGTCTGCATGGAGCCGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGACCCGTGATCCGGGTTGCGGACa > 2:2892577/1‑149 (MQ=11) ataGCTTGACACTGAACATTGAGCCTTGATGTGTAGAAAACGTTGGACGCATTTAAATTTAAACTACAGTCTGCATGGAGCCGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGACCCGTGATCCGGGTTGCGGACAg > 2:2892574/1‑149 (MQ=11) taGCTTGACACTGTACTTTGAGCCTTGATGTGTAGGATAGGTGGGAGGCTTTGAAGTGTGGAGGCCAGTCGGCATGGAGCCGACCTTGAAATACGTCCCTTTAATGTTTGATGTTCTAACGTGGACCCGTGATCCGGGTTGCGGACAgt < 1:2892511/149‑1 (MQ=11) taGCTTGACACTGAACATTGAGCCTTGATGTGTAGGATAGGTGGGAGGCTTTGAAGTGTGGACGCCAGTATGCATGGAGCCGACCTTGAAATACCCCCCTTTAATGTTTGATGTTCTAACGTGGACCCGTGATCCGGGTTGCGGACAg > 1:2892559/1‑148 (MQ=21) aGCTTGACACTGAACATTGAGCCTTGATGTGTAGGATAGGTGGGAGGCTTTGAAGTGTGGACGCCAGTCTGCATGGAGCCGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGACCCGTGATCCGGGTTGCGGACAgtg < 2:2892513/149‑1 (MQ=34) cTTGACACTGAACATTGAGCCTTGATGTGTAGGATAGGTGGGAGGCTTTGAAGTGTGGACGCCAGTCTGCATGGAGCCGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGACCCGTGATCCGGGTTGCGGACAGTGTc > 1:2892627/1‑149 (MQ=34) ttGACACTGAACATTGAGCCTTGATGTGTAGGATAGGTGGGAGGCTTTGAAGTGTGGACGCCAGTCTGCATGGAGCCGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGACCCGTGATCCGGGTTGCGGACAGTGTCt < 2:2892515/149‑1 (MQ=34) ttGACACTGAACATTGAGCCTTGATGTGTAGGATAGGTGGGAGGCTTTGAAGTGTGGACGCCAGTCTGCATGGAGCCGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGACCCGTGATCCGGGTTGCGGACAGTGTCt < 2:2892514/149‑1 (MQ=34) tGACACTGAACATTGAGCCTTGATGTGTAGAAAACGTTGGACGCATTAAAATTAAAAAAACAGTCTGCATGGAGCCGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGACCCGTGATCCGGGTTGCGGACAGTGTCtg > 2:2892566/1‑149 (MQ=11) cacTGAACATTGAGCCTTGATGTGTAGGATAGGTGGGAGGCTTTGAAGTGTGGACGCCAGTCTGCATGGAGCCGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGACCCGTGATCCGGGTTGCGGACAGTGTCtggtg < 1:2892516/149‑1 (MQ=34) cTGAACATTGAGCCTTGATGTGTAGGATAGGTGGGAGGCTTTGAAGTGTGGACGCCAGTCTGCATGGAGCCGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGACCCGTGATCCGGGTTGCGGACAGTGTCTGGTggg < 1:2892518/149‑1 (MQ=34) tGAACATTGAGCCTTGATGTGTAGGATAGGTGGGAGGCTTTGAAGTGTGGACGCCAGTCTGCATGGAGCCGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGACCCGTGATCCGGGTTGCGGACAGTGTCTGGTGGGt < 2:2892519/149‑1 (MQ=34) aaCATTGAGCCTTGATGTGTAGGATAGGTGGGAGGCTTTGAAGTGTGGACGCCAGTCTGCATGGAGCCGACCCTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGACCCGTGATCCGGGTTGCGGACAGTGTCTGGTGGGTa > 1:2892552/1‑148 (MQ=21) aCATTGAGCCTTGATGTGTAGGATAGGTGGGAGGCTTTGAAGTGTGGACGCCAGTCTGCATGGAGCCGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGACCCGTGATCCGGGTTGCGGACAGTGTCTGGTGGGTAGt > 1:2892541/1‑149 (MQ=34) cATTGAGCCTTGATGTGTAGGATAGGTGGGAGGCTTTGAAGTGTGGACGCCAGTCTGCATGGAGCCGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGACCCGTGATCCGGGTTGCGGACAGTGTCTGGTGGGTAGtt < 2:2892520/149‑1 (MQ=34) aTTGAGCCTTGATGTGTAGGATAGGTGGGAGGCTTTGAAGTGTGGACGCCAGTCTGCATGGAGCCGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGACCCGTGATCCGGGTTGCGGACAGTGTCTGGTGGGTAGttt < 1:2892521/149‑1 (MQ=34) ttGAGCCTTGATGTGTAGGATAGGTGGGAGGCTTTGAAGTGTGGACGCCAGTCTGCATGGAGCCGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGACCCGTGATCCGGGTTGCGGACAGTGTCTGGTGGGTAGTTTg > 1:2892638/1‑149 (MQ=34) ggcgtGATGGCTAGGATAGGTGGGAGGCTTTGAAGTCTTGACGCAGGTCTGCAGGGAGCCGACCTTGAAAAACCCCCCTTTCATGTTTGCAGATCGCACGTGGCCCCGTGATCCGGGTTGAGGACAGTGACTGCTGGGTAGATTGACTg > 1:2892654/5‑149 (MQ=2) ccTTGATGTGTAGGATAGGTGGGAGGCTTTGAAGTGTGGTCGCCAGGCTGCTTGGAGCCGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACTTGGTCCCGTGATCCGGGTTGCGGACAGTGTCTGGTGGGTAGTTTGACTgg < 2:2892522/149‑1 (MQ=11) ccTTGATGTGTAGGATAGGTGGGAGGCTTTGAAGTGTGGACGCCAGTCTGCATGGAGCCGACCTTGAAATACCCCCCTTTAATGTTTGATGTTCTAACGTGGACCCGTGATCCGGGTTGCGGACAGTGTCTGGTGGGTAGTTTGAc > 1:2892592/1‑146 (MQ=21) ccTTGATGTGTAGGATAGGTGGGAGGCTTTGAAGTGTGGACGCCAGTCTGCATGGAGCCGACCTGGAAATACAACCCTTTAATGTTTGATGTTCTAACGTGGACCCGTGATCCGGGTTGCGGACAGTGTCTGGTGGGTAGTTTGACTgg > 1:2892635/1‑149 (MQ=21) cTTGATGTGTAGGATAGGTGGGAGGCTTTGAAGTGTGGACGCCAGTCTGCATGGAGCCGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGACCCGTGATCCGGGTTGCGGACAGTGTCTGGTGGGTAGTTTGACTggg < 1:2892523/149‑1 (MQ=34) cTTGATGTGTAGGATAGGTGGGAGGCTTTGAAGTGTGGACGCCAGTCTGCATGGAGCCGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGACCCGTGATCCGGGTTGCGGACAGTGTCTGGTGGGTAGTTTGACTggg > 1:2892584/1‑149 (MQ=34) tgtgTAGGATAGGTGGGAGGCTTTGAAGTGTGGACGCCAGTCTGCATGGAGCCGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGACCCGTGATCCGGGTTGCGGACAGTGTCTGGTGGGTAGTTTGACTGGGGCGGt > 1:2892588/1‑149 (MQ=34) aGGATAGGGGGGAGGCTTTGAAGTGGGGACGCCAGGCGGCATGGAGCCGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGACCCGTGAGCCGGGTTGCGGACAGTGGCTGGTGGGTAGTTTGACTGGGGCGGtctc < 1:2892525/147‑1 (MQ=11) gATAGGTGGGTGGCTTTGAAGTGTGGACGCCAGTCTGGATGGAGCCTACCTTGAAATACCACTCTTTAATGTTTGTTGTTCTAACGTGGACCCGTGATCCGGGTTGCGGTCAGTGTCTGGTGGGTAGTTTGACTGGGGCGGTCtcctcc < 2:2892526/149‑1 (MQ=11) | GGAAAGACCCCGTGAACCTTTACTATAGCTTGACACTGAACATTGAGCCTTGATGTGTAGGATAGGTGGGAGGCTTTGAAGTGTGGACGCCAGTCTGCATGGAGCCGACCTTGAAATACCACCCTTTAATGTTTGATGTTCTAACGTGGACCCGTGATCCGGGTTGCGGACAGTGTCTGGTGGGTAGTTTGACTGGGGCGGTCTCCTCC > U00096/3945759‑3945967 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |