Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 323,669 T→C 7.0% N27N (AAT→AAC htpG → molecular chaperone HSP90 family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE323,6690TC 7.0% 208.7 / 6.7 71N27N (AAT→AAChtpGmolecular chaperone HSP90 family
Reads supporting (aligned to +/- strand):  ref base T (36/30);  new base C (3/2);  total (39/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.79e-01

AAACTCGTGGTTTTCAGTCAGAAGTGAAACAGCTTCTGCACCTGATGATCCATTCTCTCTATTCCAATAAAGAAATCTTCCTGCGTGAGCTTATCTCTAACGCCTCCGATGCGGCGGACAAGCTGCGTTTCCGTGCG  >  minE/323602‑323738
                                                                   |                                                                     
aaaCTCGTGGTTTTCAGTCAGAAGTGAAACAGCTTCTGCACCTGATGATCCATTCTCTCTATTCCAATa                                                                      >  1:2144113/1‑69 (MQ=255)
aaaCTCGTGGTTTTCAGTCAGAAGTGAAACAGCTTCTGCACCTGATGATCCATTCTCTCTATTCCAATa                                                                      >  1:672075/1‑69 (MQ=255)
aaaCTCGTGGTTTTCAGTCAGAAGTGAAACAGCTTCTGCACCTGATGATCCATTCTCTCTATTCCAATa                                                                      >  1:1760882/1‑69 (MQ=255)
 aaCTCGTGGTTTTCAGTCAGAAGTGAAACAGCTTCTGCACCTGATGATCCATTCTCTCTATTCCAATaaa                                                                    <  1:1804515/70‑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:3861034/71‑1 (MQ=255)
                 tcagAAGTGAAACAGCTTCTGCACCTGATGATCCATTCTCTCTATTCCAATAAAGAAATCTTCCTGCGTGa                                                   <  1:816421/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:3326170/45‑1 (MQ=255)
                             cAGCTTCTGCACCTGATGATCCATTCTCTCTATTCCAATAAAGaa                                                                 <  1:1192059/45‑1 (MQ=255)
                             cAGCTTCTGCACCTGATGATCCATTCTCTCTATTCCAATAAAGaa                                                                 <  1:2806911/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:380800/1‑44 (MQ=255)
                                     gCACCTGATGATCCATTCTCTCTATTCCAATAAAGAAATCTTcc                                                          >  1:1323488/1‑44 (MQ=255)
                                     gCACCTGATGATCCATTCTCTCTATTCCAATAAAGAAATCTTCCTGCGTGAGCTTATctct                                         >  1:2331874/1‑61 (MQ=255)
                                     gCACCTGATGATCCATTCTCTCTATTCCAATAAAGAAATCTTCCTGCGTGAGCTTATctct                                         >  1:3042005/1‑61 (MQ=255)
                                     gCACCTGATGATCCATTCTCTCTATTCCAATAAAGAAATCTTCCTGCGTGAGCTTATctct                                         >  1:2294259/1‑61 (MQ=255)
                                      cACCTGATGATTCATTCTCTCTATTCGAACAAGGAAATCTTCCTCCGCGAACTGAtttcc                                         <  1:352644/60‑2 (MQ=255)
                                      cACCTGATGATTCATTCTCTCTATTCGAACAAGGAAATCTTCCTCCGCGAACTGAtttcc                                         <  1:578528/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:1450887/1‑71 (MQ=255)
                                            atgatCCATTCTCTCTATTCCAATAAAGAAATCTTCCTGCGTGAGCTTATCTCTAACGCCTCCGATgcggc                        >  1:182401/1‑71 (MQ=255)
                                            atgatCCATTCTCTCTATTCCAATAAAGAAATCTTCCTGCGTGAGCTTATCTCTAACGCCTCCGATgcg                          <  1:559793/69‑1 (MQ=255)
                                             tgatTCATTCTCTCTATTCGAACAAGGAAATCTTCCTCCGCGAACTGATTTCCAACGCCTCGGAcgccgcc                       >  1:1492779/1‑70 (MQ=255)
                                             tgatTCATTCTCTCTATTCGAACAAGGAAATCTTCCTCCGCGAACTGATTTCCAACGCCTCGGAcgccgcc                       >  1:3132299/1‑70 (MQ=255)
                                             tgatTCATTCTCTCTATTCGAACAAGGAAATCTTCCTCCGCGAACTGATTTCCAACGCCTCGGAcgccgcc                       >  1:1239532/1‑70 (MQ=255)
                                             tgatCCATTCTCTCTATTCCAATAAAGAAATCTTCCTGCGTGAGCTTAt                                             >  1:2667968/1‑49 (MQ=255)
                                             tgatCCATTCTCTCTATTCCAATAAAGAAATCTTCCTGCGTGAGCTTAt                                             >  1:1109411/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)
                                                  cATTCTCTCTATTCCAATAAAGAAATCTTCCTGCGTGAGCTTATCTCTAACGCCTCCGATGCGGCGGACa                   <  1:2168675/70‑1 (MQ=255)
                                                    ttCTCTCTATTCCAATAAAGAAATCTTCCTGCGTGAGCTTATCTCTAACGCCTCCGATGCGGCGGACAAGc                <  1:2886575/71‑1 (MQ=255)
                                                      ctctctATTCCAATAAAGAAATCTTCCTGCGTGAGCTTAt                                             >  1:2573957/1‑40 (MQ=255)
                                                         tctATTCCAATAAAGAAATCTTCCTGCGTGAGCTTATCTCTAACGCCTCCGATGCGGCGGACAAGCTGCGt           >  1:3369438/1‑71 (MQ=255)
                                                         tctATTCCAATAAAGAAATCTTCCTGCGTGAGCTTATCTCTAACGCCTCCGATGCGGCGGACAAGCTGCGt           >  1:3346866/1‑71 (MQ=255)
                                                          ctATTCCAATAAAGAAATCTTCCTGCGTGAGCTTATCTCTAACGCCTCCGATGCGGCGGACAAGCTGCGtt          >  1:300829/1‑71 (MQ=255)
                                                            aTTCCAATAAAGAAATCTTCCTGCGTGAGCTTATCTCTAACGCCTCCGATGCGGCGGa                     <  1:778748/58‑1 (MQ=255)
                                                             ttCCAATAAAGAAATCTTCCTGCGTGAGCTTATCTCTAACGCCTCCGATGCGGCGGa                     <  1:1535752/57‑1 (MQ=255)
                                                             ttCCAATAAAGAAATCTTCCTGCGTGAGCTTATCTCTAACGCCTCCGATGCGGCGGACAAGCTGCGTTTcc       <  1:754836/71‑1 (MQ=255)
                                                               ccAATAAAGAAATCTTCCTGCGTGAGCTTATCTCTAACGCCTCCGATGCGGCGGACAAGCTGCg            >  1:706672/1‑64 (MQ=255)
                                                                  aTAAAGAAATCTTCCTGCGTGAGCTTATCTCTAACGCCTCCGATgcggcg                       >  1:1476747/1‑50 (MQ=255)
                                                                  aTAAAGAAATCTTCCTGCGTGAGCTTATCTCTAACGCCTCCGATGCGGCGGa                     >  1:550939/1‑52 (MQ=255)
                                                                  aTAAAGAAATCTTCCTGCGTGAGCTTATCTCTAACGCCTCCGATGCGGCGGa                     >  1:3017473/1‑52 (MQ=255)
                                                                  aTAAAGAAATCTTCCTGCGTGAGCTTATCTCTAACGCCTCCGATGCGGCGGa                     >  1:2623514/1‑52 (MQ=255)
                                                                  aTAAAGAAATCTTCCTGCGTGAGCTTATCTCTAACGCCTCCGATGCGGCGGACAAGCTGCGTTTc        <  1:2891562/65‑1 (MQ=255)
                                                                  aTAAAGAAATCTTCCTGCGTGAGCTTATCTCTAACGCCTCCGATGCGGCGGACAAGCTGCGTTTCCGTgcg  >  1:2467666/1‑71 (MQ=255)
                                                                  aTAAAGAAATCTTCCTGCGTGAGCTTATCTCTAACGCCTCCGATGCGGCGGACAAGCTGCGTTTCCGTgcg  >  1:917190/1‑71 (MQ=255)
                                                                  aTAAAGAAATCTTCCTGCGTGAGCTTATCTCTAACGCCTCCGATGCGGCGGACAAGCTGCGTTTCCGTgcg  >  1:3418423/1‑71 (MQ=255)
                                                                   tAAAGAAATCTTCCTGCGTGAGCTTATCTCTAACGCCTCCGATGCGGCGGACAAg                 <  1:2636582/55‑1 (MQ=255)
                                                                   tAAAGAAATCTTCCTGCGTGAGCTTATCTCTAACGCCTCCGATGCGGCGGACAAg                 <  1:402544/55‑1 (MQ=255)
                                                                   |                                                                     
AAACTCGTGGTTTTCAGTCAGAAGTGAAACAGCTTCTGCACCTGATGATCCATTCTCTCTATTCCAATAAAGAAATCTTCCTGCGTGAGCTTATCTCTAACGCCTCCGATGCGGCGGACAAGCTGCGTTTCCGTGCG  >  minE/323602‑323738

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: