Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,682,416 C→T 14.4% E432E (GAG→GAA yhgF ← predicted transcriptional accessory protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,682,4160CT14.4% 213.3 / 23.9 83E432E (GAG→GAAyhgFpredicted transcriptional accessory protein
Reads supporting (aligned to +/- strand):  ref base C (41/30);  new base T (8/4);  total (49/34)
Fisher's exact test for biased strand distribution p-value = 7.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

GCCAGCGGATCCTGCAAACGGCGGGCGATAGACACCGCGCCACGCAGCGAAACGTCGAGATCCGGGAACTCCTGTGCAGCCAGCTCGGAAGCCGAGTAAACCGACGCGCCTGCTTCGCTGACGATTACTTTCTGCGC  >  minE/2682348‑2682484
                                                                    |                                                                    
gCCAGCGGATCCTGCAAACGGCGGGCGATAGACACCGCGCCACGCAGCGAAACGTCGAGATCCGGGAACTc                                                                    <  1:4993876/71‑1 (MQ=255)
gCCAGCGGATCCTGCAAACGGCGGGCGATAGACACCGCGCCACGCAGCGAAACGTCGAGATCCGGGAACTc                                                                    <  1:629792/71‑1 (MQ=255)
gCCAGCGGATCCTGCAAACGGCGGGCGATAGACACCGCGCCACGCAGCGAAACGTCGAGATCCGGGAACTc                                                                    <  1:1883291/71‑1 (MQ=255)
gCCAGCGGATCCTGCAAACGGCAGGCGATAGACACCGCGCCACGCAGCGAAACGTCGAGATCCGGGAACTc                                                                    <  1:3701922/71‑1 (MQ=255)
       gATCCTGCAAACGGCGGGCGATAGACACCGCGCCACGCAGCGAAACGTCGAGATCCGGGAACTCCTGTGca                                                             <  1:2608787/71‑1 (MQ=255)
       gATCCTGCAAACGGCGGGCGATAGACACCGCGCCACGCAGCGAAACGTCGAGATCCGGGAACTCCTGTGca                                                             <  1:4375506/71‑1 (MQ=255)
       gATCCTGCAAACGGCGGGCGATAGACACCGCGCCACGCAGCGAAACGTCGAGATCCGGGAACTCCTGTGca                                                             <  1:4288474/71‑1 (MQ=255)
       gATCCTGCAAACGGCGGGCGATAGACACCGCGCCACGCAGCGAAACGTCGAGATCCGGGAACTCCTGTGca                                                             <  1:3855165/71‑1 (MQ=255)
           cTGCAAACGGCGGGCGATAGACACCGCGCCACGCAGCGAAACGTCGAGATCCGGGAACTCCTGTGcat                                                            <  1:4551221/68‑2 (MQ=255)
            tGCAAACGGCGGGCGATAGACACCGCGCCACGCAGCGAAACGTCGAGATCCGGGAACTcc                                                                   <  1:3384354/60‑1 (MQ=255)
             gCAAACGGCGGGCGATAGACACCGCGCCACGCAGCGAAACGTCGAGATCCGGGAACTCCTGTGcagccagc                                                       >  1:3874243/1‑71 (MQ=255)
             gCAAACGGCGGGCGATAGACACCGCGCCACGCAGCGAAACGTCGAGATCCGGGAACTCCTGTGcagccagc                                                       >  1:1422299/1‑71 (MQ=255)
             gCAAACGGCGGGCGATAGACACCGCGCCACGCAGCGAAACGTCGAGATCCGGGAACTCCTGTGcagccagc                                                       >  1:1285544/1‑71 (MQ=255)
               aaaCGGCGGGCGATAGACACCGCGCCACGCAGCGAAACGTCGAGATCCGGGAACTCCTGTGCAGCCAGCt                                                      >  1:4725042/1‑70 (MQ=255)
                       ggcgATAGACACCGCGCCACGCAGCGAAACGTCGAGATCCGGGAACTCCTGTGCAGCCAGCTCGGAAGcc                                              >  1:2407914/1‑70 (MQ=255)
                         cgATAGACACCGCGCCACGCAGCGAAACGTCGAGATCCGGGAACTCCTGTGCAGCCAGCTCGGAAGCCGAg                                           >  1:5190736/1‑71 (MQ=255)
                         cgATAGACACCGCGCCACGCAGCGAAACGTCGAGATCCGGGAACTCCTGTGCAGCCAGCTCGGAAGCCGAg                                           >  1:621769/1‑71 (MQ=255)
                           aTAGACACCGCGCCACGCAGCGAAACGTCGAGATCCGGGAACTCCTGTGCAGCCAGCTCGGAAGCCGAGt                                          <  1:1510798/70‑1 (MQ=255)
                           aTAGACACCGCGCCACGCAGCGAAACGTCGAGATCCGGGAACTCCTGTGCAGCCAGCTCGGAAGCCGAGt                                          <  1:1950570/70‑1 (MQ=255)
                           aTAGACACCGCGCCACGCAGCGAAACGTCGAGATCCGGGAACTCCTGTGCAGCCAGCTCGGAAGCCGAGTa                                         <  1:5055462/71‑1 (MQ=255)
                             aGACACCGCGCCACGCAGCGTAACGTCGAGATCCGGGAACTCCTGTGCAGCCAGCTCGGAAGCCGAGTa                                         >  1:2054650/1‑69 (MQ=255)
                             aGACACCGCGCCACGCAGCGAAACGTCGAGATCCGGGAACTCCTGTGCAGCCAGCTCGGAAGCCGAGTaaa                                       >  1:4538035/1‑71 (MQ=255)
                             aGACACCGCGCCACGCAGCGAAACGTCGAGATCCGGGAACTCCTGTGCAGCCAGCTCGGAAGCCGAGTa                                         >  1:3807170/1‑69 (MQ=255)
                               acacCGCGCCACGCAGCGAAACGTCGAGATCCGGGAACTCCTGTGCAGCCAGCTCGGAAGCCGAGTAAAcc                                     <  1:2372462/71‑1 (MQ=255)
                               acacCGCGCCACGCAGCGAAACGTCGAGATCCGGGAACTCCTGTGCAGCCAGCTCGGAAGCCGAGTAAAcc                                     <  1:3422261/71‑1 (MQ=255)
                               acacCGCGCCACGCAGCGAAACGTCGAGATCCGGGAACTCCTGTGCAGCCAGCTCGGAAGCCGAGTAAAcc                                     <  1:78670/71‑1 (MQ=255)
                                   cgcgCCACGCAGCGAAACGTCGAGATCCGGGAACTCCt                                                                  >  1:918809/1‑38 (MQ=255)
                                   cgcgCCACGCAGCGAAACGTCGAGATCCGGGAACTCCTGTGCAGCCAGCTCGGAAGc                                               >  1:781550/1‑57 (MQ=255)
                                   cgcgCCACGCAGCGAAACGTCGAGATCCGGGAACTCCTGTGCAGCCAGCTCGGAAGc                                               >  1:1153278/1‑57 (MQ=255)
                                   cgcgCCACGCAGCGAAACGTCGAGATCCGGGAACTCCTGTGCAGCCAGCTCGGAAGc                                               >  1:1307376/1‑57 (MQ=255)
                                   cgcgCCACGCAGCGAAACGTCGAGATCCGGGAACTCCTGTGCAGCCAGCTCGGAAGc                                               >  1:1369766/1‑57 (MQ=255)
                                   cgcgCCACGCAGCGAAACGTCGAGATCCGGGAACTCCTGTGCAGCCAGCTCGGAAGc                                               >  1:4558343/1‑57 (MQ=255)
                                   cgcgCCACGCAGCGAAACGTCGAGATCCGGGAACTCCTGTGCAGCCAGCTCGGAAGc                                               >  1:2557391/1‑57 (MQ=255)
                                   cgcgCCACGCAGCGAAACGTCGAGATCCGGGAACTCCTGTGCAGCCAGCTCGGAAGc                                               >  1:197291/1‑57 (MQ=255)
                                     cgcCACGCAGCGAAACGTCGAGATCCGGGAACTCCTGTGCAGCCAGCTCGGAAGCCGAGTAAACCGAcgcg                               <  1:2625832/71‑1 (MQ=255)
                                           gCAGCGAAACGTCGAGATCCGGGAACTCCTGTGCAGCCAGCTCg                                                    <  1:2864473/44‑1 (MQ=255)
                                           gCAGCGAAACGTCGAGATCCGGGAACTCCTGTGCAGCCAGCTCg                                                    <  1:3452871/44‑1 (MQ=255)
                                           gCAGCGAAACGTCGAGATCCGGGAACTCCTGTGCAGCCAGCTCGGAAGCCGAGTAAACCGa                                   >  1:729798/1‑61 (MQ=255)
                                           gCAGCGAAACGTCGAGATCCGGGAACTCCTGTGCAGCCAGCTCGGAAGCCGAGTAAACCGa                                   >  1:399412/1‑61 (MQ=255)
                                           gCAGCGAAACGTCGAGATCCGGAAATTCCTGTGCTGCCAGCTCGGa                                                  <  1:4808405/46‑1 (MQ=38)
                                           gCAGCGAAACGTCGAGATCCGGAAATTCCTGTGCTGCCAGCTCGGa                                                  <  1:2332545/46‑1 (MQ=38)
                                           gCAGCGAAACGTCGAGATCCGGAAATTCCTGTGCTGCCAGCTCGGa                                                  <  1:1251784/46‑1 (MQ=38)
                                           gCAGCGAAACGTCGAGATCCGGAAATTCCTGTGCTGCCAGCTCGGa                                                  <  1:3673476/46‑1 (MQ=38)
                                            cAGCGAAACGTCGAGATCCGGGAACTCCTGTGCAGCCAGCTCGGAAGCCGAGTAAACCGACGCGCCTGCtt                        >  1:4575444/1‑71 (MQ=255)
                                            cAGCGAAACGTCGAGATCCGGGAACTCCTGTGCAGCCAGCTCGGAAGCCGAGTAAACCGACGCGCCTGCtt                        >  1:557806/1‑71 (MQ=255)
                                             aGCGAAACGTCGAGATCCGGGAACTCCTGTGCAGCCAGCTCGGAAGCCGAGTAAACCGACGCGCCTGCTTc                       <  1:2755289/71‑1 (MQ=255)
                                              gCGAAACGTCGAGATCCGGGAACTCCTGTGCAGCCAGCTCGGAAGCCGAGTAAACCGACGCGCCTGCTTCg                      <  1:2126403/71‑1 (MQ=255)
                                              gCGAAACGTCGAGATCCGGGAACTCCTGTGCAGCCAGCTCGGAAGCCGAGTAAACCGACGCGCCTGCTTCg                      <  1:1465475/71‑1 (MQ=255)
                                                gAAACGTCGAGATCCGGGAACTCCTGTGCAGCCAGCTCGGAAGCCGAGTAAACCGACGCGCCTGCTTCGCt                    >  1:326606/1‑71 (MQ=255)
                                                gAAACGTCGAGATCCGGGAACTCCTGTGCAGCCAGCTCGGAAGCCGAGTAAACCGACGCGCCTGCTTCGCt                    >  1:5074013/1‑71 (MQ=255)
                                                 aaaCGTCGAGATCCGGGAACTCCTGTGCAGCCAGCTCGGAAGCCGAGTAAACCGACGCGCCTGCTTCGCt                    >  1:3593262/1‑70 (MQ=255)
                                                 aaaCGTCGAGATCCGGGAACTCCTGTGCAGCCAGCTCGGAAGCCGAGTAAACCGACGCGCCTGCTTCGCt                    >  1:421967/1‑70 (MQ=255)
                                                 aaaCGTCGAGATCCGGGAACTCCTGTGCAGCCAGCTCGGAAGCCGAGTAAACCGACGCGCCTGCTTCGCt                    >  1:1850092/1‑70 (MQ=255)
                                                 aaaCGTCGAGATCCGGGAACTCCTGTGCAGCCAGCTCGGAAGCCGAGTAAACCGACGCGCCTGCTTCGCt                    >  1:1597236/1‑70 (MQ=255)
                                                   aCGTCGAGATCCGGGAACTCCTGTGCAGCCAGCTCGGAAGCCGAGTAAACCGa                                   >  1:4936372/1‑53 (MQ=255)
                                                   aCGTCGAGATCCGGGAACTCCTGTGCAGCCAGCTCGGAAGCCGAGTAAACCGa                                   >  1:4904284/1‑53 (MQ=255)
                                                   aCGTCGAGATCCGGGAACTCCTGTGCAGCCAGCTCGGAAGCCGAGTAAACCGa                                   >  1:5064502/1‑53 (MQ=255)
                                                   aCGTCGAGATCCGGGAACTCCTGTGCAGCCAGCTCGGAAGCCGAGTAAACCGa                                   >  1:543794/1‑53 (MQ=255)
                                                   aCGTCGAGATCCGGGAACTCCTGTGCAGCCAGCTCGGAAGCCGAGTAAACCGa                                   >  1:970074/1‑53 (MQ=255)
                                                   aCGTCGAGATCCGGGAACTCCTGTGCAGCCAGCTCGGAAGCCGAGTAAACCGa                                   >  1:3829915/1‑53 (MQ=255)
                                                   aCGTCGAGATCCGGGAACTCCTGTGCAGCCAGCTCGGAAGCCGAGTAAACCGa                                   >  1:2728664/1‑53 (MQ=255)
                                                   aCGTCGAGATCCGGGAACTCCTGTGCAGCCAGCTCGGAAGCCGAGTAAACCGa                                   >  1:3382563/1‑53 (MQ=255)
                                                   aCGTCGAGATCCGGAAATTCCTGTGCTGCCAGCTCGGaa                                                 >  1:1725561/1‑39 (MQ=25)
                                                   aCGTCGAGATCCGGAAATTCCTGTGCTGCCAGCTCGGAAg                                                >  1:536324/1‑40 (MQ=37)
                                                   aCGTCGAGATCCGGAAATTCCTGTGCTGCCAGCTCGGAAg                                                >  1:1185155/1‑40 (MQ=37)
                                                   aCGTCGAGATCCGGAAATTCCTGTGCTGCCAGCTCGGAAg                                                >  1:1634558/1‑40 (MQ=37)
                                                   aCGTCGAGATCCGGAAATTCCTGTGCTGCCAGCTCGGAAg                                                >  1:1786186/1‑40 (MQ=37)
                                                   aCGTCGAGATCCGGAAATTCCTGTGCTGCCAGCTCGGAAg                                                >  1:2524166/1‑40 (MQ=37)
                                                   aCGTCGAGATCCGGAAATTCCTGTGCTGCCAGCTCGGAAg                                                >  1:2827550/1‑40 (MQ=37)
                                                   aCGTCGAGATCCGGAAATTCCTGTGCTGCCAGCTCGGAAg                                                >  1:4770040/1‑40 (MQ=37)
                                                          gaTCCGGGAACTCCTGTGCAGCCAGCTCGGAAGCCGAGTAAACCGACGCGCCTGCTTCGCTGACGATTACt          <  1:4782023/71‑1 (MQ=255)
                                                           aTCCGTGAACTCCTGTGCAGCCAGCTCGGAAGCCGAGTAAACCGACGCGCCTGCTTCGCTGACGATTACtt         <  1:3438560/71‑1 (MQ=255)
                                                           aTCCGGGAACTCCTGTGCAGCCAGCTCGGAAGCCGAGTAAACCGACGCGCCTGCTTCGCTGACGATTACt          <  1:4451179/70‑1 (MQ=255)
                                                           aTCCGGGAACTCCTGTGCAGCCAGCTCGGAAGCCGAGTAAACCGACGCGCCTGCTTCGCTGACGATTACt          <  1:548699/70‑1 (MQ=255)
                                                           aTCCGGGAACTCCTGTGCAGCCAGCTCGGAAGCCGAGTAAACCGACGCGCCTGCTTCGCTGACGATTACt          <  1:3990698/70‑1 (MQ=255)
                                                           aTCCGGGAACTCCTGTGCAGCCAGCTCGGAAGCCGAGTAAACCGACGCGCCTGCTTCGCTGACGATTACt          <  1:118950/70‑1 (MQ=255)
                                                           aTCCGGGAACTCCTGTGCAGCCAGCTCGGAAGCCGAGTAAACCGACGCGCCTGCTTCGCTGACGATTACt          <  1:4725795/70‑1 (MQ=255)
                                                            tCCGGTAACTCCTGTGCAGCCAGCTCGGAAGCCGAGTAAACCGACGCGCCTGCTTCGCTGACGATTACt          <  1:3129299/69‑1 (MQ=255)
                                                                  aaCTCCTGTGCAGCCAGCTCGGAAGCCGAGTAAACCGACGCGCCTGCTTCGCTGACGATTACTTTCTgcgc  >  1:3871161/1‑71 (MQ=255)
                                                                  aaCTCCTGTGCAGCCAGCTCGGAAGCCGAGTAAACCGACGCGCCTGCTTCGCTGACGATTACTTTCTgcgc  >  1:1982520/1‑71 (MQ=255)
                                                                  aaCTCCTGTGCAGCCAGCTCGGAAGCCGAGTAAACCGACGCGCCTGCTTCGCTGACGATTACTTTCTgcgc  >  1:690492/1‑71 (MQ=255)
                                                                   aCTCCTGTGCAGCCAGCTCGGAAGCCGAGTAAACCGACGCGCCTGCTTCGCTGACGATTACTTTCTgc    >  1:2935125/1‑68 (MQ=255)
                                                                   aCTCCTGTGCAGCCAGCTCGGAAGCCGAGTAAACCGACGCGCCTGCTTCGCTGACGATTACTTTCTgc    >  1:3151731/1‑68 (MQ=255)
                                                                    |                                                                    
GCCAGCGGATCCTGCAAACGGCGGGCGATAGACACCGCGCCACGCAGCGAAACGTCGAGATCCGGGAACTCCTGTGCAGCCAGCTCGGAAGCCGAGTAAACCGACGCGCCTGCTTCGCTGACGATTACTTTCTGCGC  >  minE/2682348‑2682484

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: