Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 982,629 C→T 9.3% T25T (ACC→ACT)  ydfA → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE982,6290CT 9.3% 291.3 / 15.4 97T25T (ACC→ACT) ydfAhypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (40/48);  new base T (2/7);  total (42/55)
Fisher's exact test for biased strand distribution p-value = 2.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.61e-01

CTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAAC  >  minE/982562‑982698
                                                                   |                                                                     
cTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGc                                                                     <  1:991560/70‑1 (MQ=255)
cTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCa                                                                    <  1:2137121/71‑1 (MQ=255)
aaaTCGACCTCGGCAACAACGAATCCCTGGTGTGCGGCGTGGTCCCCAACCAGGATGGAACGTTCACTGcc                                                                    <  1:587623/69‑2 (MQ=255)
 tATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCaa                                                                   >  1:617461/1‑71 (MQ=255)
  aTCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCaa                                                                   >  1:3538746/1‑70 (MQ=255)
  aTCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAAt                                                                  >  1:1185113/1‑71 (MQ=255)
   tCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAAt                                                                  <  1:360413/70‑1 (MQ=255)
   tCGACCTCGGCAACAACGAATCCCTGGTGTGCGGCGTGTTCCCCAACCAGGATGGAACGTTCACTGCCAt                                                                  >  1:368021/1‑70 (MQ=255)
     gATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGAc                                                               >  1:2597968/1‑71 (MQ=255)
     gATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGAc                                                               >  1:2418339/1‑71 (MQ=255)
        cTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGta                                                            >  1:690803/1‑71 (MQ=255)
         ttGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGAc                                                               >  1:1072219/1‑67 (MQ=255)
         ttGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGta                                                            >  1:3167414/1‑70 (MQ=255)
           ggCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACg                                                              <  1:1254250/66‑1 (MQ=255)
           ggCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGta                                                            <  1:3724591/68‑1 (MQ=255)
           ggCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGta                                                            <  1:3304546/68‑1 (MQ=255)
           ggCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGta                                                            <  1:2625073/68‑1 (MQ=255)
           ggCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGt                                                             <  1:64645/67‑1 (MQ=255)
            gCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATAcc                                                        <  1:1647967/71‑1 (MQ=255)
               acaacGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCaa                                                      >  1:3407328/1‑70 (MQ=255)
                  acGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACg                                                              <  1:1649955/59‑1 (MQ=255)
                   cGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCa                                                 <  1:4121568/71‑1 (MQ=255)
                    gAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGc                                                                     >  1:1921473/1‑50 (MQ=255)
                    gAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAAt                                                                  <  1:2877200/53‑1 (MQ=255)
                    gAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAAt                                                                  <  1:1924328/53‑1 (MQ=255)
                    gAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCaa                                                <  1:2559127/71‑1 (MQ=255)
                        ctctGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGc                                                                     <  1:2868212/46‑1 (MQ=255)
                        ctctGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACg                                            >  1:79402/1‑71 (MQ=255)
                          ctGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCaa                                                                   >  1:3436327/1‑46 (MQ=255)
                           tGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACg                                            <  1:2976596/68‑1 (MQ=255)
                             gtgtACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTaa                                       <  1:3445543/71‑1 (MQ=255)
                             gtgtACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTa                                        >  1:2074359/1‑70 (MQ=255)
                             gtgtACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTa                                        >  1:1961839/1‑70 (MQ=255)
                             gtgtACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTa                                        >  1:1749456/1‑70 (MQ=255)
                                tACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGt                                           <  1:2939933/64‑1 (MQ=255)
                                tACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTaa                                       <  1:4012211/68‑1 (MQ=255)
                                tACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTaa                                       <  1:4132219/68‑1 (MQ=255)
                                  cGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCaaa                                               >  1:2091483/1‑58 (MQ=255)
                                   ggCGTGTTTCCAAACCAGGACTGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGa                                 <  1:902362/71‑1 (MQ=255)
                                      gtgtTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGaaaa                              <  1:3696420/71‑1 (MQ=255)
                                      gtgtTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGaaaa                              <  1:3465116/71‑1 (MQ=255)
                                      gtgtTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGaaaa                              <  1:3312984/71‑1 (MQ=255)
                                      gtgtTTCCAAACAAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGaaaa                              <  1:2902657/71‑1 (MQ=255)
                                      gtgtTCCCCAACCAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTAAAACCGaaac                              <  1:3423917/71‑2 (MQ=255)
                                       tgtTTCCAAACCAGGACGGCACGTTCACCGCAATGACGta                                                            <  1:3185363/40‑1 (MQ=255)
                                       tgtTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCaa                                                      <  1:1403516/46‑1 (MQ=255)
                                       tgtTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCaa                                                      <  1:411044/46‑1 (MQ=255)
                                       tgtTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAAt                             <  1:2794418/71‑1 (MQ=255)
                                       tgtTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAAt                             <  1:1176852/71‑1 (MQ=255)
                                       tgtTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAAt                             <  1:2966535/71‑1 (MQ=255)
                                       tgtTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAAt                             <  1:1925809/71‑1 (MQ=255)
                                       tgtTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAAt                             <  1:3095275/71‑1 (MQ=255)
                                       tgtTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAAt                             <  1:857174/71‑1 (MQ=255)
                                       tgtTCCCCAACCAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAAcat                                           <  1:2597139/57‑3 (MQ=255)
                                       tgtTCCCCAACCAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAAcat                                           <  1:3212676/57‑3 (MQ=255)
                                        gtTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGaa                                <  1:660993/67‑1 (MQ=255)
                                         tttCCAAACCAGGACGGCACGTTCACCGCAATGACg                                                              <  1:113154/36‑1 (MQ=255)
                                         tttCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGaa                                <  1:2608187/66‑1 (MQ=255)
                                         tttCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATg                            <  1:3325922/70‑1 (MQ=255)
                                         tttCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACCTTTAAAACCGaaaa                              <  1:3151060/68‑1 (MQ=255)
                                          ttCCAAACCAGGACGGCACGTTCACCGCAATGACg                                                              <  1:3806118/35‑1 (MQ=255)
                                          ttCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCaa                                                >  1:1860149/1‑49 (MQ=255)
                                          ttCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCaa                                                >  1:3195584/1‑49 (MQ=255)
                                          ttCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGt                          >  1:1701101/1‑71 (MQ=255)
                                           tCCAAACCAGGACGGCACGTTCACCGCAATGACGta                                                            >  1:928099/1‑36 (MQ=255)
                                           tCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCaaaa                                              >  1:2743020/1‑50 (MQ=255)
                                           tCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCaa                                                >  1:2187549/1‑48 (MQ=255)
                                           tCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATgg                           >  1:2943252/1‑69 (MQ=255)
                                           tCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATgg                           >  1:2938365/1‑69 (MQ=255)
                                             ccaaCCAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTAAAACCGAAACTGGTGcg                       <  1:2021901/69‑2 (MQ=255)
                                             ccaaCCAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTAAAACCGAAACTGGTGcg                       <  1:3405802/69‑2 (MQ=255)
                                             cAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAAt                             <  1:1902389/65‑1 (MQ=255)
                                             cAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAAt                             <  1:2615246/65‑1 (MQ=255)
                                             cAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGcc                       >  1:3640925/1‑71 (MQ=255)
                                             cAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGcc                       <  1:3279710/71‑1 (MQ=255)
                                             aaaaCCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGcc                       >  1:2602351/2‑71 (MQ=255)
                                               aaCCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCg                     >  1:2528531/1‑71 (MQ=255)
                                                aCCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGtt                                          >  1:2418317/1‑49 (MQ=255)
                                                aCCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGt                    >  1:4142994/1‑71 (MQ=255)
                                                 ccAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTAAAACCGAAACTGGTGcgcgcc                   <  1:1982902/71‑3 (MQ=255)
                                                 ccAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTc                   >  1:2620667/1‑71 (MQ=255)
                                                 ccAGGACGGCACGTTCACCGCAATGACGTATACAAAAAGCAAAACg                                            >  1:1177457/1‑46 (MQ=255)
                                                    ggACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCgct                <  1:1588817/71‑1 (MQ=255)
                                                     gACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCgctg               <  1:1125718/71‑1 (MQ=255)
                                                     gACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCgctg               >  1:3832337/1‑71 (MQ=255)
                                                        ggCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAAt                             <  1:3469862/54‑1 (MQ=255)
                                                         gaaCGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTAAAACCGAAACTGGTGCGCGCCgatggttg           >  1:1052884/3‑71 (MQ=255)
                                                          cACGTTCACCGCAATGACGTATACCAAAAGCaaaa                                              >  1:1684374/1‑35 (MQ=255)
                                                          cACGTTCACCGCAATGACGTATACCAAAAGCaaaa                                              >  1:740026/1‑35 (MQ=255)
                                                          cACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATgg                           >  1:1660293/1‑54 (MQ=255)
                                                          cACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCgctggctg           >  1:3328170/1‑70 (MQ=255)
                                                           aCGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGga         <  1:2083453/71‑1 (MQ=255)
                                                            cGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGgaa        <  1:1497759/71‑1 (MQ=255)
                                                             gTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATgg                           >  1:1326102/1‑51 (MQ=255)
                                                                cACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAAc                                    >  1:2146386/1‑39 (MQ=255)
                                                                 aCCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGgaaagaaa   >  1:708267/1‑71 (MQ=255)
                                                                 aCCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGgaaagaaa   >  1:3224748/1‑71 (MQ=255)
                                                                  ccGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGCTGCCCGTCGCTGGCTGGAAAGAAAc  >  1:907464/1‑71 (MQ=255)
                                                                   |                                                                     
CTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGGCTGGAAAGAAAC  >  minE/982562‑982698

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑