Predicted mutation
evidence seq id position mutation annotation gene description
RA minE 1,252,050 T→C K6K (AAA→AAG yebC ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,252,0500TC100.0% 222.9 / NA 63K6K (AAA→AAGyebCconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (29/34);  total (29/34)

ATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTAC  >  minE/1251982‑1252113
                                                                    |                                                               
aTGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCActt                                                               <  1:2117102/71‑1 (MQ=255)
 tGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCACTTa                                                              >  1:1670783/1‑71 (MQ=255)
   aTTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCACTTACt                                                            <  1:1027056/71‑1 (MQ=255)
   aTTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCACTTACt                                                            <  1:1424727/71‑1 (MQ=255)
    ttttAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCACTTACTa                                                           >  1:2895269/1‑71 (MQ=255)
    ttttAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCACTTACTa                                                           >  1:2877071/1‑71 (MQ=255)
       tAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCACTTACTATGa                                                        >  1:2633796/1‑71 (MQ=255)
       tAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCACTTACTATGa                                                        >  1:1286770/1‑71 (MQ=255)
         gTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCACTTACTATGAcc                                                      <  1:1024823/71‑1 (MQ=255)
            aaGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCACTTACTATGACCTGc                                                   >  1:909784/1‑71 (MQ=255)
              gATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCACTTACTATg                                                         <  1:1347073/63‑1 (MQ=255)
              gATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCACTTACTATg                                                         <  1:1924819/63‑1 (MQ=255)
              gATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCACTTACTATGACCTGCCa                                                 >  1:1657697/1‑71 (MQ=255)
              gATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCACTTACTATGACCTGCCa                                                 >  1:1736791/1‑71 (MQ=255)
                ttttACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCACTTACTATGACCTGCCAt                                                >  1:2479775/1‑70 (MQ=255)
                 tttACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCACTTACTATGACCTg                                                    >  1:1310831/1‑65 (MQ=255)
                     ccGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCACTTACTATGACCTGCCATaaaaa                                           <  1:2120824/70‑1 (MQ=255)
                      cgcgCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCACTTACTATGACCTGCCATAAAAata                                         <  1:1606802/71‑1 (MQ=255)
                          cTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCACTTACTATGACCTGCCATAAAAATAtctc                                     <  1:2485448/71‑1 (MQ=255)
                          cTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGCCCCACTTACTATGACCTGCCATAAAAATAtctc                                     <  1:2105050/71‑1 (MQ=255)
                           ttAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCACTTACTATGACCTGCCATAAAAATAtct                                      <  1:1252698/69‑1 (MQ=255)
                             aGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCACTTACTATGACCTGCCATAAAAATATCTcc                                    <  1:1504848/69‑1 (MQ=255)
                             aGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCACTTACTATGACCTGCCATAAAAATATCTcc                                    <  1:2425423/69‑1 (MQ=255)
                             aGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCACTTACTATGACCTGCCATAAAAATATCTCCa                                   <  1:2592948/70‑1 (MQ=255)
                             aGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCACTTACTATGACCTGCCATAAAAATATCTCCa                                   <  1:608283/70‑1 (MQ=255)
                              gCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCACTTACTATGACCTGCCATAAAAATAtctc                                     >  1:1511908/1‑67 (MQ=255)
                              gCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCACTTACTATGACCTGCCATAAAAATATCTCCAg                                  >  1:1885861/1‑70 (MQ=255)
                              gCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCACTTACTATGACCTGCCATAAAAATATCTCCAGa                                 <  1:2759517/71‑1 (MQ=255)
                              gCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCACTTACTATGACCTGCCATAAAAATATCTCCAGa                                 <  1:1076168/71‑1 (MQ=255)
                               cATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCACTTACTATGAc                                                       >  1:479241/1‑48 (MQ=255)
                                aTCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCACTTACTATGACCTGCCATAAAAATAtctc                                     <  1:1684142/65‑1 (MQ=255)
                                  ccTGCGCAGCTTTACGATGTCTGGTGTTGGCCCACTTACTATGACCTGCCATAAAAATATCTCCAGATAGc                             >  1:1264804/1‑71 (MQ=255)
                                   cTGCGCAGCTTTACGATGTCTGGTGTTGGCCCACTTACTATGACCTGCCATaaaaa                                           >  1:1715972/1‑56 (MQ=255)
                                   cTGCGCAGCTTTACGATGTCTGGTGTTGGCCCACTTACTATGACCTGCCATaaaaa                                           >  1:1606987/1‑56 (MQ=255)
                                    tGCGCAGCTTTACGATGTCTGGTGTTGGCCCACTTACTATg                                                         <  1:1917737/41‑1 (MQ=255)
                                    tGCGCAGCTTTACGATGTCTGGTGTTAGCCCACTTACTATg                                                         <  1:1869938/41‑1 (MQ=255)
                                     gcgcAGCTTTACGATGTCTGGTGTTGGCCCACTTACTATGACCTGCCATAAAAat                                          >  1:2152474/1‑55 (MQ=255)
                                      cgcAGCTTTACGATGTCTGGTGTTGGCCCACTTACTa                                                           <  1:1653548/37‑1 (MQ=255)
                                      cgcAGCTTTACGATGTCTGGTGTTGGCCCACTTACTATGAc                                                       >  1:1662249/1‑41 (MQ=255)
                                        cAGCTTTACGATGTCTGGTGTTGGCCCACTTACTATGACCTGCCATaa                                              <  1:2130672/48‑1 (MQ=255)
                                           cTTTACGATGTCTGGTGTTGGCCCACTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGt                    <  1:1281609/71‑1 (MQ=255)
                                             ttACGATGTCTGGTGTTGGCCCACTTACTATGACCTGCCATAAAAat                                          <  1:1040918/47‑1 (MQ=255)
                                             ttACGATGTCTGGTGTTGGCCCACTTACTATGACCTGCCATAAAAATATCTCCAGATAGcc                            >  1:2511285/1‑61 (MQ=255)
                                             ttACGATGTCTGGTGTTGGCCCACTTACTATGACCTGCCATAAAAATATCTCCAGATAGcc                            >  1:2562626/1‑61 (MQ=255)
                                             ttACGATGTCTGGTGTTGGCCCACTTACTATGACCTGCCATAAAAATATATCCAGATAGcc                            >  1:287290/1‑61 (MQ=255)
                                             ttACGATGTCTGGTGTTGGCCCACTTACTATGAACTGCCATAAAAATATCTCCAGATAGCCCTGCCTGtt                   >  1:817994/1‑70 (MQ=255)
                                              tACGATGTCTGGTGTTGGCCCACTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGtt                   >  1:347271/1‑69 (MQ=255)
                                              tACGATGTCTGGTGTTGGCCCACTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTc                  <  1:1031373/70‑1 (MQ=255)
                                                cGATGTCTGGTGTTGGCCCACTTACTATGACCTGCCATAAAAATATCTCCAGATAGCcct                          <  1:777730/60‑1 (MQ=255)
                                                cGATGTCTGGTGTTGGCCCACTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGtt                   <  1:1474489/67‑1 (MQ=255)
                                                cGATGTCTGGTGTTGGCCCACTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAgg               <  1:1970851/71‑1 (MQ=255)
                                                cGATGTCTGGTGTTGGCCCACTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAgg               <  1:2966693/71‑1 (MQ=255)
                                                cGATGTCTGGTGTTGGCCCACTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAgg               >  1:1382138/1‑71 (MQ=255)
                                                  aTGTCTGGTGTTGGCCCACTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCa             <  1:1629225/71‑1 (MQ=255)
                                                   tGTCTGGTGTTGGCCCACTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAg            >  1:2394891/1‑71 (MQ=255)
                                                   tGTCTGGTGTTGGCCCACTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAg            <  1:2153297/71‑1 (MQ=255)
                                                    gTCTGGTGTTGGCCCACTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGc           >  1:2924571/1‑71 (MQ=255)
                                                       tGGTGTTGGCCCACTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGt         >  1:1374858/1‑70 (MQ=255)
                                                        ggTGTTGGCCCACTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGtt        >  1:2866766/1‑70 (MQ=255)
                                                            ttGGCCCACTTACTATGACCTGCCATAAAAATATCTCCAGAt                                <  1:2291244/42‑1 (MQ=255)
                                                            ttGGCCCACTTACTATGACCTGCCATAAAAATATCTCCAGAt                                <  1:1376904/42‑1 (MQ=255)
                                                            ttGGCCCACTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGtt        <  1:1392698/66‑1 (MQ=255)
                                                              ggCCCACTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTAc  >  1:3036276/1‑70 (MQ=255)
                                                                cccACTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGtt                   <  1:2313116/51‑1 (MQ=255)
                                                                    |                                                               
ATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTAC  >  minE/1251982‑1252113

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: