Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,913,579	G→A	18.2%	A413V (GCG→GTG)	prc ←	tail‑specific protease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,913,579	0	G	A	18.2%	30.4 / 3.2	23	A413V (GCG→GTG)	prc	tail‑specific protease
Reads supporting (aligned to +/- strand):  ref base G (8/10);  new base A (2/2);  total (10/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.83e-01

AAGATTTCTGAAGCCGAAGCACTGAAGCGGTCAACCAGCACCACCAGCGGGCCTTTATAGAAAACCTGTCCGTCGGTATCGCTATCTTCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATAAACAGACCGGAGAGCGATACGGCTTCAGTTAACGCCCCACCGCCATTGCTACGCAGGTCGATGATGACGCTGCTGACATTCTGTTTTTCCAGTTTCTGCAGTTGCACTTTGACATCGTCTGTCAAACCCACATA  >  NC_000913/1913447‑1913721                                                                                                                                     |                                                                                                                                               aaGATTTCTGAAGCCGAAGCAGTGAAGCGGTCAACCAGCACCACCAGTGGGTCTTTATAGAAAACCGGTCCGTCGGTATCGCTATCTTCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATaaa                                                                                                                                      <  1:208513/143‑1 (MQ=255)     tttCTGGAGCCTGCGCACTGCAGCGGTCAACCAGCACCACCCGCGGGCCTGTATAGAAAACCTGTCCGTCGGTATCGCTATCTTCACGAGCCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATAAACAGa                                                                                                                                  <  2:235388/143‑1 (MQ=255)          gAAGCCGAAGCACTGAAGCGGTCAACCAGAACCACCAGCGGGCCTTTATAGAAAACCTGTCCGTCGGTATCGCTGTCTTCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATAAACAGAACggc                                                                                                                             >  2:372262/1‑142 (MQ=255)                  aGCACTGAAGCGGTCAACCAGCACCACCAGCGGGCCTTTATAGAAAACCTGTCCGTCGGTATCGCTATCTTCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATAAACAGACCGGAGAGCGATa                                                                                                                     <  2:68964/143‑1 (MQ=255)                    cACTGAAGCGGTCAACCAGCACCACCAGCGGGCCATAATAGAAAACCTGTCCGTCGGTAGAGCTATCAGCACGCACCTTGACGTTGTTATCGGGGACCTGAACAATTGGACACGCAGGGAAAAACAGACGGTTGAGCGATaag                                                                                                                   >  2:215515/1‑141 (MQ=255)                      cTGAAGCGGTCAACCAGCCCCCCCAGCGGGCCTTTATAGAAAGCCTGTCCGTCGGTATCGCTATCTTCACGACCCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATAAACAGACCGGAGAGCGATACGGc                                                                                                                 <  2:388910/143‑1 (MQ=255)                                 aaCCAGCACCACCAGCGGGCCTTTATAGAAAACCTGTCCGTCGGTATCGCTATCTTCACGAACCTTGCAGTTGTTATCACGGACCTGAACAATGGGACCCGCAGGAATAAACAGACCGGAGAGCGATACGGCTTCAGTTAACg                                                                                                      >  1:305704/1‑143 (MQ=255)                                      gCACCACCAGCGTGCCTTTATTGATAACCTGTCCGTCGGTATCTCTATCTTCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATAAACAGACCGGAGAGCGATACGGCTTCAGTTAACGCCCCa                                                                                                 <  1:182653/143‑1 (MQ=255)                                      gCACCACCAGCGGGCCTTTATAGAAAACCTGTCCGTCGGTATCGCTATCTTCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATAAACAGACCGGAGAGCGAAACGGCCACAGTCAACGCCCCa                                                                                                 >  2:414718/1‑143 (MQ=255)                                          caccaGCGGGCCTTTATAGAAAACCTGTCCGTCGGTATCGCTATCTTCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACTCGGCGAAATAACCAGCCCGGAGAGCGATACGGCTTCAATCAACGCCCCAccac                                                                                             >  2:48886/1‑141 (MQ=255)                                          caccaGCGGGCCTTTATAGAAAACCTGTCCGTCGGTATCGCTATCTTCACGAACCTAGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATAAACAGACCGGAGACCGATACGACTTCAGTTAACGCCCCACCGc                                                                                             >  2:180716/1‑143 (MQ=255)                                                               gacgtgtctgtcGGTATCGCGTTCTTGGCGATCCGGGTCGTTGTTATCGCGGGCCTGAGGAATGGGACCCTGAGGAATAAACAGACCGGAGAGCGATACGGCTTCAGTTAACGCCCCACCGCCATTGCTACGCAGGTCgatga                                                                        <  2:459341/139‑1 (MQ=255)                                                                  cTGTCCGTCGGTATCGCTATCTTCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCCGGAATAAACAGACAGGAGAGCGATACGGCTTCAGTTAACGCCCCACCGCCATTGCTACGCAGGTCGAtgatga                                                                     >  1:316528/1‑143 (MQ=255)                                                                                                cTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATAAACAGACCGGAGAGCGATACGGCTTCAGTTAACGCCCCACCGCCATTGCTACGCAGGTCGATGATGACGCTGCTGACATTCTGTTTTTCCAGTTTCt                                       <  2:259438/143‑1 (MQ=255)                                                                                                       ttgttATCGCGGGCCTGGACAGTGGGACCCGCAGGAATAAACAGACCGGAGAGCGATACGGCTTCAGTTAACGCCCCACCGCCATTGCTACGCAGGTCGATGATGACGCTGCTGACATTCTGTTTTTCCAGTTTCTGCAGTTg                                <  2:129915/143‑1 (MQ=255)                                                                                                                cgGACCTGAACAATGGGACCCACAGGAATAAACAGACCGGAGAGCGATACGGCTTCAGTTAACGCCCCACCGCCATTGCTACGCAGGTCGATGATGACGCTGCTGACATTCTGTTTTTCCAGTTTCTGCAGTTGCACTTTGAc                       >  2:8925/1‑143 (MQ=255)                                                                                                                cgGACCTGAACAATGGGACCCACAGGAATAAACAGACCGGAGAGCGATACGGCTTCAGTTAACGCCCCACCGCCATTGCTACGCAGGTCGATGATGACGCTGCTGACATTCTGTTTTTCCAGTTTCTGCAGTTGCACTTTGAc                       >  2:123840/1‑143 (MQ=255)                                                                                                                   aCCTGAGCGATTGGTCCCGCAGGTATATACAGTCTGGAGAGCGTTTCGGCTTCAGTTATCGCCCCACCGCCATTGGTACGCAGGTCGATGATGACGCGGCTGACATTCTGTTGTTCCAGTTTCTGCAGTTGCACTTTGACATc                    <  1:478924/143‑1 (MQ=255)                                                                                                                       gAACAATGGTACCCGCAGGAATAAACAGACCGTAGAGCGATACGGCTTCAGTTAACGCCTTACCGTTATTGCTACGCAGGTCGATGATGACGCGGCTGACACTCGGTTTTTCCAGTTTCTGCAGTTGCACTTTGACATCGTCt                <  1:378430/143‑1 (MQ=255)                                                                                                                         aCAATGGGACCCACAGGAATAAACAGACCGGAGAGCGATACGGCTTCAGTTAACGCCCCACCGCCATTGCTACGCAGGTCGATGATGACGCTGCTGACATTCTGTTTTTCCAGTTTCTGCAGTTGCACTTTGACATCGTCTGt              <  1:8925/143‑1 (MQ=255)                                                                                                                         aCAATGGGACCCACAGGAATAAACAGACCGGAGAGCGATACGGCTTCAGTTAACGCCCCACCGCCATTGCTACGCAGGTCGATGATGACGCTGCTGACATTCTGTTTTTCCAGTTTCTGCAGTTGCACTTTGACATCGTCTGt              <  1:123840/143‑1 (MQ=255)                                                                                                                             tGGGACCCGCAGCCATAACCCGCCCGGAGAGCGCTACGGCTTCAGTTAACGCCCCACCGCCATTGCTACGCAGGTCGATGCTGACGCTGCTGACCTTCTGTTTTTCCAGTTTCTGCAGTTGCACTTTGACATCGTCTGTCaaa          <  2:34257/143‑1 (MQ=255)                                                                                                                                   ccGCAGGAATAAACAGACCGGAGAGCGATACGGCTTCAGTTAACGCCCCACCGCCATTGCTACGCAGGTCGATGATGACGCTGCTGACATTCTTTTTTTCCAGTTTCTGCAGTTACACTTTGACATCGTCTGTCAAACCcaca    >  1:457252/1‑143 (MQ=255)                                                                                                                                     gCAGGAATAAACAGACCGGAGAGCGATACGGCTTCAGTTAACGCCCCACCGCCATTGCTACGCAGGTCGATGATGACGCTGCTGACATTCTGTTTTTCCAGTTTATGAAGATGAAATTTAACATAGTCTGTCAAACCCACATa  >  1:165629/1‑143 (MQ=255)                                                                                                                                     |                                                                                                                                               AAGATTTCTGAAGCCGAAGCACTGAAGCGGTCAACCAGCACCACCAGCGGGCCTTTATAGAAAACCTGTCCGTCGGTATCGCTATCTTCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATAAACAGACCGGAGAGCGATACGGCTTCAGTTAACGCCCCACCGCCATTGCTACGCAGGTCGATGATGACGCTGCTGACATTCTGTTTTTCCAGTTTCTGCAGTTGCACTTTGACATCGTCTGTCAAACCCACATA  >  NC_000913/1913447‑1913721

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