Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 96,401 T→C 19.3% E62E (GAA→GAG yacG ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE96,4010TC19.3% 111.2 / 17.6 77E62E (GAA→GAGyacGconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (24/38);  new base C (11/4);  total (35/42)
Fisher's exact test for biased strand distribution p-value = 2.12e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCACTGCTGGGGATTCGTTTTTCTTCAGCAGCC  >  minE/96337‑96467
                                                                 |                                                                  
gATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGttctt                                                                <  1:1213536/69‑1 (MQ=255)
gATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGttctt                                                                <  1:3188855/69‑1 (MQ=255)
   gCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGc                                                             <  1:801026/69‑1 (MQ=255)
   gCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGc                                                             <  1:1351010/69‑1 (MQ=255)
   gCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGc                                                             <  1:2807026/69‑1 (MQ=255)
   gCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGc                                                             <  1:571375/69‑1 (MQ=255)
   gCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGc                                                             >  1:379922/1‑69 (MQ=255)
   gCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGc                                                             <  1:354349/69‑1 (MQ=255)
          tGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTc                                                      <  1:2591389/69‑1 (MQ=255)
          tGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTc                                                      <  1:2164946/69‑1 (MQ=255)
                 gttttATCAGGCCAGCGTCGCATCAGATAGAATGTCACTGCTTTGGCTCTTCTCTCCAGTCATCGCttt                                              >  1:2708900/4‑69 (MQ=255)
                 gttttATCAGGCCAGCGTCGCATCAGATAGAATGTCACTGCTTTGGCTCTTCGCTCCAGTCATCGCttt                                              >  1:3331762/4‑69 (MQ=255)
                 gttttATCAGGCCAGCGTCGCATCAGATAGAATGTCACTGCTTTGGCTCTTCGCTCCAGTCATCGCttt                                              >  1:2765166/4‑69 (MQ=255)
                 gttttATCAGGCCAGCGTCGCATCAGATAGAATGTCACTGCTTTGGCTCTTCGCTCCAGTCATCGCttt                                              >  1:3122026/4‑69 (MQ=255)
                 gttttATCAGGCCAGCGTCGCATCAGATAGAATGTCACTGCTTTGGCTCTTCGCTCCAGTCATCGCttt                                              >  1:1779954/4‑69 (MQ=255)
                 gttttATCAGGCCAGCGTCGCATCAGATAGAATGTCACTGCTTTGGCTCTTCGCTCCAGTCATCGCttt                                              >  1:1718436/4‑69 (MQ=255)
                 gttttATCAGGCCAGCGTCGCATCAGATAGAATGTCACTGCTTTGGCTCTTCGCTCCAGTCATCGCttt                                              >  1:95025/4‑69 (MQ=255)
                 gttttATCAGGCCAGCGTCGCATCAGATAGAATGTCACTGCTTTGGCTCTTCGCTCCAGTCATCGCttt                                              >  1:1683262/4‑69 (MQ=255)
                 gttttATCAGGCCAGCGTCGCATCAGATAGAATGTCACTGCTTTGGCTCTTCGCTCCAGTCATCGCttt                                              >  1:90766/4‑69 (MQ=255)
                 gttttATCAGGCCAGCGTCGCATCAGATAGAATGTCACTGCTTTGGCTCTTCGCTCCAGTCATCGCttt                                              >  1:690778/4‑69 (MQ=255)
                 gttttATCAGGCCAGCGTCGCATCAGATAGAATGTCACTGCTTTGGCTCTTCGCTCCAGTCATCGCtt                                               >  1:3186642/4‑68 (MQ=255)
                    ttATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTc                                             >  1:3213804/1‑68 (MQ=255)
                    ttATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTc                                             >  1:1810512/1‑68 (MQ=255)
                    ttATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTc                                             >  1:3021363/1‑68 (MQ=255)
                    ttATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTc                                             >  1:1332660/1‑68 (MQ=255)
                         aGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAg                                       <  1:598360/69‑1 (MQ=255)
                         aGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAg                                       <  1:1381549/69‑1 (MQ=255)
                         aGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAg                                       <  1:2710715/69‑1 (MQ=255)
                          ggCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAg                                       <  1:1387124/68‑1 (MQ=255)
                           gCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTTCAGTCATCGCTTTCGGAAAGa                                      <  1:2595754/68‑1 (MQ=255)
                           gCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGa                                      <  1:953421/68‑1 (MQ=255)
                           gCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGa                                      <  1:2533516/68‑1 (MQ=255)
                           gCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGa                                      <  1:2865212/68‑1 (MQ=255)
                           gCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGa                                      <  1:2849763/68‑1 (MQ=255)
                           gCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGa                                      <  1:26114/68‑1 (MQ=255)
                           gCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGa                                      <  1:1930476/68‑1 (MQ=255)
                           gCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGa                                      <  1:2020797/68‑1 (MQ=255)
                           gCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGa                                      <  1:1638629/68‑1 (MQ=255)
                              tACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGc                                  <  1:1832667/69‑1 (MQ=255)
                              tACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGc                                  <  1:943669/69‑1 (MQ=255)
                              tACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGc                                  <  1:1487542/69‑1 (MQ=255)
                                cGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGc                                  <  1:3068425/67‑1 (MQ=255)
                                cGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGc                                  <  1:1779710/67‑1 (MQ=255)
                                cGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGc                                  <  1:1771821/67‑1 (MQ=255)
                                cGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGc                                  <  1:2775442/67‑1 (MQ=255)
                                cGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGc                                  <  1:952900/67‑1 (MQ=255)
                                cGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGc                                  <  1:1227634/67‑1 (MQ=255)
                                cGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCa                                >  1:3318904/1‑69 (MQ=255)
                                cGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCa                                >  1:1273635/1‑69 (MQ=255)
                                cGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCa                                >  1:2377729/1‑69 (MQ=255)
                                cGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCa                                >  1:1939544/1‑69 (MQ=255)
                                  tGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCAc                               <  1:1565024/68‑1 (MQ=255)
                                  tGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCAc                               <  1:1341671/68‑1 (MQ=255)
                                  tGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCAc                               <  1:958084/68‑1 (MQ=255)
                                  tGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCAc                               <  1:2029419/68‑1 (MQ=255)
                                        aTCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGa                                      <  1:1371958/55‑1 (MQ=255)
                                              tagaaTGTCACTGCTTTGGCTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCTACTGCTGGGGAt                     <  1:2532670/64‑1 (MQ=255)
                                              tagaaTGTCACTGCTTTGGCTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCTACTGCTGGGGAt                     <  1:2527406/64‑1 (MQ=255)
                                              tagaaTGTCACTGCTTTGGCTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCTACTGCTGGGGAt                     <  1:606488/64‑1 (MQ=255)
                                              cAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCActg                             >  1:3413677/1‑58 (MQ=255)
                                              cAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCACTGCTGGGGATTCg                  >  1:316700/1‑69 (MQ=255)
                                              cAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCACTGCTGGGGATTCg                  >  1:834692/1‑69 (MQ=255)
                                              cAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCACTGCTGGGGATTCg                  >  1:3040619/1‑69 (MQ=255)
                                              cAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCACTGCTGGGGATTCg                  >  1:3029163/1‑69 (MQ=255)
                                              cAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCACTGCTGGGGATTCg                  >  1:3028254/1‑69 (MQ=255)
                                                     tCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCACTGCTGGGGATTCGTTTttc            <  1:3341218/69‑1 (MQ=255)
                                                     tCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCACTGCTGGGGATTCGTTTttc            <  1:2132340/69‑1 (MQ=255)
                                                        cTGCTTTGGCTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCTACTGCTGGGGAt                     <  1:1629906/57‑1 (MQ=255)
                                                          gCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCActg                             >  1:2068440/1‑47 (MQ=255)
                                                          gCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCACTGCTGGGGATTCGTTTttcttc         >  1:795229/1‑67 (MQ=255)
                                                          gCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCACTGCTGGGGATTCGTTTttcttc         >  1:2419652/1‑67 (MQ=255)
                                                          gCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCACTGCTGGGGATTCGTTTttcttc         >  1:1004233/1‑67 (MQ=255)
                                                          gCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCACTGCTGGGGATTCGTTTttcttc         >  1:1298946/1‑67 (MQ=255)
                                                          gCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCACTGCTGGGGATTCGTTTttcttc         >  1:2087564/1‑67 (MQ=255)
                                                          gCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCACTGCTGGGGATTCGTTTttcttc         >  1:499951/1‑67 (MQ=255)
                                                          gCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCACTGCTGGGGATTCGTTTttcttc         >  1:2072446/1‑67 (MQ=255)
                                                          gCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCACTGCTGGGGATTCGTTTttctt          >  1:2703215/1‑66 (MQ=255)
                                                                gTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCACTGCTGGGGATTCGTTTTTCTTCAGCAGcc  <  1:1998053/68‑1 (MQ=255)
                                                                gTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCACTGCTGGGGATTCGTTTTTCTTCAGCAGcc  <  1:2726020/68‑1 (MQ=255)
                                                                 |                                                                  
GATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCACTGCTGGGGATTCGTTTTTCTTCAGCAGCC  >  minE/96337‑96467

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