Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 323,651 C→T 7.0% I21I (ATC→ATT htpG → molecular chaperone HSP90 family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE323,6510CT 7.0% 232.2 / 6.0 71I21I (ATC→ATThtpGmolecular chaperone HSP90 family
Reads supporting (aligned to +/- strand):  ref base C (38/28);  new base T (3/2);  total (41/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

CCTACATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGAAGTGAAACAGCTTCTGCACCTGATGATCCATTCTCTCTATTCCAATAAAGAAATCTTCCTGCGTGAGCTTATCTCTAACGCCTCCGATGCGGCGGAC  >  minE/323584‑323720
                                                                   |                                                                     
ccTACATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGAAGTGAAACAGCTTCTGCACCTGATGATCCa                                                                     <  1:2255946/70‑1 (MQ=255)
ccTACATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGAAGTGAAACAGCTTCTGCACCTGATGATCCa                                                                     <  1:1574740/70‑1 (MQ=255)
 cTACATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGAAGTGAAACAGCTTCTGCACCTGATGATCCAtt                                                                   >  1:2882928/1‑71 (MQ=255)
    cATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGAAGTGAAACAGCTTCTGCACCTGATGATCCATtctc                                                                >  1:3190739/1‑71 (MQ=255)
    cATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGAAGTGAAACAGCTTCTGCACCTGATGATCCATtctc                                                                >  1:2463670/1‑71 (MQ=255)
    cATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGAAGTGAAACAGCTTCTGCACCTGATGATCCATtctc                                                                >  1:1330465/1‑71 (MQ=255)
    cATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGAAGTGAAACAGCTTCTGCACCTGATGATCCATtct                                                                 >  1:2099096/1‑70 (MQ=255)
          aGGACAAGAAACTCGTGGTTTTCAGTCAGAAGTGAAACAGCTTCTGCACCTGATGATCCa                                                                     >  1:2066955/1‑60 (MQ=255)
             aCAAGAAACTCGTGGTTTTCAGTCAGAAGTGAAACAGCTTCTGCACCTGATGATCCATTCTCTCTATTcc                                                        >  1:387031/1‑70 (MQ=255)
             aCAAGAAACTCGTGGTTTTCAGTCAGAAGTGAAAAAGCTTCTGCACCTGATGATCCa                                                                     >  1:2923948/1‑57 (MQ=255)
              cAAGAAACTCGTGGTTTTCAGTCAGAAGTGAAACAGCTTCTGCACCTGATGATCCATtct                                                                 >  1:3414076/1‑60 (MQ=255)
               aaGAAACTCGTGGTTTTCAGTCAGAAGTGAAACAGCTTCTGCACCTGATGATCCATTCTCTCTATTCCaa                                                      >  1:925145/1‑70 (MQ=255)
               aaGAAACTCGTGGTTTTCAGTCAGAAGTGAAACAGCTTCTGCACCTGATGATCCATTCTCTCTATTCCaa                                                      >  1:2949882/1‑70 (MQ=255)
               aaGAAACTCGTGGTTTTCAGTCAGAAGTGAAACAGCTTCTGCACCTGATGATCCATTCTCTCTATTCCaa                                                      >  1:2600398/1‑70 (MQ=255)
               aaGAAACTCGTGGTTTTCAGTCAGAAGTGAAACAGCTTCTGCACCTGATGATCCATTCTCTCTATTCCaa                                                      >  1:3326708/1‑70 (MQ=255)
               aaGAAACTCGTGGTTTTCAGTCAGAAGTGAAACAGCTTCTGCACCTGATGATCCATTCTCTCTATTCCaa                                                      >  1:1329571/1‑70 (MQ=255)
                  aaaCTCGTGGTTTTCAGTCAGAAGTGAAACAGCTTCTGCACCTGATGATCCATTCTCTCTATTCCAATa                                                    >  1:1760882/1‑69 (MQ=255)
                  aaaCTCGTGGTTTTCAGTCAGAAGTGAAACAGCTTCTGCACCTGATGATCCATTCTCTCTATTCCAATa                                                    >  1:2144113/1‑69 (MQ=255)
                  aaaCTCGTGGTTTTCAGTCAGAAGTGAAACAGCTTCTGCACCTGATGATCCATTCTCTCTATTCCAATa                                                    >  1:672075/1‑69 (MQ=255)
                   aaCTCGTGGTTTTCAGTCAGAAGTGAAACAGCTTCTGCACCTGATGATCCATTCTCTCTATTCCAATaaa                                                  <  1:1804515/70‑1 (MQ=255)
                           gTTTTCAGTCAGAAGTGAAACAGCTTCTGCACCTGATGATCCATtctctc                                                              <  1:604701/50‑1 (MQ=255)
                           gTTTTCAGTCAGAAGTGAAACAGCTTCTGCACCTGATGATCCATtctctc                                                              <  1:3569052/50‑1 (MQ=255)
                           gTTTTCAGTCAGAAGTGAAACAGCTTCTGCACCTGATGATCCATTCTCTCTATTCCAATAAAGAAATCTTc                                         <  1:1233004/71‑1 (MQ=255)
                             tttCAGTCAGAAGTGAAACAGCTTCTGCACCTGATGATCCATTCTCTCTATTCCAATAAAGAAAt                                             >  1:610590/1‑65 (MQ=255)
                             tttCAGTCAGAAGTGAAACAGCTTCTGCACCTGATGATCCATTCTCTCTATTCCAATAAAGAAATCTTCCt                                       >  1:371222/1‑71 (MQ=255)
                                cagtcagAAGTGAAACAGCTTCTGCACCTGATGATCCATTCTCTCTATTCCAATAAAGAAATCTTCCTGc                                     <  1:1369884/70‑1 (MQ=255)
                                  gtcagAAGTGAAACAGCTTCTGCACCTGATGATCCATTCTCTCTATTCCAATAAAGAAATc                                            <  1:3238967/61‑1 (MQ=255)
                                   tcagAAGTGAAACAGCTTCTGCACCTGATGATCCATTCTCTCTATTCCAATAAAGAAATCTTCCTGCGTGa                                 <  1:816421/71‑1 (MQ=255)
                                   tcagAAGTGAAACAGCTTCTGCACCTGATGATCCATTCTCTCTATTCCAATAAAGAAATCTTCCTGCGTGa                                 <  1:3861034/71‑1 (MQ=255)
                                     agAAGTGAAACAGCTTCTGCACCTGATGATCCATTCTCTCTATTCCAATAAAGAAATCTTCCTGCGTGAGc                               >  1:2990543/1‑71 (MQ=255)
                                       aaGTGAAACAGCTTCTGCACCTGATGATCCATTCTCTCTATTCCAATAAAGAAATCTTCCTGCGTGAGCtt                             <  1:2270439/71‑1 (MQ=255)
                                       aaGTGAAACAGCTTCTGCACCTGATGATCCATTCTCTCTATTCCAATAAAGAAATCTTCCTGCGTGAGCt                              >  1:1903994/1‑70 (MQ=255)
                                        aGTGAAACAGCTTCTGCACCTGATGATCCATTCTCTCTATTCCAATAAAGAAATCTTCCTGCGTGAGCTTa                            <  1:3550891/71‑1 (MQ=255)
                                        aGTGAAACAGCTTCTGCACCTGATGATCCATTCTCTCTATACCAATAAAGAAATCTTCCTGCGTGAGCTTa                            <  1:32092/71‑1 (MQ=255)
                                         gTGAAACAGCTTCTGCACCTGATGATCCATTCTCTCTATTCCAATAAAGaaa                                              <  1:999422/52‑1 (MQ=255)
                                         gTGAAACAGCTTCTGCACCTGATGATCCATTCTCTCTATTCCAATAAAGAAATCTTCCTGCGTg                                  <  1:2645261/64‑1 (MQ=255)
                                            aaaCAGCTTCTGCACCTGATGATCCATTCTCTCTATTCCAATAAAGAAATCTTCCTGCGTGAGCTTAtctc                        >  1:3263640/1‑71 (MQ=255)
                                             aaCAGCTTCTGCACCTGATGATCCATTCTCTCTATTCCAATAAAGAAATCTTCCTGCGTGAGCTTATctct                       <  1:2622752/71‑1 (MQ=255)
                                              aCAGCTTCTGCACCTGATGATCCATTCTCTCTATTCCAATAAAGAAATCTTCCTGCGTGAGCTTATCTCTa                      <  1:3053736/71‑1 (MQ=255)
                                               cAGCTTCTGCACCTGATGATCCATTCTCTCTATTCCAATAAAGaa                                               <  1:2806911/45‑1 (MQ=255)
                                               cAGCTTCTGCACCTGATGATCCATTCTCTCTATTCCAATAAAGaa                                               <  1:1192059/45‑1 (MQ=255)
                                               cAGCTTCTGCACCTGATGATCCATTCTCTCTATTCCAATAAAGaa                                               <  1:3326170/45‑1 (MQ=255)
                                                aGCTTCTGCACCTGATGATCCATTCTCTCTATTCCAATAAAGAAATCTTCCTGCGTGAGCTTATCTCTAAc                    >  1:2346884/1‑71 (MQ=255)
                                                   ttCTGCACCTGATGATCCATTCTCTCTATTCCAATAAAGAAATCTTCCTGCGTGAGCTTATCTCTAACGcc                 >  1:3372737/1‑71 (MQ=255)
                                                    tCTGCACCTGATGATCCATTCTCTCTATTCCAATAAAGAAATCTTCCTGCGTGAGCTTATCTCTAACGCCt                >  1:3089595/1‑71 (MQ=255)
                                                     cTGCACCTGATGATCCATTCTCTCTATTCCAATAAAGAAAt                                             >  1:3427458/1‑41 (MQ=255)
                                                     cTGCACCTGATGATCCATTCTCTCTATTCCAATAAAGAAATCTTCCTGCGTGAGCTTATCTCTAACGCCTc               >  1:2425837/1‑71 (MQ=255)
                                                     cTGCACCTGATGATCCATTCTCTCTATTCCAATAAAGAAATCTTCCTGCGTGAGCTTATCTCTAACGCCTc               >  1:2029569/1‑71 (MQ=255)
                                                       gCACCTGATGCTCCATTCTCTCTATTCCAATAAAGAAATCTTcc                                        >  1:1367647/1‑44 (MQ=255)
                                                       gCACCTGATGATCCATTCTCTCTATTCCAATAAAGAAATCTTcc                                        >  1:1323488/1‑44 (MQ=255)
                                                       gCACCTGATGATCCATTCTCTCTATTCCAATAAAGAAATCTTcc                                        >  1:380800/1‑44 (MQ=255)
                                                       gCACCTGATGATCCATTCTCTCTATTCCAATAAAGAAATCTTCCTGCGTGAGCTTATctct                       >  1:3042005/1‑61 (MQ=255)
                                                       gCACCTGATGATCCATTCTCTCTATTCCAATAAAGAAATCTTCCTGCGTGAGCTTATctct                       >  1:2331874/1‑61 (MQ=255)
                                                       gCACCTGATGATCCATTCTCTCTATTCCAATAAAGAAATCTTCCTGCGTGAGCTTATctct                       >  1:2294259/1‑61 (MQ=255)
                                                        cACCTGATGATTCATTCTCTCTATTCGAACAAGGAAATCTTCCTCCGCGAACTGAtttcc                       <  1:578528/60‑2 (MQ=255)
                                                        cACCTGATGATTCATTCTCTCTATTCGAACAAGGAAATCTTCCTCCGCGAACTGAtttcc                       <  1:352644/60‑2 (MQ=255)
                                                           cTGATGATCCATTCTCTCTATTCCAATAAAGAAATCTTCCTGCg                                    <  1:831061/44‑1 (MQ=255)
                                                            tgatgaTCCATTCTCTCTATTCCAATAAAGAAATCTTCCTGCGTGAGCTTATctct                       <  1:2318751/56‑1 (MQ=255)
                                                            tgatgaTCCATTCTCTCTATTCCAATAAAGAAATCTTCCTGCGTGAGCTTATCTCTAACGCCTcc              <  1:329765/65‑1 (MQ=255)
                                                             gatgatCCATTCTCTCTATTCCAATAAAGAAATCTTCCTGCGTGAGCTTATctct                       <  1:3105474/55‑1 (MQ=255)
                                                              atgatCCATTCTCTCTATTCCAATAAAGAAATCTTCCTGCGTGAGCTTATCTCTAACGCCTCCGATgcggc      >  1:182401/1‑71 (MQ=255)
                                                              atgatCCATTCTCTCTATTCCAATAAAGAAATCTTCCTGCGTGAGCTTATCTCTAACGCCTCCGATgcggc      >  1:1450887/1‑71 (MQ=255)
                                                              atgatCCATTCTCTCTATTCCAATAAAGAAATCTTCCTGCGTGAGCTTATCTCTAACGCCTCCGATgcg        <  1:559793/69‑1 (MQ=255)
                                                               tgatTCATTCTCTCTATTCGAACAAGGAAATCTTCCTCCGCGAACTGATTTCCAACGCCTCGGAcgccgcc     >  1:3132299/1‑70 (MQ=255)
                                                               tgatTCATTCTCTCTATTCGAACAAGGAAATCTTCCTCCGCGAACTGATTTCCAACGCCTCGGAcgccgcc     >  1:1492779/1‑70 (MQ=255)
                                                               tgatTCATTCTCTCTATTCGAACAAGGAAATCTTCCTCCGCGAACTGATTTCCAACGCCTCGGAcgccgcc     >  1:1239532/1‑70 (MQ=255)
                                                               tgatCCATTCTCTCTATTCCAATAAAGAAATCTTCCTGCGTGAGCTTAt                           >  1:1109411/1‑49 (MQ=255)
                                                               tgatCCATTCTCTCTATTCCAATAAAGAAATCTTCCTGCGTGAGCTTAt                           >  1:2667968/1‑49 (MQ=255)
                                                                 atCCATTCTCTCTATTCCAATAAAGAAATCTTCCTGCGTGAGCTTATCTCTAACGCCTCCGATGcggcgg    <  1:789668/70‑1 (MQ=255)
                                                                  tCCATTCTCTCTATTCCAATAAAGAAATCTTCCTGCGTGAGCTTATCTCTAACGCCTCCGATGCGGCGGAc  <  1:3252897/71‑1 (MQ=255)
                                                                  tCCATTCTCTCTATTCCAATAAAGAAATCTTCCTGCGTGAGCTTATCTCTAACGCCTCCGATGCGGCGGAc  <  1:2694386/71‑1 (MQ=255)
                                                                   |                                                                     
CCTACATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGAAGTGAAACAGCTTCTGCACCTGATGATCCATTCTCTCTATTCCAATAAAGAAATCTTCCTGCGTGAGCTTATCTCTAACGCCTCCGATGCGGCGGAC  >  minE/323584‑323720

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: