Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,491,691 C→T 60.0% T163I (ACC→ATC)  gyrB → DNA gyrase, subunit B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,491,6910CT60.0% 32.2 / 62.7 55T163I (ACC→ATC) gyrBDNA gyrase, subunit B
Reads supporting (aligned to +/- strand):  ref base C (8/14);  new base T (24/9);  total (32/23)
Fisher's exact test for biased strand distribution p-value = 1.18e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCGTCAGATCTACGAACACGGTGTACCGCAGGCCCCGCTGGCGGTTACCGGCGAGACTGAAAAAACCGGCACCATGGTGCGTTTCTGGCCCAGCCTCGAAACCTTCACCAATGTGACCGAGTTCGAATATGAAA  >  minE/2491625‑2491759
                                                                  |                                                                    
aCCGTCAGATCTACGAACACGGTGTACCGCAGGCCCCGCTGGCGGTTACCGGCGAGACTGAAAAAATCgg                                                                   <  1:3540583/70‑1 (MQ=255)
  cGTCAGATCTACGAACACGGTGTACCGCAGGCCCCGCTGGCGGTTACCGGCGAGACTGAAAAAATCGGCAc                                                                >  1:745300/1‑71 (MQ=255)
     cAGATCTACGAACACGCTGTACCGCAGGCCCCGCTGGCGGTTACCGGCGAGACTGAAAAAATCGGCAc                                                                >  1:1621672/1‑68 (MQ=255)
      aGATCTACGAACACGGTGTACCGCAGGCCCCGCTGGCGGTTACCGGCGAGACTGAAAAAATCGGCACCATg                                                            <  1:2761716/71‑1 (MQ=255)
        aTCTACGAACACGGTGTACCGCAGGCCCCGCTGGCGGTTACCGGCGAGACTGAAAAAATCGGCACCa                                                              <  1:4005455/67‑1 (MQ=255)
          cTACGAACACGGTGTACCGCAGGCCCCGCTGGCGGTTACCGGCGAGACTGAAAAAATCGGCACCATGGTGc                                                        >  1:1027645/1‑71 (MQ=255)
          cTACGAACACGGTGTACCGCAGGCCCCGCTGGCGGTTACCGGCGAGACTGAAAAAATCGGCACCATGGTGc                                                        >  1:2036667/1‑71 (MQ=255)
             cGAACACGGTGTACCGCAGGCCCCGCTGGCGGTTACCGGCGAGACTGAAAAAACCGGCACCATGGTg                                                         <  1:3895171/67‑1 (MQ=255)
              gAACACGGTGTACCGCAGGCCCCGCTGGCGGTTACCGGCGAGACTGAAAAAACCGGCACCATGGTGCGt                                                      <  1:1153765/69‑1 (MQ=255)
                acacGGTGTACCGCAGGCCCCGCTGGCTGTTACCGGCGAGACTGAAAAAACCGGCACCATGGTGCGTTTCt                                                  >  1:791873/1‑71 (MQ=255)
                  acGGTGTACCGCAGGCCCCGCTGGCGGTTACCGGCGAGACTGAAAAAATCGGCAc                                                                <  1:3032079/55‑1 (MQ=255)
                   cGGTGTACCGCAGGCCCCGCTGGCGGTTACCGGCGAGACTGAAAAAATCGGCACCATGGTGCGTTTCTGGc                                               >  1:1492510/1‑71 (MQ=255)
                      tgtACCGCAGGCCCCGCTGGCGGTTACCGGCGAGACTGAAAAAATCGGCACCATGGTGCGTTTCTGGcc                                              >  1:77926/1‑69 (MQ=255)
                        tACCGCAGGCCCCGCTGGCGGTTACCGGCGAGACTGAAAAAATCGGCACCATGGTGCGTTTCTGGcc                                              <  1:421933/67‑1 (MQ=255)
                        tACCGCAGGCCCCGCTGGCGGTTACCGGCGAGACTGAAAAAATCGGCACCATGGTGCGTTTCTGGCCCAGc                                          >  1:3263137/1‑71 (MQ=255)
                        tACCGCAGGCCCCGCTGGCGGTTACCGGCGAGACTGAAAAAACCGGCACCATGGTGCGTTTCTGGcc                                              <  1:963817/67‑1 (MQ=255)
                          ccGCAGGCCCCGCTGGCGGTTACCGGCGAGACTGAAAAAACCGGCACCATGGTGCGTTTCTGGCCCAGCCt                                        >  1:151114/1‑71 (MQ=255)
                            gCAGGCCCCGCTGGCGGTTACCGGCGAGACTGAAAAAATCGGCACCATGGTGCGTTTCTGGCCCAGCCTc                                       >  1:2663975/1‑70 (MQ=255)
                             cAGGCCCCGCTGGCGGTTACCGGCGAGACTGAAAAAATCGGCACCATGGTGCGTTTCTGGccc                                             >  1:1234412/1‑63 (MQ=255)
                               ggCCCCGCTGGCGGTTACCGGCGAGACTGAAAAAATAGGCACCATgg                                                           >  1:2435114/1‑47 (MQ=39)
                                 ccccGCTGGCGGTTACCGGCGAGACTGAAAAAACCGGCACCATGGTGCGTTTCTGGCCCAGCCTCGAAAcc                                 >  1:352412/1‑71 (MQ=255)
                                  cccGCTGGCGGTTACCGGCGAGACTGAAAAAACCGGCACCa                                                              >  1:2292994/1‑41 (MQ=255)
                                   ccGCTGGCGGTTACCGGCGAGACTGAAAAAACCGGCACCATGGTGCGtt                                                     <  1:1409823/49‑1 (MQ=255)
                                     gCTGGCGGTTACCGGCGAGACTGAAAAAATCGGCACCATGGTGCg                                                       >  1:2440884/1‑45 (MQ=255)
                                     gCTGGCGGTTACCGGCGAGACTGAAAAAATCGGCACCATGGTGCg                                                       >  1:883391/1‑45 (MQ=255)
                                     gCTGGCGGTTACCGGCGAGACTGAAAAAATCGGCACCATGGTGCGTTTCTGGCCCAGCCTCGAAACCTTCa                             >  1:821472/1‑71 (MQ=255)
                                     gCTGGCGGTTACCGGCGAGACTGAAAAAATCGGCACCATGGTGCGTTTCTGGCCCAGCCTCGAAACCTTCa                             >  1:1520269/1‑71 (MQ=255)
                                     gCTGGCGGTTACCGGCGAGACTGAAAAAATCGGCACCATGGTGCGTTTCTGGCCCAGCCTCGAAACCTTCa                             >  1:3651940/1‑71 (MQ=255)
                                     gCTGGCGGTTACCGGCGAGACTGAAAAAACCGGCACCATGGTGCGt                                                      <  1:2847427/46‑1 (MQ=255)
                                      cTGGCGGTTACCGGCGAGACTGAAAAAATCGGCACCATGGTGCGTTTc                                                   >  1:1643259/1‑48 (MQ=255)
                                      cTGGCGGTTACCGGCGAGACTGAAAAAATCGGCACCATGGTGCGTTTCTGGCCCAGCCTCGAAACCTTCAc                            >  1:2068723/1‑71 (MQ=255)
                                         gCGGTTACCGGCGAGACTGAAAAAATCGGCACCATGGTGCGTTTCTGGCCCAGCCTCGAAACCTTCACCaa                         >  1:3118686/1‑71 (MQ=255)
                                          cGGTTACCGGCGAGACTGAAAAAACCGGCACCATGGTGCGTTTCTGGCCCAGCCTCGaa                                    <  1:1763356/59‑1 (MQ=255)
                                          cGGTTACCGGCGAGACTGAAAAAACCGGCACCATGGTGCGTTTCTGGCCCAGCCTCGAAACCTTCACCAAt                        <  1:2957954/71‑1 (MQ=255)
                                            gTTACCGGCGAGACTGAAAAAATCGTCACCATGGTGCGTTTCTGGCCCAGCCTCGAAACCTTCAcc                           <  1:2497756/66‑1 (MQ=255)
                                              tACCGGCGAGACTGAAAAAATCGGCACCATGGTGCGTTTCTGGCCCAGCCTCGAAACCTTCACCAAtg                       >  1:3548225/1‑68 (MQ=255)
                                               aCCGGCGAGACTGAAAAAATCGGCACCATGGTGCGTTTCTGGCCCAGCCTCGAAACCTTCACCAATGTGAc                   >  1:1795757/1‑71 (MQ=255)
                                                ccGGCGAGACTGAAAAAATCGGCACCATGGTGCGTTTCTGGCCCAGCCTCGAAACCTTCACCAATGTGAcc                  <  1:843172/71‑1 (MQ=255)
                                                ccGGCGAGACTGAAAAAATCGGCACCATGGTGCGTTTCTGGCCCAGCCTCGAAACCTTCACCAATGTGAcc                  <  1:3082293/71‑1 (MQ=255)
                                                  ggCGAGACTGAAAAAACCTGCACCATGGTGCGTTTCTGGCCCAGCCTc                                       <  1:3120276/48‑1 (MQ=255)
                                                       gaCTGAAAAAACCGGCACCATGGTGCGTTTCTGGCCCAGCCTCGAAACCTTCACCAAtg                       >  1:528610/1‑59 (MQ=255)
                                                       gaCTGAAAAAACCGGCACCATGGTGCGTTTCTGGCCCAGCCTCGAAACCTTCACCAATGTGACCGAGTTCg           <  1:2489798/71‑1 (MQ=255)
                                                          tGAAAAAATCGGCACCATGGTGCGTTTCTGGCCCAGCCTCGAAACCTTCACCa                          <  1:997113/53‑1 (MQ=255)
                                                           gAAAAAATCGGCACCATGGTGCGTTTCTGGCCCAGCCTCGAAACCTTCACCAATGTGACCGAGTTCGAata       >  1:357823/1‑71 (MQ=255)
                                                            aaaaaaTCGGCACCATGGTGCGTTTCTGGCCCAGCCTCGAAACCTTCACCAAt                        >  1:3648228/1‑53 (MQ=255)
                                                            aaaaaaCCGGCACCATGGTGCGTTTCTGGCCCAGCCTCGAAACCTTCACCAATGTGACCGAGTTCGAatat      >  1:3647949/1‑71 (MQ=255)
                                                             aaaaaTCGGCACCATGGTGCGTTTCTGGCCCAGCCTCGAAACCTTCACCAATGTGAc                   >  1:2634856/1‑57 (MQ=255)
                                                             aaaaaTCGGCACCATGGTGCGTTTCTGGCCCAGCCTCGAAACCTTCACCAATGTGAc                   >  1:2754138/1‑57 (MQ=255)
                                                             aaaaaCCGGCACCATGGTGCGTTTCTGGCCCAGCCTCGAAACCTTCACCAATGTGACCGAGTTCGAata       <  1:4127689/69‑1 (MQ=255)
                                                              aaaaCCGGCACCATGGTGCGTTTCTGGCCAAGCCTCGAAACCTTCACCAATGTGACCGAGTTCGAATATg     >  1:2626996/1‑70 (MQ=255)
                                                               aaaCCGGCACCATGGTGCGTTTCTGGCCCAGCCTCGAAACCTTCACTAAt                        <  1:3550337/50‑1 (MQ=255)
                                                               aaaCCGGCACCATGGTGCGTTTCTGGCCCAGCCTCGAAACCTTCACCAAt                        <  1:4008491/50‑1 (MQ=255)
                                                                aaCCGGCACCATGGTGCGTTTCTGGCCCAGCCTCGAAACCTTCACCAATGTGACCGAGTTCGAATATGaaa  <  1:2443081/71‑1 (MQ=255)
                                                                aaCCGGCACCATGGTGCGTTTCTGGCCCAGCCTCGAAACCTTCACCAATGTGACCGAGTTCGAATATGaaa  <  1:1444822/71‑1 (MQ=255)
                                                                 aCCGGCACCATGGTGCGTTTCTGGCCCAGCCTCGAAACCTTCACCAAtg                       >  1:1198436/1‑49 (MQ=255)
                                                                  |                                                                    
ACCGTCAGATCTACGAACACGGTGTACCGCAGGCCCCGCTGGCGGTTACCGGCGAGACTGAAAAAACCGGCACCATGGTGCGTTTCTGGCCCAGCCTCGAAACCTTCACCAATGTGACCGAGTTCGAATATGAAA  >  minE/2491625‑2491759

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: ‑