Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,634,642	C→T	G82S (GGC→AGC)	tfaQ ←	Qin prophage; putative tail fiber assembly protein TfaQ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,634,642	0	C	T	100.0%	100.7 / NA	46	G82S (GGC→AGC)	tfaQ	Qin prophage; putative tail fiber assembly protein TfaQ
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (26/20);  total (26/20)

TCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACA‑TTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTGGCAGACCT  >  NC_000913/1634510‑1634778                                                                                                                                      |                                                                                                                                         tCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACA‑TTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGcc                                                                                                                                   <  2:173137/140‑1 (MQ=17)     gAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACA‑TTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGt                                                                                                                               >  2:341407/1‑140 (MQ=17)     gAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTAAAACTCCCCTTCCGGCGATAACCAGGTAACA‑TTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGt                                                                                                                               >  2:340827/1‑140 (MQ=17)      aaCAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACA‑TTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTc                                                                                                                              >  2:456385/1‑140 (MQ=17)        cAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACA‑TTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCAt                                                                                                                            >  2:193955/1‑140 (MQ=17)                           tCCTTCACCCAGGCTGTGCCGCTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACA‑TTTTCCGGTTACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCgatg                                                                                                         <  2:786880/140‑1 (MQ=12)                                cACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACA‑TTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCtt                                                                                                    >  1:615875/1‑140 (MQ=17)                                     aGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACA‑TTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACg                                                                                               >  1:624144/1‑140 (MQ=17)                                      ggCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACA‑TTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACga                                                                                              >  2:454559/1‑140 (MQ=17)                                        cTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACA‑TTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCgatga                                                                                                        >  1:369355/1‑128 (MQ=17)                                        cTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACA‑TTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCgatga                                                                                                        <  2:369355/128‑1 (MQ=18)                                          gtgCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACA‑TTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATg                                                                                          >  2:133467/1‑140 (MQ=17)                                           tgCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACA‑TTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGc                                                                                         >  1:260483/1‑140 (MQ=17)                                            gCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACA‑TTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGcc                                                                                        <  2:260483/140‑1 (MQ=17)                                            gCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACA‑TTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGcc                                                                                        >  1:579203/1‑140 (MQ=17)                                                ttCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACA‑TTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGa                                                                                    <  2:615875/140‑1 (MQ=17)                                                 tCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACA‑TTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGAt                                                                                   >  2:378215/1‑140 (MQ=17)                                                   cACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACA‑TTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGa                                                                                 <  1:632816/140‑1 (MQ=17)                                                   cACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACA‑TTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGa                                                                                 <  2:557339/140‑1 (MQ=17)                                                     cTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACA‑TTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACt                                                                               >  1:50522/1‑140 (MQ=17)                                                     cTTCTGAAACTCCCCCTCCGGCGATAACCAGGTAACA‑TTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACt                                                                               >  1:377966/1‑140 (MQ=12)                                                          gAAACTCCCCTTCCGGCGATAACCAGGTAACA‑TTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACtcatca                                                                          <  1:193955/140‑1 (MQ=17)                                                          gAAACTCCCCTTCCGGCGATAACCAGGTAACA‑TTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACtcatca                                                                          <  1:378215/140‑1 (MQ=17)                                                          gAAACTCCCCCTCCGGCGATAACCAGGTAACA‑TTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACtcatca                                                                          <  2:377966/140‑1 (MQ=17)                                                           aaaCTCCCCTTCCGGCGATAACCAGGTAACATTTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACAATGACtcatca                                                                          <  2:470937/140‑1 (MQ=255)                                                             aCTCCCCTTCCGGCGATAACCAGGTAACA‑TTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTg                                                                       >  1:48894/1‑140 (MQ=17)                                                                  ccTTCCGGCGATAACCAGGTAACA‑TTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGaa                                                                  >  2:717028/1‑140 (MQ=17)                                                                  ccTTCCGGCGATAACCAGGTAACA‑TTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGaa                                                                  >  1:244449/1‑140 (MQ=17)                                                                                   ggTAACA‑TTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAagccag                                                 >  2:571142/1‑140 (MQ=17)                                                                                   ggTAACA‑TTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAagccag                                                 <  2:624144/140‑1 (MQ=17)                                                                                              tCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATc                                       >  2:550202/1‑140 (MQ=17)                                                                                               ccGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCt                                      >  1:452287/1‑140 (MQ=17)                                                                                                       ggACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGc                              <  2:244449/140‑1 (MQ=17)                                                                                                       ggACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGc                              <  1:550202/140‑1 (MQ=17)                                                                                                               ttCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCgg                      >  1:206025/1‑140 (MQ=17)                                                                                                                            aaCGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTgg         <  2:366512/140‑1 (MQ=17)                                                                                                                            aaCGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTgg         <  2:822485/140‑1 (MQ=17)                                                                                                                             aCGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTGGc        <  1:340827/140‑1 (MQ=17)                                                                                                                              cgcgTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTGGCa       >  2:631048/1‑140 (MQ=17)                                                                                                                              cgcgTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTGGCa       >  2:675423/1‑140 (MQ=18)                                                                                                                                cgTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTGGCAGa     >  2:455181/1‑140 (MQ=17)                                                                                                                                 gTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGttt              >  2:238676/1‑130 (MQ=14)                                                                                                                                 gTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGttt              <  1:238676/130‑1 (MQ=17)                                                                                                                                  tCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTg          <  1:793575/133‑1 (MQ=14)                                                                                                                                  tCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTg          >  2:793575/1‑133 (MQ=11)                                                                                                                                   cGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTGGCAGACCt  <  2:452287/140‑1 (MQ=11)                                                                                                                                      |                                                                                                                                         TCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACA‑TTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTGGCAGACCT  >  NC_000913/1634510‑1634778

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