| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description | 
| RA | NC_000913 | 2,079,004 | G→C | 19.0% | intergenic (+73/+28) | yoeF → / ← yeeX | pseudogene, CP4‑44 putative prophage remnant;Phage or Prophage Related/UPF0265 family protein | 
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 2,079,004 | 0 | G | C | 19.0% | 47.2 / 7.8 | 21 | intergenic (+73/+28) | yoeF/yeeX | pseudogene, CP4‑44 putative prophage remnant;Phage or Prophage Related/UPF0265 family protein | 
| Reads supporting (aligned to +/- strand): ref base G (9/8); new base C (2/2); total (11/10) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.35e-01 | |||||||||||
ATAACAATGGAGTTACAGAACTTTTTGCCTTGCTGCATACAGAAGAAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTTGATGATGTAATCGTCAACGCGATCTTCATAG  >  NC_000913/2078873‑2079149                                                                                                                                   |                                                                                                                                                 ataaCAATGGAGTTACAGAACTTTTTGCCTTGCTGCATACAGAAGAAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATTAGCAATGAGGGTTGCCGGGCa                                                                                                                                               >  4:550843/1‑136 (MQ=255)ataaCAATGGAGTTACAGAACTTTTTGACTTGCTGCATACAGAAGAAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCa                                                                                                                                               <  3:550843/136‑1 (MQ=255)                 gAACTTTTTGCCTTGCTGCATACAGAAGAAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCaa                                                                                                                                              <  4:190783/120‑1 (MQ=255)                 gAACTTTTTGCCTTGCTGCATACAGAAGAAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCaa                                                                                                                                              >  3:190783/1‑120 (MQ=255)                                      aCAGAAGAAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCCGGCAACCCTCAt                                                                                                                                       >  3:301717/1‑106 (MQ=255)                                      aCAGAAGAAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCCGGCAACCCTCAt                                                                                                                                       <  4:301717/106‑1 (MQ=255)                                        agaagaAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCg                                                                                                                                                   <  2:210055/92‑1 (MQ=255)                                        agaagaAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCg                                                                                                                                                   >  1:210055/1‑92 (MQ=255)                                                                      tctcGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAAc                                                                                                                                             <  3:138195/68‑1 (MQ=255)                                                                      tctcGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAAc                                                                                                                                             >  4:138195/1‑68 (MQ=255)                                                                      tctcGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCCGGCAACCCTCa                                                                                                                                        <  3:1119183/73‑1 (MQ=255)                                                                      tctcGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCCGGCAACCCTCa                                                                                                                                        >  4:1119183/1‑73 (MQ=255)                                                                                                 cacaGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCaa                                                                                                                                              <  1:127832/40‑1 (MQ=255)                                                                                                 cacaGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCaa                                                                                                                                              >  2:127832/1‑40 (MQ=255)                                                                                                      aTAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCACAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTtgatg                             >  3:1080149/1‑148 (MQ=255)                                                                                                             tCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCATTCTTTGGAGAGCTCGGCATTTTTGATGATGTAAt                      <  4:1080149/148‑1 (MQ=255)                                                                                                                             gTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTgg                                                                      >  3:524156/1‑84 (MQ=255)                                                                                                                             gTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTgg                                                                      <  4:524156/84‑1 (MQ=255)                                                                                                                               tGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCtttt                                                                         >  3:1011725/1‑79 (MQ=255)                                                                                                                               tGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCtttt                                                                         <  4:1011725/79‑1 (MQ=255)                                                                                                                                 ccGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTTGATGATGTAATCGTCAACGCGATCTTCATAg  >  3:919828/1‑148 (MQ=255)                                                                                                                                  cGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCAt                                                                                       >  2:172824/1‑62 (MQ=255)                                                                                                                                  cGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCAt                                                                                       <  1:172824/62‑1 (MQ=255)                                                                                                                                   |                                                                                                                                                 ATAACAATGGAGTTACAGAACTTTTTGCCTTGCTGCATACAGAAGAAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTTGATGATGTAATCGTCAACGCGATCTTCATAG  >  NC_000913/2078873‑2079149 | 
| Alignment Legend | 
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG | 
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑ |