Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,057,601 C→A 30.8% intergenic (+27/+426) yeeN → / ← asnW UPF0082 family protein/tRNA‑Asn

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,057,6010CA30.8% 21.0 / 8.7 13intergenic (+27/+426)yeeN/asnWUPF0082 family protein/tRNA‑Asn
Reads supporting (aligned to +/- strand):  ref base C (4/5);  new base A (2/2);  total (6/7)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.61e-01

TGATTGCTCAATCTGAAGTTGAGCTTTCCCCAGAAGATTTAGAAATCTTTGAAGGGCTTGTTGATGCCCTTGAAGATGACGACGATGTACAAAAAGTTTATCATAACGTCGCAAATCTCTAATTATCTTTTAAAGAAATCTGTCTTTACGGCAGATTTCTTTAATCTCATATAATTCTTATAAAAAATATAATATTCAACTCGTCATATTGATTATACCCCCCCGTTCCCAGAGAAATAATATTTATTAAAATTCCAGTTCTTCTTTTTC  >  NC_000913/2057453‑2057722
                                                                                                                                                    |                                                                                                                         
tGATTGCTCAATCTGAAGTTGAGCTTTCCCCAGAAGATTTAGAAATCTTTGAAGGGCTTGTTGATGCCCTTGAAGATGACGACGATGTACAAAAAGTTTATCATAACGTCGCAAATCTCTAATTATCTTTTAAAGAAATCTGTCTTTAc                                                                                                                           <  3:245846/149‑1 (MQ=255)
       tCAATCTGAAGTTGAGCTTTCCCCAGAAGATTTAGAAATCTTTGAAGGGCTTGTTGATGCCCTTGAAGATGACGACGATGTACAAAAAGTTTATCATAACGTCGCAAATCTCTAATTATCTTTTAAAGAAATCTGTCTTTAc                                                                                                                           >  3:361034/1‑142 (MQ=255)
       tCAATCTGAAGTTGAGCTTTCCCCAGAAGATTTAGAAATCTTTGAAGGGCTTGTTGATGCCCTTGAAGATGACGACGATGTACAAAAAGTTTATCATAACGTCGCAAATCTCTAATTATCTTTTAAAGAAATCTGTCTTTAc                                                                                                                           <  4:361034/142‑1 (MQ=255)
               aaGTTGAGCTTTCCCCAGAAGATTTAGAAATCTTTGAAGGGCTTGTTGATGCCCTTGAAGATGACGACGATGTACAAAAAGTTTATCATAACGTCGCAAATCTCTAATTATCTTTTAAAGAAATCTGTCTTTACGGCAGATTTCTTTaa                                                                                                            <  4:1713062/149‑1 (MQ=255)
                             ccAGAAGATTTAGAAATCTTTGAAGGGCTTGTTGATGCCCTTGAAGATGACGACGATGTACAAAAAGTTTATCATAACGTCGCAAATCTCTAATTATCTTTTAAAGAAATCTGCCGTAAAGACAGATTTCTTTaaaa                                                                                                          >  1:249554/1‑135 (MQ=38)
                             ccAGAAGATTTAGAAATCTTTGAAGGGCTTGTTGATGCCCTTGAAGATGACGACGATGTACAAAAAGTTTATCATAACGTCGCAAATCTCTAATTATCTTTTAAAGAAATCTGCCGTAAAGACAGATTTCTTTaaaa                                                                                                          <  2:249554/137‑3 (MQ=38)
                                                            ttgATGCCCTTGAAGATGACGACGATGTACAAAAAGTTTATCATAACGTCGCAAATCTCTAATTATCTTTTAAAGAAATCTGTCTTTACGGCAGATTTCTTTaa                                                                                                            >  1:1833878/1‑104 (MQ=255)
                                                            ttgATGCCCTTGAAGATGACGACGATGTACAAAAAGTTTATCATAACGTCGCAAATCTCTAATTATCTTTTAAAGAAATCTGTCTTTACGGCAGATTTCTTTaa                                                                                                            <  2:1833878/104‑1 (MQ=255)
                                                                     ttGAAGATGACGACGATGTACAAAAAGTTTATCATAACGTCGCAAATCTCTAATTATCTTTTAAAGAAATCTGTCTTTACGGCAGATTTCTTTAATCTCata                                                                                                     >  1:694801/1‑102 (MQ=255)
                                                                     ttGAAGATGACGACGATGTACAAAAAGTTTATCATAACGTCGCAAATCTCTAATTATCTTTTAAAGAAATCTGTCTTTACGGCAGATTTCTTTAATCTCata                                                                                                     <  2:694801/102‑1 (MQ=255)
                                                                                      gTACAAAAAGTTTATCATAACGTCGCAAATCTCTAATTATCTTTTAAAGAAATCTGCCGTAAAGACAGATTTCtt                                                                                                               >  3:962370/1‑75 (MQ=25)
                                                                                      gTACAAAAAGTTTATCATAACGTCGCAAATCTCTAATTATCTTTTAAAGAAATCTGCCGTAAAGACAGATTTCtt                                                                                                               <  4:962370/75‑1 (MQ=21)
                                                                                                    tCATAACGTCGCAAATCTCTAATTATCTTTTAAAGAAATCTGTCTTTACGGCAGATTTCtt                                                                                                               <  3:1494786/61‑1 (MQ=255)
                                                                                                    tCATAACGTCGCAAATCTCTAATTATCTTTTAAAGAAATCTGTCTTTACGGCAGATTTCtt                                                                                                               >  4:1494786/1‑61 (MQ=255)
                                                                                                         aCGTCGCAAATCTCTAATTATCTTTTAAAGAAATCTGTCTTTACGGCAGATTTCTTTAATCTCATATAATTCTTATAAAAAATATAATATTCAACTCGTCATATTGATTATACCCCCCCGTTCCCAGAGAAATAATATTTATTAAAAtt                  >  4:1358881/1‑149 (MQ=255)
                                                                                                                         aTTATCTTTTAAAGAAATCTGTCTTTACGGCAGATTTCTTTAATCTCATATAATTCTTATAAAAAATATAATATTCAACTCGTCATATTGATTATACCCCCCCGTTCCCAGAGAAATAATATTTATTAAAATTCCAGTTCTTCTTTTTc  <  3:1358881/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                         
TGATTGCTCAATCTGAAGTTGAGCTTTCCCCAGAAGATTTAGAAATCTTTGAAGGGCTTGTTGATGCCCTTGAAGATGACGACGATGTACAAAAAGTTTATCATAACGTCGCAAATCTCTAATTATCTTTTAAAGAAATCTGTCTTTACGGCAGATTTCTTTAATCTCATATAATTCTTATAAAAAATATAATATTCAACTCGTCATATTGATTATACCCCCCCGTTCCCAGAGAAATAATATTTATTAAAATTCCAGTTCTTCTTTTTC  >  NC_000913/2057453‑2057722

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: