Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,913,560	G→T	V419V (GTC→GTA)	prc ←	tail‑specific protease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,913,560	0	G	T	100.0%	91.4 / NA	28	V419V (GTC→GTA)	prc	tail‑specific protease
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (14/14);  total (14/14)

CCGTAATCCTGCATTGCCGCGGCAAAGATTTCTGAAGCCGAAGCACTGAAGCGGTCAACCAGCACCACCAGCGGGCCTTTATAGAAAACCTGTCCGTCGGTATCGCTATCTTCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATAAACAGACCGGAGAGCGATACGGCTTCAGTTAACGCCCCACCGCCATTGCTACGCAGGTCGATGATGACGCTGCTGACATTCTGTTTTTCCAGTTTCTGCAGTT  >  NC_000913/1913423‑1913690                                                                                                                                          |                                                                                                                                   ccGTAATCCTGCATTGCCGCGGCAAAGATTTCTGAAGCCGAAGCACTGAAGCGGTCAACCAGCACCACCAGCGGGCCTTTATAGAAAACCTGTCCGTCGGTATCGCTATCTTCACGAACCTTGCCGTTGTTATCGCGTAc                                                                                                                                  <  1:308332/140‑1 (MQ=255)  cGTAATCCTGCATTGCCGCGGCAAAGATTTCTGAAGCCGAAGCACTGAAGCGGTCAACCAGCACCACCAGCGGGCCTTTATAGAAAACCTGTCCGTCGGTATCGCTATCTTCACGAACCTTGCCGTTGTTATCGCGTAcc                                                                                                                                 >  2:207911/1‑140 (MQ=255)         cTGCATTGCCGCGGCAAAGATTTCTGAAGCCGAAGCACTGAAGCGGTCAACCAGCACCACCAGCGGGCCTTTATAGAAAACCTGTCCGTCGGTATCGCTATCTTCACGAACCTTGCCGTTGTTATCGCGTACCTGAACaa                                                                                                                          >  1:81067/1‑140 (MQ=255)           gCATTGCCGCGGCAAAGATTTCTGAAGCCGAAGCACTGAAGCGGTCAACCAGCACCACCAGCGGGCCTTTATAGAAAACCTGTCCGTCGGTATCGCTATCTTCACGAACCTTGCCGTTGTTATCGCGTACCTGAACAATg                                                                                                                        >  2:165772/1‑140 (MQ=255)           gCATTGCCGCGGCAAAGATTTCTGAAGCCGAAGCACTGAAGCGGTCAACCAGCACCACCAGCGGGCCTTTATAGAAAACCTGTCCGTCGGTATCGCTATCTTCACGAACCTTGCCGTTGTTATCGCGTACCTGAACAATg                                                                                                                        >  1:363350/1‑140 (MQ=255)             aTTGCCGCGGCAAAGATTTCTGAAGCCGAAGCACTGAAGCGGTCAACCAGCACCACCAGCGGGCCTTTATAGAAAACCTGTCCGTCGGTATCGCTATCTTCACGAACCTTGCCGTTGTTATCGCGTACCTGAACAATggg                                                                                                                      >  2:317517/1‑140 (MQ=255)                   gcgGCAAAGATTTCTGAAGCCGAAGCACTGAAGCGGTCAACCAGCACCACCAGCGGGCCTTTATAGAAAACCTGTCCGTCGGTATCGCTATCTTCACGAACCTTGCCGTTGTTATCGCGTACCTGAACAATGGGACCCGc                                                                                                                >  2:162002/1‑140 (MQ=255)                    cgGCAAAGATTTCTGAAGCCGAAGCACTGAAGCGGTCAACCAGCACCACCAGCGGGCCTTTATAGAAAACCTGTCCGTCGGTATCGCTATCTTCACGAACCTTGCCGTTGTTATCGCGTACCTGAACAATGGGACCCGCa                                                                                                               <  1:351551/140‑1 (MQ=255)                              ttCTGAAGCCGAAGCACTGAAGCGGTCAACCAGCACCACCAGCGGGCCTTTATAGAAAACCTGTCCGTCGGTATCGCTATCTTCACGAACCTTGCCGTTGTTATCGCGTACCTGAACAATGGGACCCGCAGGAATAAACa                                                                                                     <  2:81067/140‑1 (MQ=255)                                             aCTGAAGCGGTCAACCAGCACCACCAGCGGGCCTTTATAGAAAACCTGTCCGTCGGTATCGCTATCTTCACGAACCTTGCCGTTGTTATCGCGTACCTGAACAATGGGACCCGCAGGAATAAACAGACCGGAGAGCGATa                                                                                      <  2:57516/140‑1 (MQ=255)                                                  aGCGGTCAACCAGCACCACCAGCGGGCCTTTATAGAAAACCTGTCCGTCGGTATCGCTATCTTCACGAACCTTGCCGTTGTTATCGCGTACCTGAACAATGGGACCCGCAGGAATAAACAGACCGGAGAGCGATACGGCt                                                                                 <  2:372275/140‑1 (MQ=255)                                                   gCGGTCAACCAGCACCACCAGCGGGCCTTTATAGAAAACCTGTCCGTCGGTATCGCTATCTTCACGAACCTTGCCGTTGTTATCGCGTACCTGAACAATGGGACCCGCAGGAATAAACAGACCGGAGAGCGATACGGCtt                                                                                >  2:48263/1‑140 (MQ=255)                                                      gTCAACCAGCACCACCAGCGGGCCTTTATAGAAAACCTGTCCGTCGGTATCGCTATCTTCACGAACCTTGCCGTTGTTATCGCGTACCTGAACAATGGGACCCGCAGGAATAAACAGACCGGAGAGCGATACGGCTTCAg                                                                             >  1:278173/1‑140 (MQ=255)                                                      gTCAACCAGCACCACCAGCGGGCCTTTATAGAAAACCTGTCCGTCGGTATCGCTATCTTCACGAACCTTGCCGTTGTTATCGCGTACCTGAACAATGGGACCCGCAGGAATAAACAGACCGGAGAGCGATACGGCTTCAg                                                                             >  1:246609/1‑140 (MQ=255)                                                      gTCAACCAGCACCACCAGCGGGCCTTTATAGAAAACCTGTCCGTCGGTATCGCTATCTTCACGAACCTTGCCGTTGTTATCGCGTACCTGAACAATGGGACCCGCAGGAATAAACAGACCGGAGAGCGATACGGCTTCAg                                                                             >  2:243061/1‑140 (MQ=255)                                                            cAGCACCACCAGCGGGCCTTTATAGAAAACCTGTCCGTCGGTATCGCTATCTTCACGAACCTTGCCGTTGTTATCGCGTACCTGAACAATGGGACCCGCAGGAATAAACAGACCGGAGAGCGATACGGCTTCAGTTAACg                                                                       <  1:279904/140‑1 (MQ=255)                                                                   accaGCGGGCCTTTATAGAAAACCTGTCCGTCGGTATCGCTATCTTCACGAACCTTGCCGTTGTTATCGCGTACCTGAACAATGGGACCCGCAGGAATAAACAGACCGGAGAGCGATACGGCTTCAGTTAACGCCCCAcc                                                                >  2:171007/1‑140 (MQ=255)                                                                     caGCGGGCCTTTATAGAAAACCTGTCCGTCGGTATCGCTATCTTCACGAACCTTGCCGTTGTTATCGCGTACCTGAACAATGGGACCCGCAGGAATAAACAGACCGGAGAGCGATACGGCTTCAGTTaa                                                                         >  1:412227/1‑129 (MQ=255)                                                                     caGCGGGCCTTTATAGAAAACCTGTCCGTCGGTATCGCTATCTTCACGAACCTTGCCGTTGTTATCGCGTACCTGAACAATGGGACCCGCAGGAATAAACAGACCGGAGAGCGATACGGCTTCAGTTaa                                                                         <  2:412227/129‑1 (MQ=255)                                                                     caGCGGGCCTTTATAGAAAACCTGTCCGTCGGTATCGCTATCTTCACGAACCTTGCCGTTGTTATCGCGTACCTGAACAATGGGACCCGCAGGAATAAACAGACCGGAGAGCGATACGGCTTCAGTTAACGCCCCACCGc                                                              <  1:48263/140‑1 (MQ=255)                                                                               ttATAGAAAACCTGTCCGTCGGTATCGCTATCTTCACGAACCTTGCCGTTGTTATCGCGTACCTGAACAATGGGACCCGCAGGAATAAACAGACCGGAGAGCGATACGGCTTCAGTTAACGCCCCACCGCCATTGCTACg                                                    >  1:128253/1‑140 (MQ=255)                                                                               ttATAGAAAACCTGTCCGTCGGTATCGCTATCTTCACGAACCTTGCCGTTGTTATCGCGTACCTGAACAATGGGACCCGCAGGAATAAACAGACCGGAGAGCGATACGGCTTCAGTTAACGCCCCACCGCCATTGCTACg                                                    <  1:171007/140‑1 (MQ=255)                                                                                tataGAAAACCTGTCCGTCGGTATTGCTTTCTTCACGAACCTTGCCGTTGTTATCGCGTCCCTGAACAATGGGACCCGCAGGAATAAACAGACCGGAGAGCGATACGGCTTCAGTTAACGCCCCACCGCCATTGCTACGc                                                   <  2:176005/140‑1 (MQ=255)                                                                                tataGAAAACCTGTCCGTCGGTATCGCTATCTTCACGAACCTTGCCGTTGTTATCGCGTACCTGAACAATGGGACCCGCAGGAATAAACAGACCGGAGAGCGATACGGCTTCAGTTAACGCCCCACCGCCATTGCTACGc                                                   <  1:59677/140‑1 (MQ=255)                                                                                          cTGTCCGTCGGTATCGCTATCTTCACGAACCTTGCCGTTGTTATCGCGTACCTGAACAATGGGACCCGCAGGAATAAACAGACCGGAGAGCGATACGGCTTCAGTTAACGCCCCACCGCCATTGCTACGCAGGTCgatga                                         >  1:31504/1‑140 (MQ=255)                                                                                               cGTCGGTATCGCTATCTTCACGAACCTTGCCGTTGTTATCGCGTACCTGAACAATGGGACCCGCAGGAATAAACAGACCGGAGAGCGATACGGCTTCAGTTAACGCCCCACCGCCATTGCTACGCAGGTCGATGATGACg                                    <  1:326695/140‑1 (MQ=255)                                                                                                gTCGGTATCGCTATCTTCACGAACCTTGCCGTTGTTATCGCGTACCTGAACAATGGGACCCGCAGGAATAAACAGACCGGAGAGCGATACGGCTTCAGTTAACGCCCCACCGCCATTGCTACGCAGGTCGATGATGACgc                                   <  2:120279/140‑1 (MQ=255)                                                                                                                                tgttATCGCGTACCTGAACAATGGGACCCGCAGGAATAAACAGACCGGAGAGCGATACGGCTTCAGTTAACGCCCCGCCGCCATTGCTACGCAGGTCGATGATGACGCTGCTGACATTCTGTTTTTCCAGTTTCTGCAGt   <  2:246609/140‑1 (MQ=255)                                                                                                                                 gttATCGCGTACCTGAACAATGGGACCCGCAGGAATAAACAGACCGGAGAGCGATACGGCTTCAGTTAACGCCCCACCGCCATTGCTACGCAGGTCGATGATGACGCTGCTGACATTCTGTTTTTCCAGTTTCTGCAGtt  <  1:165772/140‑1 (MQ=255)                                                                                                                                          |                                                                                                                                   CCGTAATCCTGCATTGCCGCGGCAAAGATTTCTGAAGCCGAAGCACTGAAGCGGTCAACCAGCACCACCAGCGGGCCTTTATAGAAAACCTGTCCGTCGGTATCGCTATCTTCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATAAACAGACCGGAGAGCGATACGGCTTCAGTTAACGCCCCACCGCCATTGCTACGCAGGTCGATGATGACGCTGCTGACATTCTGTTTTTCCAGTTTCTGCAGTT  >  NC_000913/1913423‑1913690

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