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%	82.9 / NA	38	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 (18/20);  total (18/20)

CCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTGGCAGACC  >  NC_000913/1634511‑1634777                                                                                                                                    |                                                                                                                                        ccGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCa                                                                                                                                 >  1:643033/1‑140 (MQ=17)          tttttCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCACGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATaaa                                                                                                                        <  2:954174/140‑1 (MQ=16)                    tCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTAcc                                                                                                              >  1:296187/1‑140 (MQ=17)                         aTCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCgat                                                                                                         <  2:269609/140‑1 (MQ=17)                           ccTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCgatga                                                                                                       <  1:234944/140‑1 (MQ=17)                            cTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCgatgat                                                                                                      <  2:935091/140‑1 (MQ=18)                                     ggCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACga                                                                                             >  1:797527/1‑140 (MQ=17)                                          tgCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCgat                                                                                                         <  2:87069/123‑1 (MQ=18)                                          tgCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCgat                                                                                                         >  1:87069/1‑123 (MQ=17)                                               ttCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGa                                                                                   <  2:220614/140‑1 (MQ=17)                                                tCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGAt                                                                                  >  1:255990/1‑140 (MQ=17)                                                tCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGAt                                                                                  >  1:710915/1‑140 (MQ=17)                                                    cTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACt                                                                              >  1:786213/1‑140 (MQ=17)                                                         gAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACtcatca                                                                         <  2:255990/140‑1 (MQ=17)                                                         gAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACtcatca                                                                         <  2:638424/140‑1 (MQ=17)                                                               ccccTTCCGTCGATAACCAGGTATCCTTTTCCGGTAACGGACCGAGTTCAGAACTAAAGAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTg                                                                   >  2:361482/1‑140 (MQ=11)                                                               ccccTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTg                                                                   >  2:138544/1‑140 (MQ=17)                                                               ccccTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTg                                                                   >  2:860020/1‑140 (MQ=17)                                                                 ccTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGaa                                                                 >  2:23192/1‑140 (MQ=17)                                                                          gATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCAc                                                        >  1:524097/1‑140 (MQ=17)                                                                               ccAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAagc                                                   <  1:96942/140‑1 (MQ=17)                                                                               ccAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAagc                                                   <  2:747703/140‑1 (MQ=17)                                                                                   gTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAagccagc                                               <  1:779889/140‑1 (MQ=17)                                                                                     aaCATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCg                                             <  1:23192/140‑1 (MQ=17)                                                                                       cATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAACCACGAAGCCAGCCGGa                                           >  2:913773/1‑140 (MQ=17)                                                                                            tCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATc                                      >  1:690593/1‑140 (MQ=17)                                                                                             ccGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCt                                     >  1:919444/1‑140 (MQ=17)                                                                                               ggTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTgg                                   <  1:87193/140‑1 (MQ=17)                                                                                               ggTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAACCACGAAGCCAGCCGGAATATCTgg                                   <  1:913773/140‑1 (MQ=17)                                                                                                 tAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCg                                 >  1:872949/1‑140 (MQ=18)                                                                                                 tAACGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCg                                 >  1:631790/1‑140 (MQ=18)                                                                                                  aaCGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCgg                                <  2:786213/140‑1 (MQ=17)                                                                                                   aCGGACCGAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGt                               <  2:872949/140‑1 (MQ=17)                                                                                                          gAGTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAatat                        <  2:710915/140‑1 (MQ=17)                                                                                                            gTTCAGAAATAAATAACGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCg                      <  2:797527/140‑1 (MQ=17)                                                                                                                     taaataaCGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGttt             <  1:138544/140‑1 (MQ=17)                                                                                                                          aaCGCGTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTgg        <  2:631790/140‑1 (MQ=17)                                                                                                                            cgcgTCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTGGCa      >  1:371154/1‑140 (MQ=18)                                                                                                                                tCGCTGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTGGCAGAcc  >  1:264075/1‑140 (MQ=11)                                                                                                                                    |                                                                                                                                        CCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATCACTGTTGAAAACAGCCACGAAGCCAGCCGGAATATCTGGCGGTGCAATATCGGTACTGTTTGCTGGCAGACC  >  NC_000913/1634511‑1634777

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