Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 217,996 A→T 64.7% intergenic (‑88/‑18) 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,9960AT64.7% 103.6 / 104.9 102intergenic (‑88/‑18)yafU/yafFpredicted inner membrane protein/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (21/15);  new base T (35/31);  total (56/46)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.43e-01

Only 100 of 106 total aligned reads displayed.

CTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTT  >  minE/217928‑218065
                                                                    |                                                                     
ctgATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTgg                                                                     >  1:2072619/1‑71 (MQ=255)
 tgATTTAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGc                                                                    >  1:2602757/1‑71 (MQ=255)
    tttAACCGCTGAGAAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGc                                                                    <  1:410811/68‑1 (MQ=255)
    tttAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTc                                                                 <  1:3655820/71‑1 (MQ=255)
       aaCCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTg                                                               >  1:3661896/1‑70 (MQ=255)
       aaCCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTg                                                               >  1:3740812/1‑70 (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)
              gagaAATTCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTCCACTGGCGTCTGGACGTgg                                                        >  1:3931825/1‑70 (MQ=255)
              gaaaaGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGt                                                       <  1:1641344/68‑1 (MQ=255)
               agaAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTgg                                                        >  1:2536499/1‑69 (MQ=255)
               aaaaGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGt                                                       >  1:1028133/3‑70 (MQ=255)
               aaaaGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGt                                                       >  1:570190/3‑70 (MQ=255)
                  aattCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTaatg                                                   <  1:1675544/69‑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:607213/61‑1 (MQ=38)
                        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)
                        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)
                           caGCAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAgac                                          >  1:794802/1‑71 (MQ=39)
                           caGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTaatgaat                                                >  1:337863/1‑65 (MQ=255)
                           caGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTaa                                                     <  1:1757602/60‑1 (MQ=255)
                           caGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTaa                                                     <  1:329316/60‑1 (MQ=255)
                            aGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAgac                                          >  1:1256877/1‑70 (MQ=255)
                              cAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAgacga                                        <  1:569411/70‑1 (MQ=255)
                              cAATCCGAAACCACTGGCATGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAgacga                                        <  1:341886/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)
                                  ccGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAAtgaatgaa                                             >  1:3422254/1‑61 (MQ=255)
                                  ccGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAgacg                                         <  1:3721768/65‑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)
                                      aaCCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTg                                                               >  1:2709508/1‑39 (MQ=255)
                                      aaCCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAg                                            <  1:388652/58‑1 (MQ=255)
                                      aaCCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAg                                            <  1:877244/58‑1 (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:3800580/1‑69 (MQ=255)
                                             ggCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATaagaa                          >  1:3805089/1‑69 (MQ=255)
                                              gCACGTGGAGAATAAGCTGCACTGGCGGCTGGACGTGGTAATGAATGAAgac                                          <  1:2671612/52‑1 (MQ=255)
                                                atgTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAgg                      >  1:1617337/3‑70 (MQ=255)
                                                atgTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGAAGACTGCAAAATAAGAAGAGGa                     <  1:1110959/69‑1 (MQ=255)
                                                 tgTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATaag                            >  1:2648226/2‑63 (MQ=255)
                                                 tgTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGa                     >  1:1323050/2‑70 (MQ=255)
                                                  gTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGaaa                   >  1:253522/1‑71 (MQ=255)
                                                   tGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTa                                                      <  1:3335977/35‑1 (MQ=255)
                                                     gagaATAAGCTGCACAGGCGTCTGGACGTGGTAAtgaatgaa                                             >  1:2505680/1‑42 (MQ=255)
                                                      agaATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAgacgac                                       >  1:3021510/1‑47 (MQ=255)
                                                        aaTAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACtg                                     >  1:1004400/1‑46 (MQ=255)
                                                        aaTAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAgac                                          <  1:1777659/42‑1 (MQ=255)
                                                             gCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAat        <  1:2540746/71‑1 (MQ=255)
                                                              cTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAat        >  1:1681841/1‑70 (MQ=255)
                                                               tGCACTGGCGTCTGGACGTGGTAATGAATGAAgacga                                        >  1:1035807/1‑37 (MQ=255)
                                                               tGCACTGGCGTCTGGACGTGGTAATGAATGAAGCCGACTGCAAAATAAGAAGAGGAAACGCCGCAGAatta      >  1:2413801/1‑71 (MQ=255)
                                                               tGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAatta      >  1:3952469/1‑71 (MQ=255)
                                                               tGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAatt       >  1:1169153/1‑70 (MQ=255)
                                                                 cacaGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAattat     <  1:616730/70‑1 (MQ=255)
                                                                   caGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTAtttt  <  1:1367906/71‑1 (MQ=255)
                                                                    aGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGcagcag           <  1:3086400/61‑1 (MQ=255)
                                                                    aGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTAtttt  >  1:2538903/1‑70 (MQ=255)
                                                                    |                                                                     
CTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTT  >  minE/217928‑218065

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: