Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,627,710 C→T 31.4% G60G (GGC→GGT livJ → leucine/isoleucine/valine transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,627,7100CT31.4% 66.4 / 33.9 73G60G (GGC→GGTlivJleucine/isoleucine/valine transporter subunit
Reads supporting (aligned to +/- strand):  ref base C (44/6);  new base T (11/12);  total (55/18)
Fisher's exact test for biased strand distribution p-value = 7.67e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

TGCGCAGTACGGTGACCAGGAGTTTACCGGCGCAGAGCAGGCGGTTGCGGATATCAACGCTAAAGGCGGCATTAAAGGCAACAAACTGCAAATCGTAAAATATGACGATGCCTGTGACCCGAAACAGGCGGTTGC  >  minE/2627644‑2627778
                                                                  |                                                                    
tGCGCAGTACGGTGACCAGGAGTTTACCGGCGCAGAGCAGGCGGTTGCGGATATCAACGCTAAAggcgg                                                                    >  1:2333492/1‑69 (MQ=255)
     aGTACGGTGACCAGGAGTTTACCGGCGCAGAGCAGGCGGTTGCGGATATCAACGCTAAAGGCGGCATTa                                                               <  1:2478350/69‑1 (MQ=255)
     aGTACGGTGACCAGGAGTTTACCGGCGCAGAGCAGGCGGTTGCGGATATCAACGCTAAAGGCGGCATTa                                                               <  1:303625/69‑1 (MQ=255)
      gTACGGTGACCAGGAGTTTACCGGCGCAGAGCAGGCGGTTGCGGATATCAACGCTAAAGGCGGCATTaa                                                              >  1:709507/1‑69 (MQ=255)
      gTACGGTGACCAGGAGTTTACCGGCGCAGAGCAGGCGGTTGCGGATATCAACGCTAAAGGCGGCATTaa                                                              >  1:3395974/1‑69 (MQ=255)
      gTACGGTGACCAGGAGTTTACCGGCGCAGAGCAGGCGGTTGCGGATATCAACGCTAAAGGCGGCATTaa                                                              >  1:58305/1‑69 (MQ=255)
      gTACGGTGACCAGGAGTTTACCGGCGCAGAGCAGGCGGTTGCGGATATCAACGCTAAAGGCGGCATTaa                                                              >  1:635529/1‑69 (MQ=255)
      gTACGGTGACCAGGAGTTTACCGGCGCAGAGCAGGCGGTTGCGGATATCAACGCTAAAGGCGGCATTaa                                                              >  1:2293922/1‑69 (MQ=255)
      gTACGGTGACCAGGAGGTTACCGGCGCAGAGCAGGCGGTTGCGGATATCAACGCTAAAGGCGGCATTaa                                                              >  1:2490717/1‑69 (MQ=255)
        aCGGTGACCAGGAGTTTACCGGCGCAGAGCAGGCGGTTGCGGATATCAACGCTAAAGGCGGCATTaaa                                                             >  1:3393415/1‑68 (MQ=255)
         cGGTGACCAGGAGTTTACCGGCGCAGAGCAGGCGGTTGCGGATATCAACGCTAAAGGCGGCATTAAAg                                                            >  1:1500980/1‑68 (MQ=255)
         cGGTGACCAGGAGTTTACCGGCGCAGAGCAGGCGGTTGCGGATATCAACGCTAAAGGCGGCATTAAAg                                                            >  1:2466989/1‑68 (MQ=255)
         cGGTGACCAGGAGTTTACCGGCGCAGAGCAGGCGGTTGCGGATATCAACGCTAAAGGCGGCATTAAAg                                                            >  1:3058623/1‑68 (MQ=255)
         cGGTGACCAGGAGTTTACCGGCGCAGAGCAGGCGGTTGCGGATATCAACGCTAAAGGCGGCATTAAAg                                                            >  1:3205523/1‑68 (MQ=255)
           gTGACCAGGAGTTTACCGGCGCAGAGCAGGCGGTTGCGGATATCAACGCTAAAGGCGGCATTAAAGGc                                                          <  1:3127653/68‑1 (MQ=255)
              aCCAGGAGTTTACCGGCGCAGAGCAGGCGGTTGCGGATATCAACGCTAAAGGCGGCATTAAAGGcaaca                                                      <  1:993730/69‑1 (MQ=255)
                 aGGAGTTTACTGGCGCAGAGCAGGCGGTTGCGGATATCAATGCTAAAGGTGGTATCAAAGGCAACAAAc                                                   <  1:2838061/69‑1 (MQ=255)
                 aGGAGTTTACTGGCGCAGAGCAGGCGGTTGCGGATATCAATGCTAAAGGTGGTATCAAAGGCAACAAAc                                                   <  1:2206665/69‑1 (MQ=255)
                 aGGAGTTTACTGGCGCAGAGCAGGCGGTTGCGGATATCAATGCTAAAGGTGGTATCAAAGGCAACAAAc                                                   <  1:2183156/69‑1 (MQ=255)
                 aGGAGTTTACTGGCGCAGAGCAGGCGGTTGCGGATATCAATGCTAAAGGTGGTATCAAAGGCAACAAAc                                                   <  1:2855373/69‑1 (MQ=255)
                 aGGAGTTTACTGGCGCAGAGCAGGCGGTTGCGGATATCAATGCTAAAGGTGGTATCAAAGGCAACAAAc                                                   <  1:3155579/69‑1 (MQ=255)
                 aGGAGTTTACTGGCGCAGAGCAGGCGGTTGCGGATATCAATGCTAAAGGTGGTATCAAAGGCAACAAAc                                                   <  1:186553/69‑1 (MQ=255)
                 aGGAGTTTACTGGCGCAGAGCAGGCGGTTGCGGATATCAATGCTAAAGGTGGTATCAAAGGCAACAAAc                                                   <  1:674234/69‑1 (MQ=255)
                 aGGAGTTTACTGGCGCAGAGCAGGCGGTTGCGGATATCAATGCTAAAGGTGGTATCAAAGGCAACAAAc                                                   <  1:683802/69‑1 (MQ=255)
                 aGGAGTTTACTGGCGCAGAGCAGGCGGTTGCGGATATCAATGCTAAAGGTGGTATCAAAGGCAACAAAc                                                   <  1:1595870/69‑1 (MQ=255)
                 aGGAGTTTACTGGCGCAGAGCAGGCGGTTGCGGATATCAATGCTAAAGGTGGTATCAAAGGCAACAAAc                                                   <  1:1171060/69‑1 (MQ=255)
                 aGGAGTTTACTGGCGCAGAGCAGGCGGTTGCGGATATCAATGCTAAAGGTGGTATCAAAGGCAACAAAc                                                   <  1:938557/69‑1 (MQ=255)
                 aGGAGTTTACTGGCGCAGAGCAGGCGGTTGCGGATATCAATGCTAAAGGTGGTATCAAAGGCAACAAAc                                                   <  1:939094/69‑1 (MQ=255)
                     gTTTACCGGCGCAGAGCAGGCGGTTGCGGATATCAACGCTAAAGGCGGCATTAAAGGCAACAAACTGCa                                               >  1:490858/1‑69 (MQ=255)
                     gTTTACCGGCGCAGAGCAGGCGGTTGCGGATATCAACGCTAAAGGCGGCATTAAAGGCAACAAACTGCa                                               >  1:2844070/1‑69 (MQ=255)
                     gTTTACCGGCGCAGAGCAGGCGGTTGCGGATATCAACGCTAAAGGCGGCATTAAAGGCAACAAACTGCa                                               >  1:1509812/1‑69 (MQ=255)
                     gTTTACCGGCGCAGAGCAGGCGGTTGCGGATATCAACGCTAAAGGCGGCATTAAAGGCAACAAACTGCa                                               >  1:1440415/1‑69 (MQ=255)
                        tACCGGCGCAGAGCAGGCGGTTGCGGATATCAATGCTAAAGGCGGCATTAAAGGCAACAAACTGCAAAt                                            >  1:1955485/1‑69 (MQ=255)
                        tACCGGCGCAGAGCAGGCGGTTGCGGATATCAACGCTAAAGGCGGCATTAAAGGCAACAAACTGCAAAt                                            >  1:772133/1‑69 (MQ=255)
                        tACCGGCGCAGAGCAGGCGGTTGCGGATATCAACGCTAAAGGCGGCATTAAAGGCAACAAACTGCAAAt                                            >  1:2777476/1‑69 (MQ=255)
                        tACCGGCGCAGAGCAGGCGGTTGCGGATATCAACGCTAAAGGCGGCATTAAAGGCAACAAACTGCAAAt                                            >  1:2861789/1‑69 (MQ=255)
                        tACCGGCGCAGAGCAGGCGGTTGCGGATATCAACGCTAAAGGCGGCATTAAAGGCAACAAACTGCAAAt                                            >  1:842715/1‑69 (MQ=255)
                        tACCGGCGCAGAGCAGGCGGTTGCGGATATCAACGCTAAAGGCGGCATTAAAGGCAACAAACTGCAAAt                                            >  1:67877/1‑69 (MQ=255)
                        tACCGGCGCAGAGCAGGCGGTTGCGGATATCAACGCTAAAGGCGGCATTAAAGGCAACAAACTGCAAAt                                            >  1:3278630/1‑69 (MQ=255)
                        tACCGGCGCAGAGCAGGCGGTTGCGGATATCAACGCTAAAGGCGGCATTAAAGGCAACAAACTGCAAAt                                            >  1:1079488/1‑69 (MQ=255)
                        tACCGGCGCAGAGCAGGCGGTTGCGGATATCAACGCTAAAGGCGGCATTAAAGGCAACAAACTGCAAAt                                            >  1:2097566/1‑69 (MQ=255)
                        tACCGGCGCAGAGCAGGCGGTTGCGGATATCAACGCTAAAGGCGGCATTAAAGGCAACAAACTGCAAAt                                            >  1:2147211/1‑69 (MQ=255)
                        tACCGGAGCAGAGCAGGCGGTTGCGGATATCAACGCTAAAGGCGGCAtt                                                                >  1:6098/1‑49 (MQ=255)
                           tggCGCAGAGCAGGCGGTTGCGGATATCAATGCTAAAGGTGGTATCAAGGGCAACAAACTGCAAATCg                                          >  1:2125563/2‑68 (MQ=255)
                           tggCGCAGAGCAGGCGGTTGCGGATATCAATGCTAAAGGTGGTATCAAAGGCAACAAACTGCAAATCg                                          >  1:3390428/2‑68 (MQ=255)
                           tggCGCAGAGCAGGCGGTTGCGGATATCAATGCTAAAGGTGGTATCAAAGGCAACAAACTGCAAATCg                                          >  1:1700805/2‑68 (MQ=255)
                           tggCGCAGAGCAGGCGGTTGCGGATATCAATGCTAAAGGTGGTATCAAAGGCAACAAACTGCAAATCg                                          >  1:667680/2‑68 (MQ=255)
                           tggCGCAGAGCAGGCGGTTGCGGATATCAATGCTAAAGGTGGTATCAAAGGCAACAAACTGCAAATCg                                          >  1:1859017/2‑68 (MQ=255)
                           tggCGCAGAGCAGGCGGTTGCGGATATCAATGCTAAAGGTGGTATCAAAGGCAACAAACTGCAAATCg                                          >  1:522364/2‑68 (MQ=255)
                           tggCGCAGAGCAGGCGGTTGCGGATATCAATGCTAAAGGTGGTATCAAAGGCAACAAACTGCAAATCg                                          >  1:2069736/2‑68 (MQ=255)
                           tggCGCAGAGCAGGCGGTTGCGGATATCAATGCTAAAGGTGGTATCAAAGGCAACAAACTGCAAATCg                                          >  1:2759442/2‑68 (MQ=255)
                           tggCGCAGAGCAGGCGGTTGCGGATATCAATGCTAAAGGTGGTATCAAAGGCAACAAACTGCAAATCg                                          >  1:2538341/2‑68 (MQ=255)
                           tggCGCAGAGCAGGCGGTTGCGGATATCAATGCTAAAGGTGGTATCAAAGGCAACAAACTGCAAATCg                                          >  1:2869925/2‑68 (MQ=255)
                           tggCGCAGAGCAGGCGGTTGCGGATATCAATGCTAAAGGTGGTATCAAAGGAAACAAACTGCAAATCg                                          >  1:3372450/2‑68 (MQ=255)
                           cGGCGCAGAGCAGGCGGTTGCGGATATCAACGCTAAAGTCGGCATTAAAGGcaact                                                      <  1:846767/56‑2 (MQ=255)
                                     cAGGCGGTTGCGGATATCAACGCTAAAGGCGGCATTAAAGGCAACAAACTGCAAATCGTAAAATATGAc                               <  1:2545670/69‑1 (MQ=255)
                                              gCGGATATCAACGCTAAAGGCGGCATTAAAGGCAACAAACTGCAAATCGTAAAATATGTCGATGCCtgt                      >  1:362126/1‑69 (MQ=255)
                                              gCGGATATCAACGCTAAAGGCGGCATTAAAGGCAACAAACTGCAAATCGTAAAATATGACGATGCCtgt                      >  1:2080066/1‑69 (MQ=255)
                                              gCGGATATCAACGCTAAAGGCGGCATTAAAGGCAACAAACTGCAAATCGTAAAATATGACGATGCCtgt                      >  1:1046450/1‑69 (MQ=255)
                                              gCGGATATCAACGCTAAAGGCGGCATTAAAGGCAACAAACTGCAAATCGTAAAATATGACGATGCCtgt                      >  1:47390/1‑69 (MQ=255)
                                              gCGGATATCAACGCTAAAGGCGGCATTAAAGGCAACAAACTGCAAATCGTAAAATATGACGATGCCtgt                      >  1:1630321/1‑69 (MQ=255)
                                              gCGGATATCAACGCTAAAGGCGGCATTAAAGGCAACAAACTGCAAATCGTAAAATATGACGATGCCtgt                      >  1:311844/1‑69 (MQ=255)
                                              gCGGATATCAACGCTAAAGGCGGCATTAAAGGCAACAAACTGCAAATCGTAAAATATGACGATGCCtgt                      >  1:3093772/1‑69 (MQ=255)
                                              gCGGATATCAACGCTAAAGGCGGCATTAAAGGCAACAAACTGCAAATCGTAAAATATGACGATGCCtgt                      >  1:1348027/1‑69 (MQ=255)
                                              gCGGATATCAACGCTAAAGGCGGCATTAAAGGCAACAAACTCCAAATCGTAAAATATGACGATGCCtgt                      >  1:1407478/1‑69 (MQ=255)
                                                    atCAACGCTAAAGGCGGCATTAAAGGCAACAAACTGCAAATCGTAAAATATGACGATGCCTGTGAccc                 >  1:2984028/1‑68 (MQ=255)
                                                       aaCGCTAAAGGCGGCATTAAAGGCAACAAACTGCAAATCGTAAAATATGACGATGCCTGTGACCCGa               >  1:3137931/1‑67 (MQ=255)
                                                       aaCGCTAAAGGCGGCATTAAAGGCAACAAACTGCAAATCGTAAAATATGACGATGCCTGTGACCCGa               >  1:1920010/1‑67 (MQ=255)
                                                           cTAAAGGCGGCATTAAAGGCAACAAACTGCAAATCGTAAAATATGACGATGCCTGTGACCCGAAACAgg         >  1:3048615/1‑69 (MQ=255)
                                                           cTAAAGGCGGCATTAAAGGCAACAAACTGCAAATCGTAAAATATGACGATGCCTGTGACCCGAAACAgg         >  1:155557/1‑69 (MQ=255)
                                                           cTAAAGGCGGCATTAAAGGCAACAAACTGCAAATCGTAAAATATGACGATGCCTGTGACCCGAAACAgg         >  1:2588866/1‑69 (MQ=255)
                                                           cTAAAGGCGGCATTAAAGGCAACAAACTGCAAATCGTAAAATATGACGATGCCTGTGACCCGAAACAgg         >  1:2562837/1‑69 (MQ=255)
                                                           cTAAAGGCGGCATTAAAGGCAACAAACTGCAAATCGTAAAATATGACGATGCCTGTGACCCGAAACAgg         >  1:2929862/1‑69 (MQ=255)
                                                           cTAAAGGCGGCATTAAAGGCAACAAACTGCAAATCGTAAAATATGACGATGCCTGTGACCCGAAACAgg         >  1:2425867/1‑69 (MQ=255)
                                                                 gcggcATTAAAGGCAACAAACTGCAAATCGTAAAATATGACGATGCCTGTGACCCGa               <  1:2041899/57‑1 (MQ=255)
                                                                 gcggcATTAAAGGCAACAAACTGCAAATCGTAAAATATGACGATGCCTGTGACCCGa               <  1:169989/57‑1 (MQ=255)
                                                                 gcggcATTAAAGGCAACAAACTGCAAATCGTAAAATATGACGATGCCTGTGACCCGa               <  1:2687608/57‑1 (MQ=255)
                                                                  cggcATTAAAGGCAACAAACTGCAAATCGTAAAATATGACGATTCCTGTGACCCGAAACAGGCGGTTGc  <  1:1328322/69‑1 (MQ=255)
                                                                  cggcATTAAAGGCAACAAACTGCAAATCGTAAAATATGACGATGCCTGTGACCCGAAACAGGCGGTTGc  <  1:1757572/69‑1 (MQ=255)
                                                                  cggcATTAAAGGCAACAAACTGCAAATCGTAAAATATGACGATGCCTGTGACCCGAAACAGGCGGTTGc  <  1:674206/69‑1 (MQ=255)
                                                                  cggcATTAAAGGCAACAAACTGCAAATCGTAAAATATGACGATGCCTGTGACCCGAAACAGGCGGTTGc  <  1:2683963/69‑1 (MQ=255)
                                                                  cggcATTAAAGGCAACAAACTGCAAATCGTAAAATATGACGATGCCTGTGACCCGAAACAGGCGGTTGc  <  1:3203513/69‑1 (MQ=255)
                                                                  cggcATTAAAGGCAACAAACTGCAAATCGTAAAATATGACGATGCCTGTGACCCGAAACAGGCGGTTGc  <  1:1449484/69‑1 (MQ=255)
                                                                  cggcATTAAAGGCAACAAACTGCAAATCGTAAAATATGACGATGCCTGTGACCCGAAACAGGCGGTTGc  <  1:2287639/69‑1 (MQ=255)
                                                                  cggcATTAAAGGCAACAAACTGCAAATCGTAAAATATGACGATGCCTGTGACCCGAAACAGGCGGTTGc  <  1:847537/69‑1 (MQ=255)
                                                                  |                                                                    
TGCGCAGTACGGTGACCAGGAGTTTACCGGCGCAGAGCAGGCGGTTGCGGATATCAACGCTAAAGGCGGCATTAAAGGCAACAAACTGCAAATCGTAAAATATGACGATGCCTGTGACCCGAAACAGGCGGTTGC  >  minE/2627644‑2627778

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: