Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 582,597 C→T intergenic (+500/‑278) tfaD → / → tfaX putative tail fiber assembly protein TfaD/protein TfaX

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913582,5970CT100.0% 28.0 / NA 11intergenic (+500/‑278)tfaD/tfaXputative tail fiber assembly protein TfaD/protein TfaX
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (4/7);  total (4/7)

ACAGCTTTTTGTAGATATATTTCCAGGCAGGCAGTGAAGAAGGTAAGACATGGAAATGGAATCAAATTGACCATGTAACGCCGCGGGATAAGGTTCAAAAACATCATGGCTAATTTTATGTTTAATTTTTGATTCCCCAGCCCTTGTAGATGCCGCGTTCAGGCTAGCTTCGTTCAAATCCATTAAAGATATCAGACTACTCTCAGGTACGTGAGTAAGGTAAAACCCAGTTCCAACACCAATATCCAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTAAATTCTAGCCCCATCATCTGTGTTTTTTTATTAATTTCACCATGTTATAGTTTTAT  >  NC_000913/582401‑582831
                                                                                                                                                                                                    |                                                                                                                                                                                                                                          
atatcTTTTTGTAGATATATTTCCAGGCAGGCAGTGAAGAAGGTAAGAAATGGAAATGGAATCAAATTTATCATGTAACGCCGCGGGATAAGGTTCAAAAACATCATGGCTAATTTTATGTTTAATTTTTGATTCCCCAGCCCTTGTAGATGCCGCGTTCAGGCTAGCTTCGTTCAAATCCATTAAAGATATCAGATTACTCTCAGGTACGTGAGTAAGGTAAAACCCAGTTCCAACACCAATATCCAGAt                                                                                                                                                                                      <  2:63937/247‑1 (MQ=255)
acaGCTTTTTGTAGATATATTTCCAGGCAGGCAGTGAAGAAGGTAAGACATGGAAATGGAATCAAATTGACCATGTAACGCCGCGGGATAAGGTTCAAAAACATCATGGCTAATTTTATGTTTAATTTTTGATTCCCCAGCCCTTGTAGATGCCGCGTTCAGGCTAGCTTCGTTCAAATCCATTAAAGATATCAGATTACTCTCAGGTACGTGAGTAAGGTAAAACCCAGTTCCAACACCAATATCCAGAt                                                                                                                                                                                      <  1:313337/251‑1 (MQ=255)
                      ccAGGCAGGCAGTGAAGAAGGTAAGACATGGAAATGGAATCAAATTGACCATGTAACGCCGCGGGATAAGGTTCAAAAACATCATGGCTAATTTTATGTTTAATTTTTGATTCCCCAGCCCTTGTAGATGCCGCGTTCAGGCTAGCTTCGTTCAAATCCATTAAAGATATCAGATTACTCTCAGGTACGTGTGTAAGGTAAAACCCAGTTCCAACACCACTATCCAGATGGTTGTTACCTAAATGGTcc                                                                                                                                                                  >  1:210429/1‑249 (MQ=255)
                       caggcaggCAGTGAAGAAGGTAAGACATGGAAATGGAATCAAATTGACCATGTAACGCCGCGGGATAAGGTTCAAAAACATCATGGCTAATTTTATGTTTAATTTTTGATTCCCCAGCCCTTGTAGATGCCGCGTTCAGGCTAGCTTCGTTCAAATCCATTAAAGATATCAGATTACTCTCAGGTACGTGAGTAAGGTAAAACCCAGTTCCAACACCAATATCCAGATGGTTGTTACCTAAATGTTCCAGa                                                                                                                                                               <  1:67655/251‑1 (MQ=255)
                                                                       cATGTAACGCCGCGGGATAAGGTTCAAAAACATCATGGCTAATTTTATGTTTAATTTTTGATTCCCCAGACCTTGTAGATGCCGCGTTCAGGCTAGCTTCGTTCAAATCCATTAAAGATATCAGATTACTATCAGGTACGAGAGTAAGGTAAAACCCAGTTCCAAAACCAATATCCCGCTGGTTGTTACCTAGATGTTCCAAAAAGATAGGAAGAAGGAGTTCCTTTGTAGGACATCCCAATGCAAGcagc                                                                                                               >  2:278211/1‑248 (MQ=255)
                                                                                              tCAAAAACATCATGGCTAATTTTATGTTTAATTTTTGATTCCCCAGCCCTTGTAGATGCCGCGTTCAGGCTAGCTTCGTTCAAATCCATTAAAGATATCAGATTACTCTCAGGTACGTGAGTAAGGTAAAACCCAGTTCCAACACCAATATCCAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAACccacc                                                                                         <  1:41090/250‑1 (MQ=255)
                                                                                                                                       cccAGCCCTTGTAGATGCCGCGTTCAGGCTAGCTTCGTTCAAATCCATTAAAGATATCAGATTACTCTCAGGTACGTGATTAAGGTAAAACCCAGTTCCAACACCAATATCCAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACATcc                                                                                                                            <  2:34444/174‑1 (MQ=255)
                                                                                                                                       cccAGCCCTTGTAGATGCCGCGTTCAGGCTAGCTTCGTTCAAATCCATTAAAGATATCAGATTACTCTCAGGTACGTGAGTAAGGTAAAACCCAGTTCCAACACCAATATCCAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACATcc                                                                                                                            >  1:34444/1‑174 (MQ=255)
                                                                                                                                                 ttataTGCCCCGTCCACGCTAGCTTCGTTCAAATCCATTAAAGATATCAGATTATTCTCAGGTACGTTAGGAAGGTAAAACCCCGTTCCAACACCAATATCCATATGTTTGTTATCTAAATGTTCCAGAAAGTGGGGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATGTGATACTCCCAAAACCCACCAGTCATAAAGCTTTAGGGTAAGGGGTTTGTAAATTCTTGCCCTATCATCtg                                      <  2:382360/246‑1 (MQ=255)
                                                                                                                                                                               aaaTCCATTAAAGATATCAGATTACTCTCAGGTACGTGAGTAAGGTAAAACCCAGTTCCAACACCAATATCCAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGGGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTAAATTCTAGCCCCATCATCTGTGTTTTTTTATTAATTTCACCATGTTATAGt       >  1:85784/1‑251 (MQ=255)
                                                                                                                                                                                    cATTAAAGATATCAGATTACTCTCAGGTACGTGAGTAAGGTAAAACCCAGTTCCAACACCAATATCCAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTAAATTCTAGCCCCATCATCTGTGTTTTTTTATTAATTTCACCATGTTATAGTTTTAt  <  1:163789/251‑1 (MQ=255)
                                                                                                                                                                                                    |                                                                                                                                                                                                                                          
ACAGCTTTTTGTAGATATATTTCCAGGCAGGCAGTGAAGAAGGTAAGACATGGAAATGGAATCAAATTGACCATGTAACGCCGCGGGATAAGGTTCAAAAACATCATGGCTAATTTTATGTTTAATTTTTGATTCCCCAGCCCTTGTAGATGCCGCGTTCAGGCTAGCTTCGTTCAAATCCATTAAAGATATCAGACTACTCTCAGGTACGTGAGTAAGGTAAAACCCAGTTCCAACACCAATATCCAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACATCCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTAAATTCTAGCCCCATCATCTGTGTTTTTTTATTAATTTCACCATGTTATAGTTTTAT  >  NC_000913/582401‑582831

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