| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_000913 | 2,072,022 | (G)5→4 | coding (484/3120 nt) | flu → | CP4‑44 prophage; self recognizing antigen 43 (Ag43) autotransporter |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 2,072,018 | 0 | G | . | 96.7% | 119.4 / ‑0.7 | 30 | coding (480/3120 nt) | flu | CP4‑44 prophage; self recognizing antigen 43 (Ag43) autotransporter |
| Reads supporting (aligned to +/- strand): ref base G (0/1); new base . (14/15); total (14/16) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01 | |||||||||||
GGAGGCGGACAGAGCCTTCAGGGACGGGCTGTGAACACCACGCTGAATGGTGGCGAACAGTGGATGCATGAGGGGGCGATAGCCACAGGAACCGTCATTAATGATAAGGGCTGGCAGGTCGTCAAGCCCGGTACAGTGGCAACGGAT > NC_000913/2071947‑2072093 | ggAGGCGGACAGAGCCTTCAGGGACGGGCTGTGAACACCACGCTGAATGGTGTCGAACAGTGGATGCATGA‑GGGGCGATAg > 1:362218/1‑81 (MQ=255) gAGGCGGACAGAGCCTTCAGGGACGGGCTGTGAACACCACGCTGAATGGTGGCGAACAGTGGATGCATGA‑GGGGCGATAGc < 2:362168/81‑1 (MQ=255) agagCCTTCAGGGACGGGCTGTGAACACCACGCTGAATGGTGGCGAACAGTGGATGCATGA‑GGGGCGATAGCCACAGGAAc < 1:362169/81‑1 (MQ=255) gagCCTTCAGGGACGGGCTGTGAACACCACGCTGAATGGTGGCGAACAGTGGATGCATGA‑GGGGCGATAGCCACAGGAAcc < 2:362170/81‑1 (MQ=255) agCCTTCAGGGACGGGCTGTGAACACCACGCTGAATGGTGGCGAACAGTGGATGCATGA‑GGGGCGATAGCCACAGGAACCg < 2:362171/81‑1 (MQ=255) cAGGGACGGGCTGTGAACACCACGCTGAATGGTGGCGAACAGTGGATGCATGA‑GGGGCGATAGCCACAGGAACCGTCAtt > 2:362213/1‑80 (MQ=255) aGGGACGGGCTGTGAACACCACGCTGAATGGTGGCGAACAGTGGATGCATGAGGGGGCGATAGCCACAGGAACCGTCAtt < 2:362172/80‑1 (MQ=255) gACGGGCTGTGAACACCACGCTGAATGGTGGCGAACAGTGGATGCATGA‑GGGGCGATAGCCACAGGAACCGTCATTAATGa > 2:362217/1‑81 (MQ=255) cGGGCTGTGAACACCACGCTGAATGGTGGCGAACAGTGGATGCATGA‑GGGGCGATAGCCACAGGAACCGTCATTAATGAt > 1:362226/1‑80 (MQ=255) ggCTGTGAACACCACGCTGAATGGTGGCGAACAGTGGATGCATGA‑GGGGCGATAGCCACAGGAACCGTCATTAATGATAAg > 2:362212/1‑81 (MQ=255) cTGTGAACACCACGCTGAATGGTGGCGAACAGTGGATGCATGA‑GGGGCGATAGCCACAGGAACCGTCATTAATGATAAggg > 2:362224/1‑81 (MQ=255) tgtgAACACCACGCTGAATGGTGGCGAACAGTGGATGCATGA‑GGGGCGATAGCCACAGGAACCGTCATTAATGATAAGGGc > 2:362214/1‑81 (MQ=255) tgAACACCACGCTGAATGGTGGCGAACAGTGGATGCATGA‑GGGGCGATAGCCACAGGAACCGTCATTAATGATAAGGGCt < 2:362173/80‑1 (MQ=255) gAACACCACGCTGAATGGTGGCGAACAGTGGATGCATGA‑GGGGCGATAGCCACAGGAACCGTCATTACTGATAAGGGCTgg > 1:362210/1‑81 (MQ=255) aaCACCACGCTGAATGGTGGCGAACAGTGGATGCATGA‑GGGGCGATAGCCACAGGAACCGTCATTAATGATAAGGGCTgg < 2:362174/80‑1 (MQ=255) accacGCTGAATGGTGGCGAACAGTGGATGCATGA‑GGGGCGATAGCCACAGGAACCGTCATTAATGATAAGGGCTGGCAgg > 2:362228/1‑81 (MQ=255) gCTGAATGGTGGCGAACAGTGGATGCATGA‑GGGGCGATAGCCACAGGAACCGTCATTAATGATAAGGGCTGGCAGgtcgtc < 2:362176/81‑1 (MQ=255) gCTGAATGGTGGCGAACAGTGGATGCATGA‑GGGGCGATAGCCACAGGAACCGTCATTAATGATAAGGGCTGGCAGgtcgtc < 2:362177/81‑1 (MQ=255) gCTGAATGGTGGCGAACAGTGGATGCATGA‑GGGGCGATAGCCACAGGAACCGTCATTAATGATAAGGGCTGGCAGgtcgtc < 1:362175/81‑1 (MQ=255) aaTGGTGGCGAACAGTGGATGCATGA‑GGGGCGATAGCCACAGGAACCGTCATTAATGATAAGGGCTGGCAGGTCGTCAAg < 2:362178/80‑1 (MQ=255) aTGGTGGCGAACAGTGGATGCATGA‑GGGGCGATAGCTACAGGAACCGTCATTAATGATAAGGGCTGGCAGGTCGTCAAGcc > 1:362225/1‑81 (MQ=255) tggCGAACAGTGGATGCATGA‑GGGGCGATAGCCACAGGAACCGTCATTAATGATAAGGGCTGGCAGGTCGTCAAGCCCGGt < 2:362179/81‑1 (MQ=255) gAACAGTGGATGCATGA‑GGGGCGATAGCCACAGGAACCGTCATTAATGATAAGGGCTGGCAGGTCGTCAAGCCCGGTACAg < 2:362180/81‑1 (MQ=255) gAACAGTGGATGCATGA‑GGGGCGATAGCCACAGGAACCGTCATTAATGATAAGGGCTGGCAGGTCGTCAAGCCCGGTACAg > 2:362221/1‑81 (MQ=255) gTGGATGCATGA‑GGGGCGATAGCCACAGGAACCGTCATTAATGATAAGGGCTGGCAGGTCGTCAAGCCCGGTACAGTGGc > 1:362223/1‑80 (MQ=255) gATGCATGA‑GGGGCGATAGCCACAGGAACCGTCATTAATGATAAGGGCTGGCAGGTCGTCAAGCCCGGTACAGTGGCAAc < 2:362181/80‑1 (MQ=255) gATGCATGA‑GGGGCGATAGCCACAGGAACCGTCATTAATGATAAGGGCTGGCAGGTCGTCAAGCCCGGTACAGTGGCAAc < 1:362182/80‑1 (MQ=255) gATGCATGA‑GGGGCGATAGCCACAGGAACCGTCATTAATGATAAGGGCTGGCAGGTCGTCAAGCCCGGTACAGTGGCAACg > 2:362222/1‑81 (MQ=255) gCATGA‑GGGGCGATAGCCACAGGAACCGTCATTAATGATAAGGGCTGGCAGGTAGTCAAGCCCGGTACAGTGGCAACGGa > 1:362215/1‑80 (MQ=255) cATGA‑GGGGCGATAGCCACAGGAACCGTCATTAATGATAAGGGCTGGCAGGTCGTCAAGCCCGGTACAGTGGCAACGGAt < 2:362183/80‑1 (MQ=255) | GGAGGCGGACAGAGCCTTCAGGGACGGGCTGTGAACACCACGCTGAATGGTGGCGAACAGTGGATGCATGAGGGGGCGATAGCCACAGGAACCGTCATTAATGATAAGGGCTGGCAGGTCGTCAAGCCCGGTACAGTGGCAACGGAT > NC_000913/2071947‑2072093 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |