Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB122	1,313,943	A→C	45.4%	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	45.4%	6.5 / 18.3	34	L614R (CTG→CGG)	CJ019_01235	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (11/3);  new base C (4/12);  total (15/19)
Fisher's exact test for biased strand distribution p-value = 9.22e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.96e-01

ACGGGATGAAGGTGTCCGGCGCCAGCAACCCTTCCTGTGGATGCCGCCAGCGCACCAAGGCCTCGACCGCGACGATTCGCAGGTCGTCCAGGCGGTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCAGGGCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCGGCCTTGGCCTGGTACAGGGCGATGTCGGCGCAGCGGAT  >  ACB122/1313809‑1314073                                                                                                                                       |                                                                                                                                   aCGGGATGAAGGTGTCCGGCGCCAGCAACCCTTCCTGTGGATGCCGCCAGCGCACCAAGGCCTCGACCGCGACGATTCGCAGGTCGTCCAGGCGGTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCa                                                                                                                                    >  1:135704/1‑135 (MQ=255)      aTGAAGGTGTCCGGCGCCAGCAACCCTTCCTGTGAATGCCGCCAGCGCACCAAGGCCTCGACCGCGACGATTCGCAGGTCGTCCAGGCGGTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCAGGGcc                                                                                                                               >  6:43159/1‑135 (MQ=255)               tCCGGCGCCAGCAACCCTTCCTGTGGTTGCCGCCAGCGCACCAAGGCCTCGACCGCGACGATTCGATGGTCGTCCTGGCTGTAGCGGGGGTGGTAATGCAACTCCAAGTCCTCTTTGCGCAGGGCCCGGCGCAgg                                                                                                                      <  8:43216/135‑1 (MQ=255)                       cAGCAACCCTTCCTGTGGATGCCGCCAGCGCACCAAGGCCTCGACCGCGACGATTCGCAGGTCGTCCAGGCGGTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCAGGGCCCTGCGCAGGTCGAtctc                                                                                                              >  2:79272/1‑135 (MQ=255)                                    tgtgGATGCCGCCAGCGCACCAAGGCCTCGACCGCGACGATTCGCAGGTCGTCCAGGCGGTAGCGTGGCTGGTAATGCAACTCGATCTCCTCATTGCGCAGGGCCCTGCGCAGGTCGTTCTCCAGCTGgcggcg                                                                                                  >  2:344584/1‑134 (MQ=255)                                           gccgccAGCGCACCAAGGCCTCGACCGCGACGATTCGCAGGTCGTCCAGGCGGTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCAGGGCCCAGCGCAGGTGGTTCTCCAGCTGGCGGCGGGACTggg                                                                                          >  5:261381/1‑134 (MQ=255)                                           gccgccAGCGCACCAAGGCCTCGACCGCGACGATTCGCAGGTCGTCCAGGCGGTAGCGGGGCTGGTAATGCAACTCGAACTCCTCATTGCGCAGGGCCCTGCGCAGGTCGTTCTCCAGCGGGCGGCGGTACTGGa                                                                                          >  2:309505/1‑135 (MQ=255)                                              gccAGCGCACCAAGGCCTCGACCGCGACGATTCGCAGGTCGTCCAGGCGGTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCCGGGCCCGGCGCCGGTGGGTCTCCAGCTGGCGGCGGGACTGGGTcg                                                                                       >  5:48869/1‑134 (MQ=255)                                              gccAGCGCACCAAGGCCTCGACCGCGACGATTCGCAGGTCGTCCAGGCGGTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCAGGGCCCTGCGCAGGTCGTTCTGCAGCGGGCGGCGGTACTGGATct                                                                                       >  5:171921/1‑135 (MQ=255)                                              gccAGCGCACCAAGGCCTCGACCGCGACGATTCGCAGGTCGTCCAGGCGGTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCAGGGCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGGACTGGATct                                                                                       >  5:262272/1‑135 (MQ=255)                                                      aCCAAGGCCTCGACCGCGACGATTCGCAGGTCGTCCAGGCGGTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCAGGGCCCTGCGCAGGTCGTTCTCAAGCTGGCGGCGGTTCTGGTTCTGGTGGTTc                                                                               >  5:249240/1‑135 (MQ=255)                                                            gCCTCGACCGCGACGATTCGCAGGTCGTCCAGGCGGTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCCGGGCCCGGCGCCGGTCGTTCTCCAGCTGGCGGCGGTACGGGCTGTGCGGGGTCATTTCg                                                                         >  7:260881/1‑135 (MQ=255)                                                                   ccGCGACGATTCGCAGGTCGTCCAGGCGGTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCCGGGCCCTGCGCAGGTCGTTCTCCAGCGGGCGGCGGTACAGGATCTGCTGGTTAATTTCGGCAAc                                                                    >  8:63644/1‑133 (MQ=255)                                                                   ccGCGACGATTCGCAGGTCGTCCAGGCGGTAGCGGGGCTGGTAATGCAACTCGCACTCCTCATTGCGCCGGCCCCTGCGCAGGTCGTTCTCCAACTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGa                                                                  <  6:159372/135‑1 (MQ=255)                                                                      ggacgaTTCGCGGGGCGTCCAGGCGGTAGGGTGGCTGTAAATGCAACTCGAACTCCTCATGGCGCTGGGCCGGGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGt                                                               <  7:286936/134‑1 (MQ=255)                                                                      cgccgaGTCGCTGGGCGCCCTGGGGGTGGGGTGGCTGTAAATACCCCTCCACCTCCTCTGGGCGCCGGGCCCGCCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGt                                                               <  3:100050/132‑1 (MQ=255)                                                                      cgacgaTTCGCAGGTCGTCCCGGCGGTAGCGTGGGGGGTAATGCAACTCGAACTCCTCATGGCCCCGGGCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGt                                                               <  8:305885/135‑1 (MQ=255)                                                                      cgacgaTCCCCCGGGCGCCCCGGCGGTCGCGGGGCTGGAAATCCAACTCCAACGCCTCATGGCGCCGGCCCCGGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGt                                                               <  7:77344/135‑1 (MQ=255)                                                                      caacgtTTCGCTGGGGGTCCCGGGGTTAGCGTGGCGGTTAATGCAACTCGAACTCCTCTTTGCGCCGGCCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGt                                                               <  2:111342/133‑1 (MQ=255)                                                                         cgaTTCGCAGGTCGTCCAGGCGGTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCCGGGCCCTGCGCCGGTCGTTCTCCCGCGGGCGGCGGGTCGGGGTCTGCGGGTTCATTTAGGCAACGCAGTAgc                                                            >  5:58687/1‑135 (MQ=255)                                                                                cAGGTCGTCCAGGCGGTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGAAGGGCCCTGCGCAGGTAGTTCTCAAGCTGGCGGCGGTACTGGATCTGCTGTTTCATTTCAGAAACGAAGTAGCGACAtgt                                                     >  3:130777/1‑135 (MQ=255)                                                                                              gttAGGGTGGTGGGTAATGAAAGTCGAACTCCTCTTTGCGCTGGGCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCgg                                       <  2:252466/133‑1 (MQ=255)                                                                                              gttAGCGTGGCTGGTATTGATACTCGACCTCCTCATTGGGCTGGGCCGTGGGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCgg                                       <  4:139536/133‑1 (MQ=255)                                                                                              ggTCGCGTGGCTGGTATTGCAACTCGAACTCCTCATTGCGCAGGGCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCgg                                       <  6:265872/135‑1 (MQ=255)                                                                                              ggTAGCGGGGGTGTAAATCAAACTCGAACTCCTCATGGCGCCGGCCCCGGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCgg                                       <  4:169732/135‑1 (MQ=255)                                                                                              ggTAGCGGGGCTTGTAATGCAACTCGAACTCCTCATTCCGCAGGGCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCgg                                       <  6:262272/135‑1 (MQ=255)                                                                                                   gggggCTGGTAATGCACCTCCACCTCCTGCTGGCGCCGGGCCCGGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCGGCCTTg                                  <  3:348282/132‑1 (MQ=255)                                                                                                   cgggggtGGTAAGCCAACTCCAACTCCTCTTGGCCCCGGGCCCGGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCGGCCTTg                                  <  6:151326/129‑1 (MQ=255)                                                                                                   cGTGGCTGTTAAGGCACCGCCAACTCCCCTTGGCGCTGGCGCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCGGCCTTg                                  <  1:344584/135‑1 (MQ=255)                                                                                                    gTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCAGGGCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCGGCCTTgg                                 >  7:320766/1‑135 (MQ=255)                                                                                                        cTGTTACTGCAACTCCAACTCCTCCTGGCGCCGGGCCCGGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCGGCCTTGGCCTg                             <  3:105889/135‑1 (MQ=255)                                                                                                          ggTAATGCAACTCGAACTCCTCATTGCGCAGGGCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCGGCCTTGGCCTGGt                           >  8:73895/1‑135 (MQ=255)                                                                                                                                atggcgcCGGCCCCGGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCGGCCTTGGCCTGGTACAGGGCGATGTCGGCGCAGCg     <  2:41099/132‑1 (MQ=255)                                                                                                                                aTTGCGCCGGGCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCGGCCTTGGCCTGGTACAGGGCGATGTCGGCGCAGCg     <  4:5282/135‑1 (MQ=255)                                                                                                                                    cgcCGGCCCCGGCGCAGGTCGTTCTCCAGCGGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCGGCCTTGGCCTGGTACAGGGCGATGTCGGCGCAGCGGAt  <  1:299901/134‑1 (MQ=255)                                                                                                                                       |                                                                                                                                   ACGGGATGAAGGTGTCCGGCGCCAGCAACCCTTCCTGTGGATGCCGCCAGCGCACCAAGGCCTCGACCGCGACGATTCGCAGGTCGTCCAGGCGGTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCAGGGCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCGGCCTTGGCCTGGTACAGGGCGATGTCGGCGCAGCGGAT  >  ACB122/1313809‑1314073

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: ‑