Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 385,411 G→T 10.0% G207G (GGG→GGT cstA → carbon starvation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE385,4110GT10.0% 206.0 / 4.4 67G207G (GGG→GGTcstAcarbon starvation protein
Reads supporting (aligned to +/- strand):  ref base G (27/33);  new base T (2/5);  total (29/38)
Fisher's exact test for biased strand distribution p-value = 6.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGCCCTGACTCATAGCCCGTGGGGAACATATACCGTTGCGTTCACCATTCCGCTGGCGCTGTTTATGGGGATCTACCTGCGCTATCTGCGTCCGGGGCGTATTGGCGAAGTGTCGGTCATCGGTCTGGTATTCCTGATTT  >  minE/385341‑385481
                                                                      |                                                                      
aaGCCCTGACTCATAGCCCGTGGGGAACATATACCGTTGCGTTCACCATTCCGCTGGCGCTGTTTATgggg                                                                        <  1:806619/71‑1 (MQ=255)
aaGCCCTGACTCATAGCCCGTGGGGAACATATACCGTTGCGTTCACCATTCCGCTGGCGCTGTTTATgggg                                                                        <  1:1612222/71‑1 (MQ=255)
 aGCCCTGACTCATAGCCCGTGGGGAACATATACCGTTGCGTTCACCATTCCGCTGGCGCTGTTTATGGGGa                                                                       <  1:1156335/71‑1 (MQ=255)
   cccTGACTCATAGCCCGTGGGGAACATATACCGTTGCGTTCACCATTCCGCTGGCGCTGTTTATGGGGATc                                                                     <  1:4021655/71‑1 (MQ=255)
    ccTGACTCATAGCCCGTGGGGAACATATACCGTTGCGTTCACCATTCCGCTGGCGCTGTTTATGGGGATCt                                                                    <  1:3577421/71‑1 (MQ=255)
     cTGACTCATAGCCCGTGGGGAACATATACCGTTGCGTTCACCATTCCGCTGGCGCTGTTTATGGGGATCTa                                                                   >  1:1465581/1‑71 (MQ=255)
        aCTCATAGCCCGTGGGGAACATATACCGTTGCGTTCACCATTCCGCTGGCGCTGTTTATGGGGATCt                                                                    <  1:335640/67‑1 (MQ=255)
        aCTCATAGCCCGTGGGGAACATATACCGTTGCGTTCACCATTCCGCTGGCGCTGTTTATGGGGATCTAc                                                                  <  1:2109660/69‑1 (MQ=255)
           cATAGCCCGTGGGGAACATATACCGTTGCGTTCACCATTCCGCTGGCGCTGTTTATGGGGATCTAc                                                                  >  1:5350610/1‑66 (MQ=255)
           cATAGCCCGTGGGGAACATATACCGTTGCGTTCACCATTCCGCTGGCGCTGTTTATGGGGATCTAc                                                                  >  1:2360353/1‑66 (MQ=255)
            aTAGCCCGTGGGGAACATATACCGTTGCGTTCACCATTCCGCTGGCGCTGTTTATGGGGATCTAc                                                                  >  1:921728/1‑65 (MQ=255)
             tAGCCCGTGGGGAACATATACCGTTGCGTTCACCATTCCGCTGGCGCTGTTTATGGGGATCTACCTGCGCt                                                           <  1:2931820/71‑1 (MQ=255)
             tAGCCCGTGGGGAACATATACCGTTGCGTTCACCATTCCGCTGGCGCTGTTTATGGGGATCTACCTGCGCt                                                           <  1:2673099/71‑1 (MQ=255)
             tAGCCCGTGGGGAACATATACCGTTGCGTTCACCATTCCGCTGGCGCTGTTTATGGGGATCTACCTGCGCt                                                           <  1:5321417/71‑1 (MQ=255)
             tAGCCCGTGGGGAACATATACCGTTGCGTTCACCATTCCGCTGGCGCTGTTTATGGGGATCTACCTGCGCt                                                           <  1:3757597/71‑1 (MQ=255)
                   gTGGGGAACATATACCGTTGCGTTCACCATTCCGCTGGCGCTGTTTATGGGGATCTAc                                                                  >  1:3358542/1‑58 (MQ=255)
                   gTGGGGAACATATACCGTTGCGTTCACCATTCCGCTGGCGCTGTTTATGGGGATCTACCTGCGCTATCTGc                                                     >  1:4831295/1‑71 (MQ=255)
                    tGGGGAACATATACCGTTGCGTTCACCATTCCGCTGGCGCTGTTTATGGGGATCTACCTGCGCTATCTGCg                                                    >  1:2575199/1‑71 (MQ=255)
                        gAACATATACCGTTGCGTTCACCATTCCGCTGGCGCTGTTTATGGGGATCTACATGCGCTATCTGCGTCCg                                                <  1:1323712/71‑1 (MQ=255)
                          aCATATACCGTTGCGTTCACCATTCCGCTGGCGCTGTTTATGGGGATCTAc                                                                  >  1:3965388/1‑51 (MQ=255)
                          aCATATACCGTTGCGTTCACCATTCCGCTGGCGCTGTTTATGGGGATCTAc                                                                  >  1:4508118/1‑51 (MQ=255)
                          aCATATACCGTTGCGTTCACCATTCCGCTGGCGCTGTTTATGGGGATCTAc                                                                  >  1:3387429/1‑51 (MQ=255)
                          aCATATACCGTTGCGTTCACCATTCCGCTGGCGCTGTTTATGGGGATCTACCTGCGCTATCTGCGTCCggg                                              <  1:951749/71‑1 (MQ=255)
                            atatACCGTTGCGTTCACCATTCCGCTGGCGCTGTTTATgggg                                                                        <  1:4955988/43‑1 (MQ=255)
                             tataCCGTTGCGTTCACCATTCCGCTGGCGCTGTTTATGGGGATCt                                                                    <  1:3577834/46‑1 (MQ=255)
                             tataCCGTTGCGTTCACCATTCCGCTGGCGCTGTTTATGGGGATCTACCTGCGCTATCTGCGTCCGGGGCg                                           >  1:2122117/1‑71 (MQ=255)
                             tataCCGTTGCGTTCACCATTCCGCTGGCGCTGTTTATGGGGATCTACCTGCGCTATCTGCGTCCGGGGCg                                           <  1:2811345/71‑1 (MQ=255)
                              ataCCGTTGCGTTCACCATTCCGCTGGCGCTGTTTATGGGGATCTACCTGCGCTAt                                                         >  1:4006322/1‑56 (MQ=255)
                              ataCCGTTGCGTTCACCATTCCGCTGGCGCTGTTTATGGGGATCTACCTGCGCTATCTGCGTCCGGGGCg                                           >  1:59461/1‑70 (MQ=255)
                               taCCGTTGCGTTCACCATTCCGCTGGCGCTGTTTATGGGGATCTACCTGCGCt                                                           <  1:3336642/53‑1 (MQ=255)
                               taCCGTTGCGTTCACCATTCCGCTGGCGCTGTTTATGGGGATCTACCTGCGCTATCTGCGTCCGGGGCGTa                                         >  1:3992845/1‑71 (MQ=255)
                                 ccGTTGCGTTCACCATTCCGCTGGCGCTGTTTATGGGGATCTACCTGCGCTATCTGCGTCCGGGGCGTAtt                                       >  1:371778/1‑71 (MQ=255)
                                 ccGTTGCGTTCACCATTCCGCTGGCGCTGTTTATGGGGATCTACCTGCGCTATCTGCGTCCGGGGCGTAtt                                       <  1:865078/71‑1 (MQ=255)
                                 ccGTTGCGTTCACCATTCCGCTGGCGCTGTTTATGGGGATCTACCTGCGCTATCTGCGTCCGGGGCGTAtt                                       <  1:4153633/71‑1 (MQ=255)
                                   gTTGCGTTCACCATTCCGCTGGCGCTGTTTATGGGTATCt                                                                    <  1:4649858/40‑1 (MQ=255)
                                    ttGCGTTCACCATTCCGCTGGCGCTGTTTATGGGGATCTACCTGCGCTATCTGCGTCCGGGGCGTATTGGc                                    <  1:3186701/71‑1 (MQ=255)
                                       cGTTCACCATTCCGCTGGCGCTGTTTATGGGGATCTACCTGCGCTATCTGCGTCCGGGGCGTATTGGCGa                                  <  1:3754821/70‑1 (MQ=255)
                                         ttCACCATTCCGCTGGCGCTGTTTATGGGGATCTACCTGCGCTATCTGCGTCCGGGGCGt                                          <  1:1014246/60‑1 (MQ=255)
                                           cACCATTCCGCTGGCGCTGTTTATGGGGATCTACCTGCGCTATCTGCGTCCGGGGCGTATTGGCGAAgtgt                             >  1:4635427/1‑71 (MQ=255)
                                           cACCATTCCGCTGGCGCTGTTTATGGGGATCTACCTGCGCTATCTGCGTCCGGGGCGTATTGGCGAAgtgt                             >  1:107078/1‑71 (MQ=255)
                                           cACCATTCCGCTGGCGCTGTTTATGGGGATCTACCTGCGCTATCTGCGTCCGGGGCGTATTGGCGAAgtg                              <  1:1659745/70‑1 (MQ=255)
                                           cACCATTCCGCTGGCGCTGTTTATGGGGATCTACCTGCGCTATCTGCGTCCGGGGCGTATTGGCGAAgtg                              >  1:5071817/1‑70 (MQ=255)
                                            aCCATTCCGCTGGCGCTGTTTATTGGGATCTACCTGCGCTATCTGCGTCCGGGGCGTATTGGCGAAGTGTc                            >  1:4454511/1‑71 (MQ=255)
                                            aCCATTCCGCTGGCGCTGTTTATGGGGATCTACCTGCGCTATCt                                                       <  1:3700600/44‑1 (MQ=255)
                                            aCCATTCCGCTGGCGCTGTTTATGGGGATCTACCTGCGCTATCt                                                       <  1:4454013/44‑1 (MQ=255)
                                            aCCATTCCGCTGGCGCTGTTTATGGGGATCTACCTGCGCTATCt                                                       <  1:324285/44‑1 (MQ=255)
                                            aCCATTCCGCTGGCGCTGTTTATGGGGATCTACCTGCGCTATCTGCGTc                                                  >  1:295969/1‑49 (MQ=255)
                                                ttCCGCTGTCGCTGTTTATGGGGATCTACCTGCGCTATCTGCGTCCGGGGCGTATTGGCGAAgtgg                             <  1:163748/66‑2 (MQ=255)
                                                ttCCGCTGGCGCTGTTTATGGGGATCTACCTGCGCt                                                           <  1:5526565/36‑1 (MQ=255)
                                                ttCCGCTGGCGCTGTTTATGGGGATCTACCTGCGCt                                                           <  1:2132766/36‑1 (MQ=255)
                                                ttCCGCTGGCGCTGTTTATGGGGATCTACCTGCGCTATCTGCGTCCGGGGCGTATTGGCg                                   <  1:2667754/60‑1 (MQ=255)
                                                ttCCGCTGGCGCTGTTTATGGGGATCTACCCGCGCTATCTGCGTCCGGGGCGTATTGGCGAAGTGTCg                           <  1:3361179/68‑1 (MQ=255)
                                                 tCCGCTGGCGCTGTTTATGGGGATCTACCTGCGCTATCTGCGTCCGGGGCGTATTGGCGAAGTGTCGGTCa                       >  1:4031201/1‑71 (MQ=255)
                                                       tgcgcTGTTTATGGGTATCTACATGCGCTTTATTCGTCCGGGGCGTGTGGGTGAAGTCTCTGTc                        >  1:4899344/2‑64 (MQ=255)
                                                       tgcgcTGTTTATGGGTATCTACATGCGCTTTATTCGTCCGGGGCGTGTGGGTGAAGTCTCTGTc                        >  1:4975636/2‑64 (MQ=255)
                                                        gcgcTGTTTATGGGTATCTACATGCGCTTTATTCGTCCGGGGCGTGTGGGTGAAGTCTCTGTCATTGGta                 <  1:554944/70‑2 (MQ=255)
                                                        gcgcTGTTTATGGGGATCTACCTGCGCTATCTGCGTCCGGGGCGTATTGGCGAAGTGTCGGTCATCGGTc                 <  1:3730904/70‑1 (MQ=255)
                                                         cgcTGTTTATGGGTATCTACATGCGCTTTATTCGTCCGGGGCGTGTGGGTGAAGTCTCTGTCATTGGtatc               <  1:2715536/71‑4 (MQ=255)
                                                         cgcTGTTTATGGGGATCTACCTGCGCTATCTGCGTCCgggg                                             >  1:1152620/1‑41 (MQ=255)
                                                          gcTGTTTATGGGGATCTACCTGCGCTATCTGCGTCCGGGGCGTATTGGCGAAGTGTCGGTCATCGGTCTgg              <  1:2204072/71‑1 (MQ=255)
                                                           cTGTTTATGGGTATCTACATGCGCTTTATTCGTCCGGGGCGTGTGGGTGAAgtct                             <  1:4865738/55‑3 (MQ=255)
                                                             gTTTATGGGGATCTACCTGCGCTATCTGCGTCCGGGGCGTATTGGCGAAGTGTCGGTCATCGGTCTgg              <  1:4519854/68‑1 (MQ=255)
                                                             gTTTATGGGGATCTACCTGCGCTATCTGCGTCCGGGGCGTATTGGCGAAGTGTCGGTCATCGGTCTGGTAt           >  1:814768/1‑71 (MQ=255)
                                                             gTTTATGGGGATCTACATGCGCTATCTGCGTCCGGGGCGTATTGGCGAAGTGTCGGTCATCGGTCTGGTAt           >  1:1186824/1‑71 (MQ=255)
                                                              tttATGGGTATCTACATGCGCTTTATTCGTCCGGGGCGt                                          <  1:2453367/39‑1 (MQ=19)
                                                              tttATGGGGATCTACCTGCGCTATCTGCGTCCGGGGCGTATTGGCGAAGTGTCGGTCATCGGTCTGGTAtt          <  1:914747/71‑1 (MQ=255)
                                                               ttATGGGGATCTACCTGCGCTATCTGCGTCCGGGGCGTATTGGCGAAGTGTCGGTCATCGGTCTGGTa            >  1:671109/1‑68 (MQ=255)
                                                                tATGGGGATCTACCTGCGCTATCTGCGTCCGGGGCGTATTGGCGAAGTGTCGGTCATc                     >  1:2601464/1‑58 (MQ=255)
                                                                  tGGGGATCTACCTGCGCTATCTGCGTCCGGGGCGTATTGGCGAAGTGTCGGTCAt                      >  1:2479440/1‑55 (MQ=255)
                                                                  tGGGGATCTACCTGCGCTATCTGCGTCCGGGGCGTATTGGCGAAGTGTCGGTCATCGGTCTGGTATTCCTg      <  1:2598206/71‑1 (MQ=255)
                                                                    gggATCTACCTGCGCTATCTGCGTCCGGGGCGTATTGGCGAAg                                >  1:5550889/1‑43 (MQ=255)
                                                                    gggATCTACCTGCGCTATCTGCGTCCGGGGCGTATTGGCGAAg                                >  1:4471752/1‑43 (MQ=255)
                                                                    gggATCTACCTGCGCTATCTGCGTCCGGGGCGTATTGGCGAAg                                >  1:956833/1‑43 (MQ=255)
                                                                      gATCTACCTGCGCTATCTGCGTCCGGGGCGTATTGGCGAAGTGTCGGTCATCGGTCTGGTATTCCTGAttt  <  1:1930464/71‑1 (MQ=255)
                                                                      |                                                                      
AAGCCCTGACTCATAGCCCGTGGGGAACATATACCGTTGCGTTCACCATTCCGCTGGCGCTGTTTATGGGGATCTACCTGCGCTATCTGCGTCCGGGGCGTATTGGCGAAGTGTCGGTCATCGGTCTGGTATTCCTGATTT  >  minE/385341‑385481

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: