Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 217,968 C→T 14.7% intergenic (‑60/‑46) 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,9680CT14.7% 292.1 / 33.2 115intergenic (‑60/‑46)yafU/yafFpredicted inner membrane protein/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (44/54);  new base T (9/8);  total (53/62)
Fisher's exact test for biased strand distribution p-value = 6.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 116 total aligned reads displayed.

AACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAAT  >  minE/217902‑218036
                                                                  |                                                                    
aaCTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCa                                                                    <  1:2888069/69‑1 (MQ=255)
aaCTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCAc                                                                   <  1:3761061/70‑1 (MQ=255)
aaCTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCAc                                                                   <  1:2663366/70‑1 (MQ=255)
aaCTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCAc                                                                   <  1:424125/70‑1 (MQ=255)
aaCTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACt                                                                  >  1:2458313/1‑71 (MQ=255)
aaCTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACt                                                                  >  1:724093/1‑71 (MQ=255)
  cTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTgg                                                                <  1:1390659/71‑1 (MQ=255)
  cTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTgg                                                                <  1:2077113/71‑1 (MQ=255)
    cATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCAc                                                                   <  1:245371/66‑1 (MQ=255)
         taaATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTgg                                                         <  1:2942122/71‑1 (MQ=255)
                 tcagTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAg                                                 >  1:1686334/4‑71 (MQ=255)
                  tGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATa                                                   <  1:2987154/68‑1 (MQ=255)
                  cagTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGgag                                                       <  1:1418205/62‑1 (MQ=255)
                   ggTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGaa                                                     <  1:550551/65‑1 (MQ=255)
                   agTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCAtgtg                                                          >  1:1451689/2‑60 (MQ=255)
                     ttCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAAc                                                                      >  1:2779023/1‑46 (MQ=255)
                     ttCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGc                                             <  1:3013101/71‑1 (MQ=255)
                     ttCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTgg                                                                <  1:3736631/52‑1 (MQ=255)
                      tCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGca                                            >  1:2583378/1‑71 (MQ=255)
                       ctgctgATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGcac                                           >  1:1189466/1‑71 (MQ=255)
                         gctgTTTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGc                                                               <  1:2825032/49‑1 (MQ=255)
                          ctgATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGcactg                                         <  1:918029/70‑3 (MQ=255)
                          ctgATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGcactg                                         <  1:2612697/70‑3 (MQ=255)
                          ctgATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGcactg                                         <  1:696136/70‑3 (MQ=255)
                           tgATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGc                                       >  1:2602757/1‑71 (MQ=255)
                           tgATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCa                                                                    >  1:920853/1‑42 (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)
                                   ccGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGAc                               <  1:2699247/71‑1 (MQ=255)
                                   ccGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGAc                               >  1:1126576/1‑71 (MQ=255)
                                   ccGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGAc                               >  1:2536068/1‑71 (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)
                                       agagAAGTTCGCCACAGCAATCCGAAACCACTGGCACg                                                            <  1:1954758/38‑1 (MQ=255)
                                       agagAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGca                                            >  1:2096358/1‑54 (MQ=255)
                                        gagaAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTCCACTGGCGTCTGGACGTgg                           >  1:3931825/1‑70 (MQ=255)
                                        gaaaaGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGt                          <  1:1641344/68‑1 (MQ=255)
                                         agaAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTgg                           >  1:2536499/1‑69 (MQ=255)
                                         aaaaGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTg                                              <  1:2232146/48‑1 (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)
                                             attCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCt                                   >  1:2721216/2‑57 (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)
                                                  ccACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTaatgaatg                  <  1:3806731/69‑1 (MQ=39)
                                                  ccACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTaatgaat                   <  1:618958/68‑1 (MQ=39)
                                                  ccACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTaa                        <  1:1809217/63‑1 (MQ=255)
                                                   cacaGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAg               <  1:3904199/71‑1 (MQ=255)
                                                   cacaGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAg               <  1:1179097/71‑1 (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)
                                                     caGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAgac             >  1:794802/1‑71 (MQ=39)
                                                     caGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTaatgaat                   >  1:337863/1‑65 (MQ=255)
                                                     caGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTaa                        <  1:329316/60‑1 (MQ=255)
                                                        cAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAgacga           <  1:341886/70‑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)
                                                           tCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGCCGTCTGGACGTGGTaatgaat                   <  1:183286/59‑1 (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)
                                                                aaCCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAg               <  1:877244/58‑1 (MQ=255)
                                                                aaCCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAg               <  1:388652/58‑1 (MQ=255)
                                                                aaCCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACaaaa   >  1:2347406/1‑70 (MQ=255)
                                                                 aCCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACaaaa   >  1:215781/1‑69 (MQ=255)
                                                                 aCCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAAt  >  1:3581221/1‑70 (MQ=255)
                                                                  |                                                                    
AACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAAT  >  minE/217902‑218036

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: