Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,834,094 A→G 33.3% D65G (GAT→GGT)  ygbM → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,834,0940AG33.3% 72.7 / 58.3 63D65G (GAT→GGT) ygbMconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (22/20);  new base G (8/13);  total (30/33)
Fisher's exact test for biased strand distribution p-value = 4.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTGCAAATCCAAAAGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTATCCGCCCTTCCCGGACGTGAGCACGAAGCACACGCAGACA  >  minE/1834026‑1834159
                                                                    |                                                                 
ccTGCAAATCCAAAAGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGgta                                                                 <  1:1889291/71‑1 (MQ=255)
ccTGCAAATCCAAAAGCAACTGGAGCAAAACCATCTGACACTCGCGCTGTGCAATACCGCTCCCGGAGgta                                                                 <  1:1995791/71‑1 (MQ=255)
  tGCAAATCCAAAAGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGgtat                                                                >  1:1703342/1‑70 (MQ=255)
  tGCAAATCCAAAAGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGatat                                                                >  1:3992672/1‑70 (MQ=255)
  tGCAAATCCAAAAGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGatat                                                                >  1:2838729/1‑70 (MQ=255)
   gCAAATCCAAAAGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATAtt                                                               >  1:19860/1‑70 (MQ=255)
    cAAATCCAAAAGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGGTAtt                                                               <  1:1268010/69‑1 (MQ=255)
    cAAATCCAAAAGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATAtt                                                               <  1:978731/69‑1 (MQ=255)
    cAAATCCAAAAGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATAtt                                                               <  1:1857101/69‑1 (MQ=255)
    cAAATCCAAAAGCAAATGGAGCAAAACCATCTGAGACTGGCGCTGTTCAATACCGCTCCCGGAGATAtt                                                               <  1:102885/69‑1 (MQ=255)
         ccAAAAGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGgaga                                                                   >  1:176224/1‑60 (MQ=255)
         ccAAAAGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGGTATTAATGCgg                                                        >  1:3844107/1‑71 (MQ=255)
          cAAAAGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCggg                                                       >  1:3050709/1‑71 (MQ=255)
          cAAAAGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCggg                                                       >  1:3242364/1‑71 (MQ=255)
            aaaGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCggggg                                                     <  1:897529/71‑1 (MQ=255)
             aaGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCggggg                                                     >  1:60980/1‑70 (MQ=255)
              aGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGGTATTAATGCgg                                                        >  1:1764193/1‑66 (MQ=255)
              aGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGGTATTAATGCgg                                                        >  1:3502976/1‑66 (MQ=255)
              aGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCGGGGTaa                                                   >  1:2696666/1‑71 (MQ=255)
                   cTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGGTATTaa                                                             <  1:3562751/56‑1 (MQ=255)
                   cTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGGTATTaa                                                             <  1:1141773/56‑1 (MQ=255)
                   cTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGGTATTaa                                                             <  1:699248/56‑1 (MQ=255)
                   cTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGGTATTaa                                                             <  1:611859/56‑1 (MQ=255)
                       aGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTa                                          <  1:1911740/71‑1 (MQ=255)
                       aGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTa                                          <  1:2036527/71‑1 (MQ=255)
                       aGCAAAACCATCTGACACTGGCGCAGTTCAATACCGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTa                                          <  1:481653/71‑1 (MQ=255)
                        gCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGGTATTAATGCGGGGGAATGGGGa                                             <  1:2027651/67‑1 (MQ=255)
                        gCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGGTATTAATGCGGGGGAATGGGGa                                             <  1:1535504/67‑1 (MQ=255)
                        gCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGGTATTAATGCGGGGGAATGGGGa                                             <  1:4169535/67‑1 (MQ=255)
                        gCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTAt                                         <  1:694394/71‑1 (MQ=255)
                          aaaaCCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGGTATTAATGCGGGGGAATGGGGATTATcc                                       <  1:3464214/71‑1 (MQ=255)
                             aCCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGGTATTAATGCg                                                         >  1:2365559/1‑50 (MQ=255)
                                aTCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATg                                                           <  1:3889604/45‑1 (MQ=255)
                                aTCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATg                                                           <  1:3702437/45‑1 (MQ=255)
                                aTCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATg                                                           <  1:1137914/45‑1 (MQ=255)
                                  cTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCGGGGGaa                                                   >  1:2398858/1‑51 (MQ=255)
                                   tGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTATCCGCCCTTccc                              >  1:1801727/1‑71 (MQ=255)
                                   tGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTATCCGCCCTTccc                              >  1:3112606/1‑71 (MQ=255)
                                   tGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTATCCGCCCTTcc                               >  1:3161035/1‑70 (MQ=255)
                                     acacTGGCGCTGTTCAATACCGCTCCCGGAGGTATTAATGCGGGGGAATGGGGATTATCCGCACTTcc                               <  1:23513/68‑1 (MQ=255)
                                     acacTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTATCCGCCCTTcc                               <  1:3549487/68‑1 (MQ=255)
                                     acacTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTATCCGCCCTTCCCg                             <  1:1406020/70‑1 (MQ=255)
                                     acacTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTATCCGCCCTTCCCg                             <  1:1358254/70‑1 (MQ=255)
                                       tcTGGCGCTGTTCAATACCGCTCCCGGAGGTATTAATGCGGGGGAATGGGGATTATc                                        <  1:1877379/56‑1 (MQ=40)
                                           gcgcTGTTCAATACCGCTCCCGGAGGTATTAATGCGGGGGAATGGGGATTATCCGCCCTTCCCGGACGt                        <  1:2791910/69‑1 (MQ=255)
                                           gcgcTGTTCAATACCGCTCCCGGAGGTATTAATGCGGGGGAATGGGGATTATCCGCCCTTCCCGGACGt                        >  1:3223607/1‑69 (MQ=255)
                                           gcgcTGTTCAATACCGCTCCCGGAGGTATTAATGCGGGGGAATGGGGATTATCCGCCCTTCCCGGACGt                        >  1:2363006/1‑69 (MQ=255)
                                                gTTCAATACCGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTATCCg                                      <  1:2637958/50‑1 (MQ=255)
                                                gTTCAATACCGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTATCCGCCCTTCCCGGACGTGAGCAc                  >  1:3323195/1‑70 (MQ=255)
                                                   cAATACCGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTATCCGCCCTTCCCGGACGTGAGCACGAAg              >  1:401933/1‑71 (MQ=255)
                                                      tACCGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTATCCGCCCTTCCCGGACGTGAGCACGAAGcac           <  1:3688205/71‑1 (MQ=255)
                                                      tACCGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTATCCGCCCTTCCCGGACGTGAGCACGAAGcac           <  1:1400621/71‑1 (MQ=255)
                                                      tACCGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTATCCGCCCTTCCCGGACGTGAGCACGAAGcac           <  1:537609/71‑1 (MQ=255)
                                                       aCCGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTATCCGCCCTTCCCGGACGTGAGCACGAAGcaca          >  1:1116676/1‑71 (MQ=255)
                                                       aCCGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTATCCGCCCTTCCCGGACGTGAGCACGAAGcaca          >  1:2265463/1‑71 (MQ=255)
                                                        ccGCTCCCGGAGGTATTAATGCGGGGGAATGGGGATTATCCGCCCTTCCCGGACGTGAGCACGAAGCacac         >  1:1114847/1‑71 (MQ=255)
                                                        ccGCTCCCGGAGGTATTAATGCGGGGGAATGGGGATTATCCGCCCTTCCCGGACGTGAGCACGAAGCacac         >  1:3790594/1‑71 (MQ=255)
                                                        ccGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTATCCGCCCTTCCCGGACGTGAGCACGAAGcaca          >  1:4140442/1‑70 (MQ=255)
                                                        ccGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTATCCGCCCTTCCCGGACGTGAGCACGAAGcaca          >  1:217774/1‑70 (MQ=255)
                                                        ccGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTATCCGCCCTTCCCGGACGTGAGCACGAAGCacac         >  1:326595/1‑71 (MQ=255)
                                                         cGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTATCCGCCCTTCCCGGACGTGAGCACGAAGCACACg        >  1:1988692/1‑71 (MQ=255)
                                                         cGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTATCCGCCCTTCCCGGACGTGAGCACGAAGCACACg        >  1:3111596/1‑71 (MQ=255)
                                                           cTCCCGGAGGTATTAATGCGGGGGAATGGGGATTATCCGCCCTTCCCGGACGTGAGCACGAAGCACACGCa      <  1:3831183/71‑1 (MQ=255)
                                                           cTCCCGGAGATATTAATGCGGGGGAATGGGGATTATCCGCCCTTCCCGGACGTGAGCACGAAGCacac         >  1:1538280/1‑68 (MQ=255)
                                                           cTCCCGGAGATATTAATGCGGGGGAATGGGGATTATCCGCCCTTCCCGGACGTGAGCACGAAGCACACg        >  1:58409/1‑69 (MQ=255)
                                                           cTCCCGGAGATATTAATGCGGGGGAATGGGGATTATCCGCCCTTCCCGGACGTGAGCACGAAGCACACg        >  1:17023/1‑69 (MQ=255)
                                                           cTCCCGGAGATATTAATGCGGGGGAATGGGGATTATCCGCCCTTCCCGGACGTGAGCACGAAGCACACg        >  1:3326552/1‑69 (MQ=255)
                                                               cGGAGATATTAATGCGGGGGAATGGGGATTATCCGCCCTTCCCGGACGTGAGCACGAAGCACACGCAGACa  <  1:614809/71‑1 (MQ=255)
                                                               cGGAGATATTAATGCGGGGGAATGGGGATTATCCGCCCTTCCCGGACGTGAGCACGAAGCACACGCAGACa  <  1:1464251/71‑1 (MQ=255)
                                                                    |                                                                 
CCTGCAAATCCAAAAGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTATCCGCCCTTCCCGGACGTGAGCACGAAGCACACGCAGACA  >  minE/1834026‑1834159

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: