Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 217,977 C→T 21.6% intergenic (‑69/‑37) yafU ← / → yafF predicted inner membrane protein/conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE217,9770CT21.6% 224.4 / 55.8 111intergenic (‑69/‑37)yafU/yafFpredicted inner membrane protein/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (43/44);  new base T (9/15);  total (52/59)
Fisher's exact test for biased strand distribution p-value = 3.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.24e-01

Only 100 of 113 total aligned reads displayed.

ATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAA  >  minE/217910‑218041
                                                                   |                                                                
ataaATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGt                                                                <  1:2888312/70‑1 (MQ=255)
 taaATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTgg                                                              <  1:2942122/71‑1 (MQ=255)
         tcagTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAg                                                      >  1:1686334/4‑71 (MQ=255)
          cagTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGgag                                                            <  1:1418205/62‑1 (MQ=255)
           agTTCTGCTGATTTAACCGCTGAGAAATTCGCCAGAGCGATCCGAAATCACTGGCACg                                                                 <  1:3873982/57‑1 (MQ=255)
           agTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCAtgtg                                                               >  1:1451689/2‑60 (MQ=255)
             ttCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGc                                                  <  1:3013101/71‑1 (MQ=255)
              tCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGca                                                 >  1:2583378/1‑71 (MQ=255)
               ctgctgATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGcac                                                >  1:1189466/1‑71 (MQ=255)
                  ctgATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGcactg                                              <  1:918029/70‑3 (MQ=255)
                  ctgATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTgg                                             >  1:2072619/1‑71 (MQ=255)
                  ctgATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGcactg                                              <  1:2612697/70‑3 (MQ=255)
                  ctgATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGcactg                                              <  1:696136/70‑3 (MQ=255)
                   tgATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGc                                            >  1:2602757/1‑71 (MQ=255)
                    gATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGgaga                                                           <  1:1378176/55‑1 (MQ=255)
                      tttAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGc                                            <  1:410811/68‑1 (MQ=255)
                      tttAACCGCAGAGAAGTTCGCCACAGCAATCTGAAACCACTGGCATGTGGAGAATAAGCTg                                                   >  1:464516/1‑61 (MQ=255)
                      tttAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTc                                         <  1:3655820/71‑1 (MQ=255)
                       ttAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCt                                        >  1:964310/1‑71 (MQ=255)
                        tAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGga                                                             <  1:2104179/49‑1 (MQ=255)
                         aaCCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTg                                       >  1:3661896/1‑70 (MQ=255)
                         aaCCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTg                                       >  1:3740812/1‑70 (MQ=255)
                          aCCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCAc                                                                  >  1:2507046/1‑42 (MQ=255)
                          aCCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGa                                     <  1:1274479/71‑1 (MQ=255)
                          aCCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTgg                                                              >  1:2422969/1‑46 (MQ=255)
                           ccGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGAc                                    <  1:2699247/71‑1 (MQ=255)
                           ccGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGAc                                    >  1:2536068/1‑71 (MQ=255)
                           ccGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGAc                                    >  1:1126576/1‑71 (MQ=255)
                           ccGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGga                                                             >  1:1471115/1‑46 (MQ=255)
                           ccGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGAc                                    <  1:906744/71‑1 (MQ=255)
                            cGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGACg                                   <  1:3368884/71‑1 (MQ=255)
                             gCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTg                                                   <  1:2992773/54‑1 (MQ=255)
                              ctgagaAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAAt                                                         >  1:2266602/3‑47 (MQ=37)
                               agagAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTgg                                                              <  1:1674402/41‑1 (MQ=255)
                               agagAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGACGTgg                                >  1:473284/1‑71 (MQ=255)
                               agagAAGTTCGCCACAGCAATCCGAAACCACTGGCACg                                                                 <  1:1954758/38‑1 (MQ=255)
                               agagAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGca                                                 >  1:2096358/1‑54 (MQ=255)
                                gagaAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTCCACTGGCGTCTGGACGTgg                                >  1:3931825/1‑70 (MQ=255)
                                 aaaaGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTg                                                   <  1:294466/48‑1 (MQ=255)
                                 aaaaGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGt                               >  1:1028133/3‑70 (MQ=255)
                                 aaaaGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGt                               >  1:570190/3‑70 (MQ=255)
                                    aattCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTaatg                           <  1:1675544/69‑1 (MQ=255)
                                    aGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTaatg                           <  1:3521627/71‑1 (MQ=255)
                                    aGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTaatg                           <  1:1511678/71‑1 (MQ=255)
                                    aGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTaatg                           <  1:2857703/71‑1 (MQ=255)
                                     gTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTaatg                           >  1:124471/1‑70 (MQ=39)
                                     gTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTaatg                           >  1:2895760/1‑70 (MQ=255)
                                     attCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTaatg                           <  1:3447269/69‑1 (MQ=255)
                                        cGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGACGTg                                 >  1:3804863/1‑61 (MQ=38)
                                        cGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTaatgaat                        >  1:752732/1‑70 (MQ=255)
                                        cGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAg                                              >  1:1765062/1‑48 (MQ=255)
                                        cGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACg                                   <  1:3808671/59‑1 (MQ=255)
                                        cGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACg                                   <  1:1816668/59‑1 (MQ=255)
                                        cGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAACCTGCACAGGCGTCTGGACGTGGTaatgaat                        >  1:3490900/1‑70 (MQ=255)
                                         gCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTaatgaatg                       >  1:1425715/1‑70 (MQ=255)
                                         gCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAatgaatga                      <  1:522549/71‑1 (MQ=255)
                                          ccACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGt                               <  1:2608476/61‑1 (MQ=38)
                                          ccACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGt                               <  1:607213/61‑1 (MQ=38)
                                          ccACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTaatgaatg                       <  1:3806731/69‑1 (MQ=39)
                                          ccACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTaatgaat                        <  1:618958/68‑1 (MQ=39)
                                          ccACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTaa                             <  1:1809217/63‑1 (MQ=255)
                                          ccACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAAtgaatgaa                     <  1:3257089/71‑1 (MQ=255)
                                           cacaGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAg                    <  1:1179097/71‑1 (MQ=255)
                                           cacaGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAg                    <  1:3904199/71‑1 (MQ=255)
                                           cacaGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGa                                     >  1:876085/1‑54 (MQ=255)
                                            acaGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGCAGTGGTAATGAATGAAga                   >  1:3267468/1‑71 (MQ=255)
                                            acaGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGt                               <  1:2976693/59‑1 (MQ=255)
                                             caGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTgg                                             >  1:281643/1‑44 (MQ=38)
                                             caGCAATCGGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAgac                  >  1:392748/1‑71 (MQ=255)
                                             caGCAATCCGAAACCCCTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAgac                  >  1:1126792/1‑71 (MQ=255)
                                             caGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGAc                                    <  1:252004/53‑1 (MQ=38)
                                             caGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTaatgaat                        >  1:337863/1‑65 (MQ=255)
                                             caGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTaa                             <  1:329316/60‑1 (MQ=255)
                                             caGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTaa                             <  1:1757602/60‑1 (MQ=255)
                                                cAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAgacga                <  1:569411/70‑1 (MQ=255)
                                                cAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTa                              >  1:945283/1‑56 (MQ=255)
                                                   tCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAg                    <  1:522824/63‑1 (MQ=255)
                                                   tCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTg                                 >  1:1150609/1‑50 (MQ=255)
                                                   tCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCt                                        <  1:194261/43‑1 (MQ=255)
                                                    ccGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAAtgaatgaa                     >  1:3422254/1‑61 (MQ=255)
                                                     cGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTgg                                      >  1:599143/1‑43 (MQ=255)
                                                      gAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCaaa         <  1:3322932/71‑1 (MQ=39)
                                                       aaaTCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCaaaa        >  1:173946/1‑71 (MQ=255)
                                                        aaCCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCaaaa        >  1:726669/1‑70 (MQ=39)
                                                        aaCCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTg                                       >  1:2709508/1‑39 (MQ=255)
                                                        aaCCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAg                    <  1:877244/58‑1 (MQ=255)
                                                        aaCCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACaaaa        >  1:2347406/1‑70 (MQ=255)
                                                        aaCCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACaaaa        >  1:2559308/1‑70 (MQ=255)
                                                         aCCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACaaaa        >  1:215781/1‑69 (MQ=255)
                                                         aCCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAAt       >  1:3581221/1‑70 (MQ=255)
                                                            aCTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAAtgaatgaa                     <  1:2103144/53‑1 (MQ=255)
                                                            aCTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATaag    >  1:1165177/1‑70 (MQ=255)
                                                             cTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTaatgaa                         <  1:852031/48‑1 (MQ=255)
                                                             cTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTaatgaa                         <  1:3582816/48‑1 (MQ=255)
                                                              tGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTaa                             <  1:2315908/43‑1 (MQ=38)
                                                              tGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAAtgaatgaa                     <  1:129343/51‑1 (MQ=255)
                                                              tGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACaaaa        >  1:2237008/1‑64 (MQ=255)
                                                               ggCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATaagaa  >  1:3805089/1‑69 (MQ=255)
                                                               ggCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATaagaa  >  1:3800580/1‑69 (MQ=255)
                                                                gCACGTGGAGAATAAGCTGCACTGGCGGCTGGACGTGGTAATGAATGAAgac                  <  1:2671612/52‑1 (MQ=255)
                                                                   |                                                                
ATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAA  >  minE/217910‑218041

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: