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,6090CT88.9% 35.0 / ‑5.1 18P327L (CCG→CTG) yfeOputative transport protein YfeO
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (11/5);  minor base G (0/2);  total (11/7)
Fisher's exact test for biased strand distribution p-value = 1.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

GCACGCGATGATGCATCAAATGAAAAATCCGGTGCTCGTGCTGGGTATTGGCGGATTTATTCTCGGTATTCTGGGGGTTATTGGTGGACCAGTTTCGCTGTTTAAAGGGCTGGATGAGATGCAGCAGATGGTGGCAAATCAGGCTTTCAGCACCAGCGATTACTTTTTGCTGGCGGTAATTAAACTTGCCGCCCTGGTCGTTGCTGCCGCCAGTGGCTTTCGCGGTGGGCGAATCTTCCCGGCAGTGTTTGTCGGCGTGGCATTAGGGTTGATGCTGCATGAGCACGTTCCCGCCGTACCAGCGGCAATAACCGTTTCTTGCGCTATTCTCGGCATCGTGCTGGTGGTAACACGCGATGGCTGGTTAAGTCTTTTTATGGCGGCAGTCGTTGTACCCAATACCACATTGCTACCGCTGCTCTGTATCGTCATGCTTCCGGCATGGCTGTTATTAGCAGGTAAG  >  NC_000913/2510370‑2510832
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                
gCACGCGATGATGCATCAAATGAAAAATCCGGTGCTCGTGCTGGGTATTGGCGGATTTATTCTCGGTATTCTGGGGGTTATTGGTGGACCAGTTTCGCTGTTTAAAGGGCTGGATGAGATGCAGCAGATGGTGGCAAATCAGGCTTTCAGCACCAGCGATTACTTTTTGCTGGCGGTAATTAAACTTGCCGCCCTGGTCGTTGCTGCCGCCAGTGGCTTTCGCGGTGGGCGAATCTTCCTGGCATTGttt                                                                                                                                                                                                                        <  2:28939/250‑1 (MQ=255)
       atgatgCATCAAATGAAAAATCCGGTGCTCGTGCTGGGTATTGGCGGATTTATTCTCGGTATTCTGGGGGTTATTGGTGGACCAGTTTCGCTGTTTAAAGGGCTGGGTGAGATGCAGCAGATGGTGGCAAATCAGGCTTTCAGCACCAGCGATTACTTTTTGCTGGCGGTAATTAAAACTTCCCCCCTGGTCGTTGCTGACCCCAGTGGCTTTTGCGGGGGGCGGATCTTCCTGGGAGGGTTTGttcgggg                                                                                                                                                                                                                >  1:105831/1‑245 (MQ=255)
          atgCATCAAATGAAAAATCCGGTGCTCGTGCTGGGTATTGGCGGATTTATTCTCGGTATTCTGGGGGTTATTGGTGGACCAGTTTCGCTGTTTAAAGGGCTGGATGAGATGCAGCAGATGGTGGCAAATCAGGCTTTCAGCACCAGCGATTACTTTTTGCTGGCGGTTATTAAAATTTCCGCCCTGGTCGGTGCTGCCGCCCGTGGCTTTTGCGGGGGGGGAGTCTTCCTGGCCGTGGTTGTCGGCgcggg                                                                                                                                                                                                             >  1:207390/1‑247 (MQ=255)
                      aaaaaTCCGGTGCTCGTGCGGGTTATTGGCGGATTTATTCTCGGTATTCGGGGGGTTATTGGTGGACCAGTTTCGCTGTTTAAAGGGCTGGATGAGATGCAGCAGATGGTGGCAAATCAGGCTTTCAGCACCAGCGATTACTTTTTGCTGGCGGTAATTAAACTTGCCGCCCTGGTCGTTGCTGCCGCCAGTGGCGTTCGCGGTGGGCGAATCTTCCTGGCAGTGTTTGTCGGCGTGGCATTAGGGTTGAt                                                                                                                                                                                                <  2:5867/251‑1 (MQ=255)
                                    cGTGCTGGGTATTGGCGGATTTATTCTCGGTATTCTGGGGGTTATTGGTGGACCAGTTTCGCTGTTTAAAGGGCTGGATGAGATGCAGCAGATGGTGGCAAATCAGGCTTTCAGCACCAGCGATTACTTTTTGCTGGCGGTAATTAAACTTGCCGCCCTGGTCGTTGCTGCCGCCAGTGGCTTTCGCGGTGGGCGAATCTTCCTGGCAGTGTTTGTCGGCGTGGCATTAGGGTTGATGCTGCATGAGCAc                                                                                                                                                                                   <  1:22772/250‑1 (MQ=255)
                                      ggCTGGGTATTGGCGGATTTATTCTCGGTAGTCTGGGGGTTATTGGTGGACCAGTTTCGCTGTTTAAAGGGCTGGATGAGATGCAGCAGATGGTGGCAAATCAGGCTTTCAGCACCAGCGATTACTTTTTGCTGGCGGTAATTAAACTTGCCGCCCTGGTCGTTGCTGCCGCCAGTGGCTGTCGCGGTGGGCGAATCTTCCTGGCAGTGTTTGTCGGCGTGGCATTAGGGTTGATGCTGCATTAGCACGtt                                                                                                                                                                                <  2:63384/250‑1 (MQ=255)
                                                 ggCGGATTTATTCTCGGTATTCTGGGGGTTATTGGTGGACCAGTTTCGCTGTTTAAAGGGCTGGATGAGATGCAGCAGATGGTGGCAAATCAGGCTTTCAGCACCAGCGATTACTTTTTGCTGGCGGTAATTAAACTTGCCGCCCTGGTCGTTGCTGCCGCCAGTGGCTTTCGCGGTGGGCGAATCTTCCTGGCCGTGTTTGTCGGCGTGGCATTAGGGTTGATGCTGCATGAGCACGTTCCCGCCGTAcc                                                                                                                                                                     >  1:48582/1‑251 (MQ=255)
                                                                         ggggTTATTGGTGGACCAGTTTCGCTGTTTAAAGGGCTGGATGAGATGCAGCAGATGGTGGCAAATCAGGCTTTCAGCACCAGCGATTACTTTTTGCTGGCGGTAATTAAACTTGCCGCCCTGGTCGTTGCTGCCGCCAGTGGCTTTCGCGGTGGGCGAATCTTCCTGGCAGGGTTTGTCGGCGTGGCATTAGGGTTGATGCTGCATTAGCACGTTCCCGCCGTACCAGCGGCAATAACCGTTTCTTtcgc                                                                                                                                             >  1:495486/1‑251 (MQ=255)
                                                                                ttGGTGGACCAGTTTCGCTGTTTAAAGGGCTGGATGAGATGCAGCAGATGGTGGCAAATCAGGCTTTCAGCACCAGCGATTACTTTTTGCTGGCGGTAATTAAACTTGCCGCCCTGGGCGGTGCTGCCGCCAGTGGCTTTCGCGGGGGGCGCATCTTTCTGGCAGGGTTTGTCGGCGTGGCATTTGGGATGATGCTGCCTGGGCAAGTTCCCGCCGTACCAGCGGCCATAACCGTTTCTTgctgcattccc                                                                                                                                      >  2:196822/1‑242 (MQ=255)
                                                                                 gggtggACCAGTTTCGCTGTTTAAGGGGTTGGATGAGATGCAGCAGGTGGTGGAAAATCGGTCTTCCAGCACCAGGGTTTATTTTTGGGTGGGGATAATTAACCGCGCCGCCCTGGTGTTTGCTGCCGCCGGTGGCTTTCGGGGGGGGCTAATCTCCGGGGCTGTTTTTGCCGGCGGGGCATTGGGGTTGTGGCTGCAGAAGCCCTTGCCCGCCGTACCAGCGGCAATAACCGTTTCTTGCGCTATTCTCg                                                                                                                                     <  1:31548/250‑1 (MQ=255)
                                                                                        ccAGTTTCGCTGTTTAAAGGGCTGGATGAGATGCAGCAGATGGTGGCAAATCAGGCTTTCAGCACCAGCGATTACTTTTTGCTGGCGGTAATTAAACTTGCCGCCCTGGTCGTTGCTGCCGCCAGTGGCTTTCGCGGTGGGGGCATCTTCCTGGCAGTGTTTGTCGGCGTGGCATTAGGGTTGGTGCTGCATGAGCACGTTCCCGCCGTACCAGCGGCAATAACCGTTTCTTGCGCTATTCTCGGCATCgg                                                                                                                              >  1:415840/1‑250 (MQ=255)
                                                                                                  tGTTTAAAGGGCTGGATGAGATGCAGCAGATGGTGGCAAATCAGGATTTCCGCACCAGCGATTACTTTTTGCTGGCGGTAATTAAACTTGCCGCCCTGGTCGTTGCTGCCGCCAGTGGCTTTCGCGGTGGGCGAATCTTCCTGGCAGTGTTTGTCGGCGTGGCATTAGGGTTGCTGCTGCATGAGCAAGTTCCCGCCGTACCAGCGGCAATAACCGTTTCTTGCGCTATTCTCGGCATCGGGCTggtggt                                                                                                                     >  1:253576/1‑250 (MQ=255)
                                                                                                                                              gCTTTCAGCACCAGCGATTACTTTTTGCTGGCGGTAATTTAACTTGCCGCCCCGGGCGTTGCTGCCGCCAGTGGCTTTCGCGGTGGGCGAATCTTCCTGGCAGTGTTTGTCGGCGTGGCATTAGGGGTGATGCTGCATGAGCACGTTTCCGCCGTACCAGCGGCAATAACCGTTTCTTGCGCTATTCTCGGCATCGGGCTGGTGGGAACACGCGATGGCTGGGTAAGTCTTTTTATggcggcgggtgtttt                                                                        >  2:324897/1‑242 (MQ=255)
                                                                                                                                                            cGATTACTTTTTGCTGGCGGTAATTAAACTTGCCGCCCTGGTCGTTGCTGCCGCCAGTGGCTTTCGCGGTGGGCGAATCTTCCTGGCAGTGTTTGTCGGCGTGGCATTAGGGTTGATGCTGCATGAGCACGTTCCCGCCGTACCAGCGGCAATAACCGTTTCTTGCGCTATTCTCGGCATCGTGCt                                                                                                                           >  1:121668/1‑186 (MQ=255)
                                                                                                                                                                                                        gTTGCTGCCGCCAGTGGCTTTCGCGGTGGGCGAATCTTCCTGGCAGTGTTTGTCGGCGTGGCATTAGGGTTGATGCTGCATGAGCACGTTCCCGCCGTACCAGCGGCAATAACCGTTTCTTGCGCTATTCTCGGCATCGTGCTGGTGGTAACACGCGATGGCTGGTTAAGTCTTTTTATg                                                                                      >  2:189645/1‑180 (MQ=255)
                                                                                                                                                                                                         gtgcggccgctggggtttttCGGGGTGGGCGATTCTTCCTGGCTGGGTTTGGCGGCGGGGCATTGGGGTTGTTGCGGCAGGAGCCCGTTCCCGCCGTACCAGCGGCAATAACCGTTTCTTGCGCTATTCTCGGCATCGTGCTGGTGGTAACACGCGATGGCTGGTTAAGTCTTTTTATGGCGGCAGTCGTTGTACCCAATACCACATTGCTACCGCTGCTCTGTATCGTCATGCTTCCGGCATGGCTGtt               <  1:471355/233‑1 (MQ=255)
                                                                                                                                                                                                                  ccattGGCTTTCGCGGTGGGCGAATCTTCCTGGCAGTGTTTGTCGGCGTGGCATTAGGGTTGATGCTGCATGAGCACGTTCCCGCCGTACCAGCGGCAATAACCGTTTCTTGCGCTATTCTCGGCATCGTGCTGGTGGTAACACGCGATGGCTGGTTAAGTCTTTTTATGGCGGCAGTCGTTGTACCCAATACCACATTGCTACCGCTGCTCTGGCTCGTCCTTCTTCCGGCCTGGCCGttctttgccagg     >  2:501616/5‑241 (MQ=255)
                                                                                                                                                                                                                     gggtcttccgggttGGGCGAATCTCCCGGGCTGTGTTTGTCGGCTTGGATGTGGGGTTGTGGCGTCTTGCGCACTTTCCCGCCGTACCAGCGGCAATAACCGTTTCTTGCGCTATTCTCGGCATCGTGCTGGTGGTAACACGCGATGGCTGGTTAAGTCTTTTTATGGCGGCAGTCGTTGTACCCAATACCACATTGCTACCGCTGCTCTGTATCGTCATGCTTCCGGCATGGCTGTTATTAGCAGGTAAg  <  1:114219/238‑1 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                
GCACGCGATGATGCATCAAATGAAAAATCCGGTGCTCGTGCTGGGTATTGGCGGATTTATTCTCGGTATTCTGGGGGTTATTGGTGGACCAGTTTCGCTGTTTAAAGGGCTGGATGAGATGCAGCAGATGGTGGCAAATCAGGCTTTCAGCACCAGCGATTACTTTTTGCTGGCGGTAATTAAACTTGCCGCCCTGGTCGTTGCTGCCGCCAGTGGCTTTCGCGGTGGGCGAATCTTCCCGGCAGTGTTTGTCGGCGTGGCATTAGGGTTGATGCTGCATGAGCACGTTCCCGCCGTACCAGCGGCAATAACCGTTTCTTGCGCTATTCTCGGCATCGTGCTGGTGGTAACACGCGATGGCTGGTTAAGTCTTTTTATGGCGGCAGTCGTTGTACCCAATACCACATTGCTACCGCTGCTCTGTATCGTCATGCTTCCGGCATGGCTGTTATTAGCAGGTAAG  >  NC_000913/2510370‑2510832

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