Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 39,415 A→G 23.7% G65G (GGT→GGC caiT ← predicted transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE39,4150AG23.7% 136.4 / 36.3 109G65G (GGT→GGCcaiTpredicted transporter
Reads supporting (aligned to +/- strand):  ref base A (35/48);  new base G (12/14);  total (47/62)
Fisher's exact test for biased strand distribution p-value = 8.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.92e-01

Only 100 of 151 total aligned reads displayed.

GCTAAATTCTGGCGGTTCGTTACCTAAACGCTTTTTGGCATACGGGCCAAACACCAGCCAGAACCAACCGAAAAGCATCACCACCATATACCATTCAAATGCCCATCCCCATACATTGGTGACGTAACTGAA  >  minE/39349‑39480
                                                                  |                                                                 
gCTAAATTCTGGCGGTTCGTTACCTAAACGCTTTTTGGCATACGGGCCAAACACCAGCCAGaaccaac                                                                  >  1:3027751/1‑68 (MQ=255)
gCTAAATTCTGGCGGTTCGTTACCTAAACGCTTTTTGGCATACGGGCCAAACACCAGCCAGaaccaac                                                                  >  1:740061/1‑68 (MQ=255)
gCTAAATTCTGGCGGTTCGTTACCTAAACGCTTTTTGGCATACGGGCCAAACACCAGCCAGaaccaac                                                                  >  1:1432452/1‑68 (MQ=255)
gCTAAATTCTGGCGGTTCGTTACCTAAACGCTTTTTGGCATACGGGCCAAACACCAGCCAGaaccaac                                                                  >  1:3227859/1‑68 (MQ=255)
 cTAAATTCTGGCGGTTCGTTACCTAAACGCTTTTTGGCATACGGGCCAAACACCAGCCAGaaccaacc                                                                 >  1:2311138/1‑68 (MQ=255)
 cTAAATTCTGGCGGTTCGTTACCTAAACGCTTTTTGGCATACGGGCCAAACACCAGCCAGaaccaacc                                                                 >  1:816987/1‑68 (MQ=255)
 cTAAATTCTGGCGGTTCGTTACCTAAACGCTTTTTGGCATACGGGCCAAACACCAGCCAGaaccaacc                                                                 >  1:3164186/1‑68 (MQ=255)
   aaaTTCTGGTGGTTCGTTACCTAAACGCTTTTTGGCATACGGGCCAAACACCAGCCAGAACCAGCCg                                                                >  1:933840/1‑67 (MQ=255)
   aaaTTCTGGTGGTTCGTTACCTAAACGCTTTTTGGCATACGGGCCAAACACCAGCCAGAACCAGCCg                                                                >  1:2531566/1‑67 (MQ=255)
   aaaTTCTGGTGGTTCGTTACCTAAACGCTTTTTGGCATACGGGCCAAACACCAGCCAGAACCAGCCg                                                                >  1:2756526/1‑67 (MQ=255)
   aaaTTCTGGTGGTTCGTTACCTAAACGCTTTTTGGCATACGGGCCAAACACCAGCCAGAACCAGCCg                                                                >  1:3055009/1‑67 (MQ=255)
   aaaTTCTGGTGGTTCGTTACCTAAACGCTTTTTGGCATACGGGCCAAACACCAGCCAGAACCAGCCg                                                                >  1:3218113/1‑67 (MQ=255)
   aaaTTCTGGTGGTTCGTTACCTAAACGCTTTTTGGCATACGGGCCAAACACCAGCCAGAACCAGCCg                                                                >  1:919731/1‑67 (MQ=255)
   aaaTTCTGGTGGTTCGTTACCTAAACGCTTTTTGGCATACGGGCCAAACACCAGCCAGAACCAGCCg                                                                >  1:835748/1‑67 (MQ=255)
        cTGGCGGTTCGTTACCTAAACGCTTTTTGGCATACGGGCCAAACACCAGCCAGAACCAACCGAAAAGc                                                          >  1:2380840/1‑68 (MQ=255)
             ggTTCGTTACCTAAACGCTTTTTGGCATACGGGCCAAACACCAGCCAGAACCAACCGAAAAGCAta                                                       <  1:2483873/66‑2 (MQ=255)
             ggTTCGTTACCTAAACGCTTTTTGGCATACGGGCCAAACACCAGCCAGAACCAACCGAAAAGCATc                                                       <  1:3183271/66‑1 (MQ=255)
             ggTTCGTTACCTAAACGCTTTTTGGCATACGGGCCAAACACCAGCCAGAACCAACCGAAAAGCATc                                                       <  1:174380/66‑1 (MQ=255)
                tCGTTACCTAAACGCTTTTTGGCATACGGGCCAAACACCAGCCAGAACCAACCGAAAAGCATcaccac                                                  >  1:3429475/1‑68 (MQ=255)
                tCGTTACCTAAACGCTTTTTGGCATACGGGCCAAACACCAGCCAGAACCAACCGAAAAGCATcaccac                                                  >  1:2258959/1‑68 (MQ=255)
                tCGTTACCTAAACGCTTTTTGGCATACGGGCCAAACACCAGCCAGAACCAACCGAAAAGCATcaccac                                                  >  1:2837230/1‑68 (MQ=255)
                tCGTTACCTAAACGCTTTTTGGCATACGGGCCAAACACCAGCCAGAACCAACCGAAAAGCATcaccac                                                  >  1:3052808/1‑68 (MQ=255)
                tCGTTACCTAAACGCTTTTTGGCATACGGGCCAAACACCAGCCAGAACCAACCGAAAAGCATcaccac                                                  >  1:385957/1‑68 (MQ=255)
                 cGTTACCTAAACGCTTTTTGGCATACGGGCCAAACACCAGCCAGAACCAACCGAAAAGCATCAccacca                                                <  1:2983525/69‑1 (MQ=255)
                 cGTTACCTAAACGCTTTTTGGCATACGGGCCAAACACCAGCCAGAACCAACCGAAAAGCATCAccacca                                                <  1:792680/69‑1 (MQ=255)
                 cGTTACCTAAACGCTTTTTGGCATACGGGCCAAACACCAGCCAGAACCAACCGAAAAGCATCAccacca                                                <  1:1307892/69‑1 (MQ=255)
                 cGTTACCTAAACGCTTTTTGGCATACGGGCCAAACACCAGCCAGAACCAACCGAAAAGCATCAccacca                                                <  1:1443409/69‑1 (MQ=255)
                 cGTTACCTAAACGCTTTTTGGCATACGGGCCAAACACCAGCCAGAACCAACCGAAAAGCATCAccacca                                                <  1:2057062/69‑1 (MQ=255)
                         aaaCGCTTTTTGGCATACGGGCCAAACACCAGCCAGAACCAACCGAAAAGc                                                          <  1:207376/51‑1 (MQ=255)
                            cGCTTTTTGGCATACGGGCCAAACACCAGCCAGAACCAACCGAAAAGCATCACCACCATATACCATTCa                                     <  1:2461338/69‑1 (MQ=255)
                            cGCTTTTTGGCATACGGGCCAAACACCAGCCAGAACCAACCGAAAAGCATCACCACCATATACCATTCa                                     <  1:2375549/69‑1 (MQ=255)
                            cGCTTTTTGGCATACGGGCCAAACACCAGCCAGAACCAACCGAAAAGCATCACCACCATATACCATTCa                                     <  1:1017109/69‑1 (MQ=255)
                            cGCTTTTTGGCATACGGGCCAAACACCAGCCAGAACCAACCGAAAAGCATCACCACCATATACCATTCa                                     <  1:2114824/69‑1 (MQ=255)
                            cGCTTTTTGGCATACGGGCCAAACACCAGCCAGAACCAACCGAAAAGCATCACCACCATATACCATTCa                                     <  1:2024802/69‑1 (MQ=255)
                            cGCTTTTTGGCATACGGGCCAAACACCAGCCAGAACCAACCGAAAAGCATCACCACCATATACCATTCa                                     <  1:1321117/69‑1 (MQ=255)
                            cGCTTTTTGGCATACGGGCCAAACACCAGCCAGAACCAACCGAAAAGCATCACCACCATATACCATTCa                                     <  1:1477429/69‑1 (MQ=255)
                            cGCTTTTTGGCATACGGGCCAAACACCAGCCAGAACCAACCGAAAAGCATCACCACCATATACCATTCa                                     <  1:1510573/69‑1 (MQ=255)
                                 tttGGCATACGGGCCAAACACCAGCCAGAACCAACCGAAAAGCATCACCACCATATACCATTCa                                     <  1:3365821/64‑1 (MQ=255)
                                   tGGCATACGGGCCAAACACCAGCCAGAACCAACCGAAAAGCATc                                                       <  1:1945534/44‑1 (MQ=255)
                                          cGGGCCAAACACCAGCCAGAACCAACCGAAAAGCAt                                                        <  1:1864368/36‑1 (MQ=255)
                                          cGGGCCAAACACCAGCCAGAACCAACCGAAAAGCATCACCACCATATACCATTCAAATg                                 <  1:1004178/59‑1 (MQ=255)
                                          cGGGCCAAACACCAGCCAGAACCAACCGAAAAGCATCACCACCATATACCATTCAAATg                                 <  1:3091004/59‑1 (MQ=255)
                                          cGGGCCAAACACCAGCCAGAACCAACCGAAAAGCATCACCACCATATACCATTCAAATg                                 <  1:2528672/59‑1 (MQ=255)
                                          cGGGCCAAACACCAGCCAGAACCAACCGAAAAGCATCACCACCATATACCATTCAAATGCCCATCCCCa                       >  1:648406/1‑69 (MQ=255)
                                          cGGGCCAAACACCAGCCAGAACCAACCGAAAAGCATCACCACCATATACCATTCAAATGCCCATCCCCa                       >  1:3292203/1‑69 (MQ=255)
                                          cGGGCCAAACACCAGCCAGAACCAACCGAAAAGCATCACCACCATATACCATTCAAATGCCCATCCCCa                       >  1:2184741/1‑69 (MQ=255)
                                          cGGGCCAAACACCAGCCAGAACCAACCGAAAAGCATCACCACCATATACCATTCAAATGCCCATCCCCa                       >  1:1966367/1‑69 (MQ=255)
                                          cGGGCCAAACACCAGCCAGAACCAACCGAAAAGCATCACCACCATATACCATTCAAATGCCCATCCCCa                       >  1:2361105/1‑69 (MQ=255)
                                          cGGGCCAAACACCAGCCAGAACCAACCGAAAAGCATCACCACCATATACCATTCAAATGCCCATCCCCa                       >  1:508953/1‑69 (MQ=255)
                                          cGGGCCAAACACCAGCCAGAACCAACCGAAAAGCATCACCACCATATACCATTCAAATGCCCATCCCCa                       >  1:3031646/1‑69 (MQ=255)
                                          cGGGCCAAACACCAGCCAGAACCAACCGAAAAGCATCACCACCATATACCATTCAAATGCCCATCCCCa                       >  1:1570997/1‑69 (MQ=255)
                                          cGGGCCAAACACCAGCCAGAACCAACCGAAAAGCATCACCACCATATACCATTCAAATGCCCATCCCCa                       >  1:560949/1‑69 (MQ=255)
                                          cGGGCCAAACACCAGCCAGAACCAACCGAAAAGCATCACCACCATATACCATTCAAATGCCCATCCCCa                       >  1:635186/1‑69 (MQ=255)
                                          cGGGCCAAACACCAGCCAGAACCAACCGAAAAGCATCACCACCATATACCATTCAAATGCCCATCCCCa                       >  1:2964000/1‑69 (MQ=255)
                                          cGGGCCAAACACCAGCCAGAACCAACCGAAAAGCATCACCACCATATACCATTCAAATGCCCATCCCCa                       >  1:2872748/1‑69 (MQ=255)
                                          cGGGCCAAACACCAGCCAGAACCAACCGAAAAGCATCACCACCATATACCATTCAAATGCCCATCCCCa                       <  1:1629179/69‑1 (MQ=255)
                                          cGGGCCAAACACCAGCCAGAACCAACCGAAAAGCATCACCACCATATACCATTCAAATGCCCATCCCCa                       >  1:900359/1‑69 (MQ=255)
                                          cGGGCCAAACACCAGCCAGAACCAACCGAAAAGCATCACCACCATATACCATTCAAATGCCCATCCCCa                       >  1:1173503/1‑69 (MQ=255)
                                          cGGGCCAAACACCAGCCAGAACCAACCGAAAAGCATCACCACCATATACCATTCAAATGCCCATCCCCa                       >  1:808400/1‑69 (MQ=255)
                                          cGGGCCAAACACCAGCCAGAACCAACCGAAAAGCATCACCACCATATACCATTCAAATGCCCATCCCCa                       >  1:1667150/1‑69 (MQ=255)
                                          cGGGCCAAACACCAGCCAGAACCAACCGAAAAGCATCACCACCATATACCATTCAAATGCCCATCCCCa                       >  1:263803/1‑69 (MQ=255)
                                           gggCCAAACACCAGCCAGAACCAACCGAAAAGCATCACCACCATATACCATTCAAATg                                 <  1:1416777/58‑1 (MQ=255)
                                            ggCCAAACACCAGCCAGAACCAACCGAAAAGCATCACCACCATATACCATTCAAATGCCCATCCCCATa                     <  1:857618/69‑1 (MQ=255)
                                                  acacCAGCCAGAACCAGCCGAAAAGCATCACCACCATATACCATTCAAATGCCCATCCCCATACATTgg               <  1:2570525/69‑1 (MQ=255)
                                                  acacCAGCCAGAACCAGCCGAAAAGCATCACCACCATATACCATTCAAATGCCCATCCCCATACATTgg               <  1:1073571/69‑1 (MQ=255)
                                                  acacCAGCCAGAACCAGCCGAAAAGCATCACCACCATATACCATTCAAATGCCCATCCCCATACATTgg               <  1:650846/69‑1 (MQ=255)
                                                  acacCAGCCAGAACCAGCCGAAAAGCATCACCACCATATACCATTCAAATGCCCATCCCCATACATTgg               <  1:1771975/69‑1 (MQ=255)
                                                  acacCAGCCAGAACCAGCCGAAAAGCATCACCACCATATACCATTCAAATGCCCATCCCCATACATTgg               <  1:2483754/69‑1 (MQ=255)
                                                  acacCAGCCAGAACCAGCCGAAAAGCATCACCACCATATACCATTCAAATGCCCATCCCCATACATTgg               <  1:1280897/69‑1 (MQ=255)
                                                  acacCAGCCAGAACCAGCCGAAAAGCATCACCACCATATACCATTCAAATGCCCATCCCCATACATTgg               <  1:1475645/69‑1 (MQ=255)
                                                  acacCAGCCAGAACCAGCCGAAAAGCATCACCACCATATACCATTCAAATGCCCATCCCCATACATTgg               <  1:3222015/69‑1 (MQ=255)
                                                  acacCAGCCAGAACCAGCCGAAAAGCATCACCACCATATACCATTCAAATGCCCATCCCCATACATTgg               <  1:2469413/69‑1 (MQ=255)
                                                  acacCAGCCAGAACCAGCCGAAAAGCATCACCACCATATACCATTCAAATGCCCATCCCCATACATTgg               <  1:3340017/69‑1 (MQ=255)
                                                  acacCAGCCAGAACCAGCCGAAAAGCATCACCACCATATACCATTCAAATGCCCATCCCCATACATTgg               <  1:2859743/69‑1 (MQ=255)
                                                   cacCAGCCAGAACCAGCCGAAAAGCATCACCACCATATACCATTCAAATGCCCATCCCCATACATTgg               <  1:1793512/68‑1 (MQ=255)
                                                    acCAGCCAGAACCAACCGAAAATCATCACCACCATATACCATTCa                                     <  1:432358/45‑1 (MQ=255)
                                                        gccagAACCAACCGAAAAGCATCACCACCATATACCATTCAAATGCCCATCCCCATACATTGGTGACGt         <  1:3411011/69‑1 (MQ=255)
                                                             aaccaaccGAAAAGCATCACCACCATATACCATTCAAATGCCCATCCCCATACATTGGTGACGTAACTg    <  1:326952/69‑1 (MQ=255)
                                                             aaccaaccGAAAAGCATCACCACCATATACCATTCAAATGCCCATCCCCATACATTGGTGACGTAACTg    <  1:2126773/69‑1 (MQ=255)
                                                             aaccaaccGAAAAGCATCACCACCATATACCATTCAAATGCCCATCCCCATACATTGGTGACGTAACTg    <  1:2825210/69‑1 (MQ=255)
                                                               ccagccGAAAAGCATCACCACCATATACCATTCAAATGCCCATCCCCATACATTGGTGACGTAACTg    >  1:1407288/1‑67 (MQ=255)
                                                               ccagccGAAAAGCATCACCACCATATACCATTCAAATGCCCATCCCCATACATTGGTGACGTAACTg    >  1:2003149/1‑67 (MQ=255)
                                                               ccagccGAAAAGCATCACCACCATATACCATTCAAATGCCCATCCCCATACATTGGTGACGTAACTg    >  1:1073039/1‑67 (MQ=255)
                                                               ccagccGAAAAGCATCACCACCATATACCATTCAAATGCCCATCCCCATACATTGGTGACGTAACTg    >  1:848964/1‑67 (MQ=255)
                                                               ccagccGAAAAGCATCACCACCATATACCATTCAAATGCCCATCCCCATACATTGGTGACGTAACTGa   >  1:157919/1‑68 (MQ=255)
                                                               ccagccGAAAAGCATCACCACCATATACCATTCAAATGCCCATCCCCATACATTGGTGACGTAACTGa   >  1:2981018/1‑68 (MQ=255)
                                                               ccagccGAAAAGCATCACCACCATATACCATTCAAATGCCCATCCCCATACATTGGTGACGTAACTGa   >  1:644511/1‑68 (MQ=255)
                                                               ccagccGAAAAGCATCACCACCATATACCATTCAAATGCCCATCCCCATACATTGGTGACGTAACTGa   >  1:544362/1‑68 (MQ=255)
                                                               ccagccGAAAAGCATCACCACCATATACCATTCAAATGCCCATCCCCATACATTGGTGACGTAACTGa   >  1:702291/1‑68 (MQ=255)
                                                               ccagccGAAAAGCATCACCACCATATACCATTCAAATGCCCATCCCCATACATTGGTGACGTAACTGa   >  1:1258951/1‑68 (MQ=255)
                                                               ccagccGAAAAGCATCACCACCATATACCATTCAAATGCCCATCCCCATACATTGGTGACGTAACTGa   >  1:2799/1‑68 (MQ=255)
                                                               ccagccGAAAAGCATCACCACCATATACCATTCAAATGCCCATCCCCATACATTGGTGACGTAACTGa   >  1:1155595/1‑68 (MQ=255)
                                                               ccagccGAAAAGCATCACCACCATATACCATTCAAATGCCCATCCCCATACATTGGTGACGTAACTGa   >  1:2699068/1‑68 (MQ=255)
                                                               ccagccGAAAAGCATCACCACCATATACCATTCAAATGCCCATCCCCATACATTGGTGACGTAACTGa   >  1:1516434/1‑68 (MQ=255)
                                                               ccagccGAAAAGCATCACCACCATATACCATTCAAATGCCCATCCCCATACATTGGTGACGTAACTGa   >  1:2379256/1‑68 (MQ=255)
                                                               ccagccGAAAAGCATCACCACCATATACCATTCAAATGCCCATCCCCATACATTGGTGACGTAACTGa   >  1:2243326/1‑68 (MQ=255)
                                                               ccaaccGAAAAGCATCACCACCATATACCATTCAAATGCCCATCCCCATACATTGGTGACGTAACTGaa  <  1:2915501/69‑1 (MQ=255)
                                                               ccaaccGAAAAGCATCACCACCATATACCATTCAAATGCCCATCCCCATACATTGGTGACGTAACTGaa  <  1:1635760/69‑1 (MQ=255)
                                                               ccaaccGAAAAGCATCACCACCATATACCATTCAAATGCCCATCCCCATACATTGGTGACGTAACTGaa  <  1:1859347/69‑1 (MQ=255)
                                                               ccaaccGAAAAGCATCACCACCATATACCATTCAAATGCCCATCCCCATACATTGGTGACGTAACTGaa  <  1:3418640/69‑1 (MQ=255)
                                                                  |                                                                 
GCTAAATTCTGGCGGTTCGTTACCTAAACGCTTTTTGGCATACGGGCCAAACACCAGCCAGAACCAACCGAAAAGCATCACCACCATATACCATTCAAATGCCCATCCCCATACATTGGTGACGTAACTGAA  >  minE/39349‑39480

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: