Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,208,928	A→G	R64R (CGT→CGC)	tfaE ←	e14 prophage; putative tail fiber assembly protein TfaE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,208,928	0	A	G	100.0%	121.5 / NA	34	R64R (CGT→CGC)	tfaE	e14 prophage; putative tail fiber assembly protein TfaE
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (17/17);  total (17/17)

AAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGTGGAAATATATTCCCGTGTTTCACCATCGTAGTTAT  >  NC_000913/1208795‑1209058                                                                                                                                      |                                                                                                                                   aaGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAGCGGCAg                                                                                                                              <  2:171118/140‑1 (MQ=255)        gATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAGCGGCAGATTGCAt                                                                                                                       >  2:31409/1‑140 (MQ=255)           aGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAGCGGCAGATTGCATAAc                                                                                                                    <  2:67874/140‑1 (MQ=255)           aGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAGCGGCAGATTGCATAAc                                                                                                                    >  1:65402/1‑140 (MQ=255)                   tATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAGCGGCAGATTGCATaa                                                                                                                     >  1:315011/1‑131 (MQ=255)                   tATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAGCGGCAGATTGCATaa                                                                                                                     <  2:315011/131‑1 (MQ=255)                       ttCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAGCGGCAGATTGCATAACCAGCCTTAtgtg                                                                                                        >  2:384169/1‑140 (MQ=255)                           gggTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTGTGGCACATATTCCCATGAGTTAAAATCTGCAGAGCGGCAGATTGCATAACCAGCCTTATGTGTAcc                                                                                                    <  1:384169/140‑1 (MQ=255)                           gggTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAGCGGCAGATTGCATAACCAGCCTTAt                                                                                                           >  2:50596/1‑133 (MQ=255)                           gggTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAGCGGCAGATTGCATAACCAGCCTTAt                                                                                                           <  1:50596/133‑1 (MQ=255)                                    cccGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAGCGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATc                                                                                           <  1:152194/140‑1 (MQ=255)                                       ggAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAGCGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTaa                                                                                        >  1:345052/1‑140 (MQ=255)                                       ggAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAGCGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTaa                                                                                        >  1:154399/1‑140 (MQ=255)                                         aGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAGCGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAAc                                                                                      <  1:1948/140‑1 (MQ=255)                                          gCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAGCGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACa                                                                                     >  2:191400/1‑140 (MQ=255)                                            tttGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAGCGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAgg                                                                                   >  1:277765/1‑140 (MQ=255)                                             ttGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAGCGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAGGa                                                                                  <  2:345052/140‑1 (MQ=255)                                               gATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAGCGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAGGAAc                                                                                >  1:307997/1‑140 (MQ=255)                                                               ccAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAGCGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAGGAACATGCCGGAATGCCGac                                                                >  1:253819/1‑140 (MQ=255)                                                                        gTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAGCGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGc                                                       >  1:361380/1‑140 (MQ=255)                                                                           cTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAGCGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAg                                                    >  2:347312/1‑140 (MQ=255)                                                                           cTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAGCGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAg                                                    <  2:277765/140‑1 (MQ=255)                                                                                      ttCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAGCGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTg                                         <  1:347312/140‑1 (MQ=255)                                                                                        cACCGCGATGTTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAGCGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGaa                                       >  2:270163/1‑140 (MQ=255)                                                                                        cACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAGCGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGaa                                       <  2:154399/140‑1 (MQ=255)                                                                                         aCCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAGCGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAg                                      <  2:307997/140‑1 (MQ=255)                                                                                         aCCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAGCGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAg                                      <  2:65402/140‑1 (MQ=255)                                                                                         aCCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAGCGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAg                                      <  1:171320/140‑1 (MQ=255)                                                                                              gATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAGCGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGTGGaa                                 <  2:253819/140‑1 (MQ=255)                                                                                                tGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAGCGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGTGGAAat                               >  1:380072/1‑140 (MQ=255)                                                                                                       gCACATATTCCCATGAGTTAAAATCTGCAGAGCGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGTGGAAATATATTcc                        <  2:361380/140‑1 (MQ=255)                                                                                                                      aGTTAAAATCTGCAGAGCGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGTGGAAATATATTCCCGTGTTTCACCATCg         >  2:294076/1‑140 (MQ=255)                                                                                                                        ttAAAATCTGCAGAGCGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGTGGAAATATATTCCCGTGTTTCACCATCGTa       >  1:329976/1‑140 (MQ=255)                                                                                                                             aTCTGCAGAGCGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGTGGAAATATATTCCCGTGTTTCACCATCGTAGTtat  <  2:329976/140‑1 (MQ=255)                                                                                                                                      |                                                                                                                                   AAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAGGAACATGCCGGAATGCCGACACCAACGGCAAGATATTCATTTGAAGTGGAAATATATTCCCGTGTTTCACCATCGTAGTTAT  >  NC_000913/1208795‑1209058

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑