| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 699,722 | 0 | A | G | 32.3% | 31.3 / 19.7 | 31 | intergenic (+115/‑129) | PP_t06/PP_23SD | tRNA‑Ala/23S ribosomal RNA |
| Reads supporting (aligned to +/- strand): ref base A (13/8); new base G (6/4); total (19/12) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
AGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTG‑‑‑GCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGAATTT > NC_002947/699580‑699802 | aGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGATTGTTGATTTCTG‑‑‑GCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGaa > 7:168960/1‑148 (MQ=255)aGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTG‑‑‑GCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGaa > 3:119460/1‑148 (MQ=37) tCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTG‑‑‑GCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCa < 3:259587/149‑1 (MQ=37) cccGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTG‑‑‑GCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCAc > 4:58839/1‑149 (MQ=37) ccGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTG‑‑‑GCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCAct < 2:297508/149‑1 (MQ=37) cGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTG‑‑‑GCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCActg < 2:290214/149‑1 (MQ=37) gCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTG‑‑‑GCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCtggt < 5:152971/148‑1 (MQ=37) gCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTG‑‑‑GCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCtggt < 4:92683/148‑1 (MQ=37) caccacTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTG‑‑‑GCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGaa < 5:319202/125‑1 (MQ=38) caccacTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTG‑‑‑GCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGaa > 6:319202/1‑125 (MQ=38) accacTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTG‑‑‑GCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGAt > 5:361166/1‑147 (MQ=37) accacTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTG‑‑‑GCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGAt > 5:255520/1‑147 (MQ=37) cacTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTTAGCTCGAATGTTGATTTCTG‑‑‑GCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAgg < 2:129740/149‑1 (MQ=37) cTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTGATAGCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGACCACCAGTTTCAGTGATGGTGGATAAg > 1:125454/1‑149 (MQ=11) aGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTG‑‑‑GCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTtgt > 2:185668/1‑149 (MQ=37) ttACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTG‑‑‑GCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGa < 6:361166/149‑1 (MQ=25) tACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTG‑‑‑GCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGa < 4:217327/148‑1 (MQ=25) ttGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTG‑‑‑GCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATATGATGAACACCAGTTTCAATGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTc > 7:214696/1‑149 (MQ=255) aGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTG‑‑‑GCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTc > 1:404445/1‑149 (MQ=21) aGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTG‑‑‑GCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTc > 1:357288/1‑149 (MQ=34) gAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTG‑‑‑GCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTc > 8:29523/1‑148 (MQ=34) aGAAATGAACATTCAGGTCGAATGTTGATTTATG‑‑‑GCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGTTAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTATGGTAAATTTTGTGAGTTATGATCGaat > 8:66544/1‑145 (MQ=11) aGAAATGAACATTCAGCTCGAATGTTGATTTCTG‑‑‑GCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGaaa < 7:66544/146‑1 (MQ=34) aGAAATGAACATTCAGCTCGAATGTTGATTTCTG‑‑‑GCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAagag > 6:106513/1‑149 (MQ=34) gAACATTCAGCTCGAATGTTGATTTCTG‑‑‑GCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACg > 5:287740/1‑148 (MQ=34) atcATTCAGCTCGAATGTTGATTTCTG‑‑‑GCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGTCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTa < 2:83266/147‑1 (MQ=21) aaCATTCAGCTCGAATGTTGATTTCTG‑‑‑GCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGa > 7:40331/1‑126 (MQ=25) aaCATTCAGCTCGAATGTTGATTTCTG‑‑‑GCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGa < 8:40331/126‑1 (MQ=25) aTTCAGCTCGAATGTTGATTTCTG‑‑‑GCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACg > 2:154595/1‑148 (MQ=34) tCAGCTCGAATGTTGATTTCTG‑‑‑GCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGGCTGAAAACCAGTTTCACTGCTGTTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGGATCGTACGAAt > 6:262291/1‑149 (MQ=11) aGCTCGAATGTTGATTTCTG‑‑‑GCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGAAttt > 2:327824/1‑149 (MQ=32) | AGCGGTTCGATCCCGCTTGGCTCCACCACTTTCGCCGTACGCAGTTACTTGTCAGAACTTAGAAATGAACATTCAGCTCGAATGTTGATTTCTG‑‑‑GCTTTTGTCAGATCGTTCTTTAAAAATTCGGATATGTGATAGATATAGACTGAACACCAGTTTCACTGCTGGTGGATCAGGCTAAGGTAAAATTTGTGAGTTCTGCTCGAAAGAGCAACGTACGAATTT > NC_002947/699580‑699802 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |