Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 102,699 T→A 60.7% G229G (GGA→GGT nadC ← quinolinate phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE102,6990TA60.7% 32.9 / 55.6 62G229G (GGA→GGTnadCquinolinate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base T (12/12);  new base A (20/17);  total (32/30)
Fisher's exact test for biased strand distribution p-value = 7.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.55e-01

GCGTTTGACGGCTTCGCGCATCTGTTCTGTTTCGAAGTTATCCAGCATGATGATATCGGCTCCTGCTTTCAGGGCTTCATCAAGTTCTTCCAGATTCTCTACTTCGACTTCTACTGGCGCA  >  minE/102639‑102759
                                                            |                                                            
gCGTTTGACGGCTTCGCGCATCTGTTCTGTTTCGAAGTTATCCAGCATGATGATATCGGCTCCTGCttt                                                      <  1:399259/69‑1 (MQ=255)
gCGTTTGACGGCTTCGCGCATCTGTTCTGTTTCGAAGTTATCCAGCATGATGATATCGGCCCCTGCttt                                                      <  1:547297/69‑1 (MQ=255)
     tGACGGCTTCGCGCATCTGTTCTGTTTCGAAGTTATCCAGCATGATGATATCGGCTCCTGCTTTCAggg                                                 >  1:2107739/1‑69 (MQ=255)
     tGACGGCTTCGCGCATCTGTTCTGTTTCGAAGTTATCCAGCATGATGATATCGGCTCCTGCTTTCAggg                                                 >  1:1092775/1‑69 (MQ=255)
     tGACGGCTTCGCGCATCTGTTCTGTTTCGAAGTTATCCAGCATGATGATATCGGCTCCTGCTTTCAggg                                                 >  1:643251/1‑69 (MQ=255)
     tGACGGCTTCGCGCATCTGTTCTGTTTCGAAGTTATCCAGCATGATGATATCGGCTCCTGCTTTCAggg                                                 >  1:1407089/1‑69 (MQ=255)
         ggCTTCGCGCATCTGTTCTGTTTCGAAGTTATCCAGCATGATGATATCGGCACCTGCTTTCAGAGc                                                >  1:142527/1‑66 (MQ=255)
         ggCTTCGCGCATCTGTTCTGTTTCGAAGTTATCCAGCATGATGATATCGGCACCTGCTTTCAGAGc                                                >  1:3042588/1‑66 (MQ=255)
         ggCTTCGCGCATCTGTTCTGTTTCGAAGTTATCCAGCATGATGATATCGGCACCTGCTTTCAGAGc                                                >  1:1794645/1‑66 (MQ=255)
         ggCTTCGCGCATCTGTTCTGTTTCGAAGTTATCCAGCATGATGATATCGGCACCTGCTTTCAGAGc                                                >  1:264304/1‑66 (MQ=255)
         ggCTTCGCGCATCTGTTCTGTTTCGAAGTTATCCAGCATGATGATATCGGCACCTGCTTTCAGAGCt                                               >  1:3330817/1‑67 (MQ=255)
         ggCTTCGCGCATCTGTTCTGTTTCGAAGTTATCCAGCATGATGATATCGGCACCTGCTTTCAGAGCt                                               >  1:2950090/1‑67 (MQ=255)
         ggCTTCGCGCATCTGTTCTGTTTCGAAGTTATCCAGCATGATGATATCGGCACCTGCTTTCAGAGCt                                               >  1:2912752/1‑67 (MQ=255)
         ggCTTCGCGCATCTGTTCTGTTTCGAAGTTATCCAGCATGATGATATCGGCACCTGCTTTCAGAGCt                                               >  1:2499158/1‑67 (MQ=255)
         ggCTTCGCGCATCTGTTCTGTTTCGAAGTTATCCAGCATGATGATATCGGCACCTGCTTTCAGAGCt                                               >  1:2334136/1‑67 (MQ=255)
         ggCTTCGCGCATCTGTTCTGTTTCGAAGTTATCCAGCATGATGATATCGGCACCTGCTTTCAGAGCt                                               >  1:2859802/1‑67 (MQ=255)
         ggCTTCGCGCATCTGTTCTGTTTCGAAGTTATCCAGCATGATGATATCGGCACCTGCTTTCAGAGCt                                               >  1:3033016/1‑67 (MQ=255)
         ggCTTCGCGCATCTGTTCTGTTTCGAAGTTATCCAGCATGATGATATCGGCACCTGCTTTCAGAGCt                                               >  1:2285695/1‑67 (MQ=255)
         ggCTTCGCGCATCTGTTCTGTTTCGAAGTTATCCAGCATGATGATATCGGCACCTGCTTTCAGAGCt                                               >  1:2631896/1‑67 (MQ=255)
         ggCTTCGCGCATCTGTTCTGTTTCGAAGTTATCCAGCATGATGATATCGGCACCTGCTTTCAGAGCt                                               >  1:3301647/1‑67 (MQ=255)
         ggCTTCGCGCATCTGTTCTGTTTCGAAGTTATCCAGCATGATGATATCGGCACCTGCTTTCAGAGCt                                               >  1:3311228/1‑67 (MQ=255)
         ggCTTCGCGCATCTGTTCTGTTTCGAAGTTATCCAGCATGATGATATCGGCACCTGCTTTCAGAGCt                                               >  1:1401816/1‑67 (MQ=255)
         ggCTTCGCGCATCTGTTCTGTTTCGAAGTTATCCAGCATGATGATATCGGCACCTGCTTTCAGAGCt                                               >  1:1368616/1‑67 (MQ=255)
         ggCTTCGCGCATCTGTTCTGTTTCGAAGTTATCCAGCATGATGATATCGGCACCTGCTTTCAGAGCt                                               >  1:802112/1‑67 (MQ=255)
         ggCTTCGCGCATCTGTTCTGTTTCGAAGTTATCCAGCATGATGATATCGGCACCTGCTTTCAGAGCt                                               >  1:1291782/1‑67 (MQ=255)
         ggCTTCGCGCATCTGTTCTGTTTCGAAGTTATCCAGCATGATGATATCGGCACCTGCTTTCAGAGCt                                               >  1:842958/1‑67 (MQ=255)
          gCTTCGCGCATCTGTTCTGTTTCGAAGTTATCCAGCATGATGATATCGGCTCCTGCTTTCAGGGCTtc                                             >  1:339175/1‑68 (MQ=255)
          gCTTCGCGCATCTGTTCTGTTTCGAAGTTATCCAGCATGATGATATCGGCTCCTGCTTTCAGGGCTtc                                             >  1:1006572/1‑68 (MQ=255)
                               tCGAAGTTATCCAGCATGATGATATCGGCTCCTGCTTTCAGGGCTTCATCAAGTTCTTCCAGATTctct                       >  1:3260822/1‑69 (MQ=255)
                               tCGAAGTTATCCAGCATGATGATATCGGCTCCTGCTTTCAGGGCTTCATCAAGTTCTTCCAGATTctct                       >  1:1036351/1‑69 (MQ=255)
                                                 atgatATCGGCTCCTGCTTTCAGGGCTTCATCAAGttcttc                                 >  1:2073767/1‑41 (MQ=255)
                                                 atgatATCGGCTCCTGCTTTCAGGGCTTCATCAAGTTCTTCCAGATTCTCTACTTCGACTTCTACTgg      >  1:207933/1‑68 (MQ=255)
                                                 atgatATCGGCTCCTGCTTTCAGGGCTTCATCAAGTTCTTCCAGATTCTCTACTTCGACTTCTACTgg      >  1:3010607/1‑68 (MQ=255)
                                                  tgatATCGGCTCCTGCTTTCAGGGCTTCATCAAGTTCTTCCAGATtctc                        <  1:2633350/49‑1 (MQ=255)
                                                  tgatATCGGCTCCTGCTTTCAGGGCTTCATCAAGTTCTTCCAGATtctc                        <  1:3182138/49‑1 (MQ=255)
                                                  tgatATCGGCTCCTGCTTTCAGGGCTTCATCAAGTTCTTCCAGATtctc                        <  1:22194/49‑1 (MQ=255)
                                                  tgatATCGGCTCCTGCTTTCAGGGCTTCATCAAGTTCTTCCAGATtctc                        <  1:2880085/49‑1 (MQ=255)
                                                  tgatATCGGCTCCTGCTTTCAGGGCTTCATCAAGTTCTTCCAGATtctc                        <  1:2677466/49‑1 (MQ=255)
                                                  tgatATCGGCTCCTGCTTTCAGGGCTTCATCAAGTTCTTACAGATtctc                        <  1:1045994/49‑1 (MQ=255)
                                                  tgatATCGGCACCTGCTTTCAGAGCTTCATCAAGTTCTTCCAGATTCTCTACTTCGACTTCt           <  1:164415/62‑1 (MQ=255)
                                                  tgatATCGGCACCTGCTTTCAGAGCTTCATCAAGTTCTTCCAGATTCTCTACTTCGACTTCt           <  1:1728830/62‑1 (MQ=255)
                                                  tgatATCGGCACCTGCTTTCAGAGCTTCATCAAGTTCTTCCAGATTCTCTACTTCGACTTCt           <  1:351694/62‑1 (MQ=255)
                                                  tgatATCGGCACCTGCTTTCAGAGCTTCATCAAGTTCTTCCAGATTCTCTACTTCGACTTCt           <  1:1634697/62‑1 (MQ=255)
                                                  tgatATCGGCACCTGCTTTCAGAGCTTCATCAAGTTCTTCCAGATTCTCTACTTCGACTTCt           <  1:3422210/62‑1 (MQ=255)
                                                  tgatATCGGCACCTGCTTTCAGAGCTTCATCAAGTTCTTCCAGATTCTCTACTTCGACTTCt           <  1:689815/62‑1 (MQ=255)
                                                  tgatATCGGCACCTGCTTTCAGAGCTTCATCAAGTTCTTCCAGATTCTCTACTTCGACTTCt           <  1:174601/62‑1 (MQ=255)
                                                  tgatATCGGCACCTGCTTTCAGAGCTTCATCAAGTTCTTCCAGATTCTCTACTTCGACTTCt           <  1:2003641/62‑1 (MQ=255)
                                                  tgatATCGGCACCTGCTTTCAGAGCTTCATCAAGTTCTTCCAGATTCTCTACTTCGACTTCt           <  1:2632374/62‑1 (MQ=255)
                                                  tgatATCGGCACCTGCTTTCAGAGCTTCATCAAGTTCTTCCAGATTCTCTACTTCGACTTCt           <  1:3194332/62‑1 (MQ=255)
                                                  tgatATCGGCACCTGCTTTCAGAGCTTCATCAAGTTCTTCCAGATTCTCTACTTCGACTTCt           <  1:3127646/62‑1 (MQ=255)
                                                  tgatATCGGCACCTGCTTTCAGAGCTTCATCAAGTTCTTCCAGATTCTCTACTTCGACTTCt           <  1:3105546/62‑1 (MQ=255)
                                                  tgatATCGGCACCTGCTTTCAGAGCTTCATCAAGTTCTTCCAGATTCTCTACTTCGACTTCt           <  1:2417955/62‑1 (MQ=255)
                                                  tgatATCGGCACCTGCTTTCAGAGCTTCATCAAGTTCTTCCAGATTCTCTACTTCGACTTCt           <  1:2500065/62‑1 (MQ=255)
                                                  tgatATCGGCACCTGCTTTCAGAGCTTCATCAAGTTCTTCCAGATTCTCTACTTCGACTTCt           <  1:2835951/62‑1 (MQ=255)
                                                  tgatATAGGCACCTGCTTTCAGAGCTTCATCAAGTTCTTCCAGATTCTCTACTTCGACTTCt           <  1:2302890/62‑1 (MQ=255)
                                                   gatATCGGCTCCTGCTTTCAGGGCTTCATCAAGTTCTTCCAGATTCTCTACTTCGACTTCt           >  1:2482899/1‑61 (MQ=255)
                                                    atatCGGCTCCTGCTTTCAGGGCTTCATCAAGTTCTTCCAGATTCTCTACTTCGACTTCTACTGGCGCa  <  1:897189/69‑1 (MQ=255)
                                                    atatCGGCACCTGCTTTCAGAGCTTCATCAAGTTCTTCCAGATTCTCTACTTCGACTTCt           <  1:1872950/60‑1 (MQ=255)
                                                     tatCGGCTCCTGCTTTCAGGGCTTCATCAAGTTCTTCCAGATTCTCTACTTCGACTTCTACTGGCGCa  <  1:861357/68‑1 (MQ=255)
                                                        cGGCTCCTGCTTTCAGGGCTTCATCAAGTTCTTCCAGATTCTCTACTTCg                 <  1:2844792/50‑1 (MQ=255)
                                                         ggCTCCTGCTTTCAGGGCTTCATCAAGTTCTTCCAGATTCTCTa                      <  1:2870147/44‑1 (MQ=255)
                                                         ggCTCCTGCTTTCAGGGCTTCATCAAGTTCTTCCAGATTCTCTa                      <  1:129749/44‑1 (MQ=255)
                                                            |                                                            
GCGTTTGACGGCTTCGCGCATCTGTTCTGTTTCGAAGTTATCCAGCATGATGATATCGGCTCCTGCTTTCAGGGCTTCATCAAGTTCTTCCAGATTCTCTACTTCGACTTCTACTGGCGCA  >  minE/102639‑102759

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: