Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 214,931 A→C 57.1% E14A (GAA→GCA)  dnaQ → DNA polymerase III epsilon subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE214,9310AC57.1% 361.1 / inf 331E14A (GAA→GCA) dnaQDNA polymerase III epsilon subunit
Reads supporting (aligned to +/- strand):  ref base A (63/79);  new base C (102/87);  total (165/166)
Fisher's exact test for biased strand distribution p-value = 9.60e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.66e-01

Only 100 of 333 total aligned reads displayed.

AAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGTTGA  >  minE/214862‑215000
                                                                     |                                                                     
     aCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAaccac                                                                 <  1:5302775/71‑1 (MQ=255)
        aGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCg                                                               >  1:1760818/1‑70 (MQ=255)
            tGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTAt                                                           >  1:2277968/1‑70 (MQ=255)
              aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGaa                                                        >  1:6243736/1‑71 (MQ=255)
              aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGa                                                         <  1:3623842/70‑1 (MQ=255)
              aCATAAATGACCGCTATGAGCACTGCAATTACACGCAAGATCGTTCTCGATACCGAAACCACCGGTATGaa                                                        >  1:2660815/1‑71 (MQ=255)
               cATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                       >  1:1284326/1‑71 (MQ=255)
               cATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                       >  1:467753/1‑71 (MQ=255)
               cATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                       >  1:3350769/1‑71 (MQ=255)
                 tAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                       <  1:5091762/69‑1 (MQ=255)
                   aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAAc                                                       >  1:422980/1‑67 (MQ=255)
                   aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAAc                                                       >  1:5031423/1‑67 (MQ=255)
                   aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGa                                                   <  1:4474094/71‑1 (MQ=255)
                   aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGa                                                   <  1:4879302/71‑1 (MQ=255)
                    aTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAt                                                  >  1:6195944/1‑71 (MQ=255)
                    aTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGAt                                                  <  1:5176730/71‑1 (MQ=255)
                         cGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTgg                                               <  1:2995022/69‑1 (MQ=255)
                          gCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTgc                                            <  1:5257032/71‑1 (MQ=255)
                            tATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAac                                                                    >  1:353300/1‑45 (MQ=255)
                            tATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAac                                                                    >  1:1693821/1‑45 (MQ=255)
                              tGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCAc                                        <  1:1867848/71‑1 (MQ=255)
                                aGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTa                                      >  1:6937215/1‑71 (MQ=255)
                                  cACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTAt                                     >  1:4121131/1‑70 (MQ=255)
                                  cACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTAt                                     >  1:1424228/1‑70 (MQ=255)
                                   aCTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                                    >  1:2484309/1‑70 (MQ=255)
                                   aCTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGa                                   <  1:1674339/71‑1 (MQ=255)
                                   aCTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGa                                   <  1:1449202/71‑1 (MQ=255)
                                    cTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                              <  1:5275577/43‑1 (MQ=255)
                                    cTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                              <  1:3419066/43‑1 (MQ=255)
                                    cTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTa                                      >  1:4664484/1‑67 (MQ=255)
                                    cTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa                                  <  1:5546549/71‑1 (MQ=255)
                                    cTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa                                  <  1:6361685/71‑1 (MQ=255)
                                     tGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                                    <  1:6719332/68‑1 (MQ=255)
                                     tGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                                    <  1:2171014/68‑1 (MQ=255)
                                     tGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                                    <  1:1985177/68‑1 (MQ=255)
                                      gCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAgg                                >  1:1172556/1‑71 (MQ=255)
                                      gCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAgg                                >  1:1578277/1‑71 (MQ=255)
                                      gCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAg                                 >  1:6547467/1‑70 (MQ=255)
                                       cAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCa                                         >  1:2112452/1‑61 (MQ=255)
                                       cAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCa                                         >  1:3388916/1‑61 (MQ=255)
                                       cAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa                                  >  1:2160587/1‑68 (MQ=255)
                                       cAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa                                  >  1:4009000/1‑68 (MQ=255)
                                        aaTTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGAtt                                                 <  1:4725021/52‑1 (MQ=255)
                                        aaTTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGAtt                                                 <  1:3885035/52‑1 (MQ=255)
                                          ttACACGCCAGATCGTTCTAGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCca                             >  1:6833464/1‑70 (MQ=255)
                                           tACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                                    <  1:1128819/62‑1 (MQ=255)
                                           tACACGCCAGATCGTTCTCGATACCGAAACCACCGGTa                                                            >  1:6006336/1‑38 (MQ=255)
                                            acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAg                                 >  1:962961/1‑64 (MQ=255)
                                            acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                            >  1:3137606/1‑69 (MQ=255)
                                            acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                            >  1:4019956/1‑69 (MQ=255)
                                            acacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                            >  1:6760623/1‑69 (MQ=255)
                                            acacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                            >  1:858903/1‑69 (MQ=255)
                                            acacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                            >  1:4553423/1‑69 (MQ=255)
                                            acacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                            >  1:3820012/1‑69 (MQ=255)
                                             cacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACaa                          <  1:3456723/70‑1 (MQ=255)
                                              acGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAtt                                                 <  1:508090/46‑1 (MQ=39)
                                              acGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAtt                                                 <  1:230236/46‑1 (MQ=39)
                                              acGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGa                        <  1:4866128/71‑1 (MQ=255)
                                              acGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGATGCGCACTATGAAGGCCACAAGa                        <  1:5667855/71‑1 (MQ=255)
                                               cGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGAt                       <  1:4398809/71‑1 (MQ=255)
                                                 ccAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCa                     <  1:495730/71‑1 (MQ=255)
                                                 ccAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCa                     >  1:6143119/1‑71 (MQ=255)
                                                  cAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAt                    >  1:903116/1‑71 (MQ=255)
                                                   aGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                   >  1:427007/1‑71 (MQ=255)
                                                    gATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTAt                                     <  1:5386154/52‑1 (MQ=255)
                                                    gATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTAt                                     <  1:3674739/52‑1 (MQ=255)
                                                     aTCGTTCTCGATAGCGAAACCACCGGTTTGAACGAGATTGGTGCGCACTAt                                     <  1:2546508/51‑1 (MQ=38)
                                                     aTCGTTCTCGATACCGAAACCACCGGTATGAACCa                                                     <  1:604195/35‑1 (MQ=255)
                                                     aTCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTAt                                     <  1:2590450/51‑1 (MQ=255)
                                                     aTCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTAt                                     <  1:395846/51‑1 (MQ=255)
                                                     aTCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTAt                                     <  1:3698915/51‑1 (MQ=255)
                                                      tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACt                                       >  1:6408858/1‑48 (MQ=39)
                                                      tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACt                                       >  1:1006947/1‑48 (MQ=39)
                                                      tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACt                                       >  1:2750517/1‑48 (MQ=39)
                                                      tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                   >  1:271982/1‑68 (MQ=255)
                                                      tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                   >  1:2022672/1‑68 (MQ=255)
                                                       cGTTCTCGCTACCGCAACCACCTGTATGAACCAGATTGGTGCGCACTATg                                    <  1:162799/50‑1 (MQ=255)
                                                       cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                   >  1:6589243/1‑67 (MQ=255)
                                                       cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                   >  1:5837338/1‑67 (MQ=255)
                                                       cGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGGGCACTATGa                                   >  1:3199190/1‑51 (MQ=255)
                                                         ttCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAtt             <  1:6914204/71‑1 (MQ=255)
                                                         ttCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAtt             <  1:2931095/71‑1 (MQ=255)
                                                         ttCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAg                         <  1:1332545/59‑1 (MQ=255)
                                                          tctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTg            >  1:4823061/1‑71 (MQ=255)
                                                           ctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTgg           >  1:2821388/1‑71 (MQ=255)
                                                           ctcGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTgg           >  1:5998602/1‑71 (MQ=255)
                                                           ctcGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTg            >  1:5964623/1‑70 (MQ=255)
                                                           ctcGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTg            >  1:730312/1‑70 (MQ=255)
                                                             cgctACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                   <  1:351659/58‑1 (MQ=255)
                                                             cgctACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga                 >  1:2357062/4‑63 (MQ=255)
                                                               ctACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGcc       >  1:6454845/2‑71 (MQ=255)
                                                               aTACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                                    <  1:6630410/42‑1 (MQ=255)
                                                               aTACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                                    <  1:1479543/42‑1 (MQ=255)
                                                                 aCCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGt     >  1:3576803/1‑71 (MQ=255)
                                                                  ccGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGtt    >  1:4289792/1‑71 (MQ=255)
                                                                  ccGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGtt    >  1:4332702/1‑71 (MQ=255)
                                                                  ccGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGtt    <  1:5522146/71‑1 (MQ=255)
                                                                  ccGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGt     <  1:5328999/70‑1 (MQ=255)
                                                                   cGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGTTg   <  1:4595883/71‑1 (MQ=255)
                                                                    gAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGTTg   >  1:5510795/1‑70 (MQ=255)
                                                                     |                                                                     
AAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGTTGA  >  minE/214862‑215000

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: