Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 348,499 C→T 5.2% A266V (GCA→GTA)  gcl → glyoxylate carboligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE348,4990CT 5.2% 257.9 / 3.3 77A266V (GCA→GTA) gclglyoxylate carboligase
Reads supporting (aligned to +/- strand):  ref base C (37/36);  new base T (2/2);  total (39/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.86e-01

TATCCCGGACGATCATGAACTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGGTATCGGTAACCGTTTTGCTAACCGTCATACCGGCTCGGTA  >  minE/348431‑348569
                                                                    |                                                                      
tATCCCGGACGATCATGAACTTATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTaacgc                                                                        <  1:234385/69‑1 (MQ=255)
tATCCCGGACGATCATGAACTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTaacgc                                                                        <  1:2042594/69‑1 (MQ=255)
  tCCCGGACGATCATGAACTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTaacgcaacg                                                                    >  1:3182296/1‑71 (MQ=255)
  tCCCGGACGATCATGAACTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTaacgcaacg                                                                    >  1:695045/1‑71 (MQ=255)
    ccGGACGATCATGAACTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACgct                                                                  >  1:1035695/1‑71 (MQ=255)
    ccGGACGATCATGAACTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACgct                                                                  >  1:1169986/1‑71 (MQ=255)
    ccGGACGATCATGAACTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACgct                                                                  >  1:1816790/1‑71 (MQ=255)
    ccGGACGATCATGAACTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACgct                                                                  >  1:2251056/1‑71 (MQ=255)
    ccGGACGATCATGAACTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACgct                                                                  >  1:1601472/1‑71 (MQ=255)
      ggACGATCATGAACTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACgctgc                                                                >  1:3104841/1‑71 (MQ=255)
       gACGATCATGAACTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACgctgct                                                               >  1:147836/1‑71 (MQ=255)
       gACGATCATGAACTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACgctgct                                                               <  1:3326137/71‑1 (MQ=255)
       gACGATCATGAACTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACgctgct                                                               <  1:2165573/71‑1 (MQ=255)
       gACGATCATGAACTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACgctgct                                                               <  1:2311478/71‑1 (MQ=255)
          gATCATGAACTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACgctgct                                                               <  1:2203558/68‑1 (MQ=255)
             cATGAACTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTc                                                         >  1:776064/1‑71 (MQ=255)
             cATGAACTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTc                                                         >  1:936106/1‑71 (MQ=255)
               tGAACTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTg                                                       <  1:1104799/71‑1 (MQ=255)
               tGAACTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTg                                                       <  1:98013/71‑1 (MQ=255)
                gAACTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGa                                                      >  1:1883506/1‑71 (MQ=255)
                gAACTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGa                                                      >  1:3173954/1‑71 (MQ=255)
                gAACTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGa                                                      >  1:2729795/1‑71 (MQ=255)
                gAACTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGa                                                      >  1:3245549/1‑71 (MQ=255)
                gAACTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGa                                                      >  1:1801254/1‑71 (MQ=255)
                gAACTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGa                                                      >  1:1202364/1‑71 (MQ=255)
                gAACTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGa                                                      >  1:1469423/1‑71 (MQ=255)
                 aaCTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTaacgc                                                                        <  1:3172380/52‑1 (MQ=255)
                 aaCTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACgct                                                                  <  1:45620/58‑1 (MQ=255)
                 aaCTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGa                                                      <  1:2773190/70‑1 (MQ=255)
                         gCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGgtgt                                              <  1:2135835/70‑1 (MQ=255)
                             ggATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTgg                                          >  1:1355859/1‑70 (MQ=255)
                               aTGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGGTAt                                       <  1:1297335/71‑1 (MQ=255)
                               aTGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGGTAt                                       >  1:1481189/1‑71 (MQ=255)
                                 ggtggGTCTGCAAACCGCGCATCGTTACGGTAACGCAACgctgc                                                                <  1:3339410/44‑1 (MQ=255)
                                  gtggGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGGTATCg                                     <  1:2264146/70‑1 (MQ=255)
                                   tggGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGGTATCgg                                    <  1:2509070/70‑1 (MQ=255)
                                   tggGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGGTATCgg                                    <  1:2943657/70‑1 (MQ=255)
                                        cTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGGTATCGGTAAcc                               >  1:612440/1‑70 (MQ=255)
                                        cTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGGTATCGGTAACCg                              >  1:2636986/1‑71 (MQ=255)
                                        cTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGGTATCGGTAACCg                              >  1:2190320/1‑71 (MQ=255)
                                         tGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGttt                                            <  1:1927216/56‑1 (MQ=255)
                                         tGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGttt                                            <  1:3455673/56‑1 (MQ=255)
                                         tGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGttt                                            <  1:724459/56‑1 (MQ=255)
                                         tGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGttt                                            <  1:2441978/56‑1 (MQ=255)
                                           cAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGGTATCGGTAACCGttt                           >  1:2059699/1‑71 (MQ=255)
                                           cAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGGTATCGGTAACCGttt                           <  1:3142899/71‑1 (MQ=255)
                                           cAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGGTATCGGTAACCGttt                           <  1:2459885/71‑1 (MQ=255)
                                            aaaCCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGGTATCGGTAACCGttt                           <  1:2006405/70‑1 (MQ=255)
                                             aaCCGCGCATCGTTACGGTAACGCAACGCTGCTGGCTTCTGACATGGTGTTTGGt                                         <  1:2855934/55‑1 (MQ=255)
                                             aaCCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGGt                                         <  1:2460028/55‑1 (MQ=255)
                                             aaCCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGGt                                         <  1:3232526/55‑1 (MQ=255)
                                             aaCCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGGTATCGGTAACCGTTTTg                         <  1:2803975/71‑1 (MQ=255)
                                              aCCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGGTATCGGTAACCGTTTTGc                        >  1:809175/1‑71 (MQ=255)
                                                  agcAGCGTTGCGTTACCGTAACGATGCTGGCGTCTGACATGGTGTTTGGTATCGGTAACCGTTTTGCTAAc                    >  1:501409/2‑71 (MQ=21)
                                                  agcAGCGTTGCGTTACCGTAACGATGCTGGCGTCTGACATGGTGTTTGGTATCGGTAACCGTTTTGCTAAc                    >  1:1068052/2‑71 (MQ=21)
                                                   gcATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGGTATCGGTAACCGTTTTGCTAAcc                   >  1:691588/1‑71 (MQ=255)
                                                   gcATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGGTATCGGTAACCGTTTTGCTAAcc                   >  1:1688327/1‑71 (MQ=255)
                                                   gcATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGGTATCGGTAACCGTTTTGCTAAcc                   >  1:1901217/1‑71 (MQ=255)
                                                   gcATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGGTATCGGTAACCGTTTTGCTAAcc                   >  1:3214182/1‑71 (MQ=255)
                                                   gcATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGGTATCGGTAACCGTTTTGCTAAcc                   >  1:1003651/1‑71 (MQ=255)
                                                   gcATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGGTATCGGTAACCGTTTTGCTAAcc                   >  1:2096109/1‑71 (MQ=255)
                                                     agcGTTGCGTTACCGTAACGATGCTGGCGTCTGACATGGTGTTTGGTATCGGTAACCGTTTTGCTAACCGt                 <  1:2760212/69‑1 (MQ=25)
                                                     agcGTTGCGTTACCGTAACGATGCTGGCGTCTGACATGGTGTTTGGTATCGGTAACCGTTTTGCTAACCGt                 <  1:1000793/69‑1 (MQ=25)
                                                     aTCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGGTATCGGTAACCGTTTTGCTAACCGt                 <  1:3056253/71‑1 (MQ=255)
                                                     aTCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGGTATCGGTAACCGTTTTGCTAACCGt                 <  1:3257243/71‑1 (MQ=255)
                                                     aTCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGGTATCGGTAACCGTTTTGCTAACCGt                 <  1:2189034/71‑1 (MQ=255)
                                                     aTCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGGTATCGGTAACCGTTTTGCTAACCGt                 <  1:2749386/71‑1 (MQ=255)
                                                     aTCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGGTATCGGTAACCGTTTTGCTAACCGt                 <  1:774620/71‑1 (MQ=255)
                                                       cGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTTGTATCGGTAACCg                              >  1:2719987/1‑56 (MQ=255)
                                                       cGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGGTATCGGTAACCGTTTTGCTAACCGTCa               >  1:557335/1‑71 (MQ=255)
                                                       cGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGGTATCGGTAACCGTTTTGCTAACCGTCa               >  1:2705017/1‑71 (MQ=255)
                                                       cGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGGTATCGGTAACCGTTTTGCTAACCGTCa               >  1:2700980/1‑71 (MQ=255)
                                                        gTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGGTATCGGTAACCGTTTTGCTAACCGTCa               >  1:3198605/1‑70 (MQ=255)
                                                         ttACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGGTATCgg                                    <  1:585307/48‑1 (MQ=255)
                                                         ttACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGGTATCGGTAACCGTTTTGCTAACCGTCATa             <  1:872001/71‑1 (MQ=255)
                                                         ttACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGGTATCGGTAACCGTTTTGCTAACCGTCATa             <  1:1204832/71‑1 (MQ=255)
                                                         ttACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGGTATCGGTAACCGTTTTGCTAACCGTCATa             <  1:1628150/71‑1 (MQ=255)
                                                           aCGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGGTATCGGTAACCGTTTTGCTAACCGTCATAcc           >  1:2532854/1‑71 (MQ=255)
                                                           aCGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGGTATCGGTAACCGTTTTGCTAACCGTCATAcc           >  1:2184626/1‑71 (MQ=255)
                                                           aCGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGGTATCGGTAACCGTTTTGCTAACCGTCATAcc           >  1:2010621/1‑71 (MQ=255)
                                                            cGGTAACGCAACTCTGCTGGCGTCTGACATGGTGTTTGGTATCGGTAACCGTTTTGCTAACCGTCAt              <  1:1498821/67‑1 (MQ=255)
                                                             ggTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGGTATCGGTAACCGTTTTGCTAACCGTCAt              <  1:2671314/66‑1 (MQ=255)
                                                                    caacgcTGCTGGCGTCTGACATGGTGTTTGGTATCGGTAACCGTTTTGCTAACCGTCATACCGGCTCGGTa  <  1:279176/71‑1 (MQ=255)
                                                                    caacgcTGCTGGCGTCTGACATGGTGTTTGGTATCGGTAACCGTTTTGCTAACCGTCATACCGGCTCGGTa  <  1:3029294/71‑1 (MQ=255)
                                                                    |                                                                      
TATCCCGGACGATCATGAACTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGGTATCGGTAACCGTTTTGCTAACCGTCATACCGGCTCGGTA  >  minE/348431‑348569

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