Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 102,711 G→A 65.0% A225A (GCC→GCT nadC ← quinolinate phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE102,7110GA65.0% 34.5 / 31.3 57A225A (GCC→GCTnadCquinolinate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base G (9/11);  new base A (20/17);  total (29/28)
Fisher's exact test for biased strand distribution p-value = 5.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.79e-01

TGACGGCTTCGCGCATCTGTTCTGTTTCGAAGTTATCCAGCATGATGATATCGGCTCCTGCTTTCAGGGCTTCATCAAGTTCTTCCAGATTCTCTACTTCGACTTCTACTGGCGCATCCGGGTGCA  >  minE/102644‑102769
                                                                   |                                                          
tGACGGCTTCGCGCATCTGTTCTGTTTCGAAGTTATCCAGCATGATGATATCGGCTCCTGCTTTCAggg                                                           >  1:1092775/1‑69 (MQ=255)
tGACGGCTTCGCGCATCTGTTCTGTTTCGAAGTTATCCAGCATGATGATATCGGCTCCTGCTTTCAggg                                                           >  1:1407089/1‑69 (MQ=255)
tGACGGCTTCGCGCATCTGTTCTGTTTCGAAGTTATCCAGCATGATGATATCGGCTCCTGCTTTCAggg                                                           >  1:643251/1‑69 (MQ=255)
tGACGGCTTCGCGCATCTGTTCTGTTTCGAAGTTATCCAGCATGATGATATCGGCTCCTGCTTTCAggg                                                           >  1:2107739/1‑69 (MQ=255)
    ggCTTCGCGCATCTGTTCTGTTTCGAAGTTATCCAGCATGATGATATCGGCACCTGCTTTCAGAGc                                                          >  1:3042588/1‑66 (MQ=255)
    ggCTTCGCGCATCTGTTCTGTTTCGAAGTTATCCAGCATGATGATATCGGCACCTGCTTTCAGAGc                                                          >  1:142527/1‑66 (MQ=255)
    ggCTTCGCGCATCTGTTCTGTTTCGAAGTTATCCAGCATGATGATATCGGCACCTGCTTTCAGAGc                                                          >  1:1794645/1‑66 (MQ=255)
    ggCTTCGCGCATCTGTTCTGTTTCGAAGTTATCCAGCATGATGATATCGGCACCTGCTTTCAGAGc                                                          >  1:264304/1‑66 (MQ=255)
    ggCTTCGCGCATCTGTTCTGTTTCGAAGTTATCCAGCATGATGATATCGGCACCTGCTTTCAGAGCt                                                         >  1:3311228/1‑67 (MQ=255)
    ggCTTCGCGCATCTGTTCTGTTTCGAAGTTATCCAGCATGATGATATCGGCACCTGCTTTCAGAGCt                                                         >  1:2499158/1‑67 (MQ=255)
    ggCTTCGCGCATCTGTTCTGTTTCGAAGTTATCCAGCATGATGATATCGGCACCTGCTTTCAGAGCt                                                         >  1:2950090/1‑67 (MQ=255)
    ggCTTCGCGCATCTGTTCTGTTTCGAAGTTATCCAGCATGATGATATCGGCACCTGCTTTCAGAGCt                                                         >  1:3033016/1‑67 (MQ=255)
    ggCTTCGCGCATCTGTTCTGTTTCGAAGTTATCCAGCATGATGATATCGGCACCTGCTTTCAGAGCt                                                         >  1:2334136/1‑67 (MQ=255)
    ggCTTCGCGCATCTGTTCTGTTTCGAAGTTATCCAGCATGATGATATCGGCACCTGCTTTCAGAGCt                                                         >  1:2285695/1‑67 (MQ=255)
    ggCTTCGCGCATCTGTTCTGTTTCGAAGTTATCCAGCATGATGATATCGGCACCTGCTTTCAGAGCt                                                         >  1:2631896/1‑67 (MQ=255)
    ggCTTCGCGCATCTGTTCTGTTTCGAAGTTATCCAGCATGATGATATCGGCACCTGCTTTCAGAGCt                                                         >  1:2912752/1‑67 (MQ=255)
    ggCTTCGCGCATCTGTTCTGTTTCGAAGTTATCCAGCATGATGATATCGGCACCTGCTTTCAGAGCt                                                         >  1:3301647/1‑67 (MQ=255)
    ggCTTCGCGCATCTGTTCTGTTTCGAAGTTATCCAGCATGATGATATCGGCACCTGCTTTCAGAGCt                                                         >  1:3330817/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)
    ggCTTCGCGCATCTGTTCTGTTTCGAAGTTATCCAGCATGATGATATCGGCACCTGCTTTCAGAGCt                                                         >  1:2859802/1‑67 (MQ=255)
     gCTTCGCGCATCTGTTCTGTTTCGAAGTTATCCAGCATGATGATATCGGCTCCTGCTTTCAGGGCTtc                                                       >  1:1006572/1‑68 (MQ=255)
     gCTTCGCGCATCTGTTCTGTTTCGAAGTTATCCAGCATGATGATATCGGCTCCTGCTTTCAGGGCTtc                                                       >  1:339175/1‑68 (MQ=255)
                          tCGAAGTTATCCAGCATGATGATATCGGCTCCTGCTTTCAGGGCTTCATCAAGTTCTTCCAGATTctct                                 >  1:1036351/1‑69 (MQ=255)
                          tCGAAGTTATCCAGCATGATGATATCGGCTCCTGCTTTCAGGGCTTCATCAAGTTCTTCCAGATTctct                                 >  1:3260822/1‑69 (MQ=255)
                                            atgatATCGGCTCCTGCTTTCAGGGCTTCATCAAGttcttc                                           >  1:2073767/1‑41 (MQ=255)
                                            atgatATCGGCTCCTGCTTTCAGGGCTTCATCAAGTTCTTCCAGATTCTCTACTTCGACTTCTACTgg                >  1:3010607/1‑68 (MQ=255)
                                            atgatATCGGCTCCTGCTTTCAGGGCTTCATCAAGTTCTTCCAGATTCTCTACTTCGACTTCTACTgg                >  1:207933/1‑68 (MQ=255)
                                             tgatATCGGCTCCTGCTTTCAGGGCTTCATCAAGTTCTTCCAGATtctc                                  <  1:3182138/49‑1 (MQ=255)
                                             tgatATCGGCTCCTGCTTTCAGGGCTTCATCAAGTTCTTCCAGATtctc                                  <  1:2880085/49‑1 (MQ=255)
                                             tgatATCGGCTCCTGCTTTCAGGGCTTCATCAAGTTCTTCCAGATtctc                                  <  1:2633350/49‑1 (MQ=255)
                                             tgatATCGGCTCCTGCTTTCAGGGCTTCATCAAGTTCTTCCAGATtctc                                  <  1:2677466/49‑1 (MQ=255)
                                             tgatATCGGCTCCTGCTTTCAGGGCTTCATCAAGTTCTTCCAGATtctc                                  <  1:22194/49‑1 (MQ=255)
                                             tgatATCGGCTCCTGCTTTCAGGGCTTCATCAAGTTCTTACAGATtctc                                  <  1:1045994/49‑1 (MQ=255)
                                             tgatATCGGCACCTGCTTTCAGAGCTTCATCAAGTTCTTCCAGATTCTCTACTTCGACTTCt                     <  1:3105546/62‑1 (MQ=255)
                                             tgatATCGGCACCTGCTTTCAGAGCTTCATCAAGTTCTTCCAGATTCTCTACTTCGACTTCt                     <  1:351694/62‑1 (MQ=255)
                                             tgatATCGGCACCTGCTTTCAGAGCTTCATCAAGTTCTTCCAGATTCTCTACTTCGACTTCt                     <  1:3422210/62‑1 (MQ=255)
                                             tgatATCGGCACCTGCTTTCAGAGCTTCATCAAGTTCTTCCAGATTCTCTACTTCGACTTCt                     <  1:164415/62‑1 (MQ=255)
                                             tgatATCGGCACCTGCTTTCAGAGCTTCATCAAGTTCTTCCAGATTCTCTACTTCGACTTCt                     <  1:1728830/62‑1 (MQ=255)
                                             tgatATCGGCACCTGCTTTCAGAGCTTCATCAAGTTCTTCCAGATTCTCTACTTCGACTTCt                     <  1:174601/62‑1 (MQ=255)
                                             tgatATCGGCACCTGCTTTCAGAGCTTCATCAAGTTCTTCCAGATTCTCTACTTCGACTTCt                     <  1:1634697/62‑1 (MQ=255)
                                             tgatATCGGCACCTGCTTTCAGAGCTTCATCAAGTTCTTCCAGATTCTCTACTTCGACTTCt                     <  1:2003641/62‑1 (MQ=255)
                                             tgatATCGGCACCTGCTTTCAGAGCTTCATCAAGTTCTTCCAGATTCTCTACTTCGACTTCt                     <  1:3194332/62‑1 (MQ=255)
                                             tgatATCGGCACCTGCTTTCAGAGCTTCATCAAGTTCTTCCAGATTCTCTACTTCGACTTCt                     <  1:3127646/62‑1 (MQ=255)
                                             tgatATCGGCACCTGCTTTCAGAGCTTCATCAAGTTCTTCCAGATTCTCTACTTCGACTTCt                     <  1:2417955/62‑1 (MQ=255)
                                             tgatATCGGCACCTGCTTTCAGAGCTTCATCAAGTTCTTCCAGATTCTCTACTTCGACTTCt                     <  1:2500065/62‑1 (MQ=255)
                                             tgatATCGGCACCTGCTTTCAGAGCTTCATCAAGTTCTTCCAGATTCTCTACTTCGACTTCt                     <  1:2632374/62‑1 (MQ=255)
                                             tgatATCGGCACCTGCTTTCAGAGCTTCATCAAGTTCTTCCAGATTCTCTACTTCGACTTCt                     <  1:689815/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:129749/44‑1 (MQ=255)
                                                    ggCTCCTGCTTTCAGGGCTTCATCAAGTTCTTCCAGATTCTCTa                                <  1:2870147/44‑1 (MQ=255)
                                                           gCTTTCAGGGCTTCATCAAGTTCTTCCAGATTCTCTACTTCGACTTCTACTGGCGCATCCGGGTGCa  >  1:265574/1‑67 (MQ=255)
                                                                   |                                                          
TGACGGCTTCGCGCATCTGTTCTGTTTCGAAGTTATCCAGCATGATGATATCGGCTCCTGCTTTCAGGGCTTCATCAAGTTCTTCCAGATTCTCTACTTCGACTTCTACTGGCGCATCCGGGTGCA  >  minE/102644‑102769

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: