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.5% 50.7 / ‑6.1 23P327L (CCG→CTG) yfeOputative transport protein YfeO
Reads supporting (aligned to +/- strand):  ref base C (0/1);  major base T (13/8);  minor base G (1/0);  total (14/9)
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.33e-01

CACGGTTGCACGCGATGATGCATCAAATGAAAAATCCGGTGCTCGTGCTGGGTATTGGCGGATTTATTCTCGGTATTCTGGGGGTTATTGGTGGACCAGTTTCGCTGTTTAAAGGGCTGGATGAGATGCAGCAGATGGTGGCAAATCAGGCTTTCAGCACCAGCGATTACTTTTTGCTGGCGGTAATTAAACTTGCCGCCCTGGTCGTTGCTGCCGCCAGTGGC‑TTTCGCGGTGGGCGAATCTTCCCGGC‑AGTGTTTGTCGGCGTGGCATTAGGGTTGATGCTGC‑ATGAGCACGTTCCCGCCGTACCAGCGGCAATAACCGTTTCTTGCGCTATTCTCGGCATCGTGCTGGTGGTAACACGCGATGGCTGGTTAAGTCTTTTTATGGCGGCAGTC‑GTTGTACCCAATAC‑CACATTGCTACCGCTGCTCTGTATCGTCATGCTTCCGGCATGGCTGTTATTAGCAGGTAAGCCGATGATGATGGTCA  >  NC_000913/2510363‑2510848
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   
cACGGTTGCACGCGATGATGCATCAAATGAAAAATCCGGTGCTCGTGCTGGGTATTGGCGGATTTATTCTCGGTATTCTGGGGGTTATTGGTGGACCAGTTTCGCTGTTTAAAGGGCTGGATGAGATGCAGCAGATGGTGGCAAATCAGGCTTTCAGCACCAGCGATTACTTTTTGCTGGCGGGAATTAAAATTTCCGCCCCGGTCGTTGCTGCCGCCAGTGGC‑TTTTGCGGTGGGGGACTCTTCCTGGcc                                                                                                                                                                                                                                                 >  2:108391/1‑250 (MQ=255)
      tGCACGCGATGATGCATCAAATGAAAAATCCGGGGCTCGTGCTGGGTATTGGCGGATTTATTCTCGGTATTCTGGGGGTTATTGGTGGACCAGTTTCGCTGTTTAAAGGGCTGGATGAGATGCAGCAGATGGTGGCAAATCAGGCTTTCAGCACCAGCGATTACTTTTTGCTGGCGGTAATTAAACTTGCCGCCCTGGTCGTTGCTGCCGCCAGTGGC‑TTTCGCGGTGGGCGAATCTTCCTGGC‑AGTGtt                                                                                                                                                                                                                                           <  1:353536/250‑1 (MQ=255)
         aCGCGATGATGCATCAAATGAAAAATCCGGTGCTCGTGCTGGGTATTGGCGGATTTATTCTCGGTATTCTGGGGGTTATTGGTGGACCAGTTTCGCTGTTTAAAGGGCTGGATGAGATGCAGCAGATGGTGGCAAATCAGGCTTTCAGCACCAGCGATTACTTTTTGCTGGCGGTAATTAAACTTGCCGCCCTGGTCGTTGCTGCCGCCAGTGGC‑TTTCGCGGTGGGCGAATCTTCCTGGC‑AGTGTTTGTc                                                                                                                                                                                                                                       <  2:66268/251‑1 (MQ=255)
                    cgtCAAATGAAAAATCCGGTGCTCGTGCTGGGTATTGGCGGATTTATTCTCGGTATTCTGGGGGTTAGTGGTGGACCAGTTTCGCTGTTTAAAGGGCTGGATGAGATGCAGCAGATGGTGGCAAATCAGGCTTTCAGCACCAGCGATTACTTTTTGCTGGCGGTAATTAAACGTGCCGCCCTGGTCGTTGCTGCCGCCAGTGGC‑TTTCGCGGTGGGCGAATCTTCCTGGC‑AGTGTTTGTCGGCGTGGCAtt                                                                                                                                                                                                                            <  2:176831/249‑1 (MQ=255)
                       cAAATGAAAAATCCGGTGCGCGTGCTGGGTATTGGCGGATTTATTCTCGGTATTCTGGGGGTTATTGGTGGACCAGTTGCGCTGTTTAAAGGGCTGGATGAGATGCAGCAGATGGTGGCAAATCAGGCTTTCAGCACCAGCGATTACTGTTTGCTGGCGGTAATTAAACTTGCCGCCCTGGTCGTTGCTGCCGCCAGTGGC‑TTTCGCGGTGGGCGAATCTTCCTGGC‑AGTGTTTGGCGGCGTGGCATTAgg                                                                                                                                                                                                                         <  2:41963/251‑1 (MQ=255)
                                aaTCCGGTGCTCGTGCTGGGTATTGGCGGATTTATGTTCGGTATTCTGGGGGTTATTGGTGGACCAGTTTCGATGTTTAAAGGGCTGGAGGAGATGCAGCAGATGGTGGCAAATCAGGCTTTCAGCACCAGCGATTTCTTTGTGGTGGGGGTAATTAAACTGGCCGCCTTGGGCGTGGCTGCCGCCAGTGGC‑TTTCGCGGTGGGCGAATCTTCCTGGC‑AGTGTTTGTCGGCGTGGCATTAGTGTTGAtgct                                                                                                                                                                                                                <  2:381425/251‑1 (MQ=255)
                                          tCGTGCTGGGTATTGGCGGGTTTATTCCCGGCATTCTGGGGGCTATTGGTGGAAACGATTCGCTGTTTAAAGGGCTGGATGAGATGCAGCAAATGGTGGCACATCCGGCTTTCAGCACCAGAGATTACTTTTTGCTGGGGGTTAATAAACCTGGCGCCCCGGGCGATTCTGGCGCCACGGGCGTTTC‑CGGTGGGGGAATCTTTCTGGC‑CGCGTTTGTTGGCGCGGcgtttggggtgaggttgcacgtggcc                                                                                                                                                                                                      >  2:276426/1‑226 (MQ=255)
                                           cGTGCTGGGTATTGGCGGATTTATTCTCGGTATTCTGGGGGTTATTGGTGGACCAGTTTCGCTGTTTAAAGGGCTGGATGAGATGCAGCAGATGGTGGCAAATCAGGCTTTCAGCACCAGCGATTACTTTTTGCTGGCGGTAATTAAACTTGCCGCCCTGGTCGTTGCTGCCGCCAGTGGC‑TTTCGCGGTGGGCGAATCTTCCTGGC‑AGTGTTTGTCGGCGTGGCATTAGGGTTGATGCTGC‑ATGAGCAc                                                                                                                                                                                                     <  1:108391/250‑1 (MQ=255)
                                           cGTGCTGGGTATTGGCGGATTTATTCTCGGTATTCTGGGGGTTATTGGTGGACCAGTTTCGCTGTTTAAAGGGCTGGATGAGATGCAGCAGATGGTGGCAAATCAGGCTTTCAGCACCAGCGATTACTTTTTGCTGGCGGTAATTAAACTTGCCGCCCTGGTCGTTGCTGCCGCCAGTGGC‑TTTCGCGGTGGGCGAATCTTCCTGGC‑AGTGTTTGGCGGCGTGGCATTAGGGTTTATGCTGC‑ATGAGCACg                                                                                                                                                                                                    <  2:217454/251‑1 (MQ=255)
                                                                          tttgtGGGGTTTGTTGGTGTACCTGTTTCGTTGTTTAAGGGGCGGGAGAAGATGAGGCAGATGGGGGAAATTGGGGCTTTCACCCCGAGCGTTTTCTTTTTGCGGGGGGTAATAAAGCTGGCCCCCGTGGTCGTTGCCGCCGCCGGTGGC‑TTTCGGGGTGGGCGAATCT‑CCCTGCTGGTGTTTGTCGGCGGGGAATTAGGGTGGATG‑TGCCATGAGCGCGTTCCCGCCGTACCAGCGGCAATAACCGttt                                                                                                                                                                      <  2:341442/246‑1 (MQ=255)
                                                                                        tggtggACCTGTTGCGTTGTTAAAAGGGCGGGATGAGTGGCAGGAGATGGTGGAAAGTCAGGTTTGCAGCCCCAGCGATTTTTTTTTGGGGGCGGTAATTAAACTTGCCGCTTTGGTCGGTGCGGCCGCCGGTGGC‑TTGCGGGGGGGGCGAACCTCCCTGGT‑AGTGTTGGGGGTCGTGGAATTGGGTTTGATGCTGC‑GTGGGCGGGTGCCCGCCGTACCAGCGGCAATAACCGTTTCTTGCGCTATtctc                                                                                                                                                        <  1:189964/250‑1 (MQ=255)
                                                                                             gACCAGATTAGCTGTATAAAGGGCTGGATGAGATGCAGCAGTTGGTGGCAAATCAGGCTTTCAGCACCAGCGATTACTTTTTGCTGGCGGGAAATAAACTTGCCACCCTGGTCATGGCTGCCGCCAGTGGC‑TTTTGCGGTGGGCGAAACTTGAGGGC‑AGGGTTTGTCGGCGGGGCATTAGGGGTGATGCAGC‑ATGATAACGGTCGCGCAGTACACGCGGCAATAACCGtttatttacgtcttgtaggcatt                                                                                                                                                  >  2:390070/1‑231 (MQ=255)
                                                                                                    ttCGCTGTTTAAAGGGCTGGATGAGATGCAGCAGATGGTGGCAAATCAGGCTTTCAGCACCAGCGATTACTTTTTGCTGGCGGTAATTAAACTTGCCGCCCTGGGCGTTGCTGCCGCCCGTGGC‑TTTTGCGGTGGGCGGATCTTCCTGGC‑CGTGTTTGTCGGCGTGGCCTTAGGGTTGATGCTGC‑ATGAGCACGTTCCCGCCGTACCAGCGGCAATAACAGTTTCTTGCGCTATTCTCGGCATCGGGCtgg                                                                                                                                           >  2:281238/1‑251 (MQ=255)
                                                                                                    atcgctggatgaaGGGCTGGATGAGATGCAGCAGATGGTGGAAAATCAGGCTTTCAGCACCAGCGATTACTATTTGCTGGCGGCAAATAAAATTTCCACCCCGCGCGATGCACGCGCCAGTGGC‑TTTTGCGGGGGGGGAATCTTCCTGTC‑AGTGTTTGTCGACCAGTACTTAAGGTTGATGGTGA‑CTGAGCACGTTCCCGCCGTACCAGCCACAATCACCGCTTCTTGCGCCATTCTCCGCAGCATGCtgg                                                                                                                                           >  2:262832/12‑251 (MQ=255)
                                                                                                                  ggCTGGATGAGATGCAGCAGATGGTGGCAAATCAGGCTTTCAGCACCAGCGATTACTTTTTGCTGGCGGTAATTAAACTTGCCGCCCTGGTCGTTGCTGCCGCCAGTGGC‑TTTCGCGGGGGGCGAATCTTCCTGGC‑CGTGTTTGTCGGCGTGGCATTAGGGTTGATGCTGC‑ATGAGCACGTTCCCGCCGTACCAGCGGCAATAAACGTTTCTTGCGCTATTCTCGGCATCGTGCTGGTGGTAACACGCGAt                                                                                                                             >  1:300261/1‑251 (MQ=255)
                                                                                                                                cagcagATGGTGGCAAATCAGGCTTTCAGCACCAGCGATTACTTTTTGCTGGCGGAAATTAAACTTGCCGCCCTGGGCGTTGCTGCCGCCAGTGGC‑TTTCGCGGGGGGCGCATCTTCCTGGG‑AGCGGTTGTCGGCGTGGCTTTATGGGTGATGCTGC‑ATTAGCACGTTCCCGCCCTTCCAGCGGCCATAACCGGTTCTTTTCCCATTCTCCGCACCGTGCGGGGGGTAACACGCGGGTGCTGGGTAAGTCt                                                                                                               >  2:326100/1‑251 (MQ=255)
                                                                                                                                                  cAGGCTTTCAGCACCAGCGATTACTTTTTGCTGTCGGAAATTAAAATTGCCGCCCTGGGCGTTGCTGCCGCCAGTGGC‑TTTCGCGGTGGGCGAATCTTCCTGGC‑AGTGTTTGTCGGCGTGGCATTAGGGGTGATGCTGC‑ATGTGCCCGTTCCCGCCGTAACCGCGGCAATAACCGGTTCTTGCGCTATTCTCGGCATCGGGCTGGTGGTAACCCGCGATGGCTGGTTAAGTCTTTTTATggcggggggcgg                                                                                             >  1:33917/1‑244 (MQ=255)
                                                                                                                                                  cAGGCTTTCAGCACCAGCGATTACTTTTTGCTGGCGGTAATTAAACTTGACGCCCTGGTCGTTGCTGCCGCCAGTGGC‑TTTCGCGGTGGGCGACTCTTCCTGGC‑AGTGTTTGTCGGCGTGGCCTTAGGGTTGATGCTGC‑ATGAGCACGTTCCCGCCGTACCAGCGGCAATAACCGTTTCTTGCGCTATTCTCGGCATCGTGCTGGGGGGAACACGCGATGGCTGGGTAAGTCTTTTTATggcggcggacgg                                                                                             >  2:253681/1‑245 (MQ=255)
                                                                                                                                                                cAGCGATTACTTTTTGCTGGCGGTAATTAAACTTGCCGCCCTGGTCGTTGCTGCCGCCCGTGGC‑TTTCGCGGTGGGCGAATCTTCCTGGC‑AGTGTTTGTCGGCGTGGCATTAGGGTTGATGCTGC‑ATGAGCACGTTTCCGCCGTACCAGCGGCAATAACCGTTTCTTGCGCTTTTCTCGGCATCGTGCTGGTGGGAACACGCGATGGCTGGTTAAGTCTTTTTATGGCGGCAGGCTGTTTTACCCAAAAc‑c                                                                              >  1:50233/1‑251 (MQ=255)
                                                                                                                                                                                                           gTCGTTGCTGCCGCCAGTGGC‑TTTCGCGGTGGGCGAATCTTCCTGGC‑AGTGTTTGTCGGCGTGGCATTAGGGTTGATGCTGC‑ATGAGCACGTTCCCGCCGTACCAGCGGCAATAACCGTTTCTTGCGCTATTCTCGGCATCGTGCTGGTGGTAACACGCGATGGCTGGTTAAGTCTTTTTATGGCGGCAGTC‑GGTGTACCCAATAC‑CACATTGCTACACCTGCTCTGTATCGTCATGCTTCCGGCCTGGCt                                  >  1:338278/1‑251 (MQ=255)
                                                                                                                                                                                                                            tGGC‑TTTCGCGGTGGGCGAATCTTACTGGC‑AGTGTTTGTCGGCGTGGCATTAGGGTTGATGCTGC‑ATGAGCACGTTCCCGCCGTACCAGCGGCAATAACCGTTTCTTGCGCTATTCTCGGCATCGTGCTGGTGGTAACACGCGATGGCTGGTTCAGTCTTTTTATGGCGGCAGTC‑GCTGTACCCCA‑ACACACATTTATACCCATGCTCTGTTTCGTCATACGTCCCGCATGGcctttcttttagagaaa                   >  2:92605/1‑233 (MQ=255)
                                                                                                                                                                                                                                          gggCGAATCTTCCTGGC‑AGTGTTTGTCGGCGTGGCATTAGGGTTGATGCTGC‑ATGAGCACGTTCCCGCCGTACCAGCGGCAATAACCGTTTCTTGCGCTATTCTCGGCATCGTGCTGGTGGTAACACGCGATGGCTGGTTAAGTCTTTTTATGGCGGCAGTC‑GTTGTACCCAATAC‑CACATTGCTACCCCTGCTCTGTATCGTCATGCTTCCGGCATGGCTGTTATTATCAGGTAAGCCGATGATGATggg    >  1:132013/1‑250 (MQ=255)
                                                                                                                                                                                                                                             cGAATCTTCCTGGC‑AGTGTTTGTCGGCGTGGCATTAGGGTTGATGCTGC‑ATGAGCACGTTCCCGCCGTACCAGCGGCAATAACCGTTTCTTGCGCTATTCTCGGCATCGTGCTGGTGGTAACACGCGATGGCTGGTTAAGTCTTTTTATGGCGGCAGTC‑GTTGTATCCAATAC‑CACATTGCTACCGCTGCTCTGTATCGTCATGCTTCCGGCATGGCTGTTATTAGGCGGTAAGCCGATGGTGATGGTCa  >  1:276862/1‑250 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   
CACGGTTGCACGCGATGATGCATCAAATGAAAAATCCGGTGCTCGTGCTGGGTATTGGCGGATTTATTCTCGGTATTCTGGGGGTTATTGGTGGACCAGTTTCGCTGTTTAAAGGGCTGGATGAGATGCAGCAGATGGTGGCAAATCAGGCTTTCAGCACCAGCGATTACTTTTTGCTGGCGGTAATTAAACTTGCCGCCCTGGTCGTTGCTGCCGCCAGTGGC‑TTTCGCGGTGGGCGAATCTTCCCGGC‑AGTGTTTGTCGGCGTGGCATTAGGGTTGATGCTGC‑ATGAGCACGTTCCCGCCGTACCAGCGGCAATAACCGTTTCTTGCGCTATTCTCGGCATCGTGCTGGTGGTAACACGCGATGGCTGGTTAAGTCTTTTTATGGCGGCAGTC‑GTTGTACCCAATAC‑CACATTGCTACCGCTGCTCTGTATCGTCATGCTTCCGGCATGGCTGTTATTAGCAGGTAAGCCGATGATGATGGTCA  >  NC_000913/2510363‑2510848

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