Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 1,348,419 T→G 26.4% D4A (GAC→GCC)  CJ019_01261 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1221,348,4190TG26.4% 68.6 / 3.2 53D4A (GAC→GCC) CJ019_01261hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (6/31);  new base G (9/6);  total (15/38)
Fisher's exact test for biased strand distribution p-value = 5.00e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGTCGGTGTAAGGGACGACAGTGCCGTTACGCTTGATCACGCGCAGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCTCCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGGGGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGT  >  ACB122/1348290‑1348550
                                                                                                                                 |                                                                                                                                    
tcgtcgGTGTAAGGGACGACAGTGCCGTTACGCTTGATCACGCGCAGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCgg                                                                                                                                   <  4:277502/133‑1 (MQ=255)
 cggcgGTGTAAGGGACGACAGTGCCGTTACGCTTGATCACGCGCAGTTGGCCGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCAGCTGGGCCTGAGGGTTATCGCGAGTCGTGTCGGTTTg                                                                                                                               <  4:403417/132‑1 (MQ=255)
                cgacAGTGCCGTTACGCTTGATCACGCGCAGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGtc                                                                                                                 <  2:264301/135‑1 (MQ=255)
                 gacAGTGCCGTTACGCTTGATCACGCGCAGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGtct                                                                                                               <  1:264273/135‑1 (MQ=255)
                  acAGTGCGGTTACGCTTGATCACGCGCAGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGtctc                                                                                                              <  6:150848/135‑1 (MQ=255)
                  acAGTGCCGTTACGCTTGATCACGCGCAGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGTCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGtctc                                                                                                              <  5:39919/135‑1 (MQ=255)
                  acAGTGCCGTTACGCTTGATCACGCGCAGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGtctc                                                                                                              <  4:261756/135‑1 (MQ=255)
                  acAGTGCCGTTACGCTTGATCACGCGCAGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGtctc                                                                                                              <  6:231857/135‑1 (MQ=255)
                  acAGTGCCGTTACGCTTGATCACGCGCAGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCGGCGGCAGCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGtctc                                                                                                              <  7:255012/135‑1 (MQ=255)
                   ctgTGCCGTTACGCTTGATCACGCGGAGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGGTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCTcc                                                                                                             <  5:78140/133‑1 (MQ=255)
                    aGTGCCGTTACGCTTGATCACGCGCAGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTGGGGTTGGTTGGGTGGGTGTGTCca                                                                                                            >  7:295795/1‑135 (MQ=255)
                     gTGCCGTTACGCTTGATCACGCGCAGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCTCcac                                                                                                           <  3:155814/135‑1 (MQ=255)
                     gTGCCGTTACGCTTGATCACGCGCAGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGGTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCTCcac                                                                                                           <  6:452630/135‑1 (MQ=255)
                       gCCGTTACGCTTGATCACGCGCAGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCTCCACAg                                                                                                         <  7:240133/135‑1 (MQ=255)
                         cGTTACGCTTGATCACGCGCAGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGGATGGGTGGGTGTGTGCACAgg                                                                                                        >  8:239694/1‑133 (MQ=255)
                         cGTTACGCTTGATCACGCGCAGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGGGGGTGTGTTTGGGTGGGGGTCGACACAGtt                                                                                                       >  4:380102/1‑135 (MQ=255)
                                   gATCACGCGCAGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGGCGGTTTGCATGGGTGGGTGTGTGCCCAGTGTATATAtttt                                                                                             >  7:2944/1‑133 (MQ=255)
                                     tgaCGCGCAGGTGGCCGGGGGCGGTGGCAGCCAGATCCTGGGTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCTCCACAGTTTCTATATTGTTc                                                                                           <  5:92132/133‑1 (MQ=255)
                                     tCACGCGCAGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCTCCACAGTTTCTATATTGTTc                                                                                           <  6:192354/135‑1 (MQ=255)
                                     tCACGCGCAGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCTCCACAGTTTCTATATTGTTc                                                                                           <  1:126531/135‑1 (MQ=255)
                                     tCACGCGCAGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCTCCACAGTTTCTATATTGTTc                                                                                           <  7:224102/135‑1 (MQ=255)
                                     tCACGCGCAGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCTCCACAGTTTCTATATTGTTc                                                                                           <  3:59286/135‑1 (MQ=255)
                                     tCACGCGCAGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCTCCACAGTTTCTATATTGTTc                                                                                           <  1:451029/135‑1 (MQ=255)
                                     tCACGCGCAGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCTCCACAGTTTCTATATTGTTc                                                                                           <  7:380039/135‑1 (MQ=255)
                                     tCACGCGCAGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCTCCACAGTTTCTATATTGTTc                                                                                           <  7:42166/135‑1 (MQ=255)
                                           gcAGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTGGGTTTGCGTGGGTGGGTGGGTGCACCGTTTGTATATTTTTGAGgcgc                                                                                     >  4:395262/1‑135 (MQ=255)
                                           gcAGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTGGGTGTGCATGGGGGGGTGTCGCCGCAGGTTGTATATTTTTTACgcgc                                                                                     >  4:352480/1‑135 (MQ=255)
                                           gcAGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGGGGGTTTGGATGGGTGGGTGTGGGCACAGTTTGTATATTTTTCCGgcgc                                                                                     >  4:213967/1‑135 (MQ=255)
                                            cAGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGGCGGTTTGGCTGGGTGGGTGTCTCCACCGGTTGTAGATTGTTCAGGCGcc                                                                                    >  4:394109/1‑135 (MQ=255)
                                                 ggCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCTCCACAGTTTCTATATTGTTCAGGCGCCTTATg                                                                               <  3:395926/135‑1 (MQ=255)
                                                 ggCCGGGGGCGGTGGCAGCCAGAGCCTGGTTGGAAGCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGAATGGGGGGGGGTCTCAACAGTTTCTATATTGTTCAGGCGCCTTATg                                                                               <  5:355932/135‑1 (MQ=255)
                                                  gCCGGGGGCGGGGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCTCCACAGTTTCTATATTGTTCAGGCGCCTTATg                                                                               <  7:41804/134‑1 (MQ=255)
                                                    cGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGGGGGGGTGGTTGGGTGGGGGTGTCCCCAGATTCTATATTTTTATGGCGCCGTATGGgc                                                                            >  2:128816/1‑135 (MQ=255)
                                                          cGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGGGGGTTTGGCTGGGGGGGTGGCGTCCAAGTTTCTTATTTTTTAAGGCGCCGTATTGGCGCCCAc                                                                      >  2:115859/1‑135 (MQ=255)
                                                          cGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGGGGGGTGGGTTGGGTGGGTGTGTCCGCAGGTTATATTTTTTTCTAGCGGCGTTTGGGCGccccc                                                                      >  6:328661/1‑133 (MQ=255)
                                                          cGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGGGGGGGTGGTTGGGTGGGGGTGTGCACCGGTTCGATATTTTTAAGGCGGCGTATGGGCGgcccc                                                                      >  4:266241/1‑133 (MQ=255)
                                                            gTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTGGGTTTGGGTGGGTGGGTGTGTCCACAGTTTGATTATTTTTAAGGGGCCTTATGGGGGCCCCCCg                                                                    >  4:45843/1‑135 (MQ=255)
                                                              ggCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCTCCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGtt                                                                  <  8:447567/135‑1 (MQ=255)
                                                                    cAGATCCTGGTTGGAAGCGGCGGCCTGCGGCACCGGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCTCCACAGTTTATATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTcc                                                            <  1:286135/135‑1 (MQ=255)
                                                                    cAGATCCTGGTGGGAAGCGGCGGCCTGCGGCACCTGGGCCGGCGGTTTCTCGCGAGGCGTGGCGGTTTGCAGGGGGGGGGGTCACCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTcc                                                            <  5:436042/135‑1 (MQ=255)
                                                                    cAGATCCGGGTTGGAAGCGGCGGCCGGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTGGCATGGGGGGGTGGCTCCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTcc                                                            <  6:427167/135‑1 (MQ=255)
                                                                                     cggcggCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGGGGGGTTGGCTGGGTGGGGGTGGCCCCCGGTTGTTTTTTGTTGATGCGCCGTATGGGCGCCCCCCCTTCCGGCCCCGTGCTCGGTAAccac                                           >  6:82691/1‑133 (MQ=255)
                                                                                       gcggcCTGCGGCACCTGGGCCTGCGGTTTCTCGCGAGTCGTGTCGGTTTGCAGGGGTGGGGGTCTCCACAGTTTATATTTTGTTCAGGCGCCTAATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTcgc                                         <  6:134477/135‑1 (MQ=255)
                                                                                         ggcgggCGGCCCCGGGGCGTGGGGGTTCTCGCGAGTCGTGTCGGTTGGCCTGGGGGGGTGTCTCCACTGTTTATATTTTGTTCCGGCGCCTTTTGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCgg                                       <  1:226510/130‑1 (MQ=255)
                                                                                         ggcCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTGGGGTTGCATGGGTGGGTGTGGACACAGTTTGGTTATTGTTCAGGCGCCTTATGGGCGGCCACCGGTCCGGCCACCTGCTCGATAACCTCCCgg                                       >  4:356084/1‑135 (MQ=255)
                                                                                           ccTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCTCCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCgggg                                     <  7:239694/135‑1 (MQ=255)
                                                                                             tGCGGCACCTGGGCCTGCGGGTTCTCGCGAGGCGTGGCGGTTTGCATGGGTGGGTGGCTCCACAGTTTGTATATTGTGCCGGCGCCTTAGGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGGGGAt                                   <  1:123662/135‑1 (MQ=255)
                                                                                             ggCGGCACGTGGGCCGGCGGGTGCTCGCGAGGCGTGGCGGTTTGCTTGGGTGGGTGTCTCCACTGTTTATATATTGTTCCGGCCCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGGGGAt                                   <  8:166507/134‑1 (MQ=255)
                                                                                              gCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGGCGGTTTGCATGGGTGGGTGTCTCCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGGGGATg                                  <  3:394109/135‑1 (MQ=255)
                                                                                                ggCACGGGGGCCGGCGGGTTCTCGCGAGTCGTGTCGGTTTCCATGGGGGGGTGTCTCCACAGTTTCTATATTGTGCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGGGGATg                                  <  7:422379/133‑1 (MQ=255)
                                                                                                ggCACCTGGGCCTGCGGGTTCTCGCGTGTCGTGCCGGTTTGAATGGTGGGGTGTCTCCACAGTTTCTATTTTGTTCAGGCGACTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGGGGATgc                                 <  5:328661/134‑1 (MQ=255)
                                                                                                     gtGGGCCGGCGGTTTCGCGCGCGTCGGGGGGGTTCGCCGGGGTGGGGGGCCCCACTGTGTCTATATTGTTCCGGCGCCATATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGGGGATGCGCAatg                           <  6:298622/134‑1 (MQ=255)
                                                                                                                     tcGCGCGGCGGGGCGGTGTGCCTGGGGGGGGGCCTCCACCGTTGCTATATGGTGCCGGCGCCTAAGGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGGGGATGCGCAATGCATGCGCATccgccc            <  1:203472/134‑1 (MQ=255)
                                                                                                                              gggtCGGTTTGCAGGGGGGGGTGTCTCCACAGTTTCAATATTGTCCAGGCCCCTTAGGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGGGGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGt  <  3:72378/133‑1 (MQ=255)
                                                                                                                                 |                                                                                                                                    
TCGTCGGTGTAAGGGACGACAGTGCCGTTACGCTTGATCACGCGCAGTTGGCCGGGGGCGGTGGCAGCCAGATCCTGGTTGGAATCGGCGGCCTGCGGCACCTGGGCCTGCGGGTTCTCGCGAGTCGTGTCGGTTTGCATGGGTGGGTGTCTCCACAGTTTCTATATTGTTCAGGCGCCTTATGGGCGCCCACCGTTCCGTCCACTTGCTCGATAACCTCGCGGGGATGCGCAATGCATGCGCATCCGCCCGCGCGGGCGT  >  ACB122/1348290‑1348550

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: