Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,473,718 C→T 25.9% F331F (TTC→TTT glmS → L‑glutamine:D‑fructose‑6‑phosphate aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,473,7180CT25.9% 97.8 / 34.2 54F331F (TTC→TTTglmSL‑glutamine:D‑fructose‑6‑phosphate aminotransferase
Reads supporting (aligned to +/- strand):  ref base C (23/17);  new base T (8/6);  total (31/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTTCCCGCTACTGGTTTGAATCGCTAGCAGGTATTCCGTGCGACGTCGAAATCGCCTCTGAATTCCGCTATCGCAAATCTGCCGTGCGTCGTAACAGCCTGATGATCACCTTGTCACAGTCTGGCGAAACCGCG  >  minE/2473653‑2473787
                                                                 |                                                                     
gTTTCCCGCTACTGGTTTGAATCGCTAGCAGGTATTCCGTGCGACGTCGAAATCGCCTCTGAATTTCGCTa                                                                  >  1:265157/1‑71 (MQ=255)
gTTTCCCGCTACTGGTTTGAATCGCTAGCAGGTATTCCGTGCGACGTCGAAATCGCCTCTGAATTTCGCTa                                                                  <  1:1821348/71‑1 (MQ=255)
   tCCCGCTAATGGTTTGAATCGCTAGCAGGTATTCCGTGCGACGTCGAAATCGCCTCTGAATTCCGCTATCg                                                               >  1:634918/1‑71 (MQ=255)
     ccGCTACTGGTTTGAATCGCTAGCAGGTATTCCGTGCGACGTCGAAATCGCCTCTGAATTTCGCTATCGc                                                              <  1:170253/70‑1 (MQ=255)
     ccGCTACTGGTTTGAATCGCTAGCAGGTATTCCGTGCGACGTCGAAATCGCCTCTGAATTCCGCTATCGc                                                              >  1:1376055/1‑70 (MQ=255)
      cGCTACTGGTTTGAATCGCTAGCAGGTATTCCGTGCGACGTCGAAATCGCCTCTGAATTCCGCTATCGCaa                                                            <  1:1958927/71‑1 (MQ=255)
           cTGGTTTGAATCGCTAGCAGGTATTCCGTGCGACGTCGAAATCGCCTCTGAATTcc                                                                      <  1:1298440/56‑1 (MQ=255)
             ggTTTGAATCGCTAGCAGGTATTCCGTGCGACGTCGAAATCGCCTCTGAATTCCGCTATCGCaaa                                                           <  1:849720/65‑1 (MQ=255)
             ggTTTGAATCGCTAGCAGGTATTCCGTGCGACGTCGAAATCGCCTCTGAATTCCGCTATCGCAAATCTGcc                                                     >  1:1695648/1‑71 (MQ=255)
              gTTTGAATCGCTAGCAGGTATTCCGTGCGACGTCGAAATCGCCTCTGAATTTCGCTATCGCAAATCTGCCg                                                    >  1:221005/1‑71 (MQ=255)
              gTTTGAATCGCTAGCAGGTATTCCGTGCGACGTCGAAATCGCCTCTGAATTCCGCTAt                                                                 >  1:2120803/1‑58 (MQ=255)
                ttGAATCGCTAGCAGGTATTCCGTGCGACGTCGAAATCGCCTCTGAATTCCGCTATCGCAAATCTGCCGTg                                                  >  1:248361/1‑71 (MQ=255)
                 tGAATCGCTAGCAGGTATTCCGTGCGACGTCGAAATCGCCTCTGAATTCCGCTATCGCAAATCTGCCGTg                                                  <  1:1663180/70‑1 (MQ=255)
                   aaTCGCTAGCAGGTATTCCGTGCGACGTCGAAATCGCCTCTGAATTCCGCTATCGCaaa                                                           <  1:1890233/59‑1 (MQ=255)
                    aTCGCTAGCAGGTATTCCGTGCGACGTCGAAATCGCCTCTGAATTCCGCTATCGCAAATCTGCCGTGcgtc                                              >  1:618196/1‑71 (MQ=255)
                      cGCTAGCAGGTATTCCGTGCGACGTCGAAATCGCCTCTGAATTCCGCTATCGCAAATCTGCCGTGcgtcg                                             >  1:1042841/1‑70 (MQ=255)
                      cGCTAGCAGGTATTCCGTGCGACGTCGAAATCGCCTCTGAATTCCGCTATCGCAAATCTGCCGTGcgtcg                                             >  1:145761/1‑70 (MQ=255)
                       gCTAGCAGGTATTCCGTGCGACGTCGAAATCGCCTCTGAATTTCGCTATCGCAAATCTGCCGTGCGTCGTa                                           <  1:1288039/71‑1 (MQ=255)
                       gCTAGCAGGTATTCCGTGCGACGTCGAAATCGCCTCTGAATTCCGCTATCGCAAATCTGCCGTGcgtcgt                                            <  1:316043/70‑1 (MQ=255)
                           gCAGGTATTCCGTGCGACGTCGAAATCGCCTCTGAATTCCGCTATCGCAAATCTGCCGTGCGTCGTAACa                                        <  1:1789742/70‑1 (MQ=255)
                              ggTATTCCGTGCGACGTCGAAATCGCCTCTGAATTCCGCTATCGCAAATCTGCCGTGCGTCGTAACAGCCt                                    <  1:1344888/71‑1 (MQ=255)
                               gTATTCCGTGCGACGTCGAAATCGCCTCTGAATTCCGc                                                                    <  1:173402/38‑1 (MQ=255)
                                tATTCCGTGCGACGTCGAAATCGCCTCTGAATTTCGCTATCGCAAATCTGCCGTGc                                                 >  1:425150/1‑56 (MQ=255)
                                  ttCCGTGCCACGTCGAAATCGCCTCTGAATTCCGCTATCGCaa                                                            <  1:1514631/43‑1 (MQ=255)
                                    ccGTGCGACGTCGAAATCGCCTCTGAATTCCGCTa                                                                  >  1:197205/1‑35 (MQ=255)
                                    ccGTGCGACGTCGAAATCGCCTCTGAATTCCGCTATCGCAAAt                                                          >  1:1836365/1‑43 (MQ=255)
                                    ccGTGCGACGTCGAAATCGCCTCTGAATTCCGCTATCGCAAATCTGCCGTGCGTCGTAACAGCCtgatga                               >  1:512219/1‑70 (MQ=255)
                                        gCGACGTCGAAATCGCCTCTGAATTCCGCTATCGCAAATCTGCCGTGCGTCGTAACAGCCTGATGATCa                            >  1:905600/1‑69 (MQ=255)
                                          gACGTCGAAATCGCCTCTGAATTTCGCTATCGCAAATCTGCCGTGCGTCGTAACAGCCTgatgat                              >  1:866475/1‑65 (MQ=255)
                                              tCGAAATCGCCTCTGAATTTCGCTATCGCAAATCTGCCGTGCGTCGTAACAGCCtgat                                 <  1:2104701/58‑1 (MQ=255)
                                              tCGAAATCGCCTCTGAATTTCGCTATCGCAAATCTGCCGTGCGTCGTAACAGCCtgat                                 <  1:775696/58‑1 (MQ=255)
                                                gAAATCGCCTCTGAATTCCGCTATCGCAAATCTGCCGTGCGTCGTAACAGCCTGATGATCACCTTGTcaca                  <  1:800716/71‑1 (MQ=255)
                                                gAAATCGCCTCTGAATTCCGCTATCGCAAATCTGCCGTGCGTCGTAACAGCCTGATGATCACCTTGTca                    >  1:531391/1‑69 (MQ=255)
                                                gAAATCGCCTCTGAATTCCGCTATCGCAAATCTGCCGTGCGTCGTAACAGCCTGATGATCACCTTGTca                    >  1:980424/1‑69 (MQ=255)
                                                     cGCCTCTGAATTCCGCTATCGCAAATCTGCCGTGCGTCGTAACAGCCTgatgat                              >  1:32043/1‑54 (MQ=255)
                                                     cGCCTCTGAATTCCGCTATCGCAAATCTGCCGTGCGTCGTAACAGCCTGATGATCACCTTGTCACAGTCTg             >  1:1197213/1‑71 (MQ=255)
                                                      gCCTCTGAATTCCGCTATCGCAAATCTGCCGTGCGTCGTAACAGCCTGATGATCACCTTGt                      <  1:234071/61‑1 (MQ=255)
                                                      gCCTCTGAATTCCGCTATCGCAAATCTGCCGTGCGTCGTAACAGCCTGATGATCACCTTGt                      <  1:2069471/61‑1 (MQ=255)
                                                      gCCTCTGAATTCCGCTATCGCAAATCTGCCGTGCGTCGTAACAGCCTGATGATCACCTTGt                      <  1:1501136/61‑1 (MQ=255)
                                                      gCCTCTGAATTCCGCTATCGCAAATCTGCCGTGCGTCGTAACAGCCTGATGATCACCTTGt                      <  1:1867122/61‑1 (MQ=255)
                                                      gCCTCTGAATTCCGCTATCGCAAATCTGCCGTGCGTCGTAACAGCCTGATGATCACCTTGt                      <  1:777865/61‑1 (MQ=255)
                                                      gCCTCTGAATTCCGCTATCGCAAATCTGCCGTGCGTCGTAACAGCCTGATGATCACCTTGTCACAGTCTg             <  1:455781/70‑1 (MQ=255)
                                                        ctctGAATTCCGCTATCGCAAATCTGCCGTGCGTCGTAACAGCCTGATGATCACCTTGTca                    >  1:1875116/1‑61 (MQ=255)
                                                        ctctGAATTCCGCTATCGCAAATCTGCCGTGCGTCGTAACAGCCTGATGATCACCTTGTca                    >  1:778080/1‑61 (MQ=255)
                                                        ctctGAATTCCGCTATAGCAAATCTGCCGTGCGTCGTAACAGCCTGATGATCACCTTGTca                    >  1:552105/1‑61 (MQ=255)
                                                          ctGAATTTCGCTATCGCAAATCTGCCGTGCGTCGTAACAGCCTGATGATCACCTTGTCACAGTCTGGCGa         >  1:599468/1‑70 (MQ=255)
                                                          ctGAATTTCGCTATCGCAAATCTGCCGTGCGTCGTAACAGCCTGATGATCACCTTGTCACAGTCTGGCGa         >  1:65789/1‑70 (MQ=255)
                                                          ctGAATTTCGCTATCGCAAATCAGCCGTGCGTCGTAACAGCCTGATGATCACCTTGTCACAGTCTGGCg          >  1:770892/1‑69 (MQ=255)
                                                          ctGAATTCCGCTATCGCAAATCTGCCGTGCGTCGTAACAGCCTGATGATCACCTTGTCACAGTCTGGCGa         >  1:2039083/1‑70 (MQ=255)
                                                            gAATTTCGCTATCGCAAATCTGCCGTGCGTCGTAACAGCCTGATGATCACCTTGTCACAGTCTGGCGAAAc      >  1:972656/1‑71 (MQ=255)
                                                            gAATTCCGCTATCGCAAATCTGCCGTGCGTCGTAACAGCCTGATGATCACCTTGTCACAGTCTGGCGAAAc      >  1:1673159/1‑71 (MQ=255)
                                                            gAATTCCGCTATCGCAAATCTGCCGTGCGTCGTAACAGCCTGATGATCACCTTGTCACAGTCTGGCGAAAc      >  1:1661473/1‑71 (MQ=255)
                                                             aaTTCCGCTATCGCAAATCTGCCGTGCGTCGTAACAGCCTGATGATCACCTTGTc                     <  1:1645429/55‑1 (MQ=255)
                                                              aTTTCGCTATCGCAAATCTGCCTTGCGTCGTAACAGCCTGATGATCACCTTGTCACAGTCTGGCGAAAcc     <  1:232499/70‑1 (MQ=255)
                                                               ttCCGCTATCGCAAATCTGCCGTGCGTCGTAACAGCCTGATGATCACCTTGTcac                   >  1:764557/1‑55 (MQ=255)
                                                                 ccGCTATCGCAAATCTGCCGTGCGTCGTAACAGCCTGATGATCACCTTGTCACAGTCTGGCGAAACcgcg  <  1:697514/70‑1 (MQ=255)
                                                                 ccGCTATCGCAAATCTGCCGTGCGTCGTAACAGCCTGATGATCACCTTGTCACAGTCTGGCGAAACcgcg  <  1:336040/70‑1 (MQ=255)
                                                                 |                                                                     
GTTTCCCGCTACTGGTTTGAATCGCTAGCAGGTATTCCGTGCGACGTCGAAATCGCCTCTGAATTCCGCTATCGCAAATCTGCCGTGCGTCGTAACAGCCTGATGATCACCTTGTCACAGTCTGGCGAAACCGCG  >  minE/2473653‑2473787

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: