Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	1,232,353	T→G	9.1%	T135P (ACG→CCG)	ptrB ←	protease II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	1,232,353	0	T	G	9.1%	199.8 / 13.0	76	T135P (ACG→CCG)	ptrB	protease II
Reads supporting (aligned to +/- strand):  ref base T (38/31);  new base G (3/4);  total (41/35)
Fisher's exact test for biased strand distribution p-value = 6.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAAACGAATGCCGTACTGGCGTCGGGAAAGAAAATCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTGGCATCGAGCAATGTTTCC  >  Exported/1232284‑1232423                                                                      |                                                                       gAAAACGAATGCCGTACTGGCGTCGGGAAAGAAAATCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTa                                                                       <  1:1125421/71‑1 (MQ=255)   aaaCGAATGCCGTACTGGCGTCGGGAAAGAAAATCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTaa                                                                      >  1:1138190/1‑70 (MQ=255)   aaaCGAATGCCGTACTGGCGTCGGGAAAGAAAATCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAAt                                                                     >  1:2916527/1‑71 (MQ=255)   aaaCGAATGCCGTACTGGCGTCGGGAAAGAAAATCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAAt                                                                     >  1:1282257/1‑71 (MQ=255)    aaCGAATGCCGTACTGGCGTCGGGAAAGAAAATCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATc                                                                    >  1:1901630/1‑71 (MQ=255)       gAATGCCGTACTGGCGTCGGGAAAGAAAATCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGc                                                                  >  1:822838/1‑70 (MQ=255)        aaTGCCGTACTGGCGTCGGGAAAGAAAATCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATc                                                                    >  1:2973247/1‑67 (MQ=255)         aTGCCGTACTGGCGTCGGGAAAGAAAATCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCAt                                                               <  1:21634/71‑1 (MQ=255)         aTGCCGTACTGGCGTCGGGAAAGAAAATCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCAt                                                               <  1:3233465/71‑1 (MQ=255)           gCCGTACTGGCGTCGGGAAAGAAAATCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTaa                                                                      <  1:2532211/62‑1 (MQ=255)            ccGTACTGGCGTCGGGAAAGAAAATCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTcc                                                            <  1:1742057/71‑1 (MQ=255)            ccGTACTGGCGTCGGGAAAGAAAATCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTcc                                                            <  1:3652393/71‑1 (MQ=255)              gTACTGGCGTCGGGAAAGAAAATCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGc                                                          >  1:242790/1‑71 (MQ=255)              gTACTGGCGTCGGGAAAGAAAATCTTCTGCCAGCGCCATAAAGGTGTTATCGGGCGTAAt                                                                     >  1:2519753/1‑60 (MQ=255)              gTACTGGCGTCGGGAAAGAAAATCTTCTGCCAGCGCCATAAAGGTGTTATCGGGCGTAATCGCCATTCCGc                                                          >  1:2515542/1‑71 (MQ=255)                      gTCGGGAAAGAAAATCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCAtt                                                              >  1:1028786/1‑59 (MQ=255)                       tCGGGAAAGAAAATCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCAtt                                                              >  1:3405419/1‑58 (MQ=255)                         gggAAAGAAAATCTTCTGGCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATaa                                               <  1:2173542/71‑1 (MQ=255)                         gggAAAGAAAATCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATaa                                               <  1:1189974/71‑1 (MQ=255)                           gaaagaaaATCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATaaaa                                             <  1:1435691/71‑1 (MQ=255)                              agaaaATCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGc                                                          >  1:1766811/1‑55 (MQ=255)                               gaaaATCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGc                                                                  <  1:2475121/46‑1 (MQ=255)                               gaaaATCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCa                                         <  1:1310724/71‑1 (MQ=255)                               gaaaATCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCa                                         >  1:1637556/1‑71 (MQ=255)                               gaaaATCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCa                                         >  1:546977/1‑71 (MQ=255)                               gaaaATCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCa                                         <  1:2323041/71‑1 (MQ=255)                                aaaaTCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTaa                                                                      <  1:1851108/41‑1 (MQ=255)                                aaaaTCTTCTGCCAGCGACATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCAc                                        >  1:2319817/1‑71 (MQ=255)                                 aaaTCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCAc                                        >  1:2356153/1‑70 (MQ=255)                                 aaaTCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCAc                                        <  1:1746360/70‑1 (MQ=255)                                 aaaTCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACt                                       <  1:14212/71‑1 (MQ=255)                                 aaaTCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACt                                       <  1:1673387/71‑1 (MQ=255)                                   aTCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCAtt                                                              >  1:2485635/1‑46 (MQ=255)                                    tcttctGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCa                                                                >  1:2741840/1‑43 (MQ=255)                                    tcttctGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCa                                                                >  1:765386/1‑43 (MQ=255)                                    tcttctGCCAGCGCCATAATGGTGTTATCGGGCGGAATGGCGATTACGccc                                                        >  1:3135235/1‑51 (MQ=21)                                    tcttctGCCAGCGCCATAATGGTGTTATCGGGCGGAATGGCGATTACGccc                                                        >  1:2025436/1‑51 (MQ=21)                                    tcttctGCCAGCGCCATAATGGTGTTATCGGGCGGAATGGCGATTACGccc                                                        >  1:1842334/1‑51 (MQ=21)                                     cttctGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTAt                                     >  1:107123/1‑69 (MQ=255)                                       tctGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGc                                 >  1:1840121/1‑71 (MQ=255)                                        ctGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATg                                    <  1:43657/67‑1 (MQ=255)                                        ctGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATg                                    <  1:2257400/67‑1 (MQ=255)                                        ctGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATg                                    <  1:1772999/67‑1 (MQ=255)                                            cAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTgcgc                            <  1:1553584/71‑1 (MQ=255)                                              gcgcCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCa                                         <  1:3221466/56‑1 (MQ=255)                                               cgcCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCt                         >  1:329111/1‑71 (MQ=255)                                               cgcCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCt                         >  1:355919/1‑71 (MQ=255)                                                gcCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCtt                        <  1:3243234/71‑1 (MQ=255)                                                gcCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCtt                        <  1:2452894/71‑1 (MQ=255)                                                 ccATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTgcgcgc                          <  1:2305181/68‑1 (MQ=255)                                                 ccATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTa                       >  1:2349938/1‑71 (MQ=255)                                                  cATAATGGTGTTATCGGGCGTAATCGCCATTCCGCACATCGAATAAAACTCACTATGAGCTGCGCGCTTAt                      <  1:3102388/71‑1 (MQ=255)                                                   aTAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTAtt                     >  1:664892/1‑71 (MQ=255)                                                    tAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTa                                      <  1:957888/53‑1 (MQ=255)                                                    tAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTg                    >  1:1857135/1‑71 (MQ=255)                                                    tAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTg                    >  1:1004755/1‑71 (MQ=255)                                                    tAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTg                    >  1:3287211/1‑71 (MQ=255)                                                    tAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTg                    >  1:1891808/1‑71 (MQ=255)                                                    tAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTg                    >  1:3282933/1‑71 (MQ=255)                                                     aaTGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGGGCGCTTATTgg                   >  1:2246104/1‑71 (MQ=255)                                                     aaTGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTgg                   >  1:3656216/1‑71 (MQ=255)                                                     aaTGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTgg                   >  1:3619948/1‑71 (MQ=255)                                                     aaTGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTgg                   >  1:472771/1‑71 (MQ=255)                                                     aaTGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTgg                   >  1:3205421/1‑71 (MQ=255)                                                     aaTGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTgg                   >  1:2117673/1‑71 (MQ=255)                                                              aTCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTGGCATCGa             >  1:1101784/1‑68 (MQ=255)                                                               tCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTGGCATCGAGCaa         <  1:2859813/71‑1 (MQ=255)                                                               tCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTGGCATCGAGCaa         <  1:1338074/71‑1 (MQ=255)                                                                cGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTGGCATCGAGCaa         <  1:470352/70‑1 (MQ=255)                                                                 gggCGGAATGGCGATTACGCCCATCGAATAAAACTCACTATGAGc                                 <  1:113402/45‑1 (MQ=18)                                                                 gggCGGAATGGCGATTACGCCCATCGAATAAAACTCACTATGAGc                                 <  1:3324596/45‑1 (MQ=18)                                                                 gggCGGAATGGCGATTACGCCCATCGAATAAAACTCACTATGAGc                                 <  1:1601209/45‑1 (MQ=18)                                                                 gggCGGAATGGCGATTACGCCCATCGAATAAAACTCACTATGAGc                                 <  1:1042718/45‑1 (MQ=18)                                                                    cGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGa                                   <  1:1999930/40‑1 (MQ=255)                                                                     gTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTGGCATCGAGCAATGTTTc   >  1:1735470/1‑71 (MQ=255)                                                                     gTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTGGCATCGAGCAATGTTTc   <  1:2189975/71‑1 (MQ=255)                                                                      tAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGc                                 >  1:1805924/1‑40 (MQ=255)                                                                      tAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTGGCATCGAGCAATGTTTcc  <  1:2043933/71‑1 (MQ=255)                                                                      tAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTGGCATCGAGCAATGTTTcc  >  1:238452/1‑71 (MQ=255)                                                                      |                                                                       GAAAACGAATGCCGTACTGGCGTCGGGAAAGAAAATCTTCTGCCAGCGCCATAATGGTGTTATCGGGCGTAATCGCCATTCCGCCCATCGAATAAAACTCACTATGAGCTGCGCGCTTATTGGCATCGAGCAATGTTTCC  >  Exported/1232284‑1232423

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑