Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 222,967 A→G 11.0% intergenic (+22/‑158) lpcA → / → prfH D‑sedoheptulose 7‑phosphate isomerase/predicted peptide chain release factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE222,9670AG11.0% 206.3 / 14.2 82intergenic (+22/‑158)lpcA/prfHD‑sedoheptulose 7‑phosphate isomerase/predicted peptide chain release factor
Reads supporting (aligned to +/- strand):  ref base A (22/51);  new base G (5/4);  total (27/55)
Fisher's exact test for biased strand distribution p-value = 1.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.84e-01

GTGATCCATATCCTGATCCAGTTGATTGAAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTA  >  minE/222898‑223035
                                                                     |                                                                    
gTGATCCATATCCTGATCCAGTTGATTGAAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGAc                                                                     <  1:1581574/71‑1 (MQ=255)
 tGATCCATATCCTGATCCAGTTGATTGAAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACg                                                                    <  1:1548591/71‑1 (MQ=255)
   aTCCATATCCTGATCCAGTTGATTGAAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGtt                                                                  <  1:3919900/71‑1 (MQ=255)
   aTCCATATCCTGATCCAGTTGATTGAAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGtt                                                                  <  1:1213214/71‑1 (MQ=255)
    tCCATATCCTGATCCAGTTGATTGAAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGtt                                                                  <  1:5485903/70‑1 (MQ=255)
    tCCATATCCTGATCCAGTTGATTGAAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGtt                                                                  <  1:1113132/70‑1 (MQ=255)
    tCCATATCCTGATCCAGTTGATTGAAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGtt                                                                  >  1:4519338/1‑70 (MQ=255)
     ccATATCCTGATCCAGTTGATTGAAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGtt                                                                  <  1:4063601/69‑1 (MQ=255)
     ccATATCCTGATCCAGTTGATTGAAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGtt                                                                  <  1:4873823/69‑1 (MQ=255)
     ccATATCCTGATCCAGTTGATTGAAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGtt                                                                  <  1:2842043/69‑1 (MQ=255)
      cATATCCTGATCCAGTTGATTGAAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTAc                                                               <  1:2910442/71‑1 (MQ=255)
        tatCCTGATCCAGTTGATTGAAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGttt                                                                 <  1:4165898/67‑1 (MQ=255)
        tatCCTGATCCAGTTGATTGAAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGttt                                                                 <  1:4884608/67‑1 (MQ=255)
        tatCCTGATCCAGTTGATTGAAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTActc                                                             >  1:2535177/1‑69 (MQ=255)
          tCCTGATCCAGTTGATTGAAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTa                                                                >  1:4807687/1‑66 (MQ=255)
          tCCTGATCCAGTTGATTGAAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTAcaga                                                            >  1:1883630/1‑67 (MQ=255)
            cTGATCCAGTTGATTGAAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCa                                                         <  1:1607679/71‑1 (MQ=255)
            cTGATCCAGTTGATTGAAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGCATAAGACGTTTACGCCGCa                                                         <  1:767501/71‑1 (MQ=255)
                  cAGTTGATTGAAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGc                                                   <  1:2904642/71‑1 (MQ=255)
                  cAGTTGATTGAAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGc                                                   <  1:1044458/71‑1 (MQ=255)
                  cAGTTGATTGAAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGc                                                   <  1:1881737/71‑1 (MQ=255)
                  cAGTTGATTGAAAAAGAGATGGGTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGc                                                   <  1:3162118/71‑1 (MQ=255)
                    gTTGATTGAAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCa                                                         <  1:554089/63‑1 (MQ=255)
                    gTTGATTGAAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCAt                                                 <  1:3064617/71‑1 (MQ=255)
                    gTTGATTGAAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCAt                                                 <  1:2530753/71‑1 (MQ=255)
                     ttgattgaAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCTGCAtt                                                <  1:41700/71‑1 (MQ=255)
                     ttgattgaAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCAtt                                                <  1:570566/71‑1 (MQ=255)
                         ttgaAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGtt                                                                  <  1:2418716/49‑1 (MQ=255)
                         ttgaAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGtt                                                                  <  1:67338/49‑1 (MQ=255)
                         ttgaAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGtt                                                                  <  1:3322679/49‑1 (MQ=255)
                         ttgaAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGtt                                                                  <  1:4414501/49‑1 (MQ=255)
                             aaaaGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGct                                        <  1:1090087/71‑1 (MQ=255)
                             aaaaGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGct                                        <  1:3333644/71‑1 (MQ=255)
                             aaaaGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGct                                        >  1:155804/1‑71 (MQ=255)
                              aaaGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGtt                                                                  <  1:44626/44‑1 (MQ=255)
                              aaaGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGtt                                                                  <  1:846465/44‑1 (MQ=255)
                              aaaGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGtt                                                                  <  1:3634691/44‑1 (MQ=255)
                              aaaGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGtt                                                                  <  1:3771522/44‑1 (MQ=255)
                              aaaGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGctg                                       <  1:764669/71‑1 (MQ=255)
                               aaGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGtt                                                                  <  1:2893171/43‑1 (MQ=255)
                               aaGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGtt                                                                  <  1:1809685/43‑1 (MQ=255)
                               aaGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGtt                                                                  <  1:1649470/43‑1 (MQ=255)
                               aaGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGtt                                                                  <  1:1316004/43‑1 (MQ=255)
                               aaGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTa                                                                >  1:1095810/1‑45 (MQ=255)
                               aaGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTa                                                                >  1:1610853/1‑45 (MQ=255)
                                agagATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGctgat                                     >  1:1296964/1‑71 (MQ=255)
                                agagATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGctgat                                     >  1:3852053/1‑71 (MQ=255)
                                  agaTGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTgcg                                          >  1:1077818/1‑64 (MQ=255)
                                     tGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGctgatgcctg                                <  1:1196182/71‑1 (MQ=255)
                                       gTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTAcgccg                                                           >  1:659022/1‑42 (MQ=255)
                                       gTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCAttt                                               >  1:3970016/1‑54 (MQ=255)
                                           agtaagtCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCg                           <  1:3505789/70‑1 (MQ=255)
                                               agtCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCt                      <  1:1655151/71‑1 (MQ=255)
                                                  cTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGAc                   >  1:434380/1‑71 (MQ=255)
                                                      tgTAGGCCGGACAAGACGTTTA‑‑TCGCATCCGGCAGTCGTGTGcagatgcct                                 <  1:2467046/50‑1 (MQ=9)
                                                      tgTAGGCCGGACAAGACGTTTA‑‑TCGCATCCGGCAGTCGTGTGCAGATGCCTGATGCGACGCTAACGCGTCt             <  1:3991590/70‑1 (MQ=9)
                                                      tgTAGGCCGGACAAGACGTTTA‑‑TCGCATCCGGCAGTCGTGTGCAGATGCCTGATGCGACGCTAACGCGTCt             <  1:4325091/70‑1 (MQ=9)
                                                          ggCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTgcgc                                         >  1:1101933/1‑41 (MQ=255)
                                                          ggCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTa           >  1:4039942/1‑71 (MQ=255)
                                                          ggCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTa           <  1:3968490/71‑1 (MQ=255)
                                                          ggCCGGACAAGACGTTTA‑‑TCGCATCCGGCAGTCGTGTGCAGATGCCTGATGCGACGCTAACGCGTCTTAt          <  1:3118110/70‑1 (MQ=9)
                                                          ggCCGGACAAGACGTTTA‑‑TCGCATCCGGCAGTCGTGTGCAGATGCCTGATGCGACGCTAACGCGTCTTAt          <  1:4582333/70‑1 (MQ=9)
                                                          ggCCGGACAAGACGTTTA‑‑TCGCATCCGGCAGTCGTGTGCAGATGCCTGATGCGACGCTAACGCGTCTTAt          <  1:4683535/70‑1 (MQ=9)
                                                          ggCCGGACAAGACGTTTA‑‑TCGCATCCGGCAGTCGTGTGCAGATGCCTGATGCGACGCTAACGCGTCTTAt          >  1:2766014/1‑70 (MQ=9)
                                                          ggCCGGACAAGACGTTTA‑‑TCGAATCCGGCAGTCGTGTGCAGATGCCTGATGCGACGCTAACGCGTCTTAt          <  1:834604/70‑1 (MQ=2)
                                                           gCCGGATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTa           >  1:4619478/1‑70 (MQ=11)
                                                             cGGATAAGGCGTTTACGCCGCATCCGGCAGTCGTGCGCAGATGCCTGATGCGACGCTAACGCGt               <  1:2226868/64‑1 (MQ=14)
                                                              ggATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGcagatgcctgatgc                            <  1:422928/50‑1 (MQ=11)
                                                              ggATAAGGCGTTTATGCCGCATCCGGCAGTTATGCGcagatgcctgatgc                            <  1:1756742/50‑1 (MQ=11)
                                                              ggATAAGGCGTTTACGCCGCATCCGGCAGTCGTGCGCAGATGCCTGATGCGACGCTAACGCGTCTTAt          >  1:1063915/1‑68 (MQ=14)
                                                              ggACAAGACGTTTA‑‑TCGCATCCGGCAGTCGTGTGCAGATGCCTGATGCGACGCTAACGCGTCTTATc         <  1:3573976/67‑1 (MQ=9)
                                                               gATAAGGCGTTTACGCCGCATCCGGCAGTCGTGCGCAGATGCCTGATGCGACGCTAACGCGTCTTATCagg      >  1:4350560/1‑69 (MQ=11)
                                                               gATAAGGCGTTTACGCCGCATCCGGCAGTCGTGCGCAGATGCCTGATGCGACGCTAACGCGTCTTATCagg      >  1:5554197/1‑69 (MQ=11)
                                                               gATAAGGCGTTTACGCCGCATCCGGCAGTCGTGCGCAGATGCCTGATGCGACGCTAACGCGTCTTATCagg      >  1:824945/1‑69 (MQ=11)
                                                               gATAAGACGTTTACGCCGCATCCGGCATTTGTGCGctgatgcct                                 <  1:1701772/44‑1 (MQ=255)
                                                                aTAAGACGTTTACGCCGCATCCGGCATTTGTGCGctgatgcctga                               >  1:3616525/1‑45 (MQ=255)
                                                                aTAAGACGTTTACGCCGCATCCGGCATTTGTGCGctgatgcctga                               >  1:2609340/1‑45 (MQ=255)
                                                                aTAAGACGTTTACGCCGCATCCGGCATTTGTGCGctgatgcctga                               >  1:1805980/1‑45 (MQ=255)
                                                                 tAAGGCGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGcc    <  1:3814930/71‑1 (MQ=11)
                                                                 tAAGACGTTTACGCCGCATCCGGCATTTGTGCGctga                                      >  1:4290447/1‑37 (MQ=255)
                                                                   aGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTa  >  1:3666333/1‑71 (MQ=255)
                                                                   aGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTa  >  1:989136/1‑71 (MQ=255)
                                                                     |                                                                    
GTGATCCATATCCTGATCCAGTTGATTGAAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTA  >  minE/222898‑223035

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: