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%	109.8 / NA	49	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 (23/26);  total (23/26)

CCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTGGCAGAC  >  NC_000913/1634511‑1634776                                                                                                                                    |                                                                                                                                       ccGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCa                                                                                                                                <  1:141724/140‑1 (MQ=18)                     ccGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTAccc                                                                                                            <  2:300680/140‑1 (MQ=18)                              tCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCt                                                                                                   <  2:232838/140‑1 (MQ=17)                               cACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGttt                                                                                                                 <  2:891681/125‑1 (MQ=18)                               cACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGttt                                                                                                                 >  1:891681/1‑125 (MQ=17)                               cACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCtt                                                                                                  >  1:14225/1‑140 (MQ=17)                               cACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCtt                                                                                                  >  1:709046/1‑140 (MQ=17)                                      gCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACgag                                                                                           >  1:938101/1‑140 (MQ=17)                                      gCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACgag                                                                                           >  1:335064/1‑140 (MQ=17)                                      gCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACgag                                                                                           <  1:832284/140‑1 (MQ=17)                                          tgCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCgat                                                                                                        >  2:902143/1‑123 (MQ=17)                                          tgCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCgat                                                                                                        <  1:902143/123‑1 (MQ=17)                                          tgCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGc                                                                                       <  2:938101/140‑1 (MQ=17)                                          tgCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGc                                                                                       >  2:286659/1‑140 (MQ=17)                                          tgCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGc                                                                                       >  1:391263/1‑140 (MQ=17)                                           gCCGTTCCACTTCTGAAACTCCCCTACCGGCGATATCAAGGTAACCTTTTTCGGTAGCATACGGAGTTCAGAAATAAAGAACGCGTCGATGGAAGAACCGTCATAAAGCGTTTTACCCCGATGACCTTCAGCGAGATGcg                                                                                      <  2:697480/140‑2 (MQ=2)                                            ccGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCa                                                                                     >  1:894457/1‑140 (MQ=17)                                                tCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGAt                                                                                 >  2:870823/1‑140 (MQ=17)                                                       cTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACtcat                                                                          <  2:709046/140‑1 (MQ=17)                                                       cTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACtcat                                                                          <  2:736038/140‑1 (MQ=17)                                                         gAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACtcatca                                                                        <  1:592779/140‑1 (MQ=17)                                                              tCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGtt                                                                   <  2:894457/140‑1 (MQ=17)                                                               ccccTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACtcatca                                                                        <  1:169120/134‑1 (MQ=17)                                                               ccccTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACtcatca                                                                        >  2:169120/1‑134 (MQ=17)                                                                 ccTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGa                                                                               <  2:346613/125‑1 (MQ=17)                                                                 ccTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGa                                                                               >  1:346613/1‑125 (MQ=17)                                                                  cTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACtcat                                                                          >  2:869706/1‑129 (MQ=17)                                                                  cTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACtcat                                                                          <  1:869706/129‑1 (MQ=17)                                                                     ccGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCCCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACa                                                            >  2:28417/1‑140 (MQ=17)                                                                              aCCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAag                                                   <  1:606566/140‑1 (MQ=17)                                                                                 aGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAagcca                                                <  2:391263/140‑1 (MQ=17)                                                                                  ggTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAagccag                                               >  1:734085/1‑140 (MQ=17)                                                                                   gTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAagccagc                                              <  2:14225/140‑1 (MQ=17)                                                                                      aCATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCCTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCgg                                           >  2:526264/1‑140 (MQ=17)                                                                                             ccGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCt                                    >  1:103704/1‑140 (MQ=17)                                                                                               ggTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTgg                                  <  1:870823/140‑1 (MQ=17)                                                                                                  aaCGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCgg                               >  2:562598/1‑140 (MQ=17)                                                                                                   aCGGACCGAGTTCAGAAATAAATAACGCGTCCCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGt                              <  1:28417/140‑1 (MQ=17)                                                                                                      gACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCg                                            <  2:157955/123‑1 (MQ=17)                                                                                                      gACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCg                                            >  1:157955/1‑123 (MQ=18)                                                                                                        ccGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAat                         <  2:335064/140‑1 (MQ=17)                                                                                                           aGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAatat                       >  1:250195/1‑139 (MQ=18)                                                                                                           aGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAatat                       <  2:250195/139‑1 (MQ=17)                                                                                                             ttCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCaa                          <  2:588081/134‑1 (MQ=17)                                                                                                             ttCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCaa                          >  1:588081/1‑134 (MQ=18)                                                                                                                     taaataaCGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGttt            >  2:792146/1‑140 (MQ=18)                                                                                                                          aaCGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATATTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTgg       <  1:277623/140‑1 (MQ=17)                                                                                                                              cgTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTAc                 >  1:272589/1‑126 (MQ=17)                                                                                                                              cgTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTAc                 <  2:272589/126‑1 (MQ=17)                                                                                                                               gTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTGGCAGAc  <  1:562598/140‑1 (MQ=17)                                                                                                                                    |                                                                                                                                       CCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTGGCAGAC  >  NC_000913/1634511‑1634776

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