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
*ACB1221,313,9430AC45.4% 6.5 / 18.3 34L614R (CTG→CGG) CJ019_01235hypothetical 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: