Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,927,241	(C)6→7	coding (1599/2061 nt)	ptrB ←	protease II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,927,235	1	.	C	98.1%	188.9 / ‑3.5	52	coding (1605/2061 nt)	ptrB	protease II
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base C (26/25);  total (27/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

ACCCCACTCTTCAAACTCACCAGTGGTAAGAGGAATTGATTCATCAAGCATCGTTGTTACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATG‑CC‑CCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCATCGCAGGCATCAAGATAATCATTAAACGTATTTTTCTTCTTCAGAAATTTTCCGTCTTCGTA  >  NC_000913/1927097‑1927363                                                                                                                                            |                                                                                                                                  aCCCCACTCTTCAAACTCACCAGTGGTAAGAGGAATTGATTCATCAAGCATCGTTGTTACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATg                                                                                                                                    >  1:768143/1‑139 (MQ=255)         tttCAAACTCACCAGTGGTAAGAGGAATTGATTCATCAAGCATCGTTGTTACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATG‑cc‑cccccg                                                                                                                          <  1:227322/138‑3 (MQ=255)         ctTCAAACTCACCAGTGGTAAGAGGAATTGATTCATCAAGCATCGTTGTTACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATG‑cc‑cccccg                                                                                                                          >  1:181916/1‑137 (MQ=255)           tCAAACTCACCAGTGGTAAGAGGAATTGATTCATCAAGCATCGTTGTTACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATG‑cc‑cccccgca                                                                                                                        <  2:605515/139‑5 (MQ=255)               aCTCACCAGTGGTAAGAGGAATTGATTCATCAAGCATCGTTGTTACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCC‑CCCCGCACTcc                                                                                                                     >  1:359404/1‑139 (MQ=255)                 tCACCAGTGGTAAGAGGAATTGATTCATCAAGCATCGTTGTTACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCC‑CCCCGCACTCCCg                                                                                                                   >  1:237716/1‑139 (MQ=255)                  cACCAGTGGTAAGAGGAATTGATTCATCAAGCATCGTTGTTACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCC‑CCCCGCACTCCCGc                                                                                                                  >  2:743443/1‑139 (MQ=255)                  cACCAGTGGTAAGAGGAATTGATTCATCAAGCATCGTTGTTACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCC‑CCCCGCACTCCCGc                                                                                                                  >  1:16089/1‑139 (MQ=255)                   aCCAGTGGTAAGAGGAATTGATTCATCAAGCATCGTTGTTACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCC‑CCCCGCACTCCCGcc                                                                                                                 >  2:739938/1‑139 (MQ=255)                   aCCAGTGGTAAGAGGAATTGATTCATCAAGCATCGTTGTTACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCC‑CCCCGCACTCCCGcc                                                                                                                 <  2:575118/139‑1 (MQ=255)                       gTGGTAAGAGGAATTGATTCATCAAGCATCGTTGTTACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCC‑CCCCGCACTCCCGCCCATc                                                                                                             <  2:16089/139‑1 (MQ=255)                          gTAAGAGGAATTGATTCATCAAGCATCGTTGTTACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCC‑CCCCGCACTCCCGCCCATCGCa                                                                                                          <  2:73786/139‑1 (MQ=255)                             agagGAATTGATTCATCAAGCATCGTTGTTACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCC‑CCCCGCACTCCCGCCCATCGCATaa                                                                                                       >  2:1041134/1‑139 (MQ=255)                                  aaTTGATTCATCAAGCATCGTTGTTCCAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCC‑CCCCGCACTCCCGCCCATc                                                                                                             <  2:704273/128‑1 (MQ=255)                                  aaTTGATTCATCAAGCATCGTTGTTACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCAGCAACATGCCC‑CCCCGCACTCCCGCCCATc                                                                                                             >  1:704273/1‑128 (MQ=255)                                   aTTGATTCATCAAGCATCGTTGTTACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCC‑CCCCGCACTCCCGCCCATCGCATAACAAAGc                                                                                                 >  1:594462/1‑139 (MQ=255)                                   aTTGATTCATCAAGCATCGTTGTTACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATG‑cc‑ccccc                                                                                                                           <  1:452916/112‑2 (MQ=255)                                   aTTGATTCATCAAGCATCGTTGTTACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATG‑cc‑ccccc                                                                                                                           >  2:452916/1‑111 (MQ=255)                                    ttGATTCATCAAGCATCGTTGTTACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCC‑CCCCGCACTCCCGCCCATCGCATAACAAAGCg                                                                                                <  1:647013/139‑1 (MQ=255)                                            tcaAGCATCGTTGTTACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCC‑CCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGagag                                                                                        >  1:82836/1‑139 (MQ=255)                                            tcaAGCATCGTTGTTACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCC‑CCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGagag                                                                                        >  1:1095798/1‑139 (MQ=255)                                              aaGCATCGTTGTTACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCC‑CCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGcc                                                                                      <  1:1041134/139‑1 (MQ=255)                                                 cATCGTTGTTACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCC‑CCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATa                                                                                   >  1:46632/1‑139 (MQ=255)                                                  aTCGTTGTTACAACATCAAC‑AACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCC‑CCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGc                                                                                 <  1:119569/139‑1 (MQ=255)                                                   tCGTTGTTACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCC‑CCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGc                                                                                 <  2:181916/139‑1 (MQ=255)                                                   tCGTTGTTACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCC‑CCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGc                                                                                 <  2:46632/139‑1 (MQ=255)                                                     gttgttACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCC‑CCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGccc                                                                               >  1:576134/1‑139 (MQ=255)                                                      ttgttACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCC‑CCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCa                                                                              <  1:739938/139‑1 (MQ=255)                                                      ttgttACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCC‑CCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCa                                                                              >  2:357382/1‑139 (MQ=255)                                                          tACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCC‑CCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGttt                                                                          <  2:768143/139‑1 (MQ=255)                                                           acaacaTCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATG‑CCTCCCCGCACTCCCGCCCATCGCCTAACAAAGCGAAGGAGAGCCATAGCCCAGtttt                                                                         >  1:336737/1‑139 (MQ=255)                                                                   aacaaACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCC‑CCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGc                                                                 >  1:549075/1‑139 (MQ=255)                                                                           gTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCC‑CCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCATCGCAgg                                                         <  2:1095798/139‑1 (MQ=255)                                                                            tACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCC‑CCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCATCg                                                             >  1:1072648/1‑134 (MQ=255)                                                                            tACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCC‑CCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCATCg                                                             <  2:1072648/134‑1 (MQ=255)                                                                            tACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCC‑CCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCATCGCAGGc                                                        <  2:594462/139‑1 (MQ=255)                                                                                  ggCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATg                                                                                                                                    >  1:587110/1‑58 (MQ=255)                                                                                  ggCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATg                                                                                                                                    <  2:587110/58‑1 (MQ=255)                                                                                            ccGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCC‑CCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGc                                                                                 <  2:813533/98‑1 (MQ=255)                                                                                            ccGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCC‑CCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGc                                                                                 >  1:813533/1‑98 (MQ=255)                                                                                            ccGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCC‑CCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCATCGCAGGCATCAAGATAATCATTa                                        >  1:420688/1‑139 (MQ=255)                                                                                                ggAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCC‑CCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCATCGCAGGCATCAAGATAATCATTAAACg                                    <  2:359404/139‑1 (MQ=255)                                                                                                            ggCGTTGATTAATTGCAACGCCCATCAACATGCCC‑CCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCATCGCAGGCATCAAGATAATCATTAAACGTATTTttcttct                        >  2:894610/1‑139 (MQ=255)                                                                                                            ggCGTTGATTAATTGCAACGCCCATCAACATGCCC‑CCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCATCGCAGGCATCAAGATAATCATTAAACGTATTTttcttct                        >  1:948019/1‑139 (MQ=255)                                                                                                                ttGATTAATTGCAACGCCCATCAACATGCCC‑CCCCGCACTCCCGCCCATCGCATAACAAACCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCATCGCAGGCATCAAGATAATCATTAAACGTATTTTTCTTCTTCAg                    >  1:755393/1‑139 (MQ=255)                                                                                                                 tGATTAATTGCAACGCCCATCAACATGCCC‑CCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCATCGCAGGCATCAAGATAATCATTAAACGTATTTttcttctt                       <  2:1120512/135‑1 (MQ=255)                                                                                                                 tGATTAATTGCAACGCCCATCAACATGCCC‑CCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCATCGCAGGCATCAAGATAATCATTAAACGTATTTttcttctt                       >  1:1120512/1‑135 (MQ=255)                                                                                                                         tGCAACGCCCATCAACATGCCC‑CCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCATCGCAGGCATCAAGATAATCATTAAACGTATTTTTCTTCTTCAGAAATTTTcc           <  2:237716/139‑1 (MQ=255)                                                                                                                           cAACGCCCATCAACATGCCC‑CCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCATCGCAGGCATCAAGATAATCATTAAACGTATTTTTCTTCTTCAGAAATTTTCCGt         <  2:615338/139‑1 (MQ=255)                                                                                                                           cAACGCCCATCAACATGCCC‑CCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCATCGCAGGCATCAAGATAATCATTAAACGTATTTTTCTTCTTCAGAAATTTTCCGt         <  1:357382/139‑1 (MQ=255)                                                                                                                             aCGCCCATCAACATGCCC‑CCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGt                                                                            <  2:608436/70‑1 (MQ=255)                                                                                                                             aCGCCCATCAACATGCCC‑CCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGt                                                                            >  1:608436/1‑70 (MQ=255)                                                                                                                             aCGCCCATCAACATGCCC‑CCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCATCGCAGGCATCAAGATAATCATTa                                        <  2:466863/106‑1 (MQ=255)                                                                                                                             aCGCCCATCAACATGCCC‑CCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCATCGCAGGCATCAAGATAATCATTa                                        >  1:466863/1‑106 (MQ=255)                                                                                                                             aCGCCCATCAACATGCCC‑CCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCATCGCAGGCATCAAGATAATCATTAAACGTATTTTTCTTCTTCAGAAATTTTcc           <  1:787378/135‑1 (MQ=255)                                                                                                                             aCGCCCATCAACATGCCC‑CCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCATCGCAGGCATCAAGATAATCATTAAACGTATTTTTCTTCTTCAGAAATTTTcc           >  2:787378/1‑135 (MQ=255)                                                                                                                              cGCCCATCAACATGCCC‑CCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCATCGCa                                                           <  2:222403/86‑1 (MQ=255)                                                                                                                              cGCCCATCAACATGCCC‑CCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCATCGCa                                                           >  1:222403/1‑86 (MQ=255)                                                                                                                                   aTCAACATGCCC‑CCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCATCGCAGGCATCAAGATAATCATTAAACGTATTTTTCTTCTTCAGAAATTTTCCGTCTTCGTa  <  1:967176/138‑1 (MQ=255)                                                                                                                                   aTCAACATGCCC‑CCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCATCGCAGGCATCAAGATAATCATTAAACGTATTTTTCTTCTTCAGAAATTTTCCGTCTTCGTa  >  2:967176/1‑138 (MQ=255)                                                                                                                                            |                                                                                                                                  ACCCCACTCTTCAAACTCACCAGTGGTAAGAGGAATTGATTCATCAAGCATCGTTGTTACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATG‑CC‑CCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCATCGCAGGCATCAAGATAATCATTAAACGTATTTTTCTTCTTCAGAAATTTTCCGTCTTCGTA  >  NC_000913/1927097‑1927363

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