Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 1,313,943 A→C 36.0% 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,9430AC36.0% 40.8 / 4.2 41L614R (CTG→CGG) CJ019_01235hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (19/6);  new base C (2/14);  total (21/20)
Fisher's exact test for biased strand distribution p-value = 9.29e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGGGATGAAGGTGTCCGGCGCCAGCAACCCTTCCTGTGGATGCCGCCAGCGCACCAAGGCCTCGACCGCGACGATTCGCAGGTCGTCCAGGCGGTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCAGGGCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCGGCCTTGGCCTGGTACAGGGCGATGTCGGCGCAGCGGA  >  ACB122/1313809‑1314072
                                                                                                                                      |                                                                                                                                 
aCGGGATGAAGGTGTCCGGCGCCAGCAACCCTTCCTGTGGATGCCGCCAGCGCACCAAGGCCTCGACCGCGACGATTCGCAGGTCGTCCAGGCGGTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCAg                                                                                                                                  >  3:296911/1‑136 (MQ=255)
        atggTGGCCGGCGCCAGCAACCCTTCGTGTGGATGCCGCCAGCCCACCACGGCCGCCCCCGCGACTATTCGCTGGGCGTCCCGGGGGTTGGGTGGGTTTTAATCCAACTCCAACTCCTCATTGCGCAGGGcccgg                                                                                                                           <  4:389041/133‑3 (MQ=255)
                 ggCGCCAGCAACCCTTCCTGTGGATGCCGCCAGCGCACCAAGGCCTCGACCGCGACGATTCGCAGGTCGTCCAGGCGGTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCAGGGCCCTGCGCAGGTCg                                                                                                                  >  1:101148/1‑135 (MQ=255)
                 ggCGCCAGCAACCCTTCCTGTGGATGCCGCCAGCGCACCAAGGCCTCGACCGCGACGATTCGCAGGTCGTCCAGGCGGTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCAGGGCCCGGCGCAGGTCg                                                                                                                  >  5:21317/1‑135 (MQ=255)
                           aCCCTTCCTGTGGATGCCGCCAGCGCACCAAGGCCTCGACCGCGACGATTCGCAGGTCGTCCAGGCGGTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCAGGGCCCTGCGCAGGTCGTTCTCCAGCTg                                                                                                       >  7:311754/1‑136 (MQ=255)
                             ccTTCGTGGGCATGCCGCCCGCGCACCAAGGCCTCGCCCGCGACGATGCGCCGGGCGCCCTGGCGGTCGCGTGGCGGTTAATCCAACTCCACCCCCACATGGCGCCGGGCCCGGCGCAGGTCGTTCTCCAGCTgg                                                                                                      <  6:150132/135‑1 (MQ=255)
                                    gtgGATGCCGCCAGCGCACCAAGGTCTCGACCGCGACGATTCGCAGGTCGTCCAGGCGGTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCAGGGCCCAGCGCAGGTCGTTCTCCAGCAGGccgcggga                                                                                              >  4:152898/1‑134 (MQ=255)
                                          gccgccAGCGCACCAAGGCCTCGACCGCGACGATTCGCAGGTCGTCCAGGCGGTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGAGCCCGGCCCGGCGCAGGTCGTTCTACAGCTGGCGGCGGTGCTggt                                                                                         >  5:43745/1‑134 (MQ=255)
                                           ccgcccgcGCACCAAGGCCTCGACCGCGACGATTCGCAGGTCGTCCAGGCGGTAGCGTGGCTGGTAATGAAAAGATAACTCCTCATTGCGACGGGCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGAt                                                                                        <  6:313486/135‑1 (MQ=255)
                                            cgccAGCGCACCAAGGCCTCGACCGCGACGATTCGCAGGTCGTCCAGGCGGTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCAGGGCCCGGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATc                                                                                       >  5:73820/1‑135 (MQ=255)
                                               cAGCGCACCAAGGCCTCGACCGCGACGATTCGCAGGTCGTCCAGGCGGTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCAGGGCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATctgct                                                                                   >  7:93572/1‑136 (MQ=255)
                                                  cgcACCAAGGCCTCGACCGCGACGATTCGCAGGTCGTCCAGGCGGTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCAGGGCCCTGCGCAGGTCGGTCTCCAGCTGGCGGCGGTACTGGATCTGCTgg                                                                                 >  2:257116/1‑135 (MQ=255)
                                                            ccgcgcccGCGCCGTTTCGCTGGTCGGCCCGGGGGTAGCGGGGCGGGTAATCCACCTCGAACGCCTCATTGCGCCGGGCCCGGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCgg                                                                       <  6:317447/129‑1 (MQ=255)
                                                            ccTCGCCCGCGCCGTTTCGCTGGGCGGCCTGGGGGTCGCGGGGCTGGTAATCCAACTCGAACCCCTCATGGCGCCGGGCCCGGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCgg                                                                       <  2:272029/135‑1 (MQ=255)
                                                                  ccGCGACGATTCGCAGGTCGTCCAGGCGGTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCAGGGCCCTGCGCAGGGCGTTCTCCAGCTGGGGGCGGGACTGGCTCTGCTGGTTGATTTAGGCAAc                                                                   >  2:148871/1‑133 (MQ=255)
                                                                        cgaTTCGCAGGTCGTCCAGGCGGTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCCGGGCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGGTCATTTCGGCAACGAAGTAgc                                                           >  1:23561/1‑135 (MQ=255)
                                                                        cgaTTCGCAGGTCGTCCAGGCGGTAGCGTGGCTGGTAAAGCAACTCGAACTCCTCATTGCGCAGGGCCCTGCGCAGGTCGTTCTCCAGCTGGCCGCGGTACTGGATCTGCTGGTTCCTTTCGGCCACGAAGTAgcg                                                          >  8:379855/1‑136 (MQ=255)
                                                                                      tcCTGGCGGTAGCGGGGCTGGAAACGCAACTCCACCTCCTCATGGCGCCGGGCCGGGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTc                                            <  8:268709/136‑1 (MQ=255)
                                                                                      gcccggCGGTCGCGGGGGTGCTAATGCAACTCGACCGCCTCGTGGCGCCGGGCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTc                                            <  3:228558/132‑1 (MQ=255)
                                                                                       cccggcgggcgCGTGGCTGGTAAAGCACCTCCACCTCCTCATGGCGCCGGGCCCGGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTc                                            <  5:134389/132‑1 (MQ=255)
                                                                                            cGGTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCAGGGCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCgg                                      <  4:61761/136‑1 (MQ=255)
                                                                                             ggTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCAGGGCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCGGc                                     >  7:368910/1‑136 (MQ=255)
                                                                                             ggTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCAGGGCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATATCGGCAACGAAGTAGCGCCATGGGTTCTTTCCATCgg                                      >  5:214757/1‑135 (MQ=255)
                                                                                             ggTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCACAGGGCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCgg                                      <  6:9924/135‑1 (MQ=255)
                                                                                                 gggtgggtggAAATGAAACTCGACCTCCTCATTGCGCAGGGCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCGGCCTTg                                 <  8:368910/134‑1 (MQ=255)
                                                                                                 gCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCACAGGGCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCGGCCTTg                                 <  8:362161/136‑1 (MQ=255)
                                                                                                 gCGTGGCTGGTAATGCAACTCGAACTCCTCATGGCGCCGGCCCCGGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCGGCCTTg                                 <  7:335608/136‑1 (MQ=255)
                                                                                                  gggggCTGTTAACCCAACCCCACCCCCTCATGGCGCCGGGCCCGGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCGGCCTTg                                 <  5:350394/132‑1 (MQ=255)
                                                                                                  gggggCTGTAAATTCACCTAGAACACCTCATTGCCGAGGGCCGTTCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCGGCCTTg                                 <  6:214757/132‑1 (MQ=255)
                                                                                                   gTGGCTGTTATTGCACGTCTAACTCCTCATTGCGCCGGGCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCGGCCTTg                                 <  2:101148/134‑1 (MQ=255)
                                                                                                   gTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCAGGGCCCTGCGCAGGGCGGTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCAATTCGGCAACGAAGTAGCCCCCTGTGTTCTTTCCAACGGCCGTgg                                >  2:151440/1‑135 (MQ=255)
                                                                                                        tGGTAATGCAACTCGAACTCCTCATTGCGCAGGGCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCGGCCTTGGCCTgg                           >  7:365159/1‑135 (MQ=255)
                                                                                                            aaTGCAACTCGAACTCCTCATTGCGCAGGGCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCGGCCTTGGCCTGGTACa                       >  2:362037/1‑135 (MQ=255)
                                                                                                             aggCAACTCGAACTCCTCATTGCGCCGGGCCCGGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCGGCCTTGGCCTGGTACAg                      <  4:268902/133‑1 (MQ=255)
                                                                                                                        aCTCCTCATTGCGCAGGGCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCGATGTGTTCTTTCCATAGGCCTTGGCCTGGTACAGGGCGATGTCGgc          >  7:196637/1‑136 (MQ=255)
                                                                                                                          tccgcATGGCGCCGGGCCCGGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCGGCCTTGGCCTGGTACAGGGCGATGTCGgcg         <  6:21317/135‑1 (MQ=255)
                                                                                                                               cgcgcgcCGGCCCCGGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCGGCCTTGGCCTGGTACAGGGCGATGTCGGCGCAGCg    <  5:390341/132‑1 (MQ=255)
                                                                                                                               atgccgcCGGCCCCTGCGCAGGGCGTTCTCCCGCTGGCGGCGGTACTGGATCTGCTGGGTCATTTCGGCAACGAAGTAGCGCCATGTGGTCTTTCCATCGGCCTTGGCCTGGTACAGGGCGATGTCGGCGCAGCg    <  6:15666/131‑1 (MQ=255)
                                                                                                                                ttGCGCAGGGCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCGGCCTTGGCCTGGTACAGGGCGATGTCGGCGCAGCgg   >  2:102168/1‑135 (MQ=255)
                                                                                                                                ttGCGAAGGGCCCTGCGCAGGTCGTTCTACAGCTGGCGGCGGTACTGGATCTGCTGGTTAATTTCGGCAACGAAGTAGCGCCAGGTGTTCTTTCACTAGGCCTTGGCCTGGTACAGGGCGATGTCGCCGCAGCgg   >  2:102171/1‑135 (MQ=255)
                                                                                                                                 tGCGCAGGGCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCGGCCTTGGCCTGGTACAGGGCGATGTCGGCGCAGCGGa  >  2:43221/1‑135 (MQ=255)
                                                                                                                                      |                                                                                                                                 
ACGGGATGAAGGTGTCCGGCGCCAGCAACCCTTCCTGTGGATGCCGCCAGCGCACCAAGGCCTCGACCGCGACGATTCGCAGGTCGTCCAGGCGGTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCAGGGCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCGGCCTTGGCCTGGTACAGGGCGATGTCGGCGCAGCGGA  >  ACB122/1313809‑1314072

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