Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,913,962 G→A 11.7% intergenic (‑45/+236) pyrL ← / ← yjgH pyrBI operon leader peptide/predicted mRNA endoribonuclease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,913,9620GA11.7% 183.7 / 12.8 68intergenic (‑45/+236)pyrL/yjgHpyrBI operon leader peptide/predicted mRNA endoribonuclease
Reads supporting (aligned to +/- strand):  ref base G (23/37);  new base A (4/4);  total (27/41)
Fisher's exact test for biased strand distribution p-value = 7.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.25e-01

ACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTAAAAGTCAAGGTATATG  >  minE/2913892‑2914032
                                                                      |                                                                      
aCAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAg                                                                        <  1:3146564/71‑1 (MQ=255)
  aaaaTGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAgcg                                                                      >  1:2879278/1‑71 (MQ=255)
  aaaaTGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAgcg                                                                      >  1:694756/1‑71 (MQ=255)
    aaTGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAgcggc                                                                    >  1:3573197/1‑71 (MQ=255)
     aTGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAgcggcg                                                                   >  1:2982647/1‑71 (MQ=255)
      tGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAACGgcgc                                                                  <  1:3664730/71‑1 (MQ=255)
        tCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCa                                                                 <  1:3054220/70‑1 (MQ=255)
          gAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAgc                                                                       <  1:2514719/62‑1 (MQ=255)
          gAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAg                                                               >  1:1499593/1‑70 (MQ=255)
          gAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAg                                                               <  1:3151015/70‑1 (MQ=255)
            acacacTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCAt                                                            >  1:2647384/1‑71 (MQ=255)
              acacTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCAttt                                                          <  1:3453244/71‑1 (MQ=255)
               cacTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTg                                                         <  1:1857983/71‑1 (MQ=255)
               cacTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTg                                                         <  1:3525306/71‑1 (MQ=255)
               cacTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTg                                                         <  1:582334/71‑1 (MQ=255)
               cacTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTg                                                         <  1:1602803/71‑1 (MQ=255)
                acTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCAttt                                                          <  1:2072378/69‑1 (MQ=255)
                  tGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATg                                                      >  1:3489867/1‑71 (MQ=255)
                     aCCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCa                                                    <  1:3206468/70‑1 (MQ=255)
                         tACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATAtt                                               <  1:1053919/71‑1 (MQ=255)
                         tACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATAtt                                               <  1:2606511/71‑1 (MQ=255)
                           cATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAACGGCGCAAGCATTTGATGCAATAttt                                              >  1:1324386/1‑70 (MQ=255)
                            aTCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGTCGCAAGCATTTGATGCAATAtttt                                             <  1:1267397/70‑1 (MQ=255)
                            aTCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGc                                                              <  1:33908/53‑1 (MQ=255)
                             tcctccCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGc                                                              <  1:1404568/52‑1 (MQ=255)
                              cctccCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCAt                                                            >  1:296942/1‑53 (MQ=255)
                              cctccCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAACGGCGCAAGCAt                                                            >  1:1894181/1‑53 (MQ=255)
                               ctccCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGGAATATTTTGt                                           >  1:3796353/1‑69 (MQ=255)
                               ctccCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATAtt                                               <  1:1561139/65‑1 (MQ=255)
                               ctccCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCa                                         <  1:2665298/71‑1 (MQ=255)
                                  ccGGCAAATTGTCCGGCATTATACTCATCGTCAGAAgc                                                                       >  1:3252011/1‑38 (MQ=255)
                                  ccGGCAAATTGTCCGGCATTATACTCATCGTCAGAAgc                                                                       >  1:2074968/1‑38 (MQ=255)
                                  ccGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCaa                                                                <  1:2298454/45‑1 (MQ=255)
                                  ccGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATAtttt                                             >  1:3744849/1‑64 (MQ=255)
                                  ccGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAg                                        <  1:3427694/69‑1 (MQ=255)
                                   cGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAgcg                                      >  1:3231001/1‑70 (MQ=255)
                                    ggCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAgcgc                                     >  1:1701879/1‑70 (MQ=255)
                                    ggCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCa                                    <  1:2963944/71‑1 (MQ=255)
                                    ggCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCa                                    <  1:1479097/71‑1 (MQ=255)
                                     gCAAATTGTCCTGCATTATACTCATCGTCAGAAGCGGCGCAAGCAttt                                                          <  1:1747888/48‑1 (MQ=255)
                                     gCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCAttt                                                          <  1:1948923/48‑1 (MQ=255)
                                     gCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCaa                                   <  1:1350167/71‑1 (MQ=255)
                                     gCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCa                                    >  1:1334736/1‑70 (MQ=255)
                                      cAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCaaa                                  <  1:124052/71‑1 (MQ=255)
                                          ttGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGt                              >  1:737975/1‑71 (MQ=255)
                                           tGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGtt                             >  1:1020872/1‑71 (MQ=255)
                                             tCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGttt                            <  1:2131560/70‑1 (MQ=255)
                                             tCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTGTCAGCGCAAACGGTTTAt                          <  1:2313813/71‑1 (MQ=255)
                                              ccGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCa                                    >  1:1291162/1‑61 (MQ=255)
                                               cGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTAt                          <  1:3822982/70‑1 (MQ=255)
                                                ggCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTAttt                        <  1:3472833/71‑1 (MQ=255)
                                                 gCATTATACTCATCGTCAGAAGCGGCGCAAGCAttt                                                          <  1:812857/36‑1 (MQ=255)
                                                 gCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGAt                                                       >  1:2089587/1‑39 (MQ=255)
                                                 gCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGAt                                                       >  1:2586340/1‑39 (MQ=255)
                                                 gCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTg                       <  1:2803884/71‑1 (MQ=255)
                                                 gCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTg                       <  1:2899928/71‑1 (MQ=255)
                                                      ataCTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACg                                >  1:894234/1‑57 (MQ=255)
                                                      ataCTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAAt                    >  1:62039/1‑69 (MQ=255)
                                                         cTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCa                                         <  1:2384699/45‑1 (MQ=255)
                                                         cTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGAGCAAACGGTTTATTTGAATTaaa                <  1:3164711/70‑1 (MQ=255)
                                                         cTCATCGTCAGAAACGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTaaa                <  1:3291289/70‑1 (MQ=255)
                                                         cTCATCGTCAGAAACGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTaaa                <  1:2783022/70‑1 (MQ=255)
                                                           cATCGTCAGAAACGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTAAAAg              >  1:3894242/1‑70 (MQ=255)
                                                            aTCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTaaaa               >  1:2074535/1‑68 (MQ=255)
                                                            aTCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTAAAAg              >  1:3762367/1‑69 (MQ=255)
                                                            aTCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTAAAAGt             >  1:1286338/1‑70 (MQ=255)
                                                            aTCGTCAGAAACGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTAAAAGt             >  1:1524397/1‑70 (MQ=255)
                                                             tCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTAAAAGTCa           <  1:211287/71‑1 (MQ=255)
                                                             tCGTCAGAAACGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTAAAAGTCa           <  1:1563525/71‑1 (MQ=255)
                                                              cGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTaaaa               >  1:2857055/1‑66 (MQ=255)
                                                               gTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTAAAAGTCAAg         >  1:3649359/1‑71 (MQ=255)
                                                                tCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTAt                          <  1:912296/53‑1 (MQ=255)
                                                                 cAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTAAAAGTCAAgg        >  1:3101462/1‑70 (MQ=255)
                                                                   gAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTAAAAGTCAAGGt       >  1:3348027/1‑69 (MQ=255)
                                                                    aaGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAAc                                 <  1:2868441/42‑1 (MQ=255)
                                                                    aaGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAAc                                 <  1:1212490/42‑1 (MQ=255)
                                                                      gcggcgCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTAAAAGTCAAGGTATATg  >  1:100914/1‑71 (MQ=255)
                                                                      |                                                                      
ACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTAAAAGTCAAGGTATATG  >  minE/2913892‑2914032

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: