Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 27,925 T→C 49.2% F199F (TTT→TTC 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,9250TC49.2% ‑5.8 / 53.4 59F199F (TTT→TTCcarAcarbamoyl phosphate synthetase small subunit, glutamine amidotransferase
Reads supporting (aligned to +/- strand):  ref base T (8/22);  new base C (13/16);  total (21/38)
Fisher's exact test for biased strand distribution p-value = 1.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGGCCTGCCAGAAGCGAAAAAAGAAGACGAGCTGCCGTTCCACGTCGTGGCTTATGATTTTGGTGCCAAGCGCAACATCCTGCGGATGCTGGTGGATAGAGGCTGTCGCCTGACCATCGTTCCGG  >  minE/27864‑27989
                                                             |                                                                
gTGGCCTGCCAGAAGCGAAAAAAGAAGACGAGCTGCCGTTCCACGTCGTGGCTTATGATTTTGGTGCCa                                                           <  1:1497018/69‑1 (MQ=255)
gTGGCCTGCCAGAAGCGAAAAAAGAAGACGAGCTGCCGTTCCACGTCGTGGCTTATGATTTTGGTGCCa                                                           <  1:1258392/69‑1 (MQ=255)
   gCCTGCCAGAAGCGAAAAAAGAAGACGAGCTGCCGTTCCACGTCGTGGCTTATGATTTTGGTGCCAAg                                                         >  1:2198421/1‑68 (MQ=255)
   gCCTGCCAGAAGCGAAAAAAGAAGACGAGCTGCCGTTCCACGTCGTGGCTTATGATTTTGGTGCCAAg                                                         >  1:1514062/1‑68 (MQ=255)
   gCCTGCCAGAAGCGAAAAAAGAAGACGAGCTGCCGTTCCACGTCGTGGCTTATGATTTTGGTGCCAAg                                                         >  1:1318820/1‑68 (MQ=255)
   gCCTGCCAGAAGCGAAAAAAGAAGACGAGCTGCCGTTCCACGTCGTGGCTTATGATTTTGGTGCCAAg                                                         >  1:165973/1‑68 (MQ=255)
        ccAGAAGCGAAAAAAGAAGACGATCTGCCGTTCCACGTCGTGGCTTATGATTTTGGTGCCAAGCGCAAc                                                   >  1:2912213/1‑69 (MQ=255)
        ccAGAAGCGAAAAAAGAAGACGAGCTGCAGTTCCACGTCGTGGCTTATGATTTTGGTGCCAAGCGCaa                                                    >  1:1910957/1‑68 (MQ=255)
          aGAAGCGAAAAAAGAAGACGATCTGCCGTTCCACGTCGTGGCTTATGATTTTGGTGCCAAGCGCAACAt                                                 >  1:335807/1‑69 (MQ=255)
          aGAAGCGAAAAAAGAAGACGAGCTGCCGTTCCACGTCGTGGCTTATGATTTTGGTGCCAAGCGCAACAt                                                 <  1:3301528/69‑1 (MQ=255)
          aGAAGCGAAAAAAGAAGACGAGCTGCCGTTCCACGTCGTGGCTTATGATTTTGGTGCCAAGCGCAACAt                                                 <  1:652123/69‑1 (MQ=255)
          aGAAGCGAAAAAAGAAGACGAGCTGCCGTTCCACGTCGTGGCTTATGATTTTGGTGCCAAGCGCAACAt                                                 <  1:1954986/69‑1 (MQ=255)
          aGAAGCGAAAAAAGAAGACGAGCTGCCGTTCCACGTCGTGGCTTATGATTTTGGTGCCAAGCGCAACAt                                                 <  1:660742/69‑1 (MQ=255)
          aGAAGCGAAAAAAGAAGACGAGCTGCCGTTCCACGTCGTGGCTTATGATTTTGGTGCCAAGCGCAACAt                                                 <  1:2224209/69‑1 (MQ=255)
          aGAAGCGAAAAAAGAAGACGAGCTGCCGTTCCACGTCGTGGCTTATGATTTTGGTGCCAAGCGCAACAt                                                 >  1:820086/1‑69 (MQ=255)
                taaaaaaGAAGACGAGCTGCCGTTCCACGTCGTGGCTTATGATTTTGGTGCCAAGCGCAACATCCTGCg                                           <  1:511457/68‑1 (MQ=255)
                gAAAAAAGAAGACGAGCTGCCGTTCCACGTCGTGGCTTATGATTTTGGTGCCAAGCGCAACATCCTGCg                                           <  1:489108/69‑1 (MQ=255)
                gAAAAAAGAAGACGAGCTGCCGTTCCACGTCGTGGCTTATGATTTTGGTGCCAAGCGCAACATCCTGCg                                           <  1:1474538/69‑1 (MQ=255)
                gAAAAAAGAAGACGAGCTGCCGTTCCACGTCGTGGCTTATGATTTTGGTGCCAAGCGCAACATCCTGCg                                           <  1:385077/69‑1 (MQ=255)
                caaaaaaGAAGACGAGCTGCCGTTCCACGTCGTGGCTTATGATTTTGGTGCCAAGCGCAACATCCTGCg                                           <  1:2088484/68‑1 (MQ=255)
                    aaaGAAGACGAGCTGCCGTTCCACGTCGTGGCTTATGATTTTGGTGCCAAGCGCAACAt                                                 <  1:773221/59‑1 (MQ=255)
                          gaCGAGCTGCCGTTCCACGTCGTGGCTTATGATTTTGGTGCCAAGCTCAACATCCTGCg                                           <  1:2924259/59‑1 (MQ=255)
                                    cGTACCACGTCGTGGCCTACGATTTCGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGc                        >  1:905167/1‑68 (MQ=255)
                                    cGTACCACGTCGTGGCCTACGATTTCGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCt                       >  1:2080112/1‑69 (MQ=255)
                                    cGTACCACGTCGTGGCCTACGATTTCGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCt                       >  1:983571/1‑69 (MQ=255)
                                    cGTACCACGTCGTGGCCTACGATTTCGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCt                       >  1:2063700/1‑69 (MQ=255)
                                    cGTACCACGTCGTGGCCTACGATTTCGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCt                       >  1:2095878/1‑69 (MQ=255)
                                    cGTACCACGTCGTGGCCTACGATTTCGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCt                       >  1:2306908/1‑69 (MQ=255)
                                    cGTACCACGTCGTGGCCTACGATTTCGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCt                       >  1:2418182/1‑69 (MQ=255)
                                    cGTACCACGTCGTGGCCTACGATTTCGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCt                       >  1:3312210/1‑69 (MQ=255)
                                    cGTACCACGTCGTGGCCTACGATTTCGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCt                       >  1:381746/1‑69 (MQ=255)
                                    cGTACCACGTCGTGGCCTACGATTTCGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCt                       >  1:1226766/1‑69 (MQ=255)
                                    cGTACCACGTCGTGGCCTACGATTTCGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCt                       >  1:1257741/1‑69 (MQ=255)
                                    cGTACCACGTCGTGGCCTACGATTTCGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCt                       >  1:1341460/1‑69 (MQ=255)
                                    cGTACCACGTCGTGGCCTACGATTTCGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCt                       >  1:118354/1‑69 (MQ=255)
                                        ccACGTCGTGGCTTATGATTTTGGTGCCAAGCGCAACATCCTGCg                                           <  1:3232360/45‑1 (MQ=255)
                                        ccACGTCGTGGCTTATGATTTTGGTGCCAAGCGCAACATCCTGCg                                           <  1:1685783/45‑1 (MQ=255)
                                        ccACGTCGTGGCTTATGATTTTGGTGCCAAGCGCAACATCCTGCg                                           <  1:923358/45‑1 (MQ=255)
                                        ccACGTCGTGGCTTATGATTTTGGTGCCAAGCGCAACATCCTGCg                                           <  1:37631/45‑1 (MQ=255)
                                        ccACGTCGTGGCTTATGATTTTGGTGCCAAGCGCAACATCCTGCg                                           <  1:706895/45‑1 (MQ=255)
                                        ccACGTCGTGGCTTATGATTTTGGTGCCAAGCGCAACATCCTGCg                                           <  1:2763343/45‑1 (MQ=255)
                                         cACGTCGTGGCTTATGATTTTGGTGCCAAGCGCAACATCCTGCg                                           <  1:1153947/44‑1 (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:1172556/67‑1 (MQ=255)
                                                         aTTTCGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCTGTCGCCTGACCATCGTTcc    <  1:2585471/67‑1 (MQ=255)
                                                         aTTTCGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCTGTCGCCTGACCATCGTTcc    <  1:1255066/67‑1 (MQ=255)
                                                         aTTTCGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCTGTCGCCTGACCATCGTTcc    <  1:2855359/67‑1 (MQ=255)
                                                         aTTTCGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCTGTCGCCTGACCATCGTTcc    <  1:1277436/67‑1 (MQ=255)
                                                         aTTTCGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCTGTCGCCTGACCATCGTTcc    <  1:2621948/67‑1 (MQ=255)
                                                         aTTTCGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCTGTCGCCTGACCATCGTTcc    <  1:2570243/67‑1 (MQ=255)
                                                         aTTTCGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCTGTCGCCTGACCATCGTTcc    <  1:2481428/67‑1 (MQ=255)
                                                         aTTTCGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCTGTCGCCTGACCATCGTTcc    <  1:241167/67‑1 (MQ=255)
                                                         aTTTCGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCTGTCGCCTGACCATCGTTcc    <  1:1547208/67‑1 (MQ=255)
                                                         aTTTCGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCTGTCGCCTGACCATCGTTcc    <  1:2221872/67‑1 (MQ=255)
                                                         aTTTCGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCTGTCGCCTGACCATCGTTcc    <  1:167098/67‑1 (MQ=255)
                                                         aTTTCGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCTGTCGCCTGACCATCGTTcc    <  1:1791927/67‑1 (MQ=255)
                                                         aTTTCGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCTGTCGCCTGACCATCGTTcc    <  1:1978140/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:1066236/1‑68 (MQ=255)
                                                          tttcGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCTGTCGCCTGACCATCGTTCCgg  >  1:937336/1‑68 (MQ=255)
                                                          tttcGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCTGTCGCCTGACCATCGTTCCgg  >  1:747106/1‑68 (MQ=255)
                                                          tttcGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCTGTCGCCTGACCATCGTTCCgg  >  1:729380/1‑68 (MQ=255)
                                                          tttcGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCTGTCGCCTGACCATCGTTCCgg  >  1:738227/1‑68 (MQ=255)
                                                          tttcGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCTGTCGCCTGACCATCGTTCCgg  >  1:305675/1‑68 (MQ=255)
                                                          tttcGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCTGTCGCCTGACCATCGTTCCgg  >  1:1052809/1‑68 (MQ=255)
                                                          tttcGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCTGTCGCCTGACCATCGTTCCgg  >  1:1107778/1‑68 (MQ=255)
                                                          tttcGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCTGTCGCCTGACCATCGTTCCgg  >  1:1326253/1‑68 (MQ=255)
                                                          tttcGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCTGTCGCCTGACCATCGTTCCgg  >  1:1332106/1‑68 (MQ=255)
                                                          tttcGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCTGTCGCCTGACCATCGTTCCgg  >  1:1866263/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)
                                                          tttcGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCTGTCGCCTGACCATCGTTCCgg  >  1:2626025/1‑68 (MQ=255)
                                                          tttcGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCTGTCGCCTGACCATCGTTCCgg  >  1:2627558/1‑68 (MQ=255)
                                                          tttcGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCTGTCGCCTGACCATCGTTCCgg  >  1:2695621/1‑68 (MQ=255)
                                                          tttcGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCTGTCGCCTGACCATCGTTCCgg  >  1:3045109/1‑68 (MQ=255)
                                                          tttcGGTGCCAAGCGCAACATCCTGCGTATGCTGGTGGATAGAGGCTGTCGCCTGACCATCGTTCCgg  >  1:3420754/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)
                                                             |                                                                
GTGGCCTGCCAGAAGCGAAAAAAGAAGACGAGCTGCCGTTCCACGTCGTGGCTTATGATTTTGGTGCCAAGCGCAACATCCTGCGGATGCTGGTGGATAGAGGCTGTCGCCTGACCATCGTTCCGG  >  minE/27864‑27989

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: