Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB143 4,640,631 A→G 44.5% V30A (GTG→GCG)  bicA ← Bicarbonate transporter BicA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1434,640,6310AG44.5% 9.7 / 75.0 54V30A (GTG→GCG) bicABicarbonate transporter BicA
Reads supporting (aligned to +/- strand):  ref base A (15/15);  new base G (10/14);  total (25/29)
Fisher's exact test for biased strand distribution p-value = 5.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.99e-01

AAGGCGACGATCACCACCGCCATCGAGCCGGCAGCGCCGGAGATCATCCCCGGCCGCCCGCCGAACAGCGCCGTCAGGGTACAAATGATGAACGCGCCATACAGGCCCATCAGCGGGTTGAGGTGGGCCACCAGGGCAAAGGCGATGCACTCGGGGACCAAGGCGAACGAGGTGGTGAGGCCGGCGAGCAGGTCGGCGCGTAGTCGGGCGGGCTTCATGGGCTTCCTGTGGTGCAGTCCTGACGGACGTTAGAAAAACGAGAGGC  >  ACB143/4640502‑4640766
                                                                                                                                 |                                                                                                                                       
aaGGCGACGATCACCACCGCCATCGAGCCGGCAGCGCCGGAGATCATCCCCGGCCGCCCGCCGAACAGCGCCGTCAGGGTACAAATGATGAACGCGCCATACAGGCCCATCAGCGGGTTGAGGTGGGCCACCAggg                                                                                                                                   <  5:197332/136‑1 (MQ=255)
 aGGCGACGATCACCACCGCCATCGAGCCGGCAGCGCCGGAGATCATCCCCGGCCGCCCGCCGAACAGCGCCGTCAGGGTACAAATGATGAACGCGCCATACAGGCCCATCAGCGGGTTGAGGTGGGCCGCCAggg                                                                                                                                   >  2:171067/1‑135 (MQ=255)
 aGGCGACGATCACCACCGCCATCGAGCCGGCAGCGCCGGAGATCATCCCCGGCCGCCCGCCGAACAGCGCCGTCAGGGTACAAATGATGAACGCGCCATACAGGCCCATCAGCGGGTTGAGGTGGGCCGCCAggg                                                                                                                                   >  3:119646/1‑135 (MQ=255)
 aGGCGACGATCACCACCGCCATCGAGCCGGCAGCGCCGGAGATCATCCCCGGCCGCCCGCCGAACAGCGCCGTCAGGGTACAAATGATGAACGCGCCATACAGGCCCATCAGCGGGTTGAGGTGGGCCACCAggg                                                                                                                                   <  8:120761/135‑1 (MQ=255)
 aGGCGACGATCACCACCGCCATCGAGCCGGCAGCGCCGGAGATCATCCCCGGCCGCCCGCCGAACAGCGCCGTCAGGGTACAAATGATGAACGCGCCATACAGGCCCATCAGCGGGTTGAGGTGGGCCACCAggg                                                                                                                                   >  1:287067/1‑135 (MQ=255)
 aGGCGACGATCACCACCGCCATCGAGCCGGCAGCGCCGGAGATCATCCCCGGCCGCCCGCCGAACAGCGCCGTCAGGGTACAAATGATGAACGCGCCATACAGGCCCATCAGCGGGTTGAGGTGGGCCACCAggg                                                                                                                                   >  7:120761/1‑135 (MQ=255)
 aGGCGACGATCACCACCGCCATCGAGCCGGCAGCGCCGGAGATCATCCCCGGCCGCCCGCCGAACAGCGCCGTCAGGGTACAAATGATGAACGCGCCATACAGGCCCATCAGCGGGTTGAGGTGGGCCACCAggg                                                                                                                                   >  2:186653/1‑135 (MQ=255)
     gacgaTCACCACCGCCATCGAGCCGGCAGCGCCGGAGATCATCCCCGGCCGCCCGCCGAACAGCGCCGTCAGGGTACAAATGATGAACGCGCCATACAGGCCCATCAGCGGGTTGAGGTGGGCCACCAGGGCAAAg                                                                                                                              <  7:260510/136‑1 (MQ=255)
     gacgaTCACCACCGCCATCGAGCCGGCAGCGCCGGAGATCATCCCCGGCCGCCCGCCGAACAGCGCCGTCAGGGTACAAATGATGAACGCGCCATACAGGCCCATCAGCGGGTTGAGGTGGGCCACCAGGGCAAAg                                                                                                                              <  3:253932/136‑1 (MQ=255)
         aTCACCACCGCCATCGAGCCGGCAGCGCCGGAGATCATCCCCGGCCGCCCGCCGAACAGCGCCGTCAGGGTACAAATGATGAACGCGCCATACAGGCCCATCAGCGGGTTGAGGTGGGCCACCAGGGCAAAGGCGa                                                                                                                          >  4:190193/1‑136 (MQ=255)
                           ccGGCAGCGCCGGAGATCATCCCCGGCCGCCCGCCGAACAGCGCCGTCAGGGTACAAATGATGAACGCGCCATACAGGCCCATCAGCGGGTTGAGGTGGGCCACCAGGGCAAAGGCGATGCACTCGGGGACCAAgg                                                                                                        <  8:11649/136‑1 (MQ=255)
                              gCAGCGCCGGAGATCATCCCCGGCCGCCCGCCGAACAGCGCCGTCAGGGTACAAATGATGAACGCGCCATACAGGCCCCTCAGCGGGTTGAGGTGGGCCGCCAGGGCAACGGCGATGCACTCGGGGACCAAGGCGa                                                                                                     >  2:247347/1‑136 (MQ=255)
                                             atcCCCGGCCGCCCGCCGAACAGCGCCGTCAGGGTACAAATGATGAACGCGCCATACAGGCCCATCAGCGGGTTGAGGTGGGCCGCCAGGGCAAAGGCGATGCACTCGGGGACCAAGGCGAACGAGGTGGTGAgg                                                                                       <  1:309551/135‑1 (MQ=255)
                                             atcCCCGGCCGCCCGCCGAACAGCGCCGTCAGGGTACAAATGATGAACGCGCCATACAGGCCCATCAGCGGGTTGAGGTGGGCCACCAGGGCAAAGGCGATGCACTCGGGGACCAAGGCGAACGAGGTGGTGAgg                                                                                       <  3:316841/135‑1 (MQ=255)
                                                 ccGGCCGCCCGCCGAACAGCGCCGTCAGGGTACAAATGATGAACGCGCCATACAGGCCCATCAGCGGGTTGAGGTGGGCCACCAGGGCAAAGGCGATGCACTCGGGGACCAAGGCGAACGAGGTGGTGAGGCCGGc                                                                                  <  1:215964/136‑1 (MQ=255)
                                                        cccgccGAACAGCGCCGTCAGGGTACAAATGATGAACGCGCCATACAGGCCCATCAGCGGGTTGAGGTGGGCCGCCCGGGCAAAGGCGATGCACTCGGGGACCAAGGCGAACGAGGTGGTGAGGCCGGCGAGCAgg                                                                           <  7:277599/136‑1 (MQ=255)
                                                            ccgAACAGCGCCGTCAGGGTACAAATGATGAACGCGCCATACAGGCCCATCAGCGGGTTGAGGTGGGCCGCCAGGGCAAAGGCGATGCACTCGGGGACCAAGGCGAACGAGGTGGTGAGGCCGGCGAGCAGGTCgg                                                                       >  5:89743/1‑136 (MQ=255)
                                                            ccgAACAGCGCCGTCAGGGTACAAATGATGAACGCGCCATACAGGCCCATCAGCGGGTTGAGGTGGGCCACCAGGGCAAAGGCGATGCACTCGGGGACCAAGGCGAACGAGGTGGTGAGGCCGGCAAGCAGGTCgg                                                                       >  3:210639/1‑136 (MQ=255)
                                                                    cgcCGTCAGGGTACAAATGATGAACGCGCCATACAGGCCCATCAGCGGGTTGAGGGGGGCCGCCAGGGCAAAGGCGATGCACTCGGGGACCAAGGCGAACGAGGTGGTGAGGCCGGCGAGCAGGTCGGCGCGTAGt                                                               <  8:267181/136‑1 (MQ=255)
                                                                     gcCGTCAGGGTACAAATGATGAACGCGCCATACAGGCCCATCAGCGGGTTGAGGTGGGCCACCAGGGCAAAGGCGATGCACTCGGGGACCAAGGCGAACGAGGTGGTGAGGCCGGCGAGCAGGTCGGCGCGTAGTc                                                              <  3:269812/136‑1 (MQ=255)
                                                                     gcCGTCAGGGTACAAATGATGAACGCGCCATACAGGCCCATCAGCGGGTTGAGGTGGGCCACCAGGGCAAAGGCGATGCACTCGGGGACCAAGGCGAACGAGGTGGTGAGGCCGGCGAGCAGGTCGGCGCGTAGTc                                                              <  2:209844/136‑1 (MQ=255)
                                                                       cGTCAGGGTACAAATGATGAACGCGCCATACAGGCCCATCAGCGGGTTGAGGTGGGCCGCCAGGGCAAAGGCGATGCACTCGGGGACCAAGGCGAACGAGGTGGTGAGGCCGGCGAGCAGGTCGGCGCGTAGTcgg                                                            <  8:219229/136‑1 (MQ=255)
                                                                       cGTCAGGGTACAAATGATGAACGCGCCATACAGGCCCATCAGCGGGTTGAGGTGGGCCACCAGGGCAAAGGCGATGCACTCGGGGACCAAGGCGAACGAGGTTGTGAGGCCGGCGAGCAGGTCGGCGCGTAGTcgg                                                            <  6:101692/136‑1 (MQ=255)
                                                                              gTACAAATGATGAACGCGCCATACAGGCCCATCAGCGGGTTGAGGTGGGCCGCCAGGGCAAAGGCGATGCACTCGGGGACCAAGGCGAACGAGGTGGTGAGGCCGGCGAGCAGGTCGGCGCGTAGTCGGGCGGGCt                                                     <  5:159549/136‑1 (MQ=255)
                                                                                 cAAATGATGAACGCGCCATACAGGCCCATCAGCGGGTTGAGGTGGGCCGCCAGGGCAAAGGCGATGCACTCGGGGACCAAGGCGAACGAGGTGGTGAGGCCGGCGAGCAGGTCGGCGCGTAGTCGGGCGGGCTTCa                                                  <  3:266617/136‑1 (MQ=255)
                                                                                    atgatgAACGCGCCATACAGGCCCATCAGCGGGTTGAGGTGGGCCACCAGGGCAAAGGCGATCCACTCGGGGACCAAGGCGAACGAGGTGGTGAGGCCGGCGAGCAGGTCGGCGCGTAGACGGGCGGGCTTCATgg                                               >  5:50791/1‑136 (MQ=255)
                                                                                        tgaACGCGCCATACAGGCCCATCAGCGGGTTGAGGTGGGCCGCCAGGGCAAAGGCGATGCACTCGGGGACCAAGGCGAACGAGGTGGTGAGGCCGGCGAGCAGGTCGGCGCGTAGTCGGGCGGGCTTCATGGGCtt                                           <  6:151937/136‑1 (MQ=255)
                                                                                          aaCGCGCCATACAGGCCCATCAGCGGGTTGAGGTGGGCCGCCAGGGCAAAGGCGATGCACTCGGGGACCAAGGCGAACGAGGTGGTGAGGCCGGCGAGCAGGTCGGCGCGTAGTCGGGCGGGCTTCATGGGCTTcc                                         >  6:12358/1‑136 (MQ=255)
                                                                                          aaCGCGCCATACAGGCCCATCAGCGGGTTGAGGTGGGCCGCCAGGGCAAAGGCGATGCACTCGGGGACCAAGGCGAACGAGGTGGTGAGGCCGGCGAGCAGGTCGGCGCGTAGTCGGGCGGGCTTCATGGGCTTc                                          <  1:259831/135‑1 (MQ=255)
                                                                                          aaCGCGCCATACAGGCCCATCAGCGGGTTGAGGTGGGCCACCAGGGCAAAGGGGATGCACTCGGGGACCAAGGCGAACGAGGTGGTGAGGCCGGCGAGCAGGTCGGCGCGTAGTCGGGCGGGCTTCATGGGCTTcc                                         >  5:63053/1‑136 (MQ=255)
                                                                                          aaCGCGCCATACAGGCCCATCAGCGGGTTGAGGTGGGCCACCAGGGCAAAGGCGATGCACTCGGGGACCAAGGCGAACGAGGTGGTGAGGCCGGCGAGCAGGTCGGCGCGTAGTCGGGCGGGCTTCATGGGCTTcc                                         >  5:209702/1‑136 (MQ=255)
                                                                                            cgcgCCATACAGGCCCATCAGCGGGTTGAGGTGGGCCGCCAGGGCAAAGGCGATGCACTCGGGGACCAAGGCGAACGAGGTGGTGAGGCCGGCGAGCAGGTCGGCGCGTAGTCGGGCGGGCTTCATGGGCTTCCtg                                       >  1:65948/1‑136 (MQ=255)
                                                                                              cgcCATACAGGCCCATCAGCGGGTTGAGGTGGGCCGCCAGGTCAAAGGCGATGCACTCGGGGACCAAGGCGAACGAGGTGGTGAGGCCGGCGAGCAGGTCGGCGCGTAGTCGGGCGGGCTTCATGGGCTTCCtgtg                                     <  7:152506/136‑1 (MQ=255)
                                                                                              cgcCATACAGGCCCATCAGCGGGTTGAGGTGGGCCGCCAGGGCAAAGGCGATGCACTCGGGGACCAAGGCGAACGAGGTGGTGAGGCCGGCGAGCAGGTCGGCGCGTAGTCGGGCGGGCTTCATGGGCTTCCtgtg                                     <  4:205481/136‑1 (MQ=255)
                                                                                              cgcCATACAGGCCCATCAGCGGGTTGAGGTGGGCCGCCAGGGCAAAGGCGATGCACTCGGGGACCAAGGCGAACGAGGTGGTGAGGCCGGCGAGCAGGTCGGCGCGTAGTCGGGCGGGCTTCATGGGCTTCCtgtg                                     <  7:215056/136‑1 (MQ=255)
                                                                                              cgcCATACAGGCCCATCAGCGGGTTGAGGTGGGCCGCCAGGGCAAAGGCGATGCACTCGGGGACCAAGGCGAACGAGGTGGTGAGGCCGGCGAGCAGGTCGGCGCGTAGTCGGGCGGGCTTCATGGGCTTCCtgtg                                     <  7:201426/136‑1 (MQ=255)
                                                                                              cgcCATACAGGCCCATCAGCGGGTTGAGGTGGGCCACCAGGGCAAAGGCGATGCACTCGGGGACCAAGGCGAACGAGGTGGTGAGGCCGGCGAGCAGGTCGGCGCGTAGTCGGGCGGGCTTCATGGGCTTCCtgtg                                     <  1:187036/136‑1 (MQ=255)
                                                                                              cgcCATACAGGCCCATCAGCGGGTTGAGGTGGGCCACCAGGGCAAAGGCGATGCACTCGGGGACCAAGGCGAACGAGGTGGTGAGGCCGGCGAGCAGGTCGGCGCGTAGTCGGGCGGGCTTCATGGGCTTCCtgtg                                     <  3:2939/136‑1 (MQ=255)
                                                                                              cgcCATACAGGCCCATCAGCGGGTTGAGGGGGGCCACCAGGGCAAAGGCGATGCACTCGGGGACCAAGGCGAACGAGGTGGTGAGGCCGGCGAGCAGGTCGGCGCGTAGTCGGGCGGGCTTCATGGGCTTCCtgtg                                     <  8:10973/136‑1 (MQ=255)
                                                                                              cgcCATACAGGCCCATCAGCGGGTTGAGGGGGGCCACCAGGGCAAAGGCGATGCACTCGGGGACCAAGGCGAACGAGGTGGTGAGGCCGGCGAGCAGGTCGGCGCGTAGTCGGGCGGGCTTCATGGGCTTCCtgtg                                     <  6:141241/136‑1 (MQ=255)
                                                                                                ccATACAGGCCCATCAGCGGGTTGAGGTGGGCCGCCAGGGAAAAGGCGATGCACTCGGGGACCAAGGCGAACGAGGTGGTGAGGCCGGCGAGCAGGTCGGCGCGTAGTCGGGCGGGCTTCATGGGCTTCCTgtggt                                   <  8:202234/136‑1 (MQ=255)
                                                                                                              cAGCGGGTTGAGGTGGGCCGCCAGGGCAAAGGCGATGCACTCGGGGACCAAGGCGAACGAGGTGGTGAGGCCGGCGAGCAGGTCGGCGCGTAGTCGGGCGGGCt                                                     <  6:82422/104‑1 (MQ=255)
                                                                                                              cAGCGGGTTGAGGTGGGCCGCCAGGGCAAAGGCGATGCACTCGGGGACCAAGGCGAACGAGGTGGTGAGGCCGGCGAGCAGGTCGGCGCGTAGTCGGGCGGGCt                                                     >  5:82422/1‑104 (MQ=255)
                                                                                                              cAGCGGGTTGAGGTGGGCCGCCAGGGCAAAGGCGATGCACTCGGGGACCAAGGCGAACGAGGTGGTGAGGCCGGCGAGCAGGTCGGCGCGTAGTCGGGCGGGCTTCATGGGCTTCCTGTGGTGCAGTACTgacgga                     >  8:171002/1‑136 (MQ=255)
                                                                                                              cAGCGGGTTGAGGTGGGCCACCAGGGCAAAGGCGATGCACTCGGGGACCAAGGCGAACGAGGTGGTGAGGCCGGCGAGCAGGTCGGCGCGTAGTcgggcg                                                         <  3:262670/100‑1 (MQ=255)
                                                                                                              cAGCGGGTTGAGGTGGGCCACCAGGGCAAAGGCGATGCACTCGGGGACCAAGGCGAACGAGGTGGTGAGGCCGGCGAGCAGGTCGGCGCGTAGTcgggcg                                                         >  4:262670/1‑100 (MQ=255)
                                                                                                              cAGCGGGTTGAGGTGGGCCACCAGGGCAAAGGCGATGCACTCGGGGACCAAGGCGAACGAGGTGGTGAGGCCGGCGAGCAGGTCGGCGCGTAGTCGGGCGGGCTTCATGGGCTTCCTGTGGTGCAGTCCTgacgga                     >  1:86905/1‑136 (MQ=255)
                                                                                                                gagggTTGAGGTGGGCCAAAATTTCAAAGTCTAACAACTCGTGGACCAAGGCGAACGAGGTGGTGAGGCCGGCGAGCAGGTCGGCGCGTAGTCGGGCGGGCTTCGTGGGCTTCCTGTGGTGCAGTCCTgacggacg                   >  8:240361/3‑136 (MQ=255)
                                                                                                                gCGGGTTGAGGTGGGCCGCCAGGGCAAAGGCGATGCACTCGGGGACCAAGGCGAACGAGGTGGTGAGGCCGGCGAGCAGGTCGGCGCGTAGTCGGGCGGGCTTCATGGGCTTCCTGTGGTGCAGTCCTgacggacg                   >  3:98110/1‑136 (MQ=255)
                                                                                                                gCGGGTTGAGGTGGGCCGCCAGGGCAAAGGCGATGCACTCGGGGACCAAGGCGAACGAGGTGGTGAGGCCGGCGAGCAGGTCGGCGCGTAGTCGGGCGGGCTTCATGGGCTTCCTGTGGTGCAGTCCTgacggacg                   >  4:173595/1‑136 (MQ=255)
                                                                                                                gCGGGTTGAGGTGGGCCACCAGGGCAAAGGCGATGCACTCGGGGACCAAGGCGAACGAGGTGGTGAGGCCGGCGAGCAGGTCGGCGCGTAGTCGGGCGGGCTTCATGGGCTTCCTGTGGTGCAGTCCTgacggacg                   >  8:276893/1‑136 (MQ=255)
                                                                                                                gCGGGTTGAGGTGGGCCACCAGGGCAAAGGCGATGCACTCGGGGACCAAGGCGAACGAGGTGGTGAGGCCGGCGAGCAGGTCGGCGCGTAGTCGGGCGGGCTTCATGGGCTTCCTGTGGTGCAGTCCTgacggacg                   >  4:74993/1‑136 (MQ=255)
                                                                                                                gCGGGTTGAGGTGGGCCACCAGGGCAAAGGCGATGCACTCGGGGACCAAGGCGAACGAGGTGGTGAGGCCGGCGAGCAGGTCGGCGCGTAGTCGGGCGGGCTTCATGGGCTTCCTGTGGTGCAGTCCTgacggacg                   >  6:142670/1‑136 (MQ=255)
                                                                                                                       gAGGTGGGCCACCAGGGCAAAGGCGATGCACTCGGGGACCAAGGCGAACGAGGTGGTGAGGCCGGCGAGCAGGTCGGCGCGTAGTCGGGCGGGCTTCATGGGCTTCCTGTGGTGCAGTCCTGACGGACGTTAGaa             >  3:151661/1‑135 (MQ=255)
                                                                                                                               ccaccaGTGCAAAGGCTATGCACTCGGGGACCAAGGCGAACGAGGTGGTGAGGCCGGCGAGCAGGTCGGCGCGTAGTCGGGCGGGCTTCATGGGCTTCCTGTGGTGCAGTCCTGACGGACGTTAGAAAAACGagag    <  6:239463/136‑1 (MQ=255)
                                                                                                                               ccaccaGGGCAAAGGCGATGCACTCGGGGACCAAGGCGAACGAGGTGGTGAGGCCGGCGAGCAGGTCGGCGCGTAGTCGGGCGGGCTTCATGGGCTTCCTGTGGTGCAGTCCTGACGGACGTTAGAAAAACGagag    <  4:102132/136‑1 (MQ=255)
                                                                                                                                 accaGGGCAAAGGCGATGCACTCGGGGACCAAGGCGAACGAGGTGGTGAGGCCGGCGAGCAGGTCGGCGCGTAGTCGGGCGGGCTTCATGGGCTTCCTGTGGTGCAGTCCTGACGGACGTTAGAAAAACGAGAGgc  >  6:92664/1‑136 (MQ=255)
                                                                                                                                 |                                                                                                                                       
AAGGCGACGATCACCACCGCCATCGAGCCGGCAGCGCCGGAGATCATCCCCGGCCGCCCGCCGAACAGCGCCGTCAGGGTACAAATGATGAACGCGCCATACAGGCCCATCAGCGGGTTGAGGTGGGCCACCAGGGCAAAGGCGATGCACTCGGGGACCAAGGCGAACGAGGTGGTGAGGCCGGCGAGCAGGTCGGCGCGTAGTCGGGCGGGCTTCATGGGCTTCCTGTGGTGCAGTCCTGACGGACGTTAGAAAAACGAGAGGC  >  ACB143/4640502‑4640766

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