Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 840,003 C→T G510D (GGC→GAC)  fiu ← iron catecholate outer membrane transporter Fiu

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913840,0030CT88.9% 15.0 / NA 10G510D (GGC→GAC) fiuiron catecholate outer membrane transporter Fiu
Reads supporting (aligned to +/- strand):  ref base C (0/1);  major base T (4/4);  minor base G (0/1);  total (4/6)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAACAGCGCGGCGGTGAGCAACAGACGTTTATCCAGAACCTGCCATTTGGTGCCAATCTCGCTGGTGTTGGCTTTTTGCGGTTTAAAATCGGTGCGGTTGGCACTGTTACCGCTGCCAGACTGCGCAAGGGCGAAGTTGTTGCCGCCCGGAGGCTGCTGGGAAACGGCATAGTTAATATAGACATTGCCGTTTTCCGTCAGGTGATACAGCGCCCCGGCTTTCCAGTTCATCAGATTGCCCGACTTGGCGGTGTCGACGGTGGTGACCGGAGAACCTTTTGCCACACCAGTTGGGCAGGTGATGGCACCGCGTCCGCTGCCGCCGCAGGCGGTGGCACTGTCATATTCAGTATGATAATTATCCAGACGGATCCCGCCGTTCAGCTCAAAATCACGGGTGATTTGCAGCGTATCAAAGGCG  >  NC_000913/839765‑840185
                                                                                                                                                                                                                                              |                                                                                                                                                                                      
gAACAGCGCGGCGGTGAGCAACAGACGTTTATCCAGAACCTGCCATTTGGTGCCAATCTCGCTGGTGTTGGCTTTTTGCGGTTTAAAATCGGAGCGGTTGGCACTGTTACCGCTGCCAGACTGTGCAAGGGCGAAGTTGTTGCCGCCCGGAGGCTGCTGGGAAACGGCATAGTTAATATAGACATTGCCGTTTTCCGTCAGGTGATACAGCGCCCCGGCTTTCCAGTTCATCAGATTGTCCGACTTGGCgg                                                                                                                                                                            >  1:241456/1‑251 (MQ=255)
                                                               ggTGTTGGCTTTTTGCGGTTTAAAATCGGTGCGGTTGGCACTGTTACCGCTGCCAGACTGTGCAAGGGCGAAGTTGTTGCCGCCCGGAGGCTGCTGGGAAACGGCATAGTTAATATAGACAGTGCCGTTGTCCGTCAGGTGATACAGCGCCCCGGCTTTCCAGTTCATCAGATTGGCCGACTGGGCGGTGTCGACGGTGGTGACCGGAGAAGCTTTTGCCACACCAGTTGGGCAGGTGATGGCACCGCGtn                                                                                                             <  1:32117/251‑2 (MQ=255)
                                                                             gCGGTTTAAAATCGGTGCGGTTGGCACTGTTACCGCTGCCAGACTGTGCCAGGGCGAAGATGTCGCCGCCCGGCGGCCGCTGGGAAAAGGCATAGTTAATATAAACATTGCCGTTTTACGTCAGGCGATACAGCGCCCCGGCTTTCCAGTTCATCAGATTGTCCGACTTGGCGGTGGAGCCGGAGGTGACCGGAGAACCTTTTGCCCCACCAGTTGGGCAGGCGATGGCACCGCGTCCGCTGCcgccgc                                                                                                 >  1:250680/1‑249 (MQ=255)
                                                                                                            aCCGCTGCCAGACTGTGCAAGGGCGAAGTTGTTGCCGCCCGGAGGCTGCTGGGAAACGGCATAGTTAATATAGACATTGCCGTTTTCCGTCAGGTGATACAGCGCCCCGGCTTTCCAGTTCATCAGATTGTCCGACTTGGCGGTGTCGACGGTGGTGACCGGAGAACCTTTTGCCACACCAGTTg                                                                                                                                  <  1:8150/185‑1 (MQ=255)
                                                                                                                                  gCGAAGTTGTTGCCGCCCGGAGGCTGCTGGGAAACGGCATAGTTAATATAGACATTGCCGTTTTCCGTCAGGTGATACAGCGCCCCGGCTTTCCAGTTCATCAGATTCTCCGACTTGGCGGCGTCGACGGTGGTGCCCggcgaccc                                                                                                                                                   >  1:133990/1‑140 (MQ=255)
                                                                                                                                      aGTTGTTGCCGCCCGGAGGCTGCTGGGAAACGGCATAGTTAATATAGACATTGCCGTTTTCCGTCAGGTGATACAGCGCCCCGGCTTTCCAGTTCATCAGATTGTCCGACTTGGCGGTGTCGACGGTGGTGACCGGAGAACCTTTTGCCACACCAGTTGGGCAGGTGATGGCACCGCGTCCGCTGCCGCCGCAGGCGGTGGCACTGTCATATTCAGTATGATAATTATCCAGACGGATCCCGCCGTTCAGc                                      >  1:84701/1‑251 (MQ=255)
                                                                                                                                         tgttgCCGCCCGGAGGCTGCTGGGAAACGGCATAGTTAATATAGACATTGCCGTTTTCCGTCAGGTGATACAGCGCCCCGGCTTTCCAGTTCATCAGATTGTCCGACTTGGCGGTGTCGACGGTGGTGACCng                                                                                                                                                         <  1:199556/133‑1 (MQ=255)
                                                                                                                                                      gggcggctgGGAAACGGCATGGTTAATATAGACATTGCCGTTTTCCGTCAGGGGATACAGCGCCCCGTCTTTCCAGTTCATCAGATTGCCCGACTTGGCGGGGGCGACGTGGGTGACCGGAGACCCTTTTGCCACACCAGTTGGGCAGGTGATGGCACCGCGGCCGCTGCCGCCGCAGGCGGTGGCTCTGTGATATTCAGTATGATAATTTTCCAGCCGGATCCCGCCGGTCAGCTCAAAATCACGGGTGa                      <  1:132460/246‑1 (MQ=255)
                                                                                                                                                                          atttAATATAGACATTGCCGTTTTCCGTCAGGTGATACAGCGCCCCGGCTTTCCAGTTCATCAGATTGTCCGACGGGGCGGTGTCGACGGTGGTGACCGGAGAACCTTTTGCCACACCAGTTGGGGAGGTGATGGCACCGCGTCCGCTGCCACCGCAGGCGGTGGCACTGGCATATTCAGTATGATAATTATCCAGACGGATCCCGCCGTTCAGCTCAAAATCACGGTTGATTTGCAGCGTATCAAAGGCg  <  1:209248/249‑1 (MQ=255)
                                                                                                                                                                                                        ggTGATACAGCGCCCCGGCTTTCCAGTTCATCAGATTGTCCGACTTGGCGGTGTCGACGGTGGGGACCGGAGAACCTTTTGCCACACCAGTTGGGCAGGTg                                                                                                                          <  1:95439/101‑1 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                      
GAACAGCGCGGCGGTGAGCAACAGACGTTTATCCAGAACCTGCCATTTGGTGCCAATCTCGCTGGTGTTGGCTTTTTGCGGTTTAAAATCGGTGCGGTTGGCACTGTTACCGCTGCCAGACTGCGCAAGGGCGAAGTTGTTGCCGCCCGGAGGCTGCTGGGAAACGGCATAGTTAATATAGACATTGCCGTTTTCCGTCAGGTGATACAGCGCCCCGGCTTTCCAGTTCATCAGATTGCCCGACTTGGCGGTGTCGACGGTGGTGACCGGAGAACCTTTTGCCACACCAGTTGGGCAGGTGATGGCACCGCGTCCGCTGCCGCCGCAGGCGGTGGCACTGTCATATTCAGTATGATAATTATCCAGACGGATCCCGCCGTTCAGCTCAAAATCACGGGTGATTTGCAGCGTATCAAAGGCG  >  NC_000913/839765‑840185

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