Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 1,720,951 T→G 35.7% D90A (GAC→GCC)  plsB ← Glycerol‑3‑phosphate acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1221,720,9510TG35.7% 52.0 / 6.3 49D90A (GAC→GCC) plsBGlycerol‑3‑phosphate acyltransferase
Reads supporting (aligned to +/- strand):  ref base T (18/12);  new base G (5/14);  total (23/26)
Fisher's exact test for biased strand distribution p-value = 3.87e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGTCTGGCCCCAGAACACGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCAGTCCGGGTCGGGGGTCAGGTAGAAGAACGCGGCGGGTTCCTGCAGTGGGCCGACGGCTACCGGCAGCACCGGGCGCGGCAGCCCCGCCTTGGTGCACTCGTGGTCGAGCACGGCCAGGTCGGTGAGCGCGGGCGAGGGC  >  ACB122/1720820‑1721082
                                                                                                                                   |                                                                                                                                   
gggTCTGGCCCCAGAACACGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCCGTCCGGGTCggg                                                                                                                                 >  7:326256/1‑136 (MQ=255)
gggTCTGGCCCCAGAACACGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCAGTCCGGGTCgg                                                                                                                                  <  6:112957/135‑1 (MQ=255)
gggTCTGGCCCCAGAACACGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCAGTCCGGGTCgg                                                                                                                                  <  6:107370/135‑1 (MQ=255)
gggTCTGGCCCCAGAACACGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCAGTCCGGGTCgg                                                                                                                                  <  1:163354/135‑1 (MQ=255)
gggTCTGGCCCCAGAACACGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCAGTCCGGGTCgg                                                                                                                                  <  2:129492/135‑1 (MQ=255)
gggTCTGGCCCCAGAACACGCCGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCAGTCCGGGTCgg                                                                                                                                  <  5:74322/135‑1 (MQ=255)
     gggCCCCAGAACACGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCAGTCCGGGTCGGGGGTc                                                                                                                             <  5:27561/134‑1 (MQ=255)
       gCCCCAGAACACGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCAGTCCGGGTCGGGGGTc                                                                                                                             <  8:3411/133‑1 (MQ=255)
         cccAGAACACGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTACTGGCGGCCGAGCCAGTCCGGGTCGGGGGTCAg                                                                                                                           >  6:157349/1‑133 (MQ=255)
           cAGAACACGCTGACAGGAATGATCTGCGCATCTTCTTCGGCAGGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGGGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCAGTCCGGGTCGGGGGTCAGGTag                                                                                                                       <  4:344853/135‑1 (MQ=255)
               acacGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGCCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCAGGCCGGGTCGGGGGTCAGGGagaaga                                                                                                                   >  5:155332/1‑135 (MQ=255)
                   gCTGACAGTATTGAGCGGCGCTTCTTGTGCGGCTGGCGGGCTGACGGCAGCCACCGGGCGCTCTTGCGGGGGGGGGGCGCCGCTTTTGTCGTGGCGGCCGAGCCTGTCGGGGGCGGGGGTCTGTTAGAATAAcgcg                                                                                                              <  8:254322/136‑1 (MQ=255)
                   gCTGACAGGAAGGAGATGCGCAGCTTGTTCGGCAGGCGGCGTGCCGCCAGCCACCAGGCGCTCGAGGGGGGGCGGGGCGCCGCTTTTGGCCTGGCGGCGGAGGCAGTCCGGGTGGGGGGTCAGGTAGAAGAAcgcg                                                                                                              <  8:10278/136‑1 (MQ=255)
                     tGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGCGCCAGTCCGGGGCGGGGGGCAGGGAGAAGAAccccgc                                                                                                            >  7:353962/1‑131 (MQ=255)
                      gACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCAGTACGGGTCGGGGGTCAGGTAGAAGAACgcggc                                                                                                            >  4:255361/1‑135 (MQ=255)
                        cAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCCGTCCGGGTCGGGGGGCAGGTAGAAGAACGcggcgg                                                                                                          >  2:367916/1‑135 (MQ=255)
                          ggAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCAGTCCGGGGCGGGGGTCAGGTAGAAGAACGCGGCGGGt                                                                                                        >  6:33736/1‑135 (MQ=255)
                                    gcgcATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCAGTCCGGGTCGGGGGTCCGGTAGAAGAACGCGGCGGGTTCCTGCAGTg                                                                                              >  3:71242/1‑135 (MQ=255)
                                     cgcATCTTCTTCGGCATGCTGGCTGCCGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCGGCTTTTGTCCTGGCGGCCGAGCCAGTCCGGGGCGGGGGTCAGGTAGAAGAACGCGGCGGGTTCCTGCAGTgg                                                                                             <  6:64940/135‑1 (MQ=255)
                                       cATCTTCTTCGGCTTGCTGGCGGACGGCAGCCACCCGGCGCTCGAGCGGGGGCGGGGGGCCGCTTTTGACCTGGCGGCCGAGCCAGTCCGGGGCGGGGGTCAGGTAGAAGAACGCGGCGGGTTCCTGCAGTGGGc                                                                                           <  2:46780/135‑1 (MQ=255)
                                       cATCTTCTGCGGCCGGCGGGCTGCCGGCAGCCCCCAGGCGCGCGCGCGGGGGGGGGGCGCCGCTTTGGCCCGGGCGGCCGAGCCAGGCCGGGGCGGGGGGCCGGAAGAAGACCGCGGCGGGTTCCTGCAGTGGGc                                                                                           <  2:133697/135‑1 (MQ=255)
                                        aTCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCAGTCCGGGTCGGGGGTCAGGTAGAAGAACGCGGCGGGTTCCTGCAGTGGGcc                                                                                          >  6:161499/1‑135 (MQ=255)
                                        aTCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCAGTCCGGGTCGGGGGTCAGGGAGAAGAACGCGGCGGGTTCCTGCAGTGGGcc                                                                                          >  1:360524/1‑135 (MQ=255)
                                             cttcGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCAGTCCGGGTCGGGGGTCAGGTAGAAGAACGCGGCGGGTTCCTGCAGTGGGCCGACgg                                                                                     >  3:313535/1‑135 (MQ=255)
                                                         ggctgACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCGTTTGTCCTGGCGGCCGCGCCAGTCCGGGGGGGGGGGCAGGTTGAAGAACGCGGCGGGTTCCTGCAGTGGGGCGCCGGCGTCCGgcagcc                                                                         >  6:364187/1‑134 (MQ=255)
                                                           ctgACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCGTTTGTCCTGGCGGCCGAGCCAGTCCGGGGCGGGGGTCAGGGAGAAGAACGCGGCGGGTTCCTGCAGTGGGCCGCCGGCTACCGGCAGCAcc                                                                       >  2:245193/1‑135 (MQ=255)
                                                               cGGCACCCCCCTGGCGCGCGCGGGGGGGCGGGGCGCCGCTTTTGCCCGGGCGCCCGCGCCCGGCCGGGGCGGGGGCCAGGTCCAAGAACGCGGCGGGTTCCTGCAGTGGGCCGACGGCTACCGGCAGCACCGGgcg                                                                  <  8:353962/136‑1 (MQ=255)
                                                                      cacccGGCGCTCGCGCGGGGGCGGGGCGCCGCTTTTGTCCGGGCGGCCGAGCCCGGCCGGGGCGGGGGGCAGGTCGAAGACCGCGGCGGGTTCCTGCAGCGGGCCGACGGCTACCGGCAGCACCGGGCGCGGCAg                                                            <  1:304818/135‑1 (MQ=255)
                                                                          aGGCGCGCGCGCGGGGGCGGGGCGCCGCTTTTGTCCTGGCGCCCGAGCCAGGCCGGGGCGGGGGCCAGGTAGAAGACCGCGGCGGGTTCCTGCAGTGGGCCGACGGCTACCGGCAGCACCGGGCGCGGCAGCCCCg                                                       <  8:6825/136‑1 (MQ=255)
                                                                           ggCGGTCGCGGGGGGGGGGGGCGCCGGTTTGGCCCGGGCGGCCGCGCCTGGCGGGGGGGGGGGCCGGGACCACCACCGCGGCGGGTTCCTGCAGTGGGCCGACGGCTACCGGCAGCACCGGGCGCGGCAGCCCCg                                                       <  5:268394/135‑1 (MQ=255)
                                                                              gcTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGCGCCAGTCCGGGTCGGGGGTCAGGTAGAAGAACGCGGCGGGTTACTGCAGTGGGCCGACGGCTACCGGCCGCACCGGGCGCGGccgccccgccg                                                    >  3:11805/1‑134 (MQ=255)
                                                                              gcTCGAGCGTGGGCGGGGCGCCGCGTTTGTCCTGGCGGCCGCGCCCGTCCGGGGCGGGGGGCAGGTAGCAGACCGCGGCGGGTTGCGGCAGTGGGCCGACGGCGACCGGCCGCCCCGGGCGCGGCCGCCCCCCCt                                                    >  3:16685/1‑135 (MQ=255)
                                                                                  gAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCAGTCCGGGTCGGGGGTCAGGTAGAAGAACGCGGCGGGTTCCTGCAGTGGGCCGACGGCTACCGGCAGCACCGGGCGCGGCAGCCCCGCCTTGGTg                                               >  7:309301/1‑136 (MQ=255)
                                                                                      gTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCAGTCCGGGTCGGGGGTCAGGTAGAAGAACGCGGCGGGTTCCTGCAGTGGGCCGACGGCTACCGGCAGCACCGGGCGCGGCAGCCCCGCCTTGGTGCACt                                           >  8:85894/1‑136 (MQ=255)
                                                                                              ggcgCCGCTTTGGTCCGGGCGGCCGAGCCCGTCGGGGGCGGGGGTCAGGAACAACACCGCGGCGGGTTCCTGCAGTGGGCCGACGGCTACCGGCAGCACCGGGCGCGGCAGCCCCGCCTTGGTGCACTCGTGGTc                                    <  2:80814/135‑1 (MQ=255)
                                                                                                 gccgcTTTTGTCCTGGCGGCCGAGCCAGTCCGGGTCGGGGGTCAGGTAGAAGAACGCGGCGGGTTCCTGCAGTGGGCCGACGGCTACCGGCAGCACCGGGCGCGGCCGCCCCGCCTTGGTGCACTCGTGGTCGa                                  >  1:261444/1‑134 (MQ=255)
                                                                                                     cTTTTGTCCTGGCGGCCGAGCCTGGCCGGGGCGGGGGGCCGGTAGAAGAACGCGGCGGGTTCCTGCAGTGGGCCGACGGCTACCGGCAGCACCGGGCGCGGCAGCCCCGCCTTGGTGCACTCGTGGTCGAGCACg                             <  2:227660/135‑1 (MQ=255)
                                                                                                             cgggcggcCGCGCCAGTCGGGGGCGGGGTGCAGTTAGAAGACCGCGGCGGGGTCCTGCAGGGGGCCGACGGCTACCGGCAGCACCGGGCGCGGCAGCCCCGCCTTGGTGCACTCGTGGGCGAGCACGGCCAGGTc                     <  2:263021/133‑1 (MQ=255)
                                                                                                             agggcggcCGCGCCCGGCCGGGGGGGGTGGCAGGTAGAACAACGCGGCGGGGTCCTGCAGTGGGCCGACGGCTACCGGCAGCACCGGGCGCGGCAGCCCCGCCTTGGTGCACTCGTGGTCGAGCACGGCCAGGTc                     <  2:268752/133‑1 (MQ=255)
                                                                                                                 cggcCGAGCCAGTCCGGGTCGGGGGTCAGGTAGAAGAACGCGGCGGGTTCCTGCAGTGGGCCGACGGCTACCGGCAGCACCGGGCGCGGCAGCCCCGCCTTGGTGCACTCGTGGTCGAGCACGGCCAGGTCGGTGa                <  8:309301/136‑1 (MQ=255)
                                                                                                                  ggcCGAGCCAGTCCGGGTCGGGGGTCAGGTAGAAGAACGCGGCGGGTTCCTGCAGTGGGCCGACGGCTACCGGCAGCACCGGGCGCGGCAGCCCCGCCTTGGTGCACTCGTGGTCGAGCACGGCCAGGTCGGTGa                <  4:63149/135‑1 (MQ=255)
                                                                                                                   gcCGAGCCAGTCCGGGTCGGGGGTCAGGTAGAAGAACGCGGCGGGTTCCTGCAGTGGGCCGACGGCTACCGGCAGCACCGGGCGCGGCAGCCCCGCCTTGGTGCACTCGTGGTCGAGCACGGCCAGGTCGGTGAgc              >  8:346871/1‑136 (MQ=255)
                                                                                                                   gcCGAGCCAGTCCGGGTCGGGGGTCAGGTAGAAGAACGCGGCGGGTTCCTGCAGTGGGCCGACGGCTACCGGCAGCACCGGGCGCGGCAGCCCCGCCTTGGTGCACTCCTGGTCGCGCACGGCCAGGTCgggtcgc              >  7:297018/1‑131 (MQ=255)
                                                                                                                   gcCGAGCCAGTCCGGGTCGGGGGTCAGGTAGAAGAACGCGGCGGGTTCCTGCAGTGGGCAGACGGCTACCGGCAGCACCGGGCGCGGCAGCCCCGCCTTGGTGCACTCGTGGTCGAGCACGGCCAGGTCGGTGAg               >  1:220616/1‑135 (MQ=255)
                                                                                                                        gCCAGTCCGGGTCGGGGGTCAGGTAGAAGAACGCGGCGGGTTCCTGCAGTGGGCCGACGGCTACCGGCAGCACCGGGCGCGGCAGCCCCGCCTTGGTGCACTCGTGGTCGAGCACGGCCAGGTCGGTGAGCGCgg          >  2:114421/1‑135 (MQ=255)
                                                                                                                             gccgggccgggGGCCTGTAAGACGCCCGCGGCGGGTGCCTGCAGTGGGCCGCCGGCAACCGGCAGCACCGGGCGCGGCAGCCCCGCCTGGGTGCACTCGTGGTCGAGCACGGCCAGGTCGGTGAGCGCGGGCGAgg    <  7:253632/135‑1 (MQ=255)
                                                                                                                              cgggggCGGGGGGCAGTTCGAACACCGCGGCGGGTTCCTGCAGGGGGCCGACGGCTACCGGCAGCACCGGGCGCGGCAGCCCCGCCTTGGTGCACTCGTGGTCGAGCACGGCCAGGTCGGTGAGCGCGGGCGAgg    <  1:75309/133‑1 (MQ=255)
                                                                                                                              ccGGGGGGGGGGTCTGGTAGACGACCGCGGCGGGTTCCTGCAGTGGGCCGACGGCTACCGGCAGCACCGGGCGCGGCAGCCCCGCCTTGGTGCACTCGTGGGCGAGCACGGCCAGGTCGGTGAGCGCGGGCGAgg    <  1:210774/135‑1 (MQ=255)
                                                                                                                              ccGGGGGGGGGGTCGGTTCGAACCCCGCGGGGGGTGCCTGCTGTGGGCCGACGGCGACCGGCAGCACCGGGCGCGGCAGCCCCGCCTTGGTGCACGCGTGGTCGAGCACGTCCAGGTCGGTGAGCGCGGGCGAgg    <  6:390352/135‑1 (MQ=255)
                                                                                                                               cGGGTGGGGGGCCCGGTACAAGACCGCGGCGGGTTCCTGCAGTGGGCCGACGGCTACCGGCAGCACCGGGCGCGGCAGCCCCGCCTTGGTGCACTCGTGGTCGAGCACGGCCAGGTCGGTGAGCGCGGGCGAgg    <  1:67950/134‑1 (MQ=255)
                                                                                                                                gggTCGGGGGTCAGGTAGAAGAACGCGGCGGGTTCCTGCAGTGGGCCGACGGCTACCGGCAGCACCGGGCGCGGCAGCCCCGCCTTGGTGCACTCGTGGTCGAGCACGGCCAGGTCGGTGAGCGCGGGCGAGGGc  >  2:273156/1‑135 (MQ=255)
                                                                                                                                   |                                                                                                                                   
GGGTCTGGCCCCAGAACACGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCAGTCCGGGTCGGGGGTCAGGTAGAAGAACGCGGCGGGTTCCTGCAGTGGGCCGACGGCTACCGGCAGCACCGGGCGCGGCAGCCCCGCCTTGGTGCACTCGTGGTCGAGCACGGCCAGGTCGGTGAGCGCGGGCGAGGGC  >  ACB122/1720820‑1721082

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