Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,347,119 T→C 7.7% G204G (GGT→GGC yegH → fused predicted membrane proteins

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,347,1190TC 7.7% 212.0 / 7.5 78G204G (GGT→GGCyegHfused predicted membrane proteins
Reads supporting (aligned to +/- strand):  ref base T (37/35);  new base C (2/4);  total (39/39)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.72e-01

GTCACCCGACGATCGTTATTCTCTGCTTAAGCTTCCTGTTGATGATTGGCTTTAGCCTGGTGGCAGAAGGTTTCGGCTTCGTCATTCCGAAAGGCTACCTGTACGCTGCCATTGGTTTCTCGGTGATGATCGAGGCGCTCA  >  minE/1347049‑1347189
                                                                      |                                                                      
gTCACCCGACGATCGTTATTCTCTGCTTAAGCTTCCTGTTGATGATTGGCTTTAGCCTGGTGGCAGAAGGt                                                                        <  1:2595673/71‑1 (MQ=255)
gTCACCCGACGATCGTTATTCTCTGCTTAAGCTTCCTGTTGATGATTGGCTTTAGCCTGGTGGCAGAAGGt                                                                        <  1:621652/71‑1 (MQ=255)
gTCACCCGACGATCGTTATTCTCTGCTTAAGCTTCCTGTTGATGATTGGCTTTAGCCTGGTGGCAGAAGGt                                                                        <  1:2925056/71‑1 (MQ=255)
 tCACCCGACGATCGTTATTCTCTGCTTAAGCTTCCTGTTGATGATTGGCTTTAGCCTGGTGGCAGAAGGtt                                                                       <  1:552057/71‑1 (MQ=255)
  cACCCGACGATCGTTATTCTCTGCTTAAGCTTCCTGTTGATGATTGGCTTTAGCCTGGTGGCAGAAGGtt                                                                       >  1:3560786/1‑70 (MQ=255)
   aCCCGACGATCGTTATTCTCTGCTTAAGCTTCCTGTTGATGATTGGCTTTAGCCTGGTGGCAGAAGGTTTc                                                                     >  1:2338604/1‑71 (MQ=255)
     ccGACGATCGTTATTCTCTGCTTAAGCTTCCTGTTGATGATTGGCTTTAGCCTGGTGGCAGAAGGTTTCg                                                                    <  1:3866394/70‑1 (MQ=255)
      cgacgaTCGTTATTCTCTGCTTAAGCTTCCTGTTGATGATTGGCTTTAGCCTGGTGGCAGAAGGTTTCGGc                                                                  >  1:3186732/1‑71 (MQ=255)
      cgacgaTCGTTATTCTCTGCTTAAGCTTCCTGTTGATGATTGGCTTTAGCCTGGTGGCAGAAGGTTTCGGc                                                                  >  1:3438124/1‑71 (MQ=255)
       gacgaTCGTTATTCTCTGCTTAAGCTTCCTGTTGATGATTGGCTTTAGCCTGGTGGCAGAAGGTTTCGGCt                                                                 >  1:2305353/1‑71 (MQ=255)
       gacgaTCGTTATTCTCTGCTTAAGCTTCCTGTTGATGATTGGCTTTAGCCTGGTGGCAGAAGGTTTCGGCt                                                                 <  1:3913319/71‑1 (MQ=255)
        acgaTCGTTATTCTCTGCTTAAGCTTCCTGTTGATGATTGGCTTTAGCCTGGTGGCAGAAGGTTTCGGCtt                                                                <  1:1173038/71‑1 (MQ=255)
         cgaTCGTTATTCTCTGCTTAAGCTTCCTGTTGATGATTGGCTTTAGCCTGGTGGCAGAAGGTTTCGGCTTc                                                               <  1:502480/71‑1 (MQ=255)
           aTCGTTATTCTCTGCTTAAGCTTCCTGTTGATGATTGGCTTTAGCCTGGTGGCAGAAGGTTTCGGCTTCGt                                                             <  1:1300383/71‑1 (MQ=255)
             cGTTATTCTCTGCTTAAGCTTCCTGTTGATGATTGGCTTTAGCCTGGTGGCAGAAGGTTTCg                                                                    <  1:2102827/62‑1 (MQ=255)
              gTTATTCTCTGCTTAAGCTTCCTGTTGATGATTGGCTTTAGCCTGGTGGCAGAAGGTTTCGGCTTCGTCa                                                           <  1:2501256/70‑1 (MQ=255)
              gTTATTCTCTGCTTAAGCTTCCTGTTGATGATTGGCTTTAGCCTGGTGGCAGAAGGTTTCGGCTTCGTCAt                                                          <  1:329443/71‑1 (MQ=255)
                tATTCTCTGCTTAAGCTTCCTGTTGATGATTGGCTTTAGCCTGGTGGCAGAAGGTTTCGGCTTCGTCATTc                                                        >  1:261382/1‑71 (MQ=255)
                 aTTCTCTGCTTAAGCTTCCTGTTGATGATTGGCTTTAGCCTGGTGGCAGAAGGTTTCGGCTTCGTCATTcc                                                       >  1:2160373/1‑71 (MQ=255)
                    ctctGCTTAAGCTTCCTGTTGATGATTGGCTTTAGCCTGGTGGCAGAAGGTTTCGGCTTCGTCATTCCGa                                                     >  1:1129357/1‑70 (MQ=255)
                     tctGCTTAAGCTTCCTGTTGATGATTGGCTTTAGCCTGGTGGCAGAAGGTTTCGGCTTc                                                               >  1:23943/1‑59 (MQ=255)
                     tctGCTTAAGCTTCCTGTTGATGATTGGCTTTAGCCTGGTGGCAGAAGGTTTCGGCTTCGTCa                                                           <  1:3707386/63‑1 (MQ=255)
                     tctGCTTAAGCTTCCTGTTGATGATTGGCTTTAGCCTGGTGGCAGAAGGTTTCGGCTTCGTCa                                                           >  1:386258/1‑63 (MQ=255)
                     tctGCTTAAGCTTCCTGTTGATGATTGGCTTTAGCCTGGTGGCAGAAGGTTTCGGCTTCGTCATTCCGa                                                     >  1:3770570/1‑69 (MQ=255)
                     tctGCTTAAGCTTCCTGTTGATGATTGGCTTTAGCCTGGTGGCAGAAGGTTTCGGCTTCGTCATTCCGa                                                     >  1:2286426/1‑69 (MQ=255)
                       tGCTTAAGCTTCCTGTTGATGATTGGCTTTAGCCTGGTGGCAGAAGGTTTCGGCTTCg                                                              <  1:1304697/58‑1 (MQ=255)
                         cTTAAGCTTCCTGTTGATGATTGGCTTTAGCCTGGTGGCAGAAGGt                                                                        <  1:356861/46‑1 (MQ=255)
                         cTTAAGCTTCCTGTTGATGATTGGCTTTAGCCTGGTGGCAGAAGGTTTCGGCTTCGTCATTCCGAAAGGCt                                               <  1:2572179/71‑1 (MQ=255)
                         cTTAAGCTTCCTGTTGATGATTGGCTTTAGCCTGGTGGC‑‑AAGGTTTCGGCTTCGTCATTCCGAAAGGCt                                               <  1:657006/69‑1 (MQ=25)
                          ttAAGCTTCCTGTTGATGATTGGCTTTAGCCTGGTGGCAGAAGGTTTCGGc                                                                  <  1:2287463/51‑1 (MQ=255)
                          ttAAGCTTCCTGTTGATGATTGGCTTTAGCCTGGTGGCAGAAGGTTTCGGc                                                                  <  1:3797905/51‑1 (MQ=255)
                            aaGCTTCCTGTTGATGATTGGCTTTAGCCTGGTGGCAGAAGGTTTCGGCTTCGTCATTCCGa                                                     >  1:3869031/1‑62 (MQ=255)
                            aaGCTTCCTGTTGATGATTGGCTTTAGCCTGGTGGCAGAAGGTTTCGGCTTCGTCATTCCGAAAGGCTAcc                                            >  1:1252503/1‑71 (MQ=255)
                                ttCCTGTTGATGATTGGCTTTAGCCTGGTGGCAGAAGGTTTCGGCTTCGTCAt                                                          <  1:3560445/53‑1 (MQ=255)
                                ttCCTGTTGATGATTGGCTTTAGCCTGGTGGCAGAAGGTTTCGGCTTCGTCATTCCGaa                                                    <  1:620814/59‑1 (MQ=255)
                                  ccTGTTGATGATTGGCTTTAGCCTGGTGGCAGAAGGTTTCGGCTTCGTCa                                                           >  1:2204512/1‑50 (MQ=255)
                                   cTGTTGATGATTGGCTTTAGCCTGGTGGCAGAAGGTTTCGGCTTCGTCATTCCGAAAGGCTACCTGTACGc                                     >  1:1284481/1‑71 (MQ=255)
                                      ttGATGATTGGCTTTAGCCTGGTGGCAGAAGGTTTCGGCTTCGTCATTCCGAAAg                                                  >  1:3196038/1‑55 (MQ=255)
                                      ttGATGATTGGCTTTAGCCTGGTGGCAGAAGGTTTCGGCTTCGTCATTCCGAAAGGCTAc                                             >  1:1356645/1‑60 (MQ=255)
                                        gatgatTGGCTTTAGCCTGGTGGCAGAAGGTTTCGGCTTCGTCATTCCGAAAGGCTACCTGTACGCTGcc                                 <  1:3757355/70‑1 (MQ=255)
                                         atgatTGGCTTTAGCCTGGTGGCAGAAGGTTTCGGCTTCGTCATTCCGAAAGGCTACCTGTACGCTGCCa                                <  1:3861128/70‑1 (MQ=255)
                                         atgatTGGCTTTAGCCTGGTGGCAGAAGGTTTCGGCTTCGTCATTCCGAAAGGCTACCTGTACGCTGCCa                                <  1:691659/70‑1 (MQ=255)
                                         atgatTGGCTTTAGCCTGGTGGCAGAAGGTTTCGGCTTCGTCATTCCGAAAGGCTACCTGTACGCTGCCAt                               >  1:2789107/1‑71 (MQ=255)
                                         atgatTGGCTTTAGCCTGGTGGCAGAAGGTTTCGGCTTCGTCATTCCGAAAGGCTACCTGTACGCTGCCAt                               >  1:1146873/1‑71 (MQ=255)
                                         atgatTGGCTTTAGCCTGGTGGCAGAAGGTTTCGGCTTCGTCATTCCGAAAGGCTAACTGTACGCTGCCa                                <  1:829506/70‑1 (MQ=255)
                                          tgatTGGCTTTAGCCTGGTGGCAGAAGGTTTCGGCTTCGTCATTCCGaa                                                    >  1:2526837/1‑49 (MQ=255)
                                          tgatTGGCTTTAGCCTGGTGGCAGAAGGTTTCGGCTTCGTCATTCCGa                                                     <  1:2977026/48‑1 (MQ=255)
                                            atTGGCTTTAGCCTGGTGGCAGAAGGTTTCGGCTTCGTCATTCCGAAAGGCTACCTGTACGCTGCCATTgg                            >  1:2878652/1‑71 (MQ=255)
                                             ttGGCTTTAGCCTGGTGGCAGAAGGTTTCGGCTTCGTCATTCCGaaa                                                   >  1:623049/1‑47 (MQ=255)
                                             ttGGCTTTAGCCTGGTGGCAGAAGGTTTCGGCTTCGTCATTCCGaaa                                                   >  1:433098/1‑47 (MQ=255)
                                             ttGGCTTTAGACTGGTGGCAGAAGGTTTCGGCTTCGTCATTCCGaaa                                                   >  1:3280768/1‑47 (MQ=255)
                                               ggCTTTAGCCTGGTGGCAGAAGGTTTCGGCTTCGTCATTCCGAAAGGCTa                                              <  1:1350606/50‑1 (MQ=255)
                                                 cTTTAGCCTGGTGGCAGAAGGTTTCGGCTTCGTCATTCCGAAAGGCTACCTGTACGCTGCCATTTGTTtct                       <  1:3901740/71‑1 (MQ=255)
                                                 cTTTAGCCTGGTGGCAGAAGGCTTCGGCTTCGTCATTCCGAAAGGCTACCTGTACGCTGACATTGGTTtct                       <  1:1358970/71‑1 (MQ=255)
                                                   ttAGCCTGGTGGCAGAAGGTTTCGGCTTCGTCATTCCGaaa                                                   >  1:1374977/1‑41 (MQ=255)
                                                   ttAGCCTGGTGGCAGAAGGTTTCGGCTTCGTCATTCCGAAAGGCTACCTGTAc                                       >  1:481390/1‑53 (MQ=255)
                                                   ttAGCCTGGTGGCAGAAGGTTTCGGCTTCGTCATTCCGAAAGGCTACCTGTACGCTGCCATTGGTTTCTCg                     >  1:525411/1‑71 (MQ=255)
                                                    tAGCCTGGTGGCAGAAGGTTTCGGCTTCGTCATTCCGAAAGGCTACCTGTACGCTGCCATTGGTTTCTCgg                    >  1:2368424/1‑71 (MQ=255)
                                                     aGCCTGGTGGCAGAAGGTTTCGGCTTCGTCATTCCGAAAGGCTACCTGTACGc                                     >  1:2170769/1‑53 (MQ=255)
                                                     aGCCTGGTGGCAGAAGGTTTCGGCTTCGTCATTCCGAAAGGCTACCTGTACGCTg                                   >  1:1810181/1‑55 (MQ=255)
                                                     aGCCTGGTGGCAGAAGGTTTCGGCTTCGTCATTCCGAAAGGCTACCTGTACGCTGCCATTGGTTTCTCGGt                   >  1:2815506/1‑71 (MQ=255)
                                                      gCCTGGTGGCAGAAGGCTTCGGCTTCGTCATTCCGAAAGGCTACCTGTACGCTGCCATTGGTTTCTCGGtg                  <  1:1174674/71‑1 (MQ=255)
                                                      gCCTGGTGGCAGAAGGCTTCGGCTTCGTCATTACGAAAGGCTACCTGTACGCTGCCATTGGTTTCTCGGtg                  <  1:3835723/71‑1 (MQ=255)
                                                       ccTGATCGCCGAAGGCCTGGGCTTCCACATCCCGAAGGGCTACCTGTACGccg                                   >  1:3201494/1‑51 (MQ=255)
                                                        cTGGTGGCAGAAGGTTTCGGCTTCGTCATTCCGAAAGGCt                                               <  1:2602111/40‑1 (MQ=255)
                                                         tggtggCAGAAGGTTTCGGCTTCGTCATTCCGAAAGGCTACCTGTACGCTGCCATTGGTTTCTCGGtgat                >  1:2271364/1‑70 (MQ=255)
                                                           gtggCAGAAGGTTTCGGCTTCGTCATTCCGAAAGGCTACCTGTACGCTGCCATTGGTTtc                        <  1:3056754/60‑1 (MQ=255)
                                                           gtggCAGAAGGTTTCGGCTTCGTCATTCCGAAAGGCTACCTGTACGCTGCCATTGGTTTCTCGGTgatgat             <  1:3412429/71‑1 (MQ=255)
                                                            tggCAGAAGGTTTCGGCTTCGTCATTCCGAAAGGCTACCTGTACGCt                                    <  1:2449305/47‑1 (MQ=255)
                                                            tggCAGAAGGTTTCGGCTTCGTCATTCCGAAAGGCTACCTGTACGCTGGCATTGGTTTCTCGGtgatg               <  1:290814/68‑1 (MQ=255)
                                                            tggCAGAAGGTTTCGGCTTCGTCATTCCGAAAGGCTACCTGTACGCTGCCATTGGTTTCTCGGtgatga              <  1:3908761/69‑1 (MQ=255)
                                                             ggCAGAAGGTTTCGGCTTCGTCATTCCGAAAGGCTACCTGTACGCTGCCATTGGTTTCTCGGTGATGATc            >  1:3325274/1‑70 (MQ=255)
                                                              gCAGAAGGTTTCGGCTTCGTCATTCCGAAAGGCTACCTGTACGCTGCCATTGGTTTCTCGGTGATGATc            >  1:1729369/1‑69 (MQ=255)
                                                              gCAGAAGGTTTCGGCTTCGTCATTCCGAAAGGCTACCTGTACGCTGCCATTGGTTTCTCGGTGATGATCg           <  1:3550940/70‑1 (MQ=255)
                                                              gCAGAAGGTTTCGGCTTCGTCATTCCGAAAGGCTACCTGTACGCTGCCATTGGTTTCTCGGTGATGATCGa          <  1:3189291/71‑1 (MQ=255)
                                                              gCAGAAGGCTTCGGCTTCGTCATTCCGAAAGGCTACCTGTACGCTGCCATTGGTTTCTCGGTGATGATCg           <  1:3249429/70‑1 (MQ=255)
                                                                aGAAGGTTTCGGCTTCGTCATTCCGAAAGGCTACCTGTACGc                                     <  1:1900937/42‑1 (MQ=255)
                                                                 gAAGGTTTCGGCTTCGTCATTCCGAAAGGCTACCTGTACGCTGCCATTGGTTTCTCGGTGATGATCGAGgc       <  1:2505729/71‑1 (MQ=255)
                                                                  aaGGTTTCGGCTTCGTCATTCCGAAAGGCTACCTGTACGCTGCCATTGGTTTCTCGGTGATGATCGAGgcg      >  1:2161285/1‑71 (MQ=255)
                                                                  aaGGTTTCGGCTTCGTCATTCCGAAAGGCTACCTGTACGCTGCCATTGGTTTCTCGGTGATGATCGAGgcg      <  1:2955391/71‑1 (MQ=255)
                                                                  aaGGTTTCGGCTTCGTCATTCCGAAAGGCTACCTGTACGCTGCCATTGGTTTCTCGGTGATGATCGAGgcg      <  1:3860722/71‑1 (MQ=255)
                                                                  aaGGTTTCGGCTTCGTCATTCCGAAAGGCTACCTGTACGCTGCCATTGGTTTCTCGGTGATGATCGAGgcg      >  1:3129736/1‑71 (MQ=255)
                                                                   aGGCTTCGGCTTCGTCATTCCGAAAGGCTACCTGTACGCTGCCATTGGTTTCTCGGTGATGATCGAGgcg      >  1:3735471/1‑70 (MQ=255)
                                                                    ggTTTCGGCTTCGTCATTCCGAAAGGCTACCTGTACGCTGCCATTGGTTTCTCGGTGATGATCGAGGCGCt    >  1:3678973/1‑71 (MQ=255)
                                                                      tttCGGCTTCGTCATTCCGAAAGGCTACCTGTACGCTGCCATTGGTTTCTCGGTGATGATCGAGGCGCTCa  <  1:1267586/71‑1 (MQ=255)
                                                                      |                                                                      
GTCACCCGACGATCGTTATTCTCTGCTTAAGCTTCCTGTTGATGATTGGCTTTAGCCTGGTGGCAGAAGGTTTCGGCTTCGTCATTCCGAAAGGCTACCTGTACGCTGCCATTGGTTTCTCGGTGATGATCGAGGCGCTCA  >  minE/1347049‑1347189

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: