Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB122	1,313,943	A→C	40.7%	L614R (CTG→CGG)	CJ019_01235 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB122	1,313,943	0	A	C	40.7%	19.3 / 3.9	35	L614R (CTG→CGG)	CJ019_01235	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (17/3);  new base C (2/12);  total (19/16)
Fisher's exact test for biased strand distribution p-value = 6.57e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGATGAAGGTGTCCGGCGCCAGCAACCCTTCCTGTGGATGCCGCCAGCGCACCAAGGCCTCGACCGCGACGATTCGCAGGTCGTCCAGGCGGTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCAGGGCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCGGCCTTGGCCTGGTACAGGGCGATGTCGGCGCAGCG  >  ACB122/1313812‑1314070                                                                                                                                    |                                                                                                                                ggATGAAGGTGTCCGGCGCCAGCAACCCTTCCTGTGGATGCCGCCAGCGCACCAAGGCCTCGACCGCGACGATTCGCAGGTCGTCCAGGCGGTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCAggg                                                                                                                              >  6:352246/1‑135 (MQ=255) ggATGAAGGTGTCCGGCGCCAGCAACCCTTCCTGTGGATGCCGCCAGCGCACCAAGGCCTCGACCGCGACGATTCGCAGGTCGTCCAGGCGGTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCAggg                                                                                                                              >  5:217290/1‑135 (MQ=255)   aTGAAGGTGTCCGGCGCCAGCAACCCTTCCTGTGGATGCCGCCAGCGCACCAAGGCCTCGACCGCGACGATTCGCAGGTCGTCCAGGCGGTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCAGGGcc                                                                                                                            >  6:293712/1‑135 (MQ=255)               ggCGCCAGCAACCCTTCCTGTGGATGCCGCCAGCGCACCAAGGCCTCGACCGCGACGATTCGCAGGTCGTCCAGGCGGTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCAGGGCCCTGCGCAGGTCg                                                                                                                >  1:79648/1‑135 (MQ=255)                      gCAACCCTTCCTGTGGATGCCGCCAGCGCACCAAGGCCTCGACCGCGACGATTCGCAGGTCGTCCAGGCGGTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCAGGGCCCTGCGCAGGTCGTTCTCCAg                                                                                                        >  8:42420/1‑136 (MQ=255)                      gCAACCCTTCCAGTGGATACAGCCAGCGCACCAAGGCCTCGACCGCGACGATTCGCAGGTCGTCCAGGCGTTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGAGCAGGGCCCTGCGCAGGTCAGTCTCCa                                                                                                         >  2:55190/1‑135 (MQ=255)                         aCCCTTCCTGTGGATGCCGCCAGCGCACCAAGGCCTCGACCGCGACGATTCGCAGGTCGTCCAGGCGGTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCCGGGCCCGGCGCAGGTCGTTCTCCAGcgg                                                                                                     >  7:206437/1‑134 (MQ=255)                                 tgtgGATGCCGCCAGCGCACCAAGGCCTCGACCGCGACGATTCGCAGGTCGTCCAGGCGGTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCAGGGCCCTGCGCAGGTCGTTCTCCAGCTGgcggcg                                                                                               >  6:396804/1‑134 (MQ=255)                                 tgtgGATGCCGCCAGCGCACCAAGGCCTCGACCGCGACGATTCGCAGGTCGTCCAGGCGGTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCAGGGCCCTGCGCAGGTCGTTCTCCAGCTGgcggcg                                                                                               >  2:91336/1‑134 (MQ=255)                                        gccgccAGCGCACCAAGGCCTCGACCGCGACGATTCGCAGGTCGTCCAGGCGGTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCCCGGCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGa                                                                                       >  3:72357/1‑135 (MQ=255)                                         ccgccAGCGCACCAAGGCCTCGACCGCGACGATTCGCAGGTCGTCCAGGCGGTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCAGGGCCCGGCGCAGGTCGTTCTCCCGCTTGAGGCGGGACTGGag                                                                                      >  5:16930/1‑134 (MQ=255)                                             cAGCGCACCAAGGCCTCGACCGCGACGATTCGCAGGTCGTCCAGGCGGTAGCGTGGCTGGTAATGCAACTCGAACTCCACATTGCGCAGGGCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTCCTGGATcggcg                                                                                 >  8:346293/1‑135 (MQ=255)                                                               aCCGCGACGATTCGCAGGTCGTCCAGGCGGTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCAGGGCCCTGCGCAGGTCGTTCTCCCGCGGGCGGCGGGAGGGGATCTGCTGGGTCTTTTGGGCGAc                                                                 >  5:347977/1‑134 (MQ=255)                                                                ccGCGACGATTCGCAGGTCGTCCAGGCGGTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCAGGGCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGGACTGGATGTGCTGGTTCATTTCGGCGACgca                                                              >  7:267661/1‑134 (MQ=255)                                                                     acgaTTCGCAGGTCGTCCGGGCGGTAGCGTGGCTGGTAATGCAACGCGACCGCCTCATTGCGCCGGGCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAgc                                                         <  8:374885/136‑1 (MQ=255)                                                                        aggcGCAGGGCGGCCCGGGGGTCGCGGGGCGGGTAATCCAACTCCAACTCCTCTTTGCGCCGGGCCCGGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAgcgc                                                       <  4:195482/132‑1 (MQ=255)                                                                           cGCAGGTCGTCCAGGCGGTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCAGGGCCCTGCGCAGGGCGGTCTCCAGCTGGCGGCGGGACTGGGTCGGCGGGTTCCTTTCGGCAACGAAGTAGCGCCaa                                                    >  6:241573/1‑134 (MQ=255)                                                                                     ccAGGCGGTCGCGTGGCTGTTAATGCAACGCGAACTCCTCATGGCGCCGGGCCCGGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTc                                          <  2:305565/135‑1 (MQ=255)                                                                                          gggTAGGGTGGGTGGTAATTCAACTCTACCCCCTCTTGGCGCCGGCCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCgg                                    <  8:139973/135‑1 (MQ=255)                                                                                          cGGTTGCGTGGCGGGTAATGCAACTCGAACTCCTCATTGCGCCGGGCCCTGCGCAGGTCGTGCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCgg                                    <  8:80164/136‑1 (MQ=255)                                                                                          cGGTAGGGTGGCTGGAAATGCAACTCGAACTCCTCATTGCGCCGGGCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCgg                                    <  8:89130/136‑1 (MQ=255)                                                                                           ggttgggTGGCTTCTACTCCAACTCCAACTCCTACTGGCGCCGGGCCGGGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCgg                                    <  5:396804/129‑1 (MQ=255)                                                                                           ggTAGGGTGGCTGGTAATGCAACTCGAACTCCTCATGGCCCCGGGCCGGGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCgg                                    <  1:79719/135‑1 (MQ=255)                                                                                           ggTAGCGTGGCTGGTAATGCACCTCGAACTCCTCATTGCTCAGGGCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCCGTACTGGCTCTTCTGGTTCATTTAGGCAACGAAGTAGCGCCATGTCTTCTTTCCATCgg                                    >  2:40899/1‑135 (MQ=255)                                                                                           ggTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCAGGGCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTTCTGGTTCATTTCGGCAACGAACTAGCGCCATGTGTTCTTTCCATCgg                                    >  6:65782/1‑135 (MQ=255)                                                                                               gggTGGATGGAATTGCAACTCTATCTCCTCATTGCGCCGGCCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTTCTGGATCTGCTTGTTCATTTCGGCTACGAAGTAGCGCCATGTGTTCTTTCCATCGGCCTTg                               <  8:42870/134‑1 (MQ=255)                                                                                               gCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCAGGGCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCAATCGGCCTTg                               >  7:159965/1‑136 (MQ=255)                                                                                               gCGGGGGTGGTAATGCAACTCGAACTCCTCATTGCGCGGGGCCCGGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCGGCCtt                                <  8:280575/135‑1 (MQ=255)                                                                                                ggTGGCTGTTAATGAATCTCTAACTCATCATTGCGCCGGGCCGTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCGGCCTTg                               <  8:243790/134‑1 (MQ=255)                                                                                                cggggCTTGTAATCCAACGCGAACACAGATATGCGCAGCGCCCGGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTGCTTTCCATCGGCCTTg                               <  2:49834/132‑1 (MQ=255)                                                                                                cggggCTGGAAATCCACCTCGACCCCCTCAGGGCGCCGGGCCCGGCCCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCGGCCTTg                               <  2:397805/132‑1 (MQ=255)                                                                                                cggggCTGGAAATCCAACCCCAACTCCTCATTGCGCCGGGCCCGGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCGGCCTTg                               <  4:161878/132‑1 (MQ=255)                                                                                                 gTGGCTGTTAATGAAACTCGAACTCCTCATTGCGCAGGGCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCTGTCCTGGATCTTATGGTTAATTTCGGCAACGAAGTATCGACATGTGTTCTTTCCATCGGCCTTgg                              >  4:345860/1‑135 (MQ=255)                                                                                                                          ctcATTGCGCAGGGCCCGGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCGGCCTTGGCCTGGTACAGGGCGATGTCGGCGCAg    <  8:170762/136‑1 (MQ=255)                                                                                                                              ttGCGCAGGGCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCGGCCTTGGCCTGGTACAGGGCGATGTCGGCGCAGCg  <  4:41507/134‑1 (MQ=255)                                                                                                                                    |                                                                                                                                GGATGAAGGTGTCCGGCGCCAGCAACCCTTCCTGTGGATGCCGCCAGCGCACCAAGGCCTCGACCGCGACGATTCGCAGGTCGTCCAGGCGGTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCAGGGCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCGGCCTTGGCCTGGTACAGGGCGATGTCGGCGCAGCG  >  ACB122/1313812‑1314070

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