Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,514,119 C→T 31.5% D64N (GAT→AAT)  fecB ← ferric citrate ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,514,1190CT31.5% 21.8 / 12.5 19D64N (GAT→AAT) fecBferric citrate ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base C (6/7);  new base T (3/3);  total (9/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 = 8.49e-01

TACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCG‑CATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGAC  >  NC_000913/4513940‑4514368
                                                                                                                                                                                   |                                                                                                                                                                                                                                                          
tACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATCCCGATCGGGATGGCGTCCACGGCCGCCAGGGCCGCCGGCGACCG‑GAGTTCCAAGCCCACAACCACTTGGGGCGttgtgtttg                                                                                                                                                               >  2:208019/1‑265 (MQ=255)
     gggcgATTTGCTGCAAGGGGATGTAAACCCCCGCATGGCGACTGCTGGCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCG‑CATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTa                                                                                                                                                          <  1:208019/270‑1 (MQ=255)
                             aaaCCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCG‑CATCGGCCAACGAGAGTTCCAGCCCCACAATCCGGTGGCGCGTTTTTTCGACTGTAAACGGGCCGGTTTCGttgtgatcc                                                                                                                                  >  2:163170/1‑264 (MQ=255)
                              aaCCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCG‑CATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCg                                                                                                                                                 <  1:163170/256‑1 (MQ=255)
                                       aTGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATTGTCGGCAATACCGATCGGGATGGCGTCCACGGCGGCCAGCG‑CATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCaa                                                                                                                        >  1:190996/1‑272 (MQ=255)
                                       aTGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCG‑CATCGGCGaa                                                                                                                                                                                                        >  1:385378/1‑192 (MQ=255)
                                        tGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCG‑CATCGGCGAAc                                                                                                                                                                                                       <  2:385378/192‑1 (MQ=255)
                                                               cAGGTCTGGTGTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCG‑CATCGGCGAACGAGAGTTCCAGCACCACAAt                                                                                                                                                                                   <  1:373943/189‑1 (MQ=255)
                                                                                               ggCTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATTGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCG‑CATCGGCGAACTAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCg                                                                                                                                                 <  1:408029/191‑1 (MQ=255)
                                                                                                gCTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAgc                                                                                                                                                                                                                    >  1:233379/1‑124 (MQ=255)
                                                                                                 cTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAgcg                                                                                                                                                                                                                   <  2:233379/124‑1 (MQ=255)
                                                                                                                                                acgcacTTCGGGCAGGATGCGTTTTGCATCGTTATTGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCG‑CATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCg                                                                                                                                                                >  1:369496/1‑127 (MQ=255)
                                                                                                                                                acgcacTTCGGGCAGGATGCGTTTTGCATCGTTATTGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCG‑CATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCaaaaa                                                                                    >  1:254144/1‑203 (MQ=255)
                                                                                                                                                acgcacTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCG‑CATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGcc                                                                                                                                                  >  1:26056/1‑141 (MQ=255)
                                                                                                                                                acgcacTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCG‑CATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGAt                                                                                                           >  1:243545/1‑180 (MQ=255)
                                                                                                                                                 cgcacTTCGGGCAGGATGCGTTTTGCATCGTTATTGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCG‑CATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGa                                                                                                                                                               <  2:369496/127‑1 (MQ=255)
                                                                                                                                                 cgcacTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCG‑CATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCg                                                                                                                                                 <  2:26056/141‑1 (MQ=255)
                                                                                                                                                 cgcacTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCG‑CATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATc                                                                                                          <  2:243545/180‑1 (MQ=255)
                                                                                                                                                              ggttggtttttGCATCGTTATTGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCG‑CATCGGCGAACGAGAGTTCCAGCACCAAAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGAc  <  1:369884/264‑1 (MQ=255)
                                                                                                                                                                                   |                                                                                                                                                                                                                                                          
TACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCG‑CATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGAC  >  NC_000913/4513940‑4514368

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