Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB143 3,887,510 G→A 7.4% R393C (CGC→TGC)  CJ019_03679 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1433,887,5100GA 7.4% 162.6 / 3.1 54R393C (CGC→TGC) CJ019_03679hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (20/30);  new base A (2/2);  total (22/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.07e-01

CGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAACAGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAG  >  ACB143/3887375‑3887643
                                                                                                                                       |                                                                                                                                     
cGATGCCCGGGCGTTCAACCTCAGTGGCATCGGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAACAGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACcgcg                                                                                                                                       <  6:251669/136‑1 (MQ=255)
 gATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAACAGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACcgcg                                                                                                                                       >  1:346969/1‑135 (MQ=255)
  aTGCCCGGGCCTTCAACCTCACTGGCATAAGGCATCTCACAGTCAATCGACTCAAGGATGTCGGTCATGGTCAACAGCCCGGTGAAGTCACCGAACTCGGTGACCACGAAGGCAATATGGGTCGACTGACCGCGCa                                                                                                                                     >  2:129245/1‑136 (MQ=255)
      ccGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAACAGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTg                                                                                                                                 >  4:3563/1‑136 (MQ=255)
      ccGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAACAGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTg                                                                                                                                 >  7:381082/1‑136 (MQ=255)
        gggCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAACAGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCt                                                                                                                               <  6:41231/136‑1 (MQ=255)
        gggCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAACAGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCt                                                                                                                               <  6:45734/136‑1 (MQ=255)
        gggCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAACAGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCACATTTGCt                                                                                                                               >  3:323343/1‑136 (MQ=255)
        gggCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAACAGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCACATTTGCt                                                                                                                               >  3:323340/1‑136 (MQ=255)
         ggCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAACAGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTc                                                                                                                              <  6:199916/136‑1 (MQ=255)
                aaCCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAACAGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCa                                                                                                                       >  3:89292/1‑136 (MQ=255)
                  ccTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAACAGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCAtt                                                                                                                     <  6:85008/136‑1 (MQ=255)
                        tGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAACAGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGAt                                                                                                               >  2:347581/1‑136 (MQ=255)
                         ggCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAACAGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATAGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATg                                                                                                              >  3:180980/1‑136 (MQ=255)
                          gCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAACAGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGc                                                                                                             <  5:306579/136‑1 (MQ=255)
                          gCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAACAGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGc                                                                                                             <  8:182763/136‑1 (MQ=255)
                          gCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAACAGCCCGGGGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATg                                                                                                              <  2:69491/135‑1 (MQ=255)
                            aTCCGGCTGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAACAGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTg                                                                                                           <  5:369750/136‑1 (MQ=255)
                                     cTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAACAGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTcc                                                                                                  <  5:169786/136‑1 (MQ=255)
                                              aaTCGACTCAAGGATGTCGGTCATGGTCAACAGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTTGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTc                                                                                         <  4:54849/136‑1 (MQ=255)
                                                  gACTCAAGGATGTCGGTCATGGTCAACAGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACg                                                                                     >  2:54492/1‑136 (MQ=255)
                                                  gACTCAAGGATGTCGGTCATGGTCAACAGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACg                                                                                     >  2:48284/1‑136 (MQ=255)
                                                   aCTCAAGGATGTCGGTCATGGTCAACAGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACgg                                                                                    <  3:3563/136‑1 (MQ=255)
                                                        aGGATGTCGGTCATGGTCAACAGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGccc                                                                               <  5:152673/136‑1 (MQ=255)
                                                          gATGTCGGTCATGGTCAACAGCCCGGTGAAGTCACCGAACACGTTGACCACGAAGGCAATATTGGTCGACTGACCGCGCATTTGATCCAGGGCAGTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCACgg                                                                             >  1:369339/1‑136 (MQ=255)
                                                                 gTCATGGTCAACAGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGc                                                                      >  8:260381/1‑136 (MQ=255)
                                                                 gTCATGGTCAACAGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGc                                                                      >  2:331965/1‑136 (MQ=255)
                                                                      ggTCAACAGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTc                                                                  >  4:33488/1‑135 (MQ=255)
                                                                          aaCAGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGCCCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCTCCCGGGCCATGCTTTCCAGGt                                                             >  5:164871/1‑136 (MQ=255)
                                                                                ccgcgtaaGTCACCTAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGTTCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAAGGGCTACCGGGCCATTCTTTCCAGGTCGggct                                                       <  7:162183/130‑1 (MQ=255)
                                                                                            ccGAACTCTATGACCACGAAGTCAATATGGTTCTAATGACCGAGAATTTGCTCCATGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGATGCCCgaga                                           <  4:234617/136‑1 (MQ=255)
                                                                                            ccGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCgaga                                           <  7:301959/136‑1 (MQ=255)
                                                                                            ccGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCgaga                                           <  2:346969/136‑1 (MQ=255)
                                                                                            ccGAACTCATTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCga                                             >  2:149427/1‑134 (MQ=255)
                                                                                              gAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGagag                                          >  5:73465/1‑135 (MQ=255)
                                                                                                 cTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGagcagt                                      >  6:367860/1‑136 (MQ=255)
                                                                                                 cTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGagcagt                                      <  1:232648/136‑1 (MQ=255)
                                                                                                 cTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGagcagt                                      <  4:343134/136‑1 (MQ=255)
                                                                                                     ttGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGagcagttcct                                  >  8:42032/1‑136 (MQ=255)
                                                                                                       gACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGagcagttcctt                                 <  6:317032/135‑1 (MQ=255)
                                                                                                        aCCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGagcagttccttca                               <  3:33754/136‑1 (MQ=255)
                                                                                                        aCCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGagcagttccttca                               <  8:381082/136‑1 (MQ=255)
                                                                                                        aCCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGagcagttccttca                               <  5:279534/136‑1 (MQ=255)
                                                                                                             gAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCAGTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGagcagttccttcagcagt                          <  3:144958/136‑1 (MQ=255)
                                                                                                               aGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGagcagttccttcagcagttc                        >  2:376295/1‑136 (MQ=255)
                                                                                                                       tGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTtgt                <  6:372384/136‑1 (MQ=255)
                                                                                                                         tgtcgtcTGACCGCACATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCGAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGc              <  4:323343/135‑1 (MQ=255)
                                                                                                                         ggTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGc              <  4:115891/136‑1 (MQ=255)
                                                                                                                         ggTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGc              <  2:260120/136‑1 (MQ=255)
                                                                                                                         ggTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGc              <  1:149427/136‑1 (MQ=255)
                                                                                                                         ggTCGACTGACCGCACATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGc              <  4:323340/136‑1 (MQ=255)
                                                                                                                             gACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGcaca           >  8:161700/1‑135 (MQ=255)
                                                                                                                               cTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCGATGCTTTCCAGGTCGGGCTGGCTGACCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTTCACAAAg        <  8:258395/136‑1 (MQ=255)
                                                                                                                                tGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGc       <  5:170412/136‑1 (MQ=255)
                                                                                                                                 gACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGcc      >  6:62066/1‑136 (MQ=255)
                                                                                                                                  aCCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGc       <  3:46567/134‑1 (MQ=255)
                                                                                                                                     gcgcATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCCTGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAg  >  7:62410/1‑136 (MQ=255)
                                                                                                                                       |                                                                                                                                     
CGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAACAGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAGCAGTTCCTTCTTGTGCACAAAGCCCAAG  >  ACB143/3887375‑3887643

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑