Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 3,986,621 A→C 33.2% D31E (GAT→GAG) ‡ CJ019_03763 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1223,986,6210AC33.2% 55.5 / 3.8 54D31E (GAT→GAG) ‡CJ019_03763hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (8/24);  new base C (12/4);  total (25/29)
Fisher's exact test for biased strand distribution p-value = 1.61e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACGCCATCAGCACCTGAGCCAGCAACTGATCAGGGTACAAGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGATCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATCCCCATGAGGTCAGACCGCGTGGGTTAAGGGTAGC  >  ACB122/3986496‑3986755
                                                                                                                             |                                                                                                                                      
gACGCCATCAGCACCTGAGCCAGCAACTGATCAGGGTACAAGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGACCACCGCGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGCTCCGGcgcg                                                                                                                               >  4:294987/1‑135 (MQ=255)
gACGCCATCAGCACCTGAGCCAGCAACTGATCAGGGTACAAGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGGGCGGGcgcg                                                                                                                               >  1:364992/1‑135 (MQ=255)
     cATCAGCACCTGAGCCAGCAACTGATCAGGGTACAAGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCCGCGGCGCCCGCGCGGGcgc                                                                                                                          >  1:54185/1‑133 (MQ=255)
      acaaGCACCTGAGCCACCAACTGAGCAGGGAACAAGGCAAGGGGGCCCAGCATCGGGTCCAGTTGCTCCGGGGAGAACCCCGGGTCCTTGGCTACGAGCGGGGCAGCGGAAGGCGCCGGACCCGGCGCGGGCACg                                                                                                                         <  2:364992/132‑1 (MQ=255)
      aTCAGCCCCTGAGCCAGCAACTGATCAGGGTACAAGCCAATGGGGGCCAGCAGCTGGTCCAGTTCCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGGGGCAGCGGAAGGCGCCGGATCCGGCGCGGGCACg                                                                                                                         <  5:232957/135‑1 (MQ=255)
      aTCAGCACCTGAGCCAGCAACTTATCTGGGTACAAGGCACTGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGTACAACCCCGGGTCGTGGGCTACGAGGGTGGCAGCGGAAGGCGCCGGATCCGGCGCGGGCACg                                                                                                                         <  8:151392/135‑1 (MQ=255)
           cACCTGAGCCCGCAACTGATCAGGGTACAAGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCCTCGGAAGGCGCCGGATCCGGCGCGCGCACGTCTGc                                                                                                                    >  6:386259/1‑135 (MQ=255)
           cACCTGAGCCAGCAACTGATCAGGGTACAAGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTCCGAGCGTGGCAGCGGAAGGCGCCCGATCCGGCGCGGGCACGTCTGc                                                                                                                    >  4:278369/1‑135 (MQ=255)
           cACCTGAGCCAGCAACTGATCAGGGTACAAGGCAATGGGGGCCAGCATCGGGGCCAGTTGCGCCTGGGACAACACCGGGGCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGATCCGGCGCGGGCACGTCTGc                                                                                                                    <  7:34767/135‑1 (MQ=255)
           cACCTGAGCCAGCAACGGATCAGGGTACAAGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGATCCGGCGCGGGCACGTCTGc                                                                                                                    <  6:307937/135‑1 (MQ=255)
                 cgCCAGCAACTGATCAGGGTACAAGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGGCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGATCCGGCGCGGGCACGTCTGCGGCCTg                                                                                                              <  1:268688/134‑1 (MQ=255)
                   ccAGCAACTGATCAGGGTACAAGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGATCCGGCGCGGGCCCGTCTGCGGCCGGCa                                                                                                            >  2:286966/1‑135 (MQ=255)
                                aGGGTACAAGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGGCGGCTCCGGCGCGGGGGCGGCGGCGGCCGGCCCCACGCCCCGGGc                                                                                               >  7:241471/1‑135 (MQ=255)
                                    tACAAGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGCTCCGGCGCGGGCCCGTCGGCGGCCGGCACCAGGCCCGGCGCcgcg                                                                                           >  3:125022/1‑135 (MQ=255)
                                     aCAAGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGCGCCCGCGCGGGGCCGGCGGGGGGCAGCGCCAGCCCCGGGGCCgcgc                                                                                          >  3:382303/1‑135 (MQ=255)
                                          gCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGCTCCGGGGCGGGCCCGTCTGCGGCCGGCCCCCGGCCCGGGGCCgcgcggcg                                                                                      >  6:310644/1‑130 (MQ=255)
                                          gCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGCTCCGGCGCGGGCACCTCTGCGGCCTGCACCAGGCCAGGGGCCGCGCGCAg                                                                                      >  8:458507/1‑134 (MQ=255)
                                            aaTGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGATCCGGCGCGGGCACGTCGGCGGCCTGCCCCCCGCCCGGGGCAGCGCGCAGGc                                                                                    >  2:119438/1‑134 (MQ=255)
                                               gggggCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGGCGGATCCGGCGCGGGCACCGCGGCGGCCTGCACCAGGGCAGGGGCAGCGCGCAGGCaaa                                                                                 >  6:244376/1‑134 (MQ=255)
                                                 gggCCAGCATCTGGTCCAGTTGCTCCTGGGAGCACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCCCCGGATCCGGAGCGGGCACGTCTGCGGCCGGCGACACGCCACCGGCAGCGCGCAGGCAAAgaa                                                                              >  8:431991/1‑133 (MQ=255)
                                                 gggCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACCCCGGGTCCTTGGCTACGAGCGTGGCAGCGGACGGCGCCGGCTCCGGCGCGGGCCCGGCTGCGGCCGGCCCCAGGCCCGGGGCAGCCCGCAGGCACAGCa                                                                              >  2:146875/1‑135 (MQ=255)
                                                 gggCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGGTCCGGCGGGGGCACGTCTGCGGCCTGCACCAGGCCCGGGGCAGCGCTCAGGCCAAGCa                                                                              >  3:66998/1‑135 (MQ=255)
                                                 gggCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGATCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCa                                                                              >  4:85419/1‑135 (MQ=255)
                                                      aGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGGTCCGGCGCGGGGCCGGCAGCGGCCCCCCCCACCCCCCGGGCCGCGCGGAGGCCAAGCAGTGcc                                                                         >  7:264747/1‑135 (MQ=255)
                                                      aGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGCTGCGGCGGGGGCCCCTCTGCGGCCGGGGCCCCGGCCGGGGCCGCGCGCCGGCCACGCAGGGcc                                                                         >  5:438779/1‑135 (MQ=255)
                                                      aGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGCTCCGGCGCGGGCCCGGCTGCGGCCGGCCCCACGCCCCGGGCCGCGACCCGGCCAAGCAGTGcc                                                                         >  4:198039/1‑135 (MQ=255)
                                                       gCATCGGGTCCAGTTGCGCCTGGGACAACACCGGGGCCTGGGCTCCGAGCGGGGCCGCGGACGGCGCCGGACCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCg                                                                        <  6:413185/135‑1 (MQ=255)
                                                          ccTGGCCCCGTTGCCCCGGGGACAACACGGGGGCCTTGGCTACGAGCGGGGCAGCGGAAGGCGCCGGAGCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACa                                                                     <  5:76001/134‑1 (MQ=255)
                                                                     tGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGATCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCaa                                                            <  6:5868/133‑1 (MQ=255)
                                                                      gCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGCTCCGGCCCGGGCCCCTCTGCGGCCAGCACCCGGCCACGGGCAGCGCGCACGCAAAGCAGAGCCGCCAGGAAAGCAAAc                                                          >  2:299040/1‑134 (MQ=255)
                                                                      gCTCCTGGGAGAACACCGGGTCCGTGGCTAAGAGCGTGGCAGCGGACGGAGCCGGCTCCGGCGCCGGCACCGCATCGGCCATCCCCACCCCCTGGGCCTCGCTCATGCAAAGACCTGCCCACAGCAAAGCAAAc                                                          >  6:27956/1‑134 (MQ=255)
                                                                           tggcacaaCCCCGGGTCCTTGGCACCGCGGGTGCCCGCTGACGGCGCCGGAGCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGc                                                    <  3:71109/129‑1 (MQ=255)
                                                                             ggAGAACACGGGGTCCTTGGCTACGAGCGGGGCAGCGGAAGGCGCCGGATCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCAt                                                  <  2:146334/135‑1 (MQ=255)
                                                                             ggAGAACACGGGGGCCTTGGCTACGAGCGGGGCAGCGGAAGGCGCCGGATCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCAt                                                  <  8:354674/135‑1 (MQ=255)
                                                                               acaaCCCCGGGGCCTTGGCTACGAGCGGGGCAGCGGACGGCGCCGGAGCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATgc                                                <  1:41144/133‑1 (MQ=255)
                                                                                      ggggTCCTTGGCTCCGGGCGTGGCCGCGGACGGCGCCGGCGCGGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGc                                         <  5:98741/134‑1 (MQ=255)
                                                                                       gggTCATTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGAGCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCt                                        <  6:319582/135‑1 (MQ=255)
                                                                                        tgTCCTTGGCTACGAGCGGTGCCGGGGAAGGGGCCGGAGCCGGGGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCtt                                       <  3:267111/134‑1 (MQ=255)
                                                                                        gggcattGGCTACGAGGGGGGCTGCGGACGGCGCGGTATCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCtt                                       <  6:53627/130‑1 (MQ=255)
                                                                                          tCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGATCCGGCGCGGGCACCTCTGCGGCCTGCACCAGGCCAGGGGCCGCGCGCAGGCAAAGCACTGCCGCCACGAAAGCAAACCATGGCAAGCGCATCGCGTAt                                     >  4:66052/1‑135 (MQ=255)
                                                                                                           ggccggggACGGCGGCGGACCGGGGGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATCCCCATGAGGTCAGAcc                    <  2:2172/129‑1 (MQ=255)
                                                                                                           ggCTGCGGCCGGCGCCGCCCCCGCCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATCCCCATGAGGTCAGAcc                    <  7:267449/135‑1 (MQ=255)
                                                                                                           ggCCGCGGCCGGCGCCGGAGCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATCCCCATGAGGTCAGAcc                    <  7:51382/135‑1 (MQ=255)
                                                                                                           ggCAGGTGACGGCGCCGGCCCGGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATCCCCATGAGGTCAGAcc                    <  3:294987/135‑1 (MQ=255)
                                                                                                           ggCAGGGGCCGGCGCCGGCGCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATCCCCATGAGGTCAGAcc                    <  1:182898/135‑1 (MQ=255)
                                                                                                           ggCAGCGGAAGGCGCCGGTGCCGCCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATCCCCATGAGGTCAGAcc                    <  6:343006/135‑1 (MQ=255)
                                                                                                           gcccggggACGGCTCCGGAGGCGGGGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATCCCCATGAGGGCAGAcc                    <  7:24644/129‑1 (MQ=255)
                                                                                                             cAGCGGAAGGCGCCGGGTCCCGCGCGGGGCCGTCTGCGGCCTGCACCCGGCCACGGGCCGCGCTCCCGGAAACCAGGGCCGCCAGGCAATCAAACAATGGCATGCGCCACGCGTATACCCATGAGGTGAGAcccc                  >  8:267449/1‑133 (MQ=255)
                                                                                                             cAGCGGAAGGCGCCGGGTCCCGCGCGGGCCCGTCTGCGGCCGGCACCCGGCCCGGGGGCGCGCGCAGGGCAACCAGTGCCGACACGAAAGCAAACAATGGCATGCGCATCGCGGATACCCCAGAGGTCAGACcgc                  >  2:19804/1‑135 (MQ=255)
                                                                                                                 gtaaGGCGCCGGATCCGGCGCGGGCACGTCTGCGGCCTGCACAAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATCCCCATGAGGTCAGACCGCGTgg              <  8:202848/133‑1 (MQ=255)
                                                                                                                 ggAGGGCGCCGGACCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATCCCCATGAGGTCAGACCGCGTgg              <  6:418273/135‑1 (MQ=255)
                                                                                                                 ggACGGCGCCGGATCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATCCCCATGAGGTCAGACCGCGTgg              <  8:241471/135‑1 (MQ=255)
                                                                                                                     ggcggcGTAGCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATCCCCATGAGGTCAGACCGCGTGGGTTa          <  3:66052/135‑1 (MQ=255)
                                                                                                                       cgcCGGAGCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGTAAAGCAAACAATGGCATGCGCATCGCTTATCCCCATGAGGTCAGACCGCGTGGGTTAAg        <  4:409171/135‑1 (MQ=255)
                                                                                                                             aTCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATCCCCATGAGGTCAGACCGCGTGGGTTAAGGGTAGc  <  4:45608/135‑1 (MQ=255)
                                                                                                                             |                                                                                                                                      
GACGCCATCAGCACCTGAGCCAGCAACTGATCAGGGTACAAGGCAATGGGGGCCAGCATCTGGTCCAGTTGCTCCTGGGAGAACACCGGGTCCTTGGCTACGAGCGTGGCAGCGGAAGGCGCCGGATCCGGCGCGGGCACGTCTGCGGCCTGCACCAGGCCAGGGGCAGCGCTCAGGCAAAGCAGTGCCGACAGGAAAGCAAACAATGGCATGCGCATCGCTTATCCCCATGAGGTCAGACCGCGTGGGTTAAGGGTAGC  >  ACB122/3986496‑3986755

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: