Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 1,313,950 T→G 33.8% R612R (AGG→CGG)  CJ019_01235 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1221,313,9500TG33.8% 37.3 / 2.3 38R612R (AGG→CGG) CJ019_01235hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (19/5);  new base G (2/12);  total (21/17)
Fisher's exact test for biased strand distribution p-value = 1.52e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTGTCCGGCGCCAGCAACCCTTCCTGTGGATGCCGCCAGCGCACCAAGGCCTCGACCGCGACGATTCGCAGGTCGTCCAGGCGGTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCAGGGCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCGGCCTTGGCCTGGTACAGGGCGATGTCGGCGCAGCGGATCAGCTCA  >  ACB122/1313819‑1314080
                                                                                                                                   |                                                                                                                                  
ggTGTCCGGCGCCAGCAACCCTTCCTGTGGATGCCGCCAGCGCACCAAGGCCTCGACCGCGACGATTCGCAGGTCGTCCAGGCGGTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCAGGGCCCTgc                                                                                                                                  >  2:289158/1‑134 (MQ=255)
    tCCGGCGCCAGCAACCCTTCCTGTGGATGCCGCCAGCGCACCAAGGCCTCGACCGCGACGATTCGCAGGTCGTCCAGGCGGTAGCGTGGCTGGTAATGCAACTCGAACGCCTCATTGCGCAGGGCCCGGCGCAgg                                                                                                                             >  3:22133/1‑135 (MQ=255)
      cGGCGCCAGCAACCCTTCCTGTGGATGCCGCCAGCGCACCAAGGCCTCGACCGCGACGATTCGCAGGTCGTCCAGGCGGTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCAGGGCCCTGCGCAGGTc                                                                                                                           >  1:177785/1‑135 (MQ=255)
       ggCGCCAGCAACCCTTCCTGTGGATGCCGCCAGCGCACCAAGGCCTCGACCGCGCCGATTCGCAGGTCGTCCAGGCGGTAGCTTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCAGGGCCCTGCGCAGGTCg                                                                                                                          >  2:262001/1‑135 (MQ=255)
       ggCGCCAGCAACCCTTCCTGTGGATGCCGCCAGCGCACCAAGGCCTCGACCGCGACGATTCGCAGGTCGTCCAGGCGGTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCAGGGCCCTGCGCAGGTCg                                                                                                                          >  1:252858/1‑135 (MQ=255)
       ggCGCCAGCAACCCTTCCTGTGGATGCCGCCAGCGCACCAAGGCCTCGACCGCGACGATTCGCAGGTCGTCCAGGCGGTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCAGGGCCCTGCGCAGGTCg                                                                                                                          >  6:219454/1‑135 (MQ=255)
            cAGCACCCCTTCCTGTGGATGCCGCCAGCGCACCAAGGCCTCGACCGCGACGATTCGCAGGTCGTCCAGGCGGTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCCGGGCCCTGCGCAGGTCGTtctc                                                                                                                     >  2:275655/1‑135 (MQ=255)
                   cctgcgtgtgGTGGCCGCCAGGTCACCACGGCCGCGACGGCGACGTTTCGCTGGTGGTCCTGGCGGTAGGGGGGCTGGTAATGCAACTCGACCACCTCATGGCGCCGGGCCCGGCGCAGGTCGTTCTCCAGCTgg                                                                                                              <  8:107402/129‑1 (MQ=255)
                   ccTTCTTGTGGCGGCCGCCCGCCCCCCACGGCCGCGCCCGCCACTATGCGCTGGGGGTCCCGGGGGTAGGGTGGCTGGAAATGCACCTCCAACACCTCATGGCCCGGGGCCGTGCCCAGGTCGTTCTCCAGCTgg                                                                                                              <  8:48834/135‑1 (MQ=255)
                   ccTTCCTGTGGATGCCGCCAGCGCACCAAGGCCTCGACCGCGACGATGCGCAGGTCGTCCAGGCGGTAGCGTGGCTGTTAATGCAACTCTAACTCCTCATTGCGCAGGGCCCGGCGCAGGTCGTTCTCCAGCTgg                                                                                                              <  4:80322/135‑1 (MQ=255)
                   ccTTCCTGGGTCTGCCGCCAGGGCACCAGGGCGTCGACCGCGACGATGCGCCGTTGGTCCAGGCGGTCGCGGGGCTGTTAATGCAACTCCAACACGTCTTTGCGCTGGGCGCGGCGCAGGTCGTTCTCCAGCTgg                                                                                                              <  2:31829/135‑1 (MQ=255)
                   ccTTCCTGGGGAGCCCGCCAGCGCACCACGGCCGCGACCGCGACTATGCGCAGGTCGACCCGGCGGTCGGGGGGCTGGAAAGCCAACTCGACCTCCTCATGGCGCCGGGCCCGGCGCAGGTCGTTCTCCAGCTgg                                                                                                              <  2:95449/135‑1 (MQ=255)
                         tgtgGATGCCGCCAGCGCACCAAGGCCTCGACCGCGACGATTCGCAGGTCGTCCAGGCGGTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCAGGGCCCTGCGCAGGTCGTTCTCCAGCTGgcggcg                                                                                                         >  2:295535/1‑134 (MQ=255)
                         tgtgGATGCCGCCAGCGCACCAAGGCCTCGACCGCGACGATTCGCAGGTCGTCCAGGCGGTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCAGGGCCCTGCGCAGGTCGTTCTCCAGCTGGcggcgg                                                                                                        >  3:191077/1‑135 (MQ=255)
                                          cACCAGGCCCTCCCCCGCGACTATGCGCTGGGGGCCCTGGCGGGTGGGGGGCTGGAAATCAAACTCTAACGCCTCATGGCGCTGGGCCCGGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGtt                                                                                       <  4:425358/135‑1 (MQ=255)
                                                  ccgcgcccGCGACGACTCGCGGGGCGGCCCGCCGTACGCGGGGCGGGTAATGCACCACGAACTCCTCATTGCCCCGGGCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCgg                                                                               <  6:217762/129‑1 (MQ=255)
                                                  ccgcGACCGCCACCCTGCGCCGGGCGGCCCGGCGGGCGGGGGGCGGGAAAACCAACGCCAACTCCTCATTGCGCAGGCCCGGGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCgg                                                                               <  3:191460/132‑1 (MQ=255)
                                                  ccccgcccGCGCCGACGCGCCGGTCGTCCCGGCGGTCGCGTGGCTGCTAACCCACCCCCACCGCCTCAGGCCGCGGGGCCCGGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCgg                                                                               <  1:44458/129‑1 (MQ=255)
                                                  ccTCCACCGCCCCGATTCGCCGGGGGTCCCGGCGGTAGCGGGGCTGGTAATCCAACGCCAACCCCTCATTGCGCCGGCCCCGGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCgg                                                                               <  4:151509/135‑1 (MQ=255)
                                                        ccGCGACGATTCGCAGGTCGTCCAGGCGGTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCAGGGCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGGACTGGATCTGCTGGTTCATTTCGGCAAc                                                                           >  6:227698/1‑133 (MQ=255)
                                                        ccGCGACGATTCGCAGGTCGTCCAGGCGGTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCAGGGCCCGGCGCAGGGCGTTCTCCAGCGGGCGGCGGTACTGGATGTGCGGGGTACTTTCgggcac                                                                           >  1:450939/1‑129 (MQ=255)
                                                              cgaTTCGCAGGTCGTCCAGGCGGTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCAGGGCCCTGCGCAGGTCGTTCTGCAGCTGGCGGCGGTACTGGATGTGGTGGTTCATTTAGGCAACGAAGTAgc                                                                   >  6:275156/1‑135 (MQ=255)
                                                                            gcccggCGTTAGGGGGGCTGGAAATCAAACTCTACCACCTCATGGCGCCGGGCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCttt                                                     <  6:337158/131‑1 (MQ=255)
                                                                             cccggCGGTCGCGGGGCTTTTAATCCACCCCGACCCCCTCCTGGCGCCGGGCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTc                                                    <  1:275655/132‑1 (MQ=255)
                                                                                   ggTCGCGGGGCTGTTAATGCACCTCGAACTCCTCTTTGCCCAGGGCCGTGGGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCgg                                              <  2:71641/135‑1 (MQ=255)
                                                                                       gCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCAGGGCCCTGCGCACGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTAGGCAACGAACTAGCGCCATGTGTTCTTTCCATCGGCCtt                                          >  1:71641/1‑135 (MQ=255)
                                                                                        gggggCTGTTAACGCAACTCTACCTCCGCTTGGCGCGGGGCCCGGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCGGCCTTg                                         <  7:455270/132‑1 (MQ=255)
                                                                                        ggTGGCTGTTAACGCAACTCGACCTCCTGTTTGGGCCGGCCCCGGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCGGCCTTg                                         <  7:99429/134‑1 (MQ=255)
                                                                                        cGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCAGGGCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATAGGCCTTg                                         >  7:50349/1‑135 (MQ=255)
                                                                                             cTGTTAATCCAACTCGACCGCCTCCTTGCGCCGGGCCGGGCGCAGGTCGTGCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCGGCCTTGGCCTg                                    <  3:459516/135‑1 (MQ=255)
                                                                                             cTGGTAATGCAACTCGAACTCCTCATTGCGCAGGGCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCGGCCTTGGCCTg                                    >  7:271439/1‑135 (MQ=255)
                                                                                              tGGTAATGCAACTCGAACTCCTCATTGCGCAGGGCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCGGCCTTGGCCTg                                    >  6:147120/1‑134 (MQ=255)
                                                                                               ggTAATGCAACTCGAACTCCTCATTGCGCAGGGCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCGGCCTTGGCCTGGt                                  >  1:30844/1‑135 (MQ=255)
                                                                                                         cTCGAACTCCTCATTGCGCAGGGCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGGACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCGGCCTTGGGCTGGTACAGGGCGAt                        >  1:451375/1‑135 (MQ=255)
                                                                                                                   tcATGGCGGTGGCCCCGGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCGGCCTTGGCCTGGTACAGGGCGATGTCGGCGCAg              <  3:124095/135‑1 (MQ=255)
                                                                                                                       tGCGCAGGGCCCTGCGCAGGTCGTTCTCAAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTAGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCGGCCTTGGCCTGGTACAGGGCGATGTCGGCGCAGCGGa          >  6:154843/1‑135 (MQ=255)
                                                                                                                          gcAGGGCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCGGCCTTGGCCTGGTACAGGGCGATGTCGGCGCAGCGGATCa       >  7:59112/1‑135 (MQ=255)
                                                                                                                               gCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCGGCCTTGGCCTGGTACAGGGCGATGTCGGCGCAGCGGATCAGCTcc  >  3:242536/1‑134 (MQ=255)
                                                                                                                                   |                                                                                                                                  
GGTGTCCGGCGCCAGCAACCCTTCCTGTGGATGCCGCCAGCGCACCAAGGCCTCGACCGCGACGATTCGCAGGTCGTCCAGGCGGTAGCGTGGCTGGTAATGCAACTCGAACTCCTCATTGCGCAGGGCCCTGCGCAGGTCGTTCTCCAGCTGGCGGCGGTACTGGATCTGCTGGTTCATTTCGGCAACGAAGTAGCGCCATGTGTTCTTTCCATCGGCCTTGGCCTGGTACAGGGCGATGTCGGCGCAGCGGATCAGCTCA  >  ACB122/1313819‑1314080

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