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	96.0%	87.2 / NA	25	G535G (GGC→GGG)	ptrB	protease II
Reads supporting (aligned to +/- strand):  ref base . (0/0);  major base C (12/12);  minor base A (1/0);  total (13/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACCAGTGGTAAGAGGAATTGATTCATCAAGCATCGTTGTTACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATG‑CCCCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCATCGCAGGCATCAAGATAATCATTAAACGTATTTTTCTTCTT  >  NC_000913/1927114‑1927342                                                                                                                           |                                                                                                            cACCAGTGGTAAGAGGAATTGATTCATCAAGCATCGTTGTTACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGc                                                                                             >  2:217308/1‑139 (MQ=255)        ggTAAGAGGAATTGATTCATCAAGCATCGTTGTTACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGc                                                                                      >  1:369045/1‑139 (MQ=255)          tAAGAGGAATTGATTCATCAAGCATCGTTGTTACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCAt                                                                                    <  2:64389/139‑1 (MQ=255)                 aaTTGATTCATCAAGCATCGTTGTTACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAg                                                                             >  2:136645/1‑139 (MQ=255)                 aaTTGATTCATCAAGCATCGTTGTTACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAg                                                                             >  1:72274/1‑139 (MQ=255)                   ttGATTCATCAAGCATCGTTGTTACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAGCg                                                                           >  1:598021/1‑139 (MQ=255)                              aGCATCGTTGTTACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCa                                                                <  1:217308/139‑1 (MQ=255)                                   cGTTGTTACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGcc                                                           <  1:505718/139‑1 (MQ=255)                                    gttgttACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAATATG‑cccccccg                                                                                                     >  1:47713/1‑93 (MQ=255)                                    gttgttACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAATATG‑cccccccg                                                                                                     <  2:47713/95‑3 (MQ=255)                                       gttACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGt                                                       <  1:99664/139‑1 (MQ=255)                                              caTCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACa                                                <  1:273515/139‑1 (MQ=255)                                               aTCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAg                                                                             <  2:552139/109‑1 (MQ=255)                                               aTCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAg                                                                             >  1:552139/1‑109 (MQ=255)                                                  aacaaACGGTACCTGGGCGATAACGCCGTGGAATAATTACGGGCGTTGATTAATTGCAACGCCCATCAACATGACCCCCCGCCCTCCCGCCCATCGCATAACAAAGCGAAGGAGCGCCATAGCCTAGTTTTAACAATGc                                            >  1:24247/1‑139 (MQ=255)                                                        cGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAATTTTAACAATGCATCGCa                                      <  1:340034/139‑1 (MQ=255)                                                           tACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACaa                                               <  2:347402/127‑1 (MQ=255)                                                           tACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACaa                                               >  1:347402/1‑127 (MQ=255)                                                              cTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATG‑cccccccg                                                                                                     >  1:276110/1‑67 (MQ=255)                                                              cTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATG‑cccccccg                                                                                                     <  2:276110/69‑3 (MQ=255)                                                                gggCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGc                                                            >  2:581183/1‑109 (MQ=255)                                                                gggCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGc                                                            <  1:581183/109‑1 (MQ=255)                                                                 ggCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGt                                                       >  2:434864/1‑113 (MQ=255)                                                                 ggCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGt                                                       <  1:434864/113‑1 (MQ=255)                                                                 ggCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATG‑cc                                                                                                           >  1:591310/1‑60 (MQ=255)                                                                 ggCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATG‑cc                                                                                                           <  2:591310/60‑1 (MQ=255)                                                                    gATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCATCGCAGGCATCAAGATa                          <  2:72274/139‑1 (MQ=255)                                                                     aTAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCATCGCAGGCATCAAGATaa                         >  1:66593/1‑139 (MQ=255)                                                                       aaCGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCATCGCAGGCATCAAGATAATc                       >  2:596333/1‑139 (MQ=255)                                                                           ccGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCATCGCAGGCATCAAGATAATCATTa                   >  1:261943/1‑139 (MQ=255)                                                                                            gCGTTGATTAATTGCAACGCCCATCAACATGCCCCCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCATCGCAGGCATCAAGATAATCATTAAACGTATTTttcttctt  <  2:369045/139‑1 (MQ=255)                                                                                                                           |                                                                                                            CACCAGTGGTAAGAGGAATTGATTCATCAAGCATCGTTGTTACAACATCAACAAACGGTACCTGGGCGATAACGCCGTGGAATAATTCCGGGCGTTGATTAATTGCAACGCCCATCAACATG‑CCCCCCGCACTCCCGCCCATCGCATAACAAAGCGAAGGAGAGCCATAGCCCAGTTTTAACAATGCATCGCAGGCATCAAGATAATCATTAAACGTATTTTTCTTCTT  >  NC_000913/1927114‑1927342

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