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_000913	582,597	0	C	T	100.0%	87.6 / NA	29	intergenic (+500/‑278)	tfaD/tfaX	putative tail fiber assembly protein TfaD/protein TfaX
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (17/12);  total (17/12)

GTTCCATCGTCAGTTAAGGCCTGCGCCGCATTGCGTATTACACAGCTTTTTGTAGATATATTTCCAGGCAGGCAGTGAAGAAGGTAAGACATGGAAATGGAATCAAATTGACCATGTAACGCCGCGGGATAAGGTTCAAAAACATCATGGCTAATTTTATGTTTAATTTTTGATTCCCCAGCCCTTGTAGATGCCGCGTTCAGGCTAGCTTCGTTCAAATCCATTAAAGATATCAGACTACTCTCAGGTACGTGAGTAAGGTAAAACCCAGTTCCAACACCAATATCCAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACAT‑CCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTAAATTCTAGCCCCATCATCTGTGTTTTTTTATTAATTTCACC  >  NC_000913/582360‑582816                                                                                                                                                                                                                                              |                                                                                                                                                                                                                             gTTCCATCGTCAGTTAAGGCCTGCGCCGCATTGCGTATTACACAGCTTTTTGTAGATATATTTCCAGGCAGGCAGTGAAGAAGGTAAGACATGGAAATGGAATCAAATTGACCATGTAACGCCGCGGGATAAGGTTCAAAAACATCATGGCTAATTTTATGTTTAATTTTTGATTCCCCAGCCCTTGTAGATGCCGCGTTCAGGCTAGCTTCGTTCAAATCCATTAAAGATATCAGATTACTCTCAggaac                                                                                                                                                                                                                 >  2:90675/1‑248 (MQ=255)      aTCGTCAGTTAAGGCCTGCGCCGCATTGCGTATTACACAGCTTTTTGTAGCTATATTTCCAGGCAGGCAGTGAAGAAGGTAAGACATGGAAATGGAATCAAATTGACCATGTAACGCCGCGGGATAAGGTTCAAAAACATCATGGCTAATTTTATGTTTAATTTTTGATTCCCCAGCCCTTTAAGATGCCGCGTTCAGGCTAGCTTCGTTCAAATACATTAAAGATATCAGATTACTCTCAGGTACGTg                                                                                                                                                                                                              >  2:206930/1‑249 (MQ=255)                                                        atatatTGCCAGGCAGGCAGTGAAGAAGGTAAGACATGGAAATTGAATCAAATTGACCATGTAACGCCGCGGGATAAGGTTCAAAAACATCATGGCTAATTTTATGTTTAAGTTTTGATTCCCCAGCCCTTGTAGATGCCGCGTTCAGGCTAGCTTCGTTCAAATCCATTAAAGATATCAGATTACTCTCAGGTACGTGAGTAAGGTAAAACCCAGTTCCAACACCAATATCCAGATGGTTGTTACCTaaa                                                                                                                                                          <  1:195750/251‑1 (MQ=255)                                                                ccAGGCAGGCAGTGAAGAAGGTAAGACATGGAAATGGAATCAAATTGACCATGTAACGCCGCGGGATAAGGTTCAAAAACATCATGGCTAATTTTATGTTTAATTTT‑GATTCCCCAGCCCTTGTAGATGCCGCGTTCAGGCTAGCTTCGTTCAAATCCATTAAAGATATCAGATTACTCTCAGGTACGTGAGTAAGGTAAAACCCAGTTCCAACACCAATATCCAGATGGTTGTTACCTAAATGTTCCAg                                                                                                                                                  <  1:6082/250‑1 (MQ=255)                                                                 caggcaggCAGAGCAGACGGTAAGACATGGAAATGGAATCAAATTGACCATGTAACGCCGCGGGATAAGGATCAAAAACATCATGGCTAATTTTATGTTTAATTTTTGATTCCCCAGCCCTTGTAGATGCCGCGTTCAGGCTAGCTTCTTTCAAATCCATTAAAGATATCAGATTACTCTCAGGTACGGTAGTAAGGTAAAACCCAGTTCCAACACCAATATCACGCTGGGTGTTACCTAAATGTTca                                                                                                                                                    >  2:327132/1‑247 (MQ=255)                                                                        gcagTGAAGAAGGTAAGACATGGAAATGGAATCAAATTGACCATGTAACGCCGCGGGATAAGGTTCAAAAACATCATGGCTAATTTTATGTTTAATTTTTGATTCCCCAGCCCTTGTAGATGCCGCGTTCAGGCTAGCTTCGTTCAAATCCATTAAAGATATCAGATTACTCTCAGGTACGTGAGTAAGGTAAAACCCAGTTCCAACACCAATATCCAGATGGTTGTTACCTAAATGTTCCAGAAAGtgtg                                                                                                                                          >  2:441896/1‑251 (MQ=255)                                                                                                   ggaaTCAAATTGCCCATGTAACGCCGCGGGATAAGGTTCAAAAACATCATGGCTAATTTTATGTTTAATTTTTGATTCCCCAGCCCTTGTAGATGCCGCGTTCAGGCTAGCTTCGTTCAAATCCATTAAAGATATCAGATTACTCTCAGGTACGTGAGTAAGGTAAAACCCAGTTCCAACACCACTATCCAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACAT‑c                                                                                                              >  2:50265/1‑251 (MQ=255)                                                                                                     aaTCAAATTGACCATGTAACGCCGCGGGATAAGGTGCAAAAACATCATGGCTAATTTTATGTTTAATTTTTGATTCCCCAGCCCTTGTAGATGCCGCGTGCATGCTAGCTTCGTTCAAATCCATTAAAGATATCAGATTACTCTCAGGTACGTGAGTAAGGTAAAACCCAGTTCCAACACCAATATCCAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCGTTGTAGGACAT‑ccc                                                                                                            <  2:172516/251‑1 (MQ=255)                                                                                                         aaaTTGACCATGTAACGCCGCGGGATAAGGTTCAAAAACATCATGGCTAATTTTATGTTTAATTTTTGATTCCCCAGCCCTTGTAGATGCCGCGTTCAGGCTAGCTTCGTTCAAATCCATTAAAGATATCAGATTACTCTCAGGTACGTGAGTAAGGTAAAACCCAGTTCCAACACCAATATCCAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACAT‑CCCCATg                                                                                                        >  1:144625/1‑251 (MQ=255)                                                                                                                       aCGCCGCGGGATAAGGTTCAAAAACATCATGGCTAATTTTATGTTTAATTTTTGATTCCCCAGCCCTTGTAGATGCCGCGTTCAGGCTAGCTTCGTTCAAATCCATTAAAGATATCAGATTACTCTCAGGTACGTGAGTAAGGTAAAACCCAGTTCCAACACCAATATCCAGATGGTTGGTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACAT‑CCCCATGCAAGCCGATTTGAt                                                                                          <  1:441896/251‑1 (MQ=255)                                                                                                                         gccgcGGGATAAGGTTCAAAAACATCATGGCTACTTTTATATTTACTTTTTGATTCCCCAGCCCTTGTAAATGCCGCGTTAAGACTAGCTTAGTTAAAATCCATTAAAGATATCAGCTTACTCTCAGGTACGTGAGTAAGGTAAAACACAGTTCCAACACCAATATCCAGATGGTTGTTAACTAAATGTTCCAGAAAGTGTGGAAGATGGTGTTCCTTTGTAGGACAT‑CCCAATGCAAGCCGTTTTGATAc                                                                                        >  1:272450/1‑251 (MQ=255)                                                                                                                                  tAAGGTTCAAAAACCTCATGGCTAATTTTATGTTTAATTTTTGATTCCCCAGCCCTTGTAGATGCCGCGTTCAGGCTAGCTTCGTTCAAATCCATTAAAGATATCAGATTACTCTCAGGTACGTGAGTAAGGTAAAACCAAGTTCCAAAACAAATATCCAGATGGTTGTTACCTAAATGTCCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGAAAT‑CCCCATGCAAGCCGATTTGATACTCCCAAAAc                                                                               >  2:366707/1‑251 (MQ=255)                                                                                                                                     ggTTCAAAAACATCATGGCTAATTTTATGTTTAATTTTTGATGCCCCAGCCCTTGTAGATGCCGCGTTCAGGCTAGCTTCGTTCAAATCCATTAAAGATATCAGATTACTCTCAGGTACGTGAGTAAGGTAAAACCCAGTTCCAACACCAATATCCAGATGGTTGTTACCTAAATTTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACAT‑CCCCATGCAAGCCGATTTGATACTCCCAAAACcca                                                                            <  1:158450/251‑1 (MQ=255)                                                                                                                                           aaaaCATCATGGCTAATTTTATGTTTAATTTTTGATTCCCCAGCCCTTGTAGATGCCGCGTTCAGGCTAGCTTCGTTCAAATCCATTAAAGATATCAGATTACTCTCAGGTACGTGAGTAAGGTAAAACCCAGTTCCAACACCAATATCCAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACAT‑CCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTc                                                                      <  1:60336/251‑1 (MQ=255)                                                                                                                                              aCATCATGGCTAATTTTATGTTTAATTTTTGATTCCCCAGCCCTTGTAGATGCCGCGTTCAGGCTAGCTTCGTTCAAATCCATTAAAGATATCAGATTACTCTCAGGTACGTGCGTAAGGTAAAACCCAGTTCCAACACCAATATCCAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGCCAT‑CCCAATGCAAGCCGATTTGATACTCCAAAAACCCACCAGTCAt                                                                    >  2:184092/1‑250 (MQ=255)                                                                                                                                               catcatGGCTAATTTTATGTTTAATTTTTGATTCCCCAGCCCTTGTAGATGCCGCGTTCAGGCTAGCTTCGTTCAAATCCATTAAAGATATCAGATTACTCTCAGGTACGTGAGTAAGGTAAAACCCAGTTCCAACACCAATATCCAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACAT‑CCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATaa                                                                  <  1:206930/251‑1 (MQ=255)                                                                                                                                                 ccatGGCTAATTTTATGTTTAATTTTTGATTCCCCAGCCCTTGTAGATGCCGCGTTCAGGCTAGCTTCGTTCAAATCCATTAAAGATATCAGATTACTCTCAGGTACGTGAGTAAGGTAAAACCCAGTTCCAACACCAATATCCAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACAG‑CCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAg                                                                <  2:102867/250‑1 (MQ=255)                                                                                                                                                     ggCTAATTTTATGTTTAATTTTTGATTCCCCAGACCTTGTAGATGCCGCGTTCAGGCTAGCTTCGTTCAAATCCATTAAAGATATCAGATTACTCTCAGGTACGTGAGTAAGGTAAAACCAAGTTCCAACACCAATATCCAGAAGGTTGTTACCTAAAAGTTCCAGAAAGTGGGGAAGAAGGTGTTCCTTTGTAGG‑CATTACCCTAGCAAGCCGATTTGATAATCCCAAAACCAACCAgacctaagacttt                                                            >  1:269350/1‑239 (MQ=255)                                                                                                                                                       cTAATTTTATGTTTAATTTTTGATTCCCCAGCCCTTGTAGATGCCGCGTTAAGGCTAGCTTCGTTCAAATCCATTAAAGATATCAGATTACTCTCAGGTACGTGAGTAAGGTAAAACCCAGTTCCAACACCAATATCCAGATGGTTGTTACCTAAATGTTCAAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACAT‑CCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCttt                                                            >  2:142967/1‑249 (MQ=255)                                                                                                                                                        tAATTTTATGTTTAATTTTTGATTCCCCAGCCCTTGTAGATGCCGCGTTCAGGCTAGCTTCGTTCAAATCCATTAAAGATATCAGATTACTCTCAGGTACGTGAGTAAGGTAAAACCCAGTTCCAACACCAATATCCAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACAT‑CCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCTTTAg                                                          <  2:32530/250‑1 (MQ=255)                                                                                                                                                         aaTTTTATGTTTAATTTTTGATTCCCCAGCCCTTGTAGATGCCGCGTTCAGGCTAGCTTCGTTCAAATCCATTAAAGATATCAGATTACTCTCAGGTACGTGAGTAAGGTAAAACCCAGTTCCAACACCAATATCCAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACAT‑CCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCTTTAggg                                                        >  1:267588/1‑251 (MQ=255)                                                                                                                                                            tttATGTTTAATTTTTGATTCCCCAGCCCTTGTAGATGCCGCGTTCAGGCTAGCTTCGTTCAAATCCATTAAAGATATCAGATTACTCTCAGGTACGTGAGTAAGGTAAAACCCAGTTCCAACACCAATATCCAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCt                                                                                                                           <  1:339278/182‑1 (MQ=255)                                                                                                                                                            tttATGTTTAATTTTTGATTCCCCAGCCCTTGTAGATGCCGCGTTCAGGCTAGCTTCGTTCAAATCCATTAAAGATATCAGATTACTCTCAGGTACGTGAGTAAGGTAAAACCCAGTTCCAACACCAATATCCAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCt                                                                                                                           >  2:339278/1‑182 (MQ=255)                                                                                                                                                                     aaTTTTTGATTCCCCAGCCCTTGTAGATGGCGCGTTCAGGCTAGCTTCGTTCAAATCCATTAAAGATATCAAATTACTATCAGGAACTTGAGTAAGGTAAAATGCAGTTCCAACACCAATATCCAGATGGTTGATACCTAAATGTTCCAGAAAGTGTGTAAGACGGTGTTCATTTGTAGGACAT‑CCCCATGCAAGCCGATTTGATACTCCCAAATCCCACCTGTCATAAAGCTTTAGGGTAAGTGgtgtgt                                            >  2:115012/1‑251 (MQ=255)                                                                                                                                                                         tttGATTCCCCAGCCCTTGTAGATGCCGCGTTCAGGCTAGCTTCGTTCAAATCCATTAAAGATATCAGATTACTCTCAGGTACGTGAGTAAGGTAAAACCCAGTTCCAACACCAATATCCAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACAT‑CCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGt                                                                       >  1:63999/1‑220 (MQ=255)                                                                                                                                                                         tttGATTCCCCAGCCCTTGTAGATGCCGCGTTCAGGCTAGCTTCGTTCAAATCCATTAAAGATATCAGATTACTCTCAGGTACGTGAGTAAGGTAAAACCCAGTTCCAACACCAATATCCAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACAT‑CCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGt                                                                       <  2:63999/220‑1 (MQ=255)                                                                                                                                                                              ttCCCCAGCCCTTGTAGATGCCGCGTTCAGGCTAGCTTCGTTCAAATCCATTAAAGATATCAGATTACTCTCAGGTACGTGAGTAAGGTAAAACCCAGTTCCAACACCAATATCCAGATGGTTGTAACCTAAATGTTCCAGAAAGTGTGGAAAAAGGTGTTCCTTTGTAGGACTT‑CCCCATGCAAGCCGATTTGATACTCCCAAACCCCACCAGCCATAAAGCTTTAGGGTAAGTGGTGAGTACATTcc                                     >  2:198273/1‑248 (MQ=255)                                                                                                                                                                                             aGATGCCGCTTTCAGGCTAGCTTCGTTCAAATCCATTAAAGATATCAGATTACTCTCAGGTACGTGAGTAAGGTAAAACCCAGTTCCAACACCAATATCCAGATGGTTGTTACCGAAATGTGCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACAT‑CCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTAAATTCTAGCCCCATCATCTgtgt                    <  1:390225/251‑1 (MQ=255)                                                                                                                                                                                                                gCTTCGTTCAAATCCATTAAAGATATCAGATTACTCTCAGGTACGTGAGTACGGTAAAACCCAGTTCCAACACCAATATCCAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACAT‑CCCCATGCAAGCCGATTTGATACTCCCAAAACACACAAGTCATAAAGCTTTAGGGTAAGTGGTGTGTAAATTCTAGCCCCATCATCTGTGTTTTTTTATTAATTTCAcc  >  1:317847/1‑250 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                             GTTCCATCGTCAGTTAAGGCCTGCGCCGCATTGCGTATTACACAGCTTTTTGTAGATATATTTCCAGGCAGGCAGTGAAGAAGGTAAGACATGGAAATGGAATCAAATTGACCATGTAACGCCGCGGGATAAGGTTCAAAAACATCATGGCTAATTTTATGTTTAATTTTTGATTCCCCAGCCCTTGTAGATGCCGCGTTCAGGCTAGCTTCGTTCAAATCCATTAAAGATATCAGACTACTCTCAGGTACGTGAGTAAGGTAAAACCCAGTTCCAACACCAATATCCAGATGGTTGTTACCTAAATGTTCCAGAAAGTGTGGAAGAAGGTGTTCCTTTGTAGGACAT‑CCCCATGCAAGCCGATTTGATACTCCCAAAACCCACCAGTCATAAAGCTTTAGGGTAAGTGGTGTGTAAATTCTAGCCCCATCATCTGTGTTTTTTTATTAATTTCACC  >  NC_000913/582360‑582816

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: ‑