Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 214,925 A→C 56.9% D12A (GAT→GCT)  dnaQ → DNA polymerase III epsilon subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE214,9250AC56.9% 356.3 / inf 315D12A (GAT→GCT) dnaQDNA polymerase III epsilon subunit
Reads supporting (aligned to +/- strand):  ref base A (67/69);  new base C (92/87);  total (159/156)
Fisher's exact test for biased strand distribution p-value = 7.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 326 total aligned reads displayed.

GTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGC  >  minE/214859‑214994
                                                                  |                                                                     
gTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATAcc                                                                   >  1:2059331/1‑71 (MQ=255)
 tAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATAcc                                                                   >  1:3219561/1‑70 (MQ=255)
   aaaCGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGa                                                                 >  1:6150938/1‑70 (MQ=255)
     aCGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAac                                                              <  1:5998510/71‑1 (MQ=255)
         cAAGTCTGACATAAATGACCG‑‑ATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                        >  1:6464617/1‑71 (MQ=255)
            gTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACAGCAACCACCgg                                                        >  1:2816281/1‑70 (MQ=255)
               tGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTAt                                                     >  1:5262300/1‑70 (MQ=255)
                 aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGaa                                                  >  1:6243736/1‑71 (MQ=255)
                 aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGa                                                   <  1:3623842/70‑1 (MQ=255)
                  cATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                 >  1:6673372/1‑71 (MQ=255)
                   aTAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGGAc                                                 <  1:3985830/70‑1 (MQ=255)
                    tAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAcc                                                <  1:1288364/70‑1 (MQ=255)
                      aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGa                                                                       >  1:6376756/1‑45 (MQ=255)
                      aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGa                                                                       >  1:4181352/1‑45 (MQ=255)
                      aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGa                                             <  1:1816423/71‑1 (MQ=255)
                       aTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAt                                            <  1:4878024/71‑1 (MQ=255)
                       aTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGAt                                            <  1:5176730/71‑1 (MQ=255)
                            cGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTgg                                         <  1:2995022/69‑1 (MQ=255)
                              cTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCgca                                                                <  1:58431/44‑3 (MQ=255)
                              cTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTg                                       >  1:2421446/1‑69 (MQ=255)
                               tATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAac                                                              >  1:1693821/1‑45 (MQ=255)
                               tATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAac                                                              >  1:353300/1‑45 (MQ=255)
                                 ttaGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCAc                                  <  1:5316751/69‑1 (MQ=255)
                                 tGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAt                                            <  1:2251587/61‑1 (MQ=255)
                                 tGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTgcgc                                    <  1:2596130/69‑1 (MQ=255)
                                  gAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCAc                                  >  1:5524798/1‑70 (MQ=255)
                                   aGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGAt                                            >  1:4754058/1‑59 (MQ=255)
                                    gCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGTACCa                                               <  1:3514138/55‑1 (MQ=255)
                                     cACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTAt                               >  1:4129109/1‑70 (MQ=255)
                                     cACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTAt                               >  1:1424228/1‑70 (MQ=255)
                                      aCTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGa                                             >  1:6728776/1‑55 (MQ=255)
                                      aCTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGa                             <  1:1449202/71‑1 (MQ=255)
                                       cTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                        <  1:3581420/43‑1 (MQ=255)
                                       cTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                        <  1:3379848/43‑1 (MQ=255)
                                       cTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTa                                >  1:4302806/1‑67 (MQ=255)
                                       cTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGt                                        >  1:4911808/1‑59 (MQ=255)
                                        tGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                              <  1:760162/68‑1 (MQ=255)
                                         gCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                              <  1:2513759/67‑1 (MQ=255)
                                         gCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAgg                          >  1:5626736/1‑71 (MQ=255)
                                         gCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAgg                          >  1:6929653/1‑71 (MQ=255)
                                         gCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAgg                          >  1:1172556/1‑71 (MQ=255)
                                         gCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAgg                          >  1:6558965/1‑71 (MQ=255)
                                         gCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAgg                          >  1:1578277/1‑71 (MQ=255)
                                          cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                         >  1:5047376/1‑71 (MQ=255)
                                          cAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGa                                             >  1:1671603/1‑51 (MQ=255)
                                          cAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa                            >  1:4009000/1‑68 (MQ=255)
                                           aaTTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGAtt                                           <  1:2280700/52‑1 (MQ=255)
                                           aaTTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGAtt                                           <  1:4261479/52‑1 (MQ=255)
                                           aaTTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGAtt                                           <  1:1870972/52‑1 (MQ=255)
                                           aaTTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGAtt                                           <  1:3885035/52‑1 (MQ=255)
                                            aTTACACGCCAGATCGTTCTCGATACCGAAACCACCgg                                                        <  1:3444969/38‑1 (MQ=255)
                                            aTTACACGCCAGATCGTTCTCGATACCGAAACCACCgg                                                        <  1:2577516/38‑1 (MQ=255)
                                             ttACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                      <  1:2149608/71‑1 (MQ=255)
                                             ttACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCca                       >  1:4564754/1‑70 (MQ=255)
                                              tACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                              <  1:5096151/62‑1 (MQ=255)
                                              tACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                              <  1:1082158/62‑1 (MQ=255)
                                              tACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                              <  1:1128819/62‑1 (MQ=255)
                                              tACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                              <  1:2161899/62‑1 (MQ=255)
                                              tACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                              <  1:3681086/62‑1 (MQ=255)
                                              tACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                      >  1:2671495/1‑70 (MQ=255)
                                               acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTa                                                      >  1:4463372/1‑37 (MQ=37)
                                               acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCAc                                  >  1:3153739/1‑57 (MQ=255)
                                               acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAg                           >  1:962961/1‑64 (MQ=255)
                                               acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAg                           >  1:6125378/1‑64 (MQ=255)
                                               acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                      >  1:633670/1‑69 (MQ=255)
                                               acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                      >  1:4693590/1‑69 (MQ=255)
                                               acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                      >  1:4019956/1‑69 (MQ=255)
                                               acacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTAt                               >  1:4501811/1‑60 (MQ=255)
                                               acacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                      >  1:3820012/1‑69 (MQ=255)
                                                 acGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAtt                                           <  1:6860921/46‑1 (MQ=39)
                                                 acGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGa                  <  1:884430/71‑1 (MQ=255)
                                                  cGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGAt                 <  1:4398809/71‑1 (MQ=255)
                                                    ccAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCa               >  1:5298352/1‑71 (MQ=255)
                                                      aGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt             <  1:3357286/71‑1 (MQ=255)
                                                       gATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTAt                               <  1:5386154/52‑1 (MQ=255)
                                                       gATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTAt                               <  1:5464595/52‑1 (MQ=255)
                                                       gATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTAt                               <  1:508269/52‑1 (MQ=255)
                                                        aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAt                                            <  1:717095/38‑1 (MQ=37)
                                                        aTCGTTCTCGATACCGAAACCACCGGTATGAACCa                                               <  1:604195/35‑1 (MQ=255)
                                                        aTCGTTCTCGATACCGAAACCACCGGTATGAACCa                                               <  1:3709698/35‑1 (MQ=255)
                                                        aTCGTTCTCGATACCGAAACCACCGGTATGAACCa                                               <  1:4641915/35‑1 (MQ=255)
                                                        aTCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTAt                               <  1:6827372/51‑1 (MQ=255)
                                                         tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACt                                 >  1:2750517/1‑48 (MQ=39)
                                                         tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACt                                 >  1:4978049/1‑48 (MQ=39)
                                                         tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACt                                 >  1:6408858/1‑48 (MQ=39)
                                                         tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACt                                 >  1:1006947/1‑48 (MQ=39)
                                                         tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt             >  1:271982/1‑68 (MQ=255)
                                                         tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt             >  1:5682331/1‑68 (MQ=255)
                                                         tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt             >  1:2022672/1‑68 (MQ=255)
                                                          cGTTCTCGCTACCGCAACCACCTGTATGAACCAGATTGGTGCGCACTATg                              <  1:162799/50‑1 (MQ=255)
                                                          cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt             >  1:467303/1‑67 (MQ=255)
                                                          cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt             >  1:6591577/1‑67 (MQ=255)
                                                          cGTTCTCGATACCGAAACCACCGGTATGAACCAGa                                             >  1:1748126/1‑35 (MQ=255)
                                                            ttCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAtt       <  1:2931095/71‑1 (MQ=255)
                                                            ttCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAg                   <  1:1332545/59‑1 (MQ=255)
                                                              ctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTgg     >  1:3392598/1‑71 (MQ=255)
                                                              ctcGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGa                  >  1:486156/1‑58 (MQ=255)
                                                              ctcGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTg      >  1:5964623/1‑70 (MQ=255)
                                                               tcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGt    >  1:2330421/1‑71 (MQ=255)
                                                                  aTACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                              <  1:1479543/42‑1 (MQ=255)
                                                                  |                                                                     
GTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGC  >  minE/214859‑214994

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: