Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 27,949 G→T 80.9% R207R (CGG→CGT carA → carbamoyl phosphate synthetase small subunit, glutamine amidotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE27,9490GT80.9% 114.8 / 23.4 68R207R (CGG→CGTcarAcarbamoyl phosphate synthetase small subunit, glutamine amidotransferase
Reads supporting (aligned to +/- strand):  ref base G (8/5);  new base T (37/18);  total (45/23)
Fisher's exact test for biased strand distribution p-value = 7.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

CGTTCCACGTCGTGGCTTATGATTTTGGTGCCAAGCGCAACATCCTGCGGATGCTGGTGGATAGAGGCTGTCGCCTGACCATCGTTCCGGCGCAAACTTCTGCGGAAGATGTG  >  minE/27900‑28012
                                                 |                                                               
cGTACCACGTCGTGGCCTACGATTTCGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGc                                               >  1:905167/1‑68 (MQ=255)
cGTACCACGTCGTGGCCTACGATTTCGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCt                                              >  1:1226766/1‑69 (MQ=255)
cGTACCACGTCGTGGCCTACGATTTCGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCt                                              >  1:1341460/1‑69 (MQ=255)
cGTACCACGTCGTGGCCTACGATTTCGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCt                                              >  1:381746/1‑69 (MQ=255)
cGTACCACGTCGTGGCCTACGATTTCGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCt                                              >  1:2306908/1‑69 (MQ=255)
cGTACCACGTCGTGGCCTACGATTTCGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCt                                              >  1:1257741/1‑69 (MQ=255)
cGTACCACGTCGTGGCCTACGATTTCGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCt                                              >  1:2418182/1‑69 (MQ=255)
cGTACCACGTCGTGGCCTACGATTTCGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCt                                              >  1:2095878/1‑69 (MQ=255)
cGTACCACGTCGTGGCCTACGATTTCGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCt                                              >  1:3312210/1‑69 (MQ=255)
cGTACCACGTCGTGGCCTACGATTTCGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCt                                              >  1:118354/1‑69 (MQ=255)
cGTACCACGTCGTGGCCTACGATTTCGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCt                                              >  1:2080112/1‑69 (MQ=255)
cGTACCACGTCGTGGCCTACGATTTCGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCt                                              >  1:2063700/1‑69 (MQ=255)
cGTACCACGTCGTGGCCTACGATTTCGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCt                                              >  1:983571/1‑69 (MQ=255)
            tGGCTTATGATTTTGGTGCCAAGCGCAACATCCTGCGGATGCTGGTGGATAGAGGCTGTCGCCTGAcc                                   <  1:1524151/68‑1 (MQ=255)
                     aTTTCGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCTGTCGCCTGACCATTGTTcc                           <  1:536223/67‑1 (MQ=255)
                     aTTTCGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCTGTCGCCTGACCATCGTTcc                           <  1:3159186/67‑1 (MQ=255)
                     aTTTCGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCTGTCGCCTGACCATCGTTcc                           <  1:167098/67‑1 (MQ=255)
                     aTTTCGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCTGTCGCCTGACCATCGTTcc                           <  1:2221872/67‑1 (MQ=255)
                     aTTTCGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCTGTCGCCTGACCATCGTTcc                           <  1:241167/67‑1 (MQ=255)
                     aTTTCGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCTGTCGCCTGACCATCGTTcc                           <  1:2855359/67‑1 (MQ=255)
                     aTTTCGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCTGTCGCCTGACCATCGTTcc                           <  1:2481428/67‑1 (MQ=255)
                     aTTTCGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCTGTCGCCTGACCATCGTTcc                           <  1:2570243/67‑1 (MQ=255)
                     aTTTCGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCTGTCGCCTGACCATCGTTcc                           <  1:2585471/67‑1 (MQ=255)
                     aTTTCGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCTGTCGCCTGACCATCGTTcc                           <  1:2621948/67‑1 (MQ=255)
                     aTTTCGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCTGTCGCCTGACCATCGTTcc                           <  1:1978140/67‑1 (MQ=255)
                     aTTTCGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCTGTCGCCTGACCATCGTTcc                           <  1:1791927/67‑1 (MQ=255)
                     aTTTCGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCTGTCGCCTGACCATCGTTcc                           <  1:1172556/67‑1 (MQ=255)
                     aTTTCGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCTGTCGCCTGACCATCGTTcc                           <  1:1547208/67‑1 (MQ=255)
                     aTTTCGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCTGTCGCCTGACCATCGTTcc                           <  1:1277436/67‑1 (MQ=255)
                     aTTTCGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCTGTCGCCTGACCATCGTTcc                           <  1:1255066/67‑1 (MQ=255)
                      tttcGGTGCCAAGCGCAACATCCTGCGTATGCTGGTTGATAGAGGCTGTCGCCTGACCATCGTTCCgg                         >  1:2035946/1‑68 (MQ=255)
                      tttcGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGAt                                                     >  1:1174494/1‑40 (MQ=255)
                      tttcGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCTGTCGCCTGACCATCGTTcc                           <  1:3174498/66‑1 (MQ=255)
                      tttcGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCTGTCGCCTGACCATCGTTCCgg                         >  1:2627558/1‑68 (MQ=255)
                      tttcGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCTGTCGCCTGACCATCGTTCCgg                         >  1:937336/1‑68 (MQ=255)
                      tttcGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCTGTCGCCTGACCATCGTTCCgg                         >  1:2695621/1‑68 (MQ=255)
                      tttcGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCTGTCGCCTGACCATCGTTCCgg                         >  1:747106/1‑68 (MQ=255)
                      tttcGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCTGTCGCCTGACCATCGTTCCgg                         >  1:738227/1‑68 (MQ=255)
                      tttcGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCTGTCGCCTGACCATCGTTCCgg                         >  1:729380/1‑68 (MQ=255)
                      tttcGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCTGTCGCCTGACCATCGTTCCgg                         >  1:3420754/1‑68 (MQ=255)
                      tttcGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCTGTCGCCTGACCATCGTTCCgg                         >  1:1052809/1‑68 (MQ=255)
                      tttcGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCTGTCGCCTGACCATCGTTCCgg                         >  1:3045109/1‑68 (MQ=255)
                      tttcGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCTGTCGCCTGACCATCGTTCCgg                         >  1:305675/1‑68 (MQ=255)
                      tttcGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCTGTCGCCTGACCATCGTTCCgg                         >  1:1866263/1‑68 (MQ=255)
                      tttcGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCTGTCGCCTGACCATCGTTCCgg                         >  1:1326253/1‑68 (MQ=255)
                      tttcGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCTGTCGCCTGACCATCGTTCCgg                         >  1:1332106/1‑68 (MQ=255)
                      tttcGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCTGTCGCCTGACCATCGTTCCgg                         >  1:1107778/1‑68 (MQ=255)
                      tttcGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCTGTCGCCTGACCATCGTTCCgg                         >  1:2626025/1‑68 (MQ=255)
                      tttcGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCTGTCGCCTGACCATCGTTCCgg                         >  1:1066236/1‑68 (MQ=255)
                      tttcGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCTGTCGCCTGACCATCGTTCCgg                         >  1:203827/1‑68 (MQ=255)
                      tttcGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCTGTCGCCTGACCATCGTTCCgg                         >  1:2042529/1‑68 (MQ=255)
                      tttcGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCTGTCGCCTGACCATCGTTCCgg                         >  1:2450550/1‑68 (MQ=255)
                      tttcGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCTGTCGCCTGACCATCGTTCAgg                         >  1:1669489/1‑68 (MQ=255)
                      tttcGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCTGTAGCCTGACCATCGTTCCgg                         >  1:987166/1‑68 (MQ=255)
                      tttcGGTGCCAAGAGCAACATCCTGCGTATGCTGGTGGATAGAGGCTGTCGCCTGACCATCGTTCCgg                         >  1:122607/1‑68 (MQ=255)
                                  gcgcAACATCCTGCGTATGCTGGTGGATAGAGGCTGTCGCCTGACCATCGTTcc                           <  1:2728613/54‑1 (MQ=255)
                                   cgcAACATCCTGCGGATGCTGGTGGATAGAGGCTGTCGCCTGACCATCGTTCCGGCGCAAACTTCTgt            >  1:1008976/1‑67 (MQ=255)
                                   cgcAACATCCTGCGGATGCTGGTGGATAGAGGCTGTCGCCTGACCATCGTTCCGGCGCAAACTTCTgt            >  1:201516/1‑67 (MQ=255)
                                   cgcAACATCCTGCGGATGCTGGTGGATAGAGGCTGTCGCCTGACCATCGTTCCGGCGCAAACTTCTgt            >  1:1451109/1‑67 (MQ=255)
                                   cgcAACATCCTGCGGATGCTGGTGGATAGAGGCTGTCGCCTGACCATCGTTCCGGCGCAAACTTCTgt            >  1:408104/1‑67 (MQ=255)
                                      aaCATCCTGCGGATGCTGGTGGATAGAGGCTGTCGCCTGACCATCGTTCCGGCGCAAACTTCTGCGGaa        <  1:1337934/69‑1 (MQ=255)
                                      aaCATCCTGCGGATGCTGGTGGATAGAGGCTGTCGCCTGACCATCGTTCCGGCGCAAACTTCTGCGGaa        <  1:432046/69‑1 (MQ=255)
                                      aaCATCCTGCGGATGCTGGTGGATAGAGGCTGTCGCCTGACCATCGTTCCGGCGCAAACTTCTGCGGaa        <  1:2547485/69‑1 (MQ=255)
                                       aCATCCTGCGGATGCTGGTGGATAGAGGCTGTCGCCTGACCATCGTTCCGGCGCAAACTTCTGCGGaa        <  1:3151604/68‑1 (MQ=255)
                                          tCCTGCGGATGCTGGTGGATAGAGGCTGTCGCCTGACCATCGTTCCGGCGCAAACTTCTGCGGAAGAtg    >  1:1237943/1‑69 (MQ=255)
                                          tCCTGCGGATGCTGGTGGATAGAGGCTGTCGCCTGACCATCGTTCCGGCGCAAACTTCTGCGGAAGAtg    >  1:685218/1‑69 (MQ=255)
                                            cTGCGGATGCTGGTGGATAGAGGCTGTCGCCTGACCATCGTTCCGGCGCAAACTTCTGCGGAAGAtgtg  >  1:699344/1‑69 (MQ=255)
                                            cTGCGGATGCTGGTGGATAGAGGCTGTCGCCTGACCATCGTTCCGGCGCAAACTTCTGCGGAAGAtgtg  >  1:3086635/1‑69 (MQ=255)
                                                 |                                                               
CGTTCCACGTCGTGGCTTATGATTTTGGTGCCAAGCGCAACATCCTGCGGATGCTGGTGGATAGAGGCTGTCGCCTGACCATCGTTCCGGCGCAAACTTCTGCGGAAGATGTG  >  minE/27900‑28012

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: