| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_000913 | 1,981,712 | C→A | 15.9% | A226S (GCA→TCA) | otsB ← | trehalose‑6‑phosphate phosphatase, biosynthetic |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,981,712 | 0 | C | A | 15.9% | 170.5 / 21.8 | 63 | A226S (GCA→TCA) | otsB | trehalose‑6‑phosphate phosphatase, biosynthetic |
| Reads supporting (aligned to +/- strand): ref base C (27/26); new base A (7/3); total (34/29) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.19e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.46e-01 | |||||||||||
GTCATCACTCCTGTTATTTTCTCTTTTTTGTTGTAATGCGGTGGTTATCATTTCAAGCCAGCTCCAGACATCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAA‑CGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGG > NC_000913/1981614‑1981811 | gTCATCACTCCTGTTATTTTCTCTTTTTTGTTGTAATGCGGTGGTTATCATTTCAAGCCAGCTCCAGACATCCGGCACACCCGCCAGTCGCCATGATGcc > 2:587745/1‑100 (MQ=255) catcaCTCCTGTTATTTTCTCTTTTTTGTTGTAATGCGGTGGTTATCATTTCAAGCCAGCTCCAGACATCCGGCACACCCGCCAGTCGCCATGATGCCTg < 1:1162356/100‑1 (MQ=255) caCTCCTGTTATTTTCTCTTTTTTGTTGTAATGCGGTGGTTATCATTTCAAGCCAGCTCCAGACATCCGGCACACCCGCCAGTCGCCATGATGCCTGAGt < 1:1541630/100‑1 (MQ=255) caCTCCTGTTATTTTCTCTTTTTTGTTGTAATGCGGTGGTTATCATTTCAAGCCAGCTCCAGACATCCGGCACACCCGCCAGTCGCCATGATGACTGAGt > 1:2241432/1‑100 (MQ=255) aCTCCTGTTATTTTCTCTTTTTTGTTGTAATGCGGTGGTTATCATTTCAAGCCAGCTCCAGACACCCGGCACACCCGCCAGTCGCCATGATGCCTGAGtt > 1:1930532/1‑100 (MQ=255) cTCCTGTTATTTTCTCTTTTTTGTTGTAATGCGGTGGTTATCATTTCAAGCCAGCTCCAGACATCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTg > 1:1377843/1‑100 (MQ=255) ccTGTTATTTTCTCTTTTTTGTTGTAATGCGGTGGTTATCATTTCAAGCCAGCTCCAGACATCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCa > 1:936308/1‑100 (MQ=255) ccTGTTATTTTCTCTTTTTTGTTGTAATGCGGTGGTTATCATTTCAAGCCAGCTCCAGACATCCGGCACACCCGCCAGTCGCCATGATGACTGAGTTGCa > 1:1546427/1‑100 (MQ=255) tGTTATTTTCTCTTTTTTATTGTAATGCGGTGGTTATCATTTCAAGCCAGATCCAGACATCCGGAACACCCGCCAGTCGCCATGATGACTGAGTTGCAcc > 2:876666/1‑100 (MQ=255) gTTATTTTCTCTTTTTTGTTGTAATGCGGTGGTTATCATTTCAAGCCAGCTCCAGACATCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCt > 2:1131360/1‑100 (MQ=255) gTTATTTTCTCTTTTTTGTTGTAATGCGGTGGTTATCATTTCAAGCCAGCTCCAGACATCCGGCACACCCGCCAGTCGCCATGATGACTGAGTTGCACCt > 2:558235/1‑100 (MQ=255) tATTTTCTCTTTTTTGTTGTAATGCGGTGGTTATCATTTCAAGCCAGCTCCAGACATCCTGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCtgt > 2:264986/1‑100 (MQ=255) aTTTTCTCTTTTTTGTTGTAATGCGGTGGTTATCATTTCAAGCCAGCTCCAGACATCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCtgtg > 2:1318471/1‑100 (MQ=255) tttttGTTGTAATGCGGTGGTTATCATTTCAAGCCAGCTCCAGACATCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAAttttt < 1:375583/100‑1 (MQ=255) ttgttgTAATGCGGTGGTTATCATTTCAAGCCAGCTCCAGACATCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACt < 1:1131370/100‑1 (MQ=255) ttgctgTAATGCGGTGGTTATCATTTCAAGCCAGCTCCAGACATCCGGCACACCCGCCAGTCGTCATGATGCTTGAGTTGCACCTGTGCCAATTTTTACt < 2:211363/100‑1 (MQ=255) tgttgtAATGCGGTGGTTATCATTTCAAGCCAGCTCCAGACATCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTg < 1:253472/100‑1 (MQ=255) gttgtAATGCGGTGGTTATCATTTCAAGCCAGCTCCAGACATCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCTCCTGTGCCAATTTTTACTGa < 1:2260434/100‑1 (MQ=255) tgtAATGCGGTGGTTATCATTTCAAGCCAGCTCCAGACATCCGGCACACCCGCCAGTCGCCATGATGACTGAGTTGCACCTGTGCCAATTTTTACTGACa > 2:538574/1‑100 (MQ=255) gtAATGCGGTGGTTATCATTTCAAGCCAGCTCCAGACATCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCTCCAGTGCCAATTTTTACTGACAt < 1:327654/100‑1 (MQ=255) gCGGTGGTTATCATTTCAAGCCAGCTCCAGACATCCGGCACACCCGCCAGTCGCCATGATGACTGAGTTGCACCTGTGCTAATTTTTACTGACATTCCGc < 2:1131162/100‑1 (MQ=255) cGGTGGTTATCATTTCAAGCCAGCTCCAGACATCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTcc > 1:1147086/1‑97 (MQ=255) cGGTGGTTATCATTTCAAGCCAGCTCCAGACATCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGc < 1:2269622/99‑1 (MQ=255) tgtttATCATTTCAAGCCAGCTCCAGACATCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAg < 2:930319/97‑1 (MQ=255) tggtTATCATTTCAAGCCAGCTCCAGACATCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGccc < 1:216413/98‑1 (MQ=255) tggtTATCATTTCAAGCCAGCTCCAGACATCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAg < 1:380046/100‑1 (MQ=255) gtTATCATTTCAAGCCAGCTCCAGACATCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTc > 1:548488/1‑100 (MQ=255) aTCATTTCAAGCCAGCTCCAGACATCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGt > 2:670738/1‑100 (MQ=255) tCATTTCAAGCCAGCTCCAGACATCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGtt < 1:203395/100‑1 (MQ=255) tCATTTCAAGCCAGCTCCAGACATCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGtt > 2:1861203/1‑100 (MQ=255) tCATTTCAAGCCAGCTCCAGACATCCGGCACACCCGCCAGTCGCCATGATGACTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGtt < 1:1502108/100‑1 (MQ=255) cATTTCAAGCCAGCTCCAGACATCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTa > 2:871897/1‑100 (MQ=255) tttCAAGCCAGCTCCAGACATCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAA‑c < 1:392288/100‑1 (MQ=255) tttCAAGCCAGCTCCAGACATCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAA‑c < 2:548485/100‑1 (MQ=255) tttCAAGCCAGCTCCAGACATCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAA‑c < 2:96791/100‑1 (MQ=255) tCAAGCCAGCTCCAGACATCCGGCACACCCGCCAGTCGCCATGATGACTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAA‑CGa > 2:1094036/1‑100 (MQ=255) cAAGCCAGCTCCAGACATCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAA‑CGAc < 1:1179804/100‑1 (MQ=255) aGCCAGCTCCAGACATCCGTCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAA‑CGACTg < 2:172343/100‑1 (MQ=255) gCCAGCTCCAGACATCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAA‑CGACTGc < 2:480118/100‑1 (MQ=255) aGCTCCAGACATCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTC‑GTTAACCGACTgggcg < 1:1875077/100‑5 (MQ=255) cTCCAGACATCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAA‑CGACTGCGAAGc > 2:1796410/1‑100 (MQ=255) tCCAGACATCCGGCACACCCGCCAGTTGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAA‑CGACTGCGAAGcc < 1:1027205/100‑1 (MQ=255) tCCAGACATCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAA‑CGACTGCGAAGcc > 1:117348/1‑100 (MQ=255) ccAGACATCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAA‑CGACTGCGAAGcc > 1:1035547/1‑99 (MQ=255) ccAGACATCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAA‑CGACTGCGAAGCCa > 2:2346107/1‑100 (MQ=255) cAGACATCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAA‑CGACTGCGAAGCCAg > 2:993382/1‑100 (MQ=255) cAGACATCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAA‑CGACTGCGAAGCCAg > 2:960402/1‑100 (MQ=255) gACATCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAA‑CGACTACGAAGCCAGa > 1:532099/1‑99 (MQ=255) cATCCGGCACACCCGCCAGTCGCCATGATGACTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAA‑CGACTGCGAAGCCAGATTc > 1:568929/1‑100 (MQ=255) ccGGCACGCCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAA‑CGACTGCGAAGCCAGATTCATc < 2:1561595/100‑1 (MQ=255) cGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAA‑CGACTGCGAAGCCAGATTCATCg > 2:1451407/1‑100 (MQ=255) cGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAA‑CGACTGCGAAGCCAGATTCATCg > 1:2237344/1‑100 (MQ=255) gCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTATTGACATTCCGCCCAGTCGGTTAA‑CGACTGCGAAGCCAGATTCATCGGt < 2:77369/100‑1 (MQ=255) cacCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAA‑CGACTGCGAAGCCAGATTCATCGGTTaa > 2:2129601/1‑100 (MQ=255) acCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAA‑CGACTGCGAAGCCAGATTCATCGGTTaaa < 2:1021580/100‑1 (MQ=255) cccGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAA‑CGACTGCGAAGCCAGATTCATCGGTTAAat < 2:434435/100‑1 (MQ=255) ccGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAA‑CGACTGCGAAGCCAGATTCATCGGTTaaa < 1:757679/98‑1 (MQ=255) ccGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAA‑CGACTGCGAAGCCAGATTCATCGGTTAAatc > 2:1509838/1‑100 (MQ=255) ccAGTCGCCATGATGACTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAA‑CGACTGCGAAGCCAGATTCATCGGTTAAatca < 1:538577/98‑1 (MQ=255) cGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAA‑CGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCa > 2:2217407/1‑100 (MQ=255) gCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTTACATTCCGCCCAGTAGGTTAA‑CGACTGCGAAGCCAAATTAATCGGTTAAATAATAGccc > 2:427290/1‑98 (MQ=255) gCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAA‑CGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAg > 1:124531/1‑100 (MQ=255) ccATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAA‑CGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGa > 1:1209811/1‑100 (MQ=255) tgCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAA‑CGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATAc > 1:847387/1‑99 (MQ=255) ccTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAA‑CGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACggg < 2:47740/100‑1 (MQ=255) ccTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAA‑CGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACggg > 1:420163/1‑100 (MQ=255) | GTCATCACTCCTGTTATTTTCTCTTTTTTGTTGTAATGCGGTGGTTATCATTTCAAGCCAGCTCCAGACATCCGGCACACCCGCCAGTCGCCATGATGCCTGAGTTGCACCTGTGCCAATTTTTACTGACATTCCGCCCAGTCGGTTAA‑CGACTGCGAAGCCAGATTCATCGGTTAAATCATCGCCCAGAAATACGGG > NC_000913/1981614‑1981811 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |