Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 385,804 G→A 34.1% V338V (GTG→GTA cstA → carbon starvation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE385,8040GA34.1% 106.1 / 86.9 105V338V (GTG→GTAcstAcarbon starvation protein
Reads supporting (aligned to +/- strand):  ref base G (33/33);  new base A (24/11);  total (61/44)
Fisher's exact test for biased strand distribution p-value = 9.28e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.10e-01

Only 100 of 110 total aligned reads displayed.

TGGCACTGGCCCGGTATGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGC  >  minE/385735‑385873
                                                                     |                                                                     
tGGCACTGGCCCGGTATGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGgtgt                                                                      <  1:381716/71‑1 (MQ=255)
tGGCACTGGCCCGGTATGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGgtg                                                                       >  1:757379/1‑70 (MQ=255)
 ggCACTGGCCCGGTATGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTc                                                                     >  1:1851656/1‑71 (MQ=255)
 ggCACTGGCCCGGTATGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTc                                                                     >  1:3297794/1‑71 (MQ=255)
  gCACTGGCCCGGTATGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTc                                                                     <  1:3043332/70‑1 (MQ=255)
  gCACTGGCCCGGTATGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCt                                                                    <  1:208726/71‑1 (MQ=255)
  gCACTGGCCCGGTATGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCt                                                                    >  1:1195516/1‑71 (MQ=255)
  gCACTGGCCCGGTATGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCt                                                                    <  1:3771173/71‑1 (MQ=255)
  gCACTGGCCCGGTATGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCt                                                                    <  1:2961331/71‑1 (MQ=255)
   cACTGGCCCGGTATGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTg                                                                   >  1:1514140/1‑71 (MQ=255)
     cTGGCCCGGTATGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGc                                                                 >  1:3721522/1‑71 (MQ=255)
     cTGGCCCGGTATGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGc                                                                 >  1:1656502/1‑71 (MQ=255)
        gCCCGGTATGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCC‑‑TGGCGCGGTGTCTGGCTTc                                                              <  1:2628831/69‑1 (MQ=255)
                tGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGc                                                                 <  1:5280871/60‑1 (MQ=255)
                  gACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTg                                                     <  1:1045122/70‑1 (MQ=255)
                  gACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTg                                                     <  1:4515735/70‑1 (MQ=255)
                   aCCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCAGTGGCGCGGTGTCTgg                                                                  <  1:3608653/56‑1 (MQ=255)
                     cGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTc                                                                     <  1:4457964/51‑1 (MQ=255)
                        tAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCtctt                                              <  1:3370116/71‑1 (MQ=255)
                        tAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCtctt                                              <  1:4003777/71‑1 (MQ=255)
                         aaCCTGTTCGCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCtcttc                                             >  1:3309006/1‑71 (MQ=255)
                         aaCCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTg                                                     >  1:1059454/1‑63 (MQ=255)
                         aaCCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCtcttc                                             >  1:5359783/1‑71 (MQ=255)
                            cTGTTCCCGTTCCTGTTTTTCACCATCGCCTGTGGCGCGGTGTCTGGCTTc                                                              <  1:754599/51‑1 (MQ=255)
                            cTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTc                                                              <  1:1546709/51‑1 (MQ=255)
                            cTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTc                                                              <  1:705027/51‑1 (MQ=255)
                            cTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTc                                                              <  1:423090/51‑1 (MQ=255)
                            cTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTgg                                          <  1:3770457/71‑1 (MQ=255)
                            cTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTgg                                          <  1:4010256/71‑1 (MQ=255)
                            cTGTTCCCGTTCCTGTTCATCACCATCGCCTGTGGTGCGGTATCTGGCTTCCACGCGCTGAt                                                   >  1:3604457/1‑62 (MQ=255)
                            cTGTTCCCGTTCCTGTTCATCACCATCGCCTGTGGTGCGGTATCTGGCTTCCACGCGCTGATCTCTTCCgg                                          >  1:3630361/1‑71 (MQ=255)
                            cTGTTCCCGTTCCTGTTCATCACCATCGCCTGTGGTGCGGTATCTGGCTTCCACGCGCTGATCTCTTCCgg                                          >  1:344535/1‑71 (MQ=255)
                             tGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTgg                                          >  1:316871/1‑70 (MQ=255)
                             tGTTCCCGTTCCTGTTCATCACCATCGCCTGTGGTGCGGTATCTGGCTTCCACGCGCTGATCtcttccg                                           <  1:379824/69‑3 (MQ=255)
                             tGTTCCCGTTCCTGTTCATCACCATCGCCTGTGGTGCGGTATCTGGCTTCCACGCGCTGATCTCTTCCgg                                          >  1:3299243/1‑70 (MQ=255)
                             tGTTCCCCTTCCTGTTCATCACCATCGCCTGCGGCGCCGTCTCGGGCTTcc                                                             >  1:3721054/1‑51 (MQ=19)
                             tGTTCCCCTTCCTGTTCATCACCATCGCCTGCGGCGCCGTCTCGGGCTTCCa                                                            >  1:204728/1‑52 (MQ=25)
                             tGTTCCCCTTCCTGTTCATCACCATCGCCTGCGGCGCCGTCTCGGGCTTCCa                                                            >  1:179118/1‑52 (MQ=25)
                             tGTTCCCCTTCCTGTTCATCACCATCGCCTGCGGCGCCGTCTCGGGCTTCCa                                                            >  1:4937681/1‑52 (MQ=25)
                               ttCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCa                                                            >  1:1626343/1‑50 (MQ=255)
                               ttCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCa                                                            >  1:2207877/1‑50 (MQ=255)
                               ttCCCGTTCCTGTTCATCACCATCGCCTGTGGTGCGGTATCTGGCTTCCACGCGCTGATCTCTTCCGGTAc                                       >  1:2880985/1‑71 (MQ=255)
                                 cccGTTCCTGTTCATCACCATCGCCTGTGGTGCGGTATCTGGCTTc                                                              >  1:1554944/1‑46 (MQ=255)
                                 cccGTTCATGTTCATCACCATCGCCTGTGGTGCGGTATCTGGCTTc                                                              >  1:1115606/1‑46 (MQ=255)
                                  ccGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCAc                                    <  1:5355415/71‑1 (MQ=255)
                                    gTTCCTGTTCATCACCATCGCCTGTGGTGCGGTATCTGGc                                                                 >  1:2734362/1‑40 (MQ=37)
                                    gTTCCTGTTCATCACCATCGCCTGTGGTGCGGTATCTGGCTTCCACGCGCTGATCTCTTCCGGTACGACGc                                  >  1:201029/1‑71 (MQ=255)
                                    gTTCCTGTTCATCACCATCGCCTGTGGTGCGGTATCTGGCTTCCACGCGCTGATCTCTTCCGGTACGACGc                                  >  1:5588902/1‑71 (MQ=255)
                                    gTTCCTGTTCATCACCATCGCCTGTGGTGCGGTATCTGGCTTCCACGCGCTGATCTCTTCCGGTACGACGc                                  >  1:4703036/1‑71 (MQ=255)
                                     ttCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGcc                                 <  1:168692/71‑1 (MQ=255)
                                     ttCCTGTTCATCACCATCGCCTGTGGTGCGGTATCTGGCTTCCACGCGCTGATCTCTTCCGGTACGACGcc                                 <  1:3677951/71‑1 (MQ=255)
                                      tCCTGTTCATCACCATCGCCTGTGGTGCGGTATCTGGCTTCCACGCGCTGATCtcttc                                             >  1:2926318/1‑58 (MQ=255)
                                       ccTGTTCATCACCATCGCCTGTGGTGCGGTATCTGGCTTCCACGCGCTGa                                                    >  1:3575214/1‑50 (MQ=255)
                                        cTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGAt                                                   >  1:1272292/1‑50 (MQ=255)
                                        cTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGAt                                                   >  1:1213713/1‑50 (MQ=255)
                                        cTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGAt                                                   >  1:577523/1‑50 (MQ=255)
                                        cTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGAt                                                   >  1:4105061/1‑50 (MQ=255)
                                        cTGTTCATCACCATCGCCTGTGGTGCGGTATCTGGCTTCCACGCGCTGATCTCTTCCGGTACGACGccaa                               <  1:345194/70‑3 (MQ=255)
                                        cTGTTCATCACCATCGCCTGTGGTGCGGTATCTGGCTTCCACGCGCTGATCTCTTCCGGTACGACGccaa                               <  1:2874781/70‑3 (MQ=255)
                                         tGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCtcttct                                            <  1:2502230/56‑1 (MQ=255)
                                           tttATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATctct                                               <  1:1289056/51‑1 (MQ=255)
                                           tttATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGa                            <  1:1576820/70‑1 (MQ=255)
                                           tttATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGa                            <  1:5116961/70‑1 (MQ=255)
                                           tttATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGa                            <  1:1661703/70‑1 (MQ=255)
                                               tCACCATCGCCTGTGGTGCGGTATCTGGCTTCCACGCGCTGATCTCTTCCGGTAcga                                     <  1:2420229/57‑3 (MQ=255)
                                                 aCCATCGCCTGTGGTGCGGTATCTGGCTTCCACGCGCTGATCTCTTCCGGTACGACGCCAAAACTGCTgg                      >  1:3700686/1‑70 (MQ=255)
                                                 aCCATCGCCTGTGGTGCGGTATCTGGCTTCCACGCGCTGATCTCTTCCGGTACGACGCCAAAACTGCTGGc                     >  1:4549970/1‑71 (MQ=255)
                                                 aCCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGc                     >  1:2658140/1‑71 (MQ=255)
                                                 aCCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGc                     >  1:5469831/1‑71 (MQ=255)
                                                  ccATCGCCTGTGGTGCGGTATCTGGCTTCCACGCGCTGATCTCTTCCGGTAcga                                     <  1:1188041/54‑3 (MQ=255)
                                                  ccATCGCCTGTGGTGCGGTATCTGGCTTCCACGCGCTGATCTCTTCCGGTACGACGCCAAAACTGCTGGct                    <  1:3375521/71‑2 (MQ=255)
                                                     tCGCCTGTGGTGCGGTATCTGGCTTCCACGCGCTGATCTCTTCCGGTACGACGCCAAAACTGCTgg                      <  1:4625930/66‑1 (MQ=255)
                                                     tCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGcgaac                 >  1:1952477/1‑71 (MQ=255)
                                                     tCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGcgaac                 >  1:5309445/1‑71 (MQ=255)
                                                      cGCCTGTGGTGCGGTATCTGGCTTCCACGCGCTGATCTCTTCCGGTACGACGCCAAAACTGCTGGctaat                 >  1:1707696/1‑69 (MQ=255)
                                                       gCCTGTGGTGCGGTATCTGGCTTCCACGCGCTGATCTCTTCCGGTACGACGCCAAAACTGCTGGctaatg                >  1:5027682/1‑68 (MQ=255)
                                                       gCCTGTGGTGCGGTATCTGGCTTCCACGCGCTGATCTCTTCCGGTACGACGCCAAAACTGCTGGctaatg                <  1:5170441/70‑3 (MQ=255)
                                                       gCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGcgaacg                <  1:2474919/70‑1 (MQ=255)
                                                       gCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGcgaacg                <  1:2238044/70‑1 (MQ=255)
                                                       gCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGcgaacg                <  1:4227084/70‑1 (MQ=255)
                                                       gCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGcgaacg                <  1:3113432/70‑1 (MQ=255)
                                                       gCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGcgaacg                <  1:133727/70‑1 (MQ=255)
                                                       gCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGcg                    >  1:3418768/1‑66 (MQ=255)
                                                         cTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAg             >  1:2125439/1‑71 (MQ=255)
                                                            tgGTGCGGTATCTGGCTTCCACGCGCTGATCTCTTCCGGTACGACGCCAAAACTGCTGGctaatgaaacc           <  1:89464/70‑4 (MQ=255)
                                                            tgGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGa                            >  1:4199975/1‑53 (MQ=255)
                                                            tgGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGa                            >  1:3948685/1‑53 (MQ=255)
                                                            tgGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGa                            >  1:2217885/1‑53 (MQ=255)
                                                            tgGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGa                            >  1:2359027/1‑53 (MQ=255)
                                                            tgGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGa                            >  1:5612503/1‑53 (MQ=255)
                                                             ggCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATg                          >  1:4360710/1‑54 (MQ=255)
                                                             ggCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATg                          >  1:5336243/1‑54 (MQ=255)
                                                             ggCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATg                          >  1:2834309/1‑54 (MQ=255)
                                                             ggCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTAGGGACCACGCCGAAGATg                          >  1:2771463/1‑53 (MQ=255)
                                                              gtgcgGTATCTGGCTTCCACGCGCTGATCTCTTCCGGTACGACGCCAAAACTGCTGGctaa                  >  1:5412112/3‑61 (MQ=255)
                                                              gtgcgGTATCTGGCTTCCACGCGCTGATCTCTTCCGGTACGACGCCAAAACTGCTGGctaa                  >  1:3497235/3‑61 (MQ=255)
                                                              gtgcgGTATCTGGCTTCCACGCGCTGATCTCTTCCGGTACGACGCCAAAACTGCTGGctaa                  >  1:36651/3‑61 (MQ=255)
                                                              gtgcgGTATCTGGCTTCCACGCGCTGATCTCTTCCGGTACGACGCCAAAACTGCTGGc                     >  1:3204708/3‑58 (MQ=255)
                                                                    tgtCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGcgaacg                <  1:4439594/57‑1 (MQ=255)
                                                                    tgtCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGc  <  1:510735/71‑1 (MQ=255)
                                                                     |                                                                     
TGGCACTGGCCCGGTATGGACCGGTAACCTGTTCCCGTTCCTGTTTATCACCATCGCCTGTGGCGCGGTGTCTGGCTTCCATGCGCTGATCTCTTCTGGGACCACGCCGAAGATGCTGGCGAACGAAGGGCAGGCGTGC  >  minE/385735‑385873

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