Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 3,308,725 A→C 27.6% L164R (CTG→CGG)  CJ019_03110 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1223,308,7250AC27.6% 41.8 / 4.5 49L164R (CTG→CGG) CJ019_03110hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (22/11);  new base C (2/11);  total (24/25)
Fisher's exact test for biased strand distribution p-value = 2.68e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCCGGCCAGCGCCTGCAGGCCCAGGTGTTGTTGCAACCAGGCGCCGAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCACCTGTTCCACGACCGCATTGATCATGCCCATGGCCAGCGTGGCCGCCAACAG  >  ACB122/3308594‑3308840
                                                                                                                                   |                                                                                                                   
tGCCGGCCAGCGCCTGCAGGCCCAGGTGTTGTTGCAACCAGGCGCCGAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCCCCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGc                                                                                                                  >  3:244979/1‑135 (MQ=255)
tGCCGGCCAGCGCCTGCAGGCCCAGGTGTTGTTGCAACCAGGCGCCGAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGc                                                                                                                  >  5:186376/1‑135 (MQ=255)
     gCCAGCGCCTGCAGGCCCAGGTGTTGTTGCAACCAGGCGCCGAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCACCTGCCCCCGCCCGCCGTTGCGCCGCTCGCCCAGCCGCCAACGCCCCAGGCCGCgg                                                                                                             >  5:395575/1‑135 (MQ=255)
                  ggCCCAGGTGTTGTTGCACCCAGCCGCCGAGGGGGGCGCCCAGCAACATGCCGGCCCAAGCCACCAACTGCACCCCCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCCGGCCGCGGCCGAAACGCGACa                                                                                                <  3:122557/135‑1 (MQ=255)
                        ggTGTTGTTGCAACCAGGCGCCGAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTCCACCCGCCAGCCGTGGCGCCGCTCGCCCAGCAGCCACCGCCCCAGGCCGCGGCCGAACCGCGACAGGgcgc                                                                                          <  6:128254/135‑1 (MQ=255)
                         gtgtTGTTGCAACCAGGCGCCGAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGgcgc                                                                                          <  8:90563/134‑1 (MQ=255)
                           gttgttGCAACCAGGCGCCGAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCCGGCCGCGGCCGAAACGCGACAGGGCGCCCg                                                                                       >  4:343241/1‑135 (MQ=255)
                              tttgCACCCAGGCGCCGCGCGGGGCGCCCAGCAACATGCCGCCCCAAGGAACAAACTGCCCCCGCCAGCCGTTGCGCCGCTCGCCCACCAGCCAACGCCCCTGGCCGCGGCCGAAACGCGACAGGGCGCCCGtt                                                                                     <  4:146629/133‑1 (MQ=255)
                                tgCAACCAGGCGCCGAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCCGCCAACGCCCCAGGCCGCGGCCGAACCGCGACAGGGCGCCCGGTAc                                                                                   >  5:240417/1‑134 (MQ=255)
                                tgCAACCAGGCGCCGAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAACCGCGACAGGGCGCCCGTTACg                                                                                  >  8:388822/1‑135 (MQ=255)
                                  cAACCAGGCGCCGAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCCGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAACCGCGACAGGGCGCCCGTTCCGt                                                                                 >  7:351152/1‑134 (MQ=255)
                                   aaCCCGGCGCGGCGCGGGGCGCCAAGAAACATGCCGGCCCATCGCACAAACTCCCCCCGCCACCGGTTGCGCCGCGCGCCCACCAGCCAACGCCCCCGGCGGCGGCCTAAACGCGACAGGGCGCCCGTTACGTaa                                                                               <  6:263923/135‑1 (MQ=255)
                                       aGGCGCCGAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCCGGCCGCGGCCGAAACGCGACAGGGCGCCCGGTACGTAAGTCa                                                                           >  2:294769/1‑135 (MQ=255)
                                       aGGCGCCGAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAg                                                                          >  7:161359/1‑136 (MQ=255)
                                        ggCGCCGCGCGGGGCGCCCACCAAATTGCCGGCCCACGCCACCAACTGCACCCGCCAGCCGTTGCGCCGCGCGCCAAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAgg                                                                         <  8:161359/136‑1 (MQ=255)
                                        ggCGCCGCGCGCGGCGCCCACAAACATGCCGGCCCAAGGCACCAACGCCACCCGCCCGCCGTGGCGCCGCGCCCCCAGCACCAAACGCCCCCGGCCGCGGCCGAACCGCGACAGGGCGCCCGTTACGTAAGTCAgg                                                                         <  8:259575/136‑1 (MQ=255)
                                            ccGAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCg                                                                      >  2:57972/1‑135 (MQ=255)
                                            ccGAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCg                                                                      >  1:94338/1‑135 (MQ=255)
                                            ccGAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCg                                                                      >  1:336555/1‑135 (MQ=255)
                                               cggggggcgCCCAGCAAATTGCGGGCCCACGGAAACAACTGCACCCGCCAGGCGGTGGGCCGCGCGCCCTGCAGCCACCGCCCCGGGCCGCGGCGGAAACGCGACAGGGCTCCCGTTACGTAAGTCAGGCCGATgg                                                                  <  8:103119/131‑1 (MQ=255)
                                               aGCGCGGCGCCCAGAAACATGCCGGCCCAAGGCACCAACTCAACCCGCCAGCGGTTGCGCCGCGCGCCCAGCACCCAACGCCCCAGGCCGGGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATgg                                                                  <  7:78972/136‑1 (MQ=255)
                                                gcgcGGCGCCCAGCACCATGCGGGCCCACTGAACCAACTCACCCCGCCACCCGGGGCGCCGCTCGCCCACCAGCCAACGCCCCAGGCCGCGGCCTAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATgg                                                                  <  2:230342/135‑1 (MQ=255)
                                                gcgcGGCGCCCAGCAACATGCCGGCCCAAGCCACCCACTGCACCCCCCAGCCGGCGCGCCCCGCGCCCAGCCGCCACCGCCCCCGGCCGCGGCCGAACCGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATgg                                                                  <  5:125347/135‑1 (MQ=255)
                                                     gcgcCCAGACCCATGCGGGCCACCGGCACCAACTGCCCCCGCCCCCGGTTGCGCCGCGCCCCCACCTGCCACCCCCCCCGGCCGGGGCCCAACCGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAgg                                                             <  1:294769/135‑1 (MQ=255)
                                                     gcgcCCACCAACATGCCGCCCAACGCCCCCAACGGCACCCCCCCGCGGTGGCGCCGCGCCCCCCGCCCCCACCGCCCCCGGCCGCGGCCCAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAgg                                                             <  6:240417/135‑1 (MQ=255)
                                                      cgcCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGtt                                                                                     >  8:36142/1‑110 (MQ=255)
                                                      cgcCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGtt                                                                                     <  7:36142/110‑1 (MQ=255)
                                                           cccAACAGGCCGGCCCACCCCACCAACTCCCCCCCCCACCCGTGGCGCCCCCCCCCCAGCCGCCCCCCCCCCCGCCCGCGGCCGAACCGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTc                                                       <  5:53999/133‑1 (MQ=255)
                                                           aGCAACATGCCGGCCCAAGGCACAAACTCCACCCTCACCCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTc                                                       <  6:186803/135‑1 (MQ=255)
                                                                     cGGCCCAAGGAACCAACTGCCCCCGCCAGCCGTTGCGCGGCTCGCCCCGCACCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCACCTGTTCCa                                             <  6:191148/135‑1 (MQ=255)
                                                                      ggCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCACCTGTTCCa                                             >  3:351719/1‑134 (MQ=255)
                                                                        cccAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCCGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCCGGCCATTCACCTTTTCCACGa                                          >  6:283322/1‑135 (MQ=255)
                                                                        cccAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCACCTGTTCCACGa                                          >  6:200136/1‑135 (MQ=255)
                                                                          cAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCCCCTGTTCCACGccc                                        >  4:227899/1‑135 (MQ=255)
                                                                                aCCAACTGCCCCCGCCAGCCTTGGCGCCGCGCGCCCAGCCCCCAAGGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCACCTGTTCCACGACCGCATTGa                                 <  8:223111/136‑1 (MQ=255)
                                                                                 caaccTGCACCCGCCCGCGGTGGCGCCGCTCGCCGAGCAGCCAACGCCCCCGGCCGCGGCCGAACCGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCACCTGTTCCACGACCGCATTGa                                 <  2:133471/131‑1 (MQ=255)
                                                                                               cAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCCCCTGTTCCACAACCGCATTGATCATGCCCATgg                     >  2:322205/1‑133 (MQ=255)
                                                                                                cgCCGTGCCGCCGCCCCCCCCGCCCCCACCCCCCCCGCCCGCGGCCCAACCGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCACCTGTTCCACGACCGCATTGATCATGCCCATGGCCa                  <  5:101478/134‑1 (MQ=255)
                                                                                                aGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGCAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCCGGCCATTCCCCTGTTCCCCGACCGCATTGCTCATGCCCATGGCCa                  >  2:192007/1‑135 (MQ=255)
                                                                                                aGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGt                                                                                 >  1:103221/1‑72 (MQ=255)
                                                                                                aGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGt                                                                                 <  2:103221/72‑1 (MQ=255)
                                                                                                aGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCACCTGTTCCACGACCGCATTGATCATGCCCATGGCCa                  >  3:342472/1‑135 (MQ=255)
                                                                                                aGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCACCTGTTCCACGACCGCATTGATCATGCCCATGGCCa                  >  5:342445/1‑135 (MQ=255)
                                                                                                        cgccgcTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCACCTGTTCCACGACCGCATTGATCATGCCCATGGCCAGCGTGgcc          >  4:231590/1‑135 (MQ=255)
                                                                                                        cgccgcTCCCCCCGCCGCCACCGCCCCGGGCGGCGCCCTAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCACCTGTTCCACGACCGCATTGATCATGCCCATGGCCAGCGTGgcc          <  3:303596/135‑1 (MQ=255)
                                                                                                          ccgcTCGCCCAGCAGCCAACGCCCCCGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCACCTGTTCCACGACCGCATTGATCATGCCCATGGCCAGCGTGgccgc        <  6:342445/135‑1 (MQ=255)
                                                                                                          ccgcTCGCCCAGCAGCCAACGCCCCCGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCACCTGTTCCACGACCGCATTGATCATGCCCATGGCCAGCGTGgccgc        <  3:343241/135‑1 (MQ=255)
                                                                                                                gCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCACCTGTTCCACGACCGCATTGATCATGCCCATGGCCAGCGTGGCCGCCAACAg  >  2:238825/1‑135 (MQ=255)
                                                                                                                 ccccGCCGCCAACGCCCCAGGCCGCGCCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCACCTGTTCCACGACCGCATTGATCATGCCCATGGCCAGCGTGGCCGCCAACa   <  4:244979/133‑1 (MQ=255)
                                                                                                                                   |                                                                                                                   
TGCCGGCCAGCGCCTGCAGGCCCAGGTGTTGTTGCAACCAGGCGCCGAGCGCGGCGCCCAGCAACATGCCGGCCCAAGGCACCAACTGCACCCGCCAGCCGTTGCGCCGCTCGCCCAGCAGCCAACGCCCCAGGCCGCGGCCGAAACGCGACAGGGCGCCCGTTACGTAAGTCAGGCCGATGGGCAGGCCATTCACCTGTTCCACGACCGCATTGATCATGCCCATGGCCAGCGTGGCCGCCAACAG  >  ACB122/3308594‑3308840

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