Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,510,609 C→T P327L (CCG→CTG)  yfeO → putative transport protein YfeO

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,510,6090CT95.0% 46.8 / ‑6.6 20P327L (CCG→CTG) yfeOputative transport protein YfeO
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (17/2);  minor base G (1/0);  total (18/2)
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.95e-01

ATCCGGTGCTCGTGCTGGGTATTGGCGGATTTATTCTCGGTATTCTGGGGGTTATTGGTGGACCAGTTTCGCTGTTTAAAGGGCTGGATGAGATGCAGCAGATGGTGGCAAATCAGGCTTTCAGCACCAGCGATTACTTTTTGCTGGCGGTAATTAAACTTGCCGCCCTGGTCGTTGCTGCCGCCAGTGGCTTTCGCGGTGGGCGAATCTTCCCGGCAGTGTTTGTCGGCGTGGCATTAGGGTTGATGCTGCATGAGCACGTTCCCGCCGTACCAGCGGCAATAACCGTTTCTTGCGCTATTCTCGGCATCGTGCTGGTGGTAACACGCGATGGCTGGTTAAGTCTTTTTATGGCGGCAGTCGTTGTACCCAATACCACATTGCTACCGCTGCTCTGTATCGTCATGCTTCCGGCATGGCTGTTATTAGCAGGTAA  >  NC_000913/2510396‑2510831
                                                                                                                                                                                                                      |                                                                                                                                                                                                                                  
aTCCGGTGCTCGTGCTGGGTATTGGCGGATTTATTCTCGGTATTCTGGGGGTTATTGGTGGACCAGTTTCGCTGTTTAAAGGGCTGGGTGAGATGCAGCAGATGGTGGCAAATCAGGCTTTCAGCACCAGCGATTACTTTTTGCTGGCGGTAATTAAACTTTGCGCCCTTGGCGTTGCTGCCGCCAGGGGCTTTTGCGGAGGGGGCATCTTCCTGGGAGTGTTTGTCGGCGTggggtttggggtggtggtt                                                                                                                                                                                               >  1:145626/1‑234 (MQ=255)
  ccGGTGCTCGTGCTGGGTATTGGCGGATTTATTCTCGGTATTCTGGGGGTTAGTGGTGGACCAGTTTCGCTGTTTAAAGGGCTGGATGAGATGCAGCAGATGGTGGCAAATCAGGCTTTCAGCACCAGCGATTACTTTTTGCTGGCGGTAATTAAACTTGCCGCCCTGGTCGTTGCTGCCGCCAGTGGCTTTCGCGGTGGGCGAATCTTCCTGGCAGTGTTTGTCGGCGTGGCATTAGGGTTGATGCTGCa                                                                                                                                                                                           <  1:526621/251‑1 (MQ=255)
                                      ggTATTCTGGGGGTTATTGGTGGACCAGTTTCGCTGTTTAAAGGGCTGGATGAGATGCAGCAGATGGTGGCAAATCAGGCTTTCAGCACCAGCGATTACTTTTTGCTGGCGGTAATTAAACTTGCCGACCTGGTCGTTGCTGCCGCCAGTGGCTTTCGCGGGGGGCGAAACTTTCTGGCCGGGTTTGTCGGCGTGGCATTAGGGTTGATGCTGAATGAGGACGTTCCCGCCGTAACAGCCGCAAAAACCgg                                                                                                                                                     >  1:84217/1‑250 (MQ=255)
                                                    tattggaagACCAGTTCCGCTGTTTAAAGGGATGGATGAGCTGAAGCAGATGGTGGCAAATCAGGCTTTCAGCACCAGCGATTACTTTTTGCTGGCGGTTATTAAACTTGACGCCCTGGGCGTTGCTGCCCCCAGGTGCGTTCGCGGCGGGGGAAACTTCCTGGGCGGTGTTGTCGGGGTGGACTTTGGGGTTGTTGTGCCTGGGTCAGTTCCCCTCGCACCAGCGGCAATGAACGTTTCTTGCGCTAttt                                                                                                                                       >  2:508028/9‑250 (MQ=255)
                                                          tggACCAGTTTCGCTGTTTAAAGGGCTGGATGAGATGCAGCAGATGGCGGCAAATCAGGCTTTAAGCACCAGCGCTTACTTTTTGCTGGCGGTAATTAAAAATGCCGCCCTGGTCGTTGCTGACGCACGTGGCTTTCGCGGGGGGGGAATCCTCCTGGCAGAGTTTGACCGGATGGGCTTAGGGTTGATGCTAACTGGGCAGGTTTCCGCCGGACCAGCAGCAATACCCGGTTCTTGCGCTCTTTTCGGc                                                                                                                                  >  1:57736/1‑250 (MQ=255)
                                                                        tGTTTAAAGGGCTGGATGAGATGCAGCAGATGGTGGCAAATCAGGCTTTCAGCACCAGCGATTACTTTTTGCTGGCGGTAATTAAACTTGCCGCCCTGGTCGTTGCTGCCGCCAGTGGCTTTCGCGGTGGGCGAATCTTCCTGGCAGTGTTTGTCGGCGTGGCCTTAGGGTTGATGCTGCATGAGCACGTTCCCGCCGTACCAGCGGCAATAACCGTTTCTTGCGCTATTCTCGGCATCGTGCTggtggt                                                                                                                    >  1:137637/1‑250 (MQ=255)
                                                                                   cTGGATGAGATGCAGCAGATGGTGGCAAATCAGGCTTTCAGCACCAGCGATTACTTTTTGCTGGCGGTAATTAAAATTGCCGCCCTGGTCGTTGCTGCCGCCAGTGGCTTTTGCGGTGGGCGCATCTTCCTGGCAGTGTTTGTCGGCGTGGCATTAGGGGTGCTGCTGCCTGAGCACGTTCCCGCCGTACCCGCGGCCATAACCGTTTCTTGCGCTATTTTCGGCATCGTGCTGGTGGGTACACGCGAttg                                                                                                        >  1:516053/1‑249 (MQ=255)
                                                                                         gagaTGCAGCAGATGGAGGCAAATCAGGCTTTCAGCACCAGCGATTACTTTTTGCTGGCGGTAATTAAACTTGCCGCCCTGGTCGTTGCTGCCGCCCATGGCTTTCGCGGGGGGAGGATCTTCCTGGAAGGGTTTGTCGGCGGGGCATTTAGGGTGGAGCTGCATGAGCACGGTCCCGCCGTACCAGAGGCAATAACACTTTCTTGCGCTATTCTCGGCATCGTGGTGGTGGGAACAAGCGgtgggtgggt                                                                                                  >  1:149263/1‑249 (MQ=255)
                                                                                               cagcagATGGTGGCAAATCAGGCTTTCAGCACCAGCGATTACTTTTTGCAGGCGGTAATTAAACTTTCCGCCCTTGTCGTTGCTGACTCCAGTGGCTTTTGCGGGGGGCGAATCTTCCTGGGAGTTTTTGTCGGCGTGGCCTTAGGGTTGATGCTGCCAGAGCCCCGTCCCGCCGTACCCACGGCCAATACCGGTTTTTGTGCTATATTCGGCACCGTGCTGGAGGTaaacagcgccggttgggttaaatc                                                                                            >  2:471427/1‑229 (MQ=255)
                                                                                                       ggtggCAAATCAGGCTTTCAGCACCAGCGATTACTTTTTGCTGGCGGTAATTAAACTTGCCGCCCTGGTCGTTGCTGCCGCCAGTGGCTTTCGCGGTGGGCGAATCTTCCTGGCAGTGTTTGTCGGCGTGGCATTAGGGTTGATGCTGCATGAGCACGTTCCCGCCGTACCAGCGGCAATAACCGTTTCTTGCGCTATTCTCGGCATCGTGCTGGTGGTAACACGCGATGGCTGGTTAAGTCTTTTTATgg                                                                                    >  1:233828/1‑251 (MQ=255)
                                                                                                                                 gCGATTACTTTTTGCTGGCGGTAATTAAACTTGCCGCCCTGGTCGTTGCTGCCGCCAGTGGCTTTCGCGGTGGGCGAATCTTCCTGGCAGTGTTTGTCGGCGTGGCATTAGGGTTGATGCTGCCTGAGCACGTTCCCGCCGTACCAGCGGCAATAACCGTTTCTTGCGCTATTCTCGGCATCGTGCTGGTGGTAACACGCGGTGGCTGGTTAAGTCTTTTTTTggcgggcggtcgtgttcccaaaaccccc                                                          >  1:67190/1‑228 (MQ=255)
                                                                                                                                                     gTAATTAAACTTGCCGCCCTGGTCGTTGCTGCCGCCAGTGGCTTTCGCGGTGGGCGAATCTTCCTGGCAGTGTTTGTCGGCGTGGCATTAGGGTTGATGCTGCATGAGCACGTTCCCGCCGTACCAGCGGCAATAACCGTTTCTTGCGCTATTCTCGGCATCGTGCTGGCGGGAACACGCGATGGCTGGTTAAGCCTTTTTATGGCGGCAGTGGTGGTATCCACAACCACCTTTCTACCGCTGctctcgtt                                      >  1:209933/1‑247 (MQ=255)
                                                                                                                                                          tAAACTTGCCGCCCTGGTCGTTGCTGCAGCCAGTGGCTTTCGCGGTGGGCGAATCTTCCTGGCAGTGTTTGTCGGCGTGGCATTAGGGTTGATGCTGCATGAGCACGTTCCCGCCGTACCAGCGGCAATAACCGTTTCTTGCGCTATTCTCGGCATCGTGCTGGTGGTAACACGCGATGGCTGGTTAAGTCTTTTTATGGCGGCCGGCGTTTTACCCCAAACCACCTTTCTAcccccgtttctgttcgtac                                 >  1:105062/1‑234 (MQ=255)
                                                                                                                                                                   cgccCTGGTCGTTGCTGCCGCCAGTGGCTTTCGCGGTGGGCGAATCTTCCTGGCAGTGTTTGTCGGCGTGGCATTAGGGTTGATGCTGCATGAGCACGTTCCCGCCGTACCAGCGGCAATAACCGTTTCTTGCGCTATTCTCGGCATCGTGCTGGTGGTAACACGCGATGGCTGGTTAAGTCTTTTTATGGCGGCAGTCGTTTTAACCCATACCCCATTTCTATCGCTTCTCTTTTTCGTcctgatttcgg                        >  1:463767/1‑241 (MQ=255)
                                                                                                                                                                   cgccCTGGTCGTTGCTGCCGCCAGTGGCTTTCGCGGTGGGCGAATCTTCCTGGCAGTGTTTGTCGGCGTGGCATTAGGGTTGATGCTGCATGAGCACGTTCCCGCCGTACCAGCGGCAATAACCGTTTCTTGCGCTATTCTCGGCATCGTGCTGGTGGTAACACGCGATGGCTGGTTAAGTCTTTTTATGGCGGCAGTCGTTGTACCCAATACCACATTTCTACCGCTGCTCTGTATCGTCATGCTTTCgg                        >  1:560514/1‑251 (MQ=255)
                                                                                                                                                                            ttttgctgccggagggggtttccgcgGTGGGCAAATCTTCCTGGATGTTTTTGTGGGGGTGGCTTTGGGGTTGATGCTGCATGAGCACGTCCCCGCCGTACCAGCGGCAATAACCGTTTCTTGCGCTATTCTCGGCATCGTGCTGGTGGTAACACGCGATGGCTGGTTAAGTCTTTTTATGGCGGCAGTCGTTGTACCCAATACCACATTGCTACCGCTGCTCTGTATCGTCATGCTTCCGGCATGGCTg                <  1:508028/228‑1 (MQ=255)
                                                                                                                                                                                gctgcCGCCAGTGGCTTTCGCGGTGGGCGAATCTTCATGGCAGTGTTTGTCGGCGTGGCATTAGGGTTGATGCTGAATGAGCACGTTCCCGCCGTACCAGCGGCAATAACCGTTTCTTGCGCTATTCTCGGCATCGTGCTGGTGGTAACACGCGATGGCt                                                                                                      >  2:507780/1‑160 (MQ=255)
                                                                                                                                                                                  tgcCGCCAGTGGCTTTCGCGGTGGGCGAATCTTCCTGGCAGTGTTTGTCGGCGTGGCATTAGGGTTGATGCTGCATGAGCACGTTCCCGCCGTACCAGCGGCAATAACCGTTTCTTGCGCTATTCTCGGCATCGTGCTGGTGGTAACACGCGATGGCTGGTTAAGTCTTTTTATGGCGGCAGTCGTTGTACCCAATACCACCTTGCTACCGCTGCTCTGTATCGTCATGCTTCCGGCATGGCTGttatt           >  1:440385/1‑249 (MQ=255)
                                                                                                                                                                                     cgccAGTGGCTTTCGCGGTGGGCGAATCTTCCTGGCAGTGTTTGTCGGCGTGGCATTAGGGTTGATGCTGCATGAGCACGTTCCCGCCGTACCAGCGGCAATAACCGTTTCTTGCGCTATTCTCGGCATCGTGCTGGTGGTAACACGCGATGGCTGGTTAAGTCTTTTTATGGCGGCAGTCGTTGTACCCAATACCCCCTTGCTACCGCTGCTCTGTATCGTCATGCTTCCGGCATGGCTGTTATTAGc        >  1:8807/1‑249 (MQ=255)
                                                                                                                                                                                          gTGGCTTTCGCGGTGGGCGAATCTTCCTGGCAGTGTTTGTCGGCGTGGCATTAGGGTTGATGCTGCATGAGCACGTTCCCGCCGTACCAGCGGCAATAACCGTTTCTTGCGCTATTCTCGGCATCGTGCTGGTGGTAACACGCGATGGCTGGTTAAGTCTTTTTATGGCGGCAGTAGTTGTACCCAATACCACATTGCTACCGCTGCTCTGTATTGTCATGCCTCCGGCATGGCTGTTTTTAGAAGGTaa  >  1:83926/1‑250 (MQ=255)
                                                                                                                                                                                                                      |                                                                                                                                                                                                                                  
ATCCGGTGCTCGTGCTGGGTATTGGCGGATTTATTCTCGGTATTCTGGGGGTTATTGGTGGACCAGTTTCGCTGTTTAAAGGGCTGGATGAGATGCAGCAGATGGTGGCAAATCAGGCTTTCAGCACCAGCGATTACTTTTTGCTGGCGGTAATTAAACTTGCCGCCCTGGTCGTTGCTGCCGCCAGTGGCTTTCGCGGTGGGCGAATCTTCCCGGCAGTGTTTGTCGGCGTGGCATTAGGGTTGATGCTGCATGAGCACGTTCCCGCCGTACCAGCGGCAATAACCGTTTCTTGCGCTATTCTCGGCATCGTGCTGGTGGTAACACGCGATGGCTGGTTAAGTCTTTTTATGGCGGCAGTCGTTGTACCCAATACCACATTGCTACCGCTGCTCTGTATCGTCATGCTTCCGGCATGGCTGTTATTAGCAGGTAA  >  NC_000913/2510396‑2510831

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