Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 351,286 A→T 31.3% G128G (GGT→GGA purK ← N5‑carboxyaminoimidazole ribonucleotide synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE351,2860AT31.3% 71.4 / 43.6 67G128G (GGT→GGApurKN5‑carboxyaminoimidazole ribonucleotide synthase
Reads supporting (aligned to +/- strand):  ref base A (11/35);  new base T (16/5);  total (27/40)
Fisher's exact test for biased strand distribution p-value = 1.02e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

TACATTCGCCGTAACACTCTGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTGACCGCGACCGTCATAACCACCAGTGCGACGCTTAACAATCGCCAGCTCACCTAAACGATCAAACACCGCAGGC  >  minE/351219‑351352
                                                                   |                                                                  
tACATTCGCCGTAACACTCTGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTGACCGCGAc                                                                   >  1:1930832/1‑69 (MQ=255)
tACATTCGCCGTAACACTCTGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTGACCGCGAc                                                                   >  1:796941/1‑69 (MQ=255)
tACATTCGCCGTAACACTCTGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTGACCGCGAc                                                                   >  1:2451658/1‑69 (MQ=255)
tACATTCGCCGTAACACTCTGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTGACCGCGAc                                                                   >  1:2507179/1‑69 (MQ=255)
tACATTCGCCGTAACACTCTGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTGACCGCGAc                                                                   >  1:182839/1‑69 (MQ=255)
tACATTCGCCGTAACACTCTGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTGACCGCGAc                                                                   >  1:279492/1‑69 (MQ=255)
tACATTCGCCGTAACACTCTGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGACATTGACCGCGAc                                                                   >  1:2051044/1‑69 (MQ=255)
tACATTCGCCGTAACACTCTGCCGGTAACAGTTCGGTTTCATTTGCGCGTAAACGCCATTGACCGCGAc                                                                   >  1:2706188/1‑69 (MQ=255)
      cGCCGTAGCATTCAGCAGGTAACTGTTCGGTTTCATCGGCGCGTAAACGCCATTGACCGCGTCCGTCAt                                                             <  1:2757272/69‑1 (MQ=255)
      cGCCGTAGCATTCAGCAGGTAACTGTTCGGTTTCATCGGCGCGTAAACGCCATTGACCGCGTCCGTCAt                                                             <  1:2780012/69‑1 (MQ=255)
      cGCCGTAGCATTCAGCAGGTAACTGTTCGGTTTCATCGGCGCGTAAACGCCATTGACCGCGTCCGTCAt                                                             <  1:1596286/69‑1 (MQ=255)
      cGCCGTAGCATTCAGCAGGTAACTGTTCGGTTTCATCGGCGCGTAAACGCCATTGACCGCGTCCGTCAt                                                             <  1:3193566/69‑1 (MQ=255)
      cGCCGTAACACTCTGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTGACCGCGACCGTCAt                                                             <  1:3239612/69‑1 (MQ=255)
       gCCGTAGCATTCAGCAGGTAACTGTTCGGTTTCATCGGCGCGTAAACGCCATTGACCGCGTCCGTCAt                                                             <  1:2293477/68‑1 (MQ=255)
        ccGTAACACTCTGCCGGTAACTGTTCGGTTTCATTTTCGCGTAAACGCCATTGACCGCGACCGTCATa                                                            <  1:2532033/68‑1 (MQ=255)
        ccGTAACACTCTGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTGACCGCGACCGTCATa                                                            <  1:3134605/68‑1 (MQ=255)
        ccGTAACACTCTGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTGACCGCGACCGTCATa                                                            <  1:3224990/68‑1 (MQ=255)
        ccGTAACACTCTGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTGACCGCGACCGTCATa                                                            <  1:3251314/68‑1 (MQ=255)
        ccGTAACACTCTGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTGACCGCGACCGTCATa                                                            <  1:3397177/68‑1 (MQ=255)
        ccGTAACACTCTGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTGACCGCGACCGTCATa                                                            <  1:429586/68‑1 (MQ=255)
        ccGTAACACTCTGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTGACCGCGACCGTCATa                                                            <  1:461115/68‑1 (MQ=255)
        ccGTAACACTCTGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTGACCGCGACCGTCATa                                                            <  1:737158/68‑1 (MQ=255)
        ccGTAACACTCTGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTGACCGCGACCGTCATa                                                            <  1:806346/68‑1 (MQ=255)
        ccGTAACACTCTGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTGACCGCGACCGTCATa                                                            <  1:879667/68‑1 (MQ=255)
        ccGTAACACTCTGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTGACCGCGACCGTCATa                                                            <  1:2465022/68‑1 (MQ=255)
        ccGTAACACTCTGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTGACCGCGACCGTCATa                                                            <  1:1085674/68‑1 (MQ=255)
        ccGTAACACTCTGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTGACCGCGACCGTCATa                                                            <  1:139280/68‑1 (MQ=255)
        ccGTAACACTCTGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTGACCGCGACCGTCATa                                                            <  1:1432436/68‑1 (MQ=255)
        ccGTAACACTCTGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTGACCGCGACCGTCATa                                                            <  1:1643377/68‑1 (MQ=255)
        ccGTAACACTCTGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTGACCGCGACCGTCATa                                                            <  1:1894580/68‑1 (MQ=255)
        ccGTAACACTCTGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTGACCGCGACCGTCATa                                                            <  1:187543/68‑1 (MQ=255)
        ccGTAACACTCTGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGACATTGACCGCGACCGTCATa                                                            <  1:614350/68‑1 (MQ=255)
        ccGTAACACTCTGCCGGTAACTGTTCGGATTCATTTGCGCGTAAACGCCATTGACCGCGACCGTCATa                                                            <  1:219233/68‑1 (MQ=255)
        ccGTAACACTCTGCCGGTAACTGTTCGGATTCATTTGCGCGTAAACGCCATTGACCGCGACCGTCATa                                                            <  1:2657158/68‑1 (MQ=255)
         cGTAACACTCTGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTGACCGCGACCGTCATa                                                            <  1:1718491/67‑1 (MQ=255)
                ctctGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTGACCGCGACCGTCATAACCACCAGt                                                   >  1:2359285/1‑69 (MQ=255)
                  ctGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTGACCGCGACCGTCATAACCACCAGTGc                                                 <  1:2847139/69‑1 (MQ=255)
                  ctGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTGACCGCGACCGTCATAACCACCAGTGc                                                 <  1:401565/69‑1 (MQ=255)
                  ctGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTGACCGCGACCGTCATAACCACCAGTGc                                                 <  1:1329288/69‑1 (MQ=255)
                  ctGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTGACCGCGACCGTCATAACCACCAGTGc                                                 <  1:2381452/69‑1 (MQ=255)
                    gCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTGACCGCGACCGTCATAACCACCAGTGc                                                 <  1:106759/67‑1 (MQ=255)
                              gTTCGGTTTCATCGGCGCGTAAACGCCATTGACCGCGTCCGt                                                                >  1:1684451/1‑42 (MQ=37)
                              gTTCGGTTTCATCGGCGCGTAAACGCCATTGACCGCGTCCGTCATAGCCACCAGTGCGGCGGTTAAc                                       >  1:2506270/1‑67 (MQ=255)
                              gTTCGGTTTCATCGGCGCGTAAACGCCATTGACCGCGTCCGTCATAGCCACCAGTGCGGCGCTTaa                                        >  1:1562334/1‑66 (MQ=255)
                              gTTCGGTTTCATCGGCGCGTAAACGCCATTGACCGCGTCCGTCATAGCCACCAGTGCGGCGCTTAAc                                       >  1:2664765/1‑67 (MQ=255)
                              gTTCGGTTTCATCGGCGCGTAAACGCCATTGACCGCGTCCGTCATAGCCACCAGTGCGGCGCTTAAc                                       >  1:3140078/1‑67 (MQ=255)
                              gTTCGGTTTCATCGGCGCGTAAACGCCATTGACCGCGTCCGTCATAGCCACCAGTGCGGCGCTTAACa                                      >  1:1272807/1‑68 (MQ=255)
                              gTTCGGTTTCATCGGCGCGTAAACGCCATTGACCGCGTCCGTCATAGCCACCAGTGCGGCGCTTAACa                                      >  1:2195853/1‑68 (MQ=255)
                              gTTCGGTTTCATCGGCGCGTAAACGCCATTGACCGCGTCCGTCATAGCCACCAGTGCGGCGCTTAACa                                      >  1:1234713/1‑68 (MQ=255)
                              gTTCGGTTTCATCGGCGCGTAAACGCCATTGACCGCGTCCGTCATAGCCACCAGTGCGGCGCTTAACa                                      >  1:2122539/1‑68 (MQ=255)
                              gTTCGGTTTCATCGGCGCGTAAACGCCATTGACCGCGTCCGTCATAGCCACCAGTGCGGCGCTTAACa                                      >  1:1259342/1‑68 (MQ=255)
                              gTTCGGTTTCATCGGCGCGTAAACGCCATTGACCGCGTCCGTCATAGCCACCAGTGCGGCGCTTAACa                                      >  1:2503493/1‑68 (MQ=255)
                              gTTCGGTTTCATCGGCGCGTAAACGCCATTGACCGCGTCCGTCATAGCCACCAGTGCGGCGCTTAACa                                      >  1:3417717/1‑68 (MQ=255)
                              gTTCGGTTTCATCGGCGCGTAAACGCCATTGACCGCGTCCGTCATAGCCACCAGTGCGGCGCTTAACa                                      >  1:1917195/1‑68 (MQ=255)
                              gTTCGGTTTCATCGGCGCGTAAACGCCATTGACCGCGTCCGTCATAGCCACCAGTGCGGCGCTTAACa                                      >  1:1833971/1‑68 (MQ=255)
                              gTTCGGTTTCATCGGCGCGTAAACGCCATTGACCGCGTCCGTCATAGCCACCAGTGCGGCGCTTAACa                                      >  1:274864/1‑68 (MQ=255)
                              gTTCGGTTTCATCGGCGCGTAAACGCCATTGACCGCGTCCGTCATAGCCACCAGTGCGGCGCTTAACa                                      >  1:1829828/1‑68 (MQ=255)
                                   gTTTCATTTGCGCGTAAACGCCATTGACCGCGACCGTCATAACCACCAGTGc                                                 <  1:2550347/52‑1 (MQ=255)
                                    tttCATTTGCGCGTAAACGCCATTGACCGCGACCGTGATa                                                            <  1:742460/40‑1 (MQ=255)
                                              gcgTAAACGCCATTGACCGCGACCGTCATAACCACCAGTGCGACGCTTAACAATCGCCAGCTCACCTa                      <  1:865784/68‑1 (MQ=255)
                                                 tAAACGCCATTGACCGCGACCGTCATAACCACCAGTGCGACGCTTAACAATCGCCAGCTCACCTa                      <  1:2294527/65‑1 (MQ=255)
                                                        cATTGACCGCGACCGTCATAACCACCAGTGCGACGCTTAACAATCGCCAGCTCACCTAAACGATCAAac           >  1:440462/1‑69 (MQ=255)
                                                            gACCGCGACCGTCATAACCACCAGTGCGACGCTTAACAATCGCAAGCTCACCTAAACGATCAAacac         >  1:1250100/1‑67 (MQ=255)
                                                                 cgACCGTCATAACCACCAGTGCGACGCTTAACAATCGCCAGCTCACCTAAACGATCAAACACCGCAGGc  <  1:2413259/69‑1 (MQ=255)
                                                                 cgACCGTCATAACCACCAGTGCGACGCTTAACAATCGCCAGCTCACCTAAACGATCAAACACCGCAGGc  <  1:603868/69‑1 (MQ=255)
                                                                 cgACCGTCATAACCACCAGTGCGACGCTTAACAATCGCCAGCTCACCTAAACGATCAAACACCGCAGGc  <  1:115451/69‑1 (MQ=255)
                                                                 cgACCGTCATAACCACCAGTGCGACGCTTAACAATCGCCAGCTCACCTAAACGATCAAACACCGCAGGc  <  1:1761892/69‑1 (MQ=255)
                                                                   |                                                                  
TACATTCGCCGTAACACTCTGCCGGTAACTGTTCGGTTTCATTTGCGCGTAAACGCCATTGACCGCGACCGTCATAACCACCAGTGCGACGCTTAACAATCGCCAGCTCACCTAAACGATCAAACACCGCAGGC  >  minE/351219‑351352

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: