Predicted mutation
evidence seq id position mutation annotation gene description
RA minE 798,949 C→T M214I (ATG→ATA treA ← periplasmic trehalase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE798,9490CT100.0% 163.9 / NA 57M214I (ATG→ATAtreAperiplasmic trehalase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (28/29);  total (28/29)

TCTTTTTGCATTTGCGGCAGGTATTGCTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACCATCAGGGCAAAGAAGGGCGGTTGCGAGCGGCTTAAATAGTAACTGCGGTTGCCGTTGGGAATATGA  >  minE/798881‑799015
                                                                    |                                                                  
tCTTTTTGCATTTGCGGCAGGTATTGCTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACTAt                                                                  >  1:2920064/1‑71 (MQ=255)
 cTTTTTGCATTTGCGGCAGGTATTGCTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACTAt                                                                  <  1:744440/70‑1 (MQ=255)
 cTTTTTGCATTTGCGGCAGGTATTGCTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACTAt                                                                  >  1:535859/1‑70 (MQ=255)
  tttttGCATTTGCGGCAGGTATTGCTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACTAt                                                                  >  1:1346239/1‑69 (MQ=255)
  tttttGCATTTGCGGCAGGTATTGCTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACTAt                                                                  <  1:1773269/69‑1 (MQ=255)
    tttGCATTTGCGGCAGGTATTGCTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACTATc                                                                 <  1:1143520/68‑1 (MQ=255)
    tttGCATTTGCGGCAGGTATTGCTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACTATCAgg                                                              >  1:2563085/1‑71 (MQ=255)
         aTTTGCGGCAGGTATTGCTTGAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACTATCAGGGCaaa                                                         <  1:1779296/71‑1 (MQ=255)
         aTTTGCGGCAGGTATTGCTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACTATCAGGGCa                                                           >  1:689471/1‑69 (MQ=255)
          tttGCGGCAGGTATTGCTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACTATCAGGGCAaag                                                        <  1:2872456/71‑1 (MQ=255)
            tGCGGCAGGTATTGCTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACTATCAGGGCAaag                                                        >  1:714812/1‑69 (MQ=255)
            tGCGGCAGGTATTGCTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACTATCAGGGCAaag                                                        >  1:2908026/1‑69 (MQ=255)
              cGGCAGGTATTGCTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACTATCAGGGCAAAGAAgg                                                    <  1:695571/71‑1 (MQ=255)
              cGGCAGGTATTGCTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACTATCAGGGCAAAGAAgg                                                    <  1:2862191/71‑1 (MQ=255)
                     tATTGCTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACTATCAGGGCAaagaag                                                     <  1:28879/63‑1 (MQ=255)
                      aTTGCTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACTATCAGGGCAaag                                                        >  1:913127/1‑59 (MQ=255)
                      aTTGCTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACTATCAGGGCAAAGAAGGGCGGTTGc                                            >  1:1789109/1‑71 (MQ=255)
                       ttGCTTCAATGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACTATCAGGGCAAAGAAGGGCGGTTGCg                                           <  1:2445896/71‑1 (MQ=255)
                         gCTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACTATCAGGGCAAAGAAGGGCGGTTGCGAg                                         >  1:1682802/1‑71 (MQ=255)
                           ttCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACTATCAGGGCAaagaag                                                     <  1:691499/57‑1 (MQ=255)
                           ttCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACTATCAGGGCAAAGAAGGGCGGTTGCg                                           <  1:1732338/67‑1 (MQ=255)
                           ttCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACTATCAGGGCAAAGAAGGGCGGTTGCGAGCg                                       >  1:432470/1‑71 (MQ=255)
                           ttCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACTATCAGGGCAAAGAAGGGCGGTTGCGAGCg                                       >  1:161083/1‑71 (MQ=255)
                            tCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACTATCAGGGCAAAGAAGGGCGGtt                                              <  1:1596911/63‑1 (MQ=255)
                             cAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACTATCAGGGCAAAGAAGGGCGGTTGCGAGCGGc                                     >  1:139052/1‑71 (MQ=255)
                              aaCGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACTATCAGGGCAAAGAAGGGCGGTTGCGAGCGGCt                                    >  1:1897303/1‑71 (MQ=255)
                               aCGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACTATCAGGGCAAAGAAGGGCg                                                 <  1:579679/57‑1 (MQ=255)
                                 gcgGCATCGCCTTCATGCTGCGCCAGTAACTCTACTATCAGGGCAAAGAAGGGCg                                                 <  1:233216/55‑1 (MQ=255)
                                 gcgGCATCGCCTTCATGCTGCGCCAGTAACTCTACTATCAGGGCAAAGAAGGGCg                                                 <  1:122578/55‑1 (MQ=255)
                                 gcgGCATCGCCTTCATGCTGCGCCAGTAACTCTACTATCAGGGCAAAGAAGGGCGCTTGCGAGCGGCTTa                                  >  1:662891/1‑70 (MQ=255)
                                  cgGCATCGCCTTCATGCTGCGCCAGTAACTCTACTATCAGGGCAAAGAAGGGCGGTTGCGAGCGTCTTaaa                                <  1:1020930/71‑1 (MQ=255)
                                  cgGCATCGCCTTCATGCTGCGCCAGTAACTCTACTATCAGGGCAAAGAAGGGCGGTTGCGAGCGGCTTaaa                                <  1:2835425/71‑1 (MQ=255)
                                   ggCATCGCCTTCATGCTGCGCCAGTAACTCTACTATCAGGGCAAAGAAGAGCGGTTGCGAGCGGCTTaaa                                >  1:1371751/1‑70 (MQ=255)
                                      aTCGCCTTCATGCTGCGCCAGTAACTCTACTATCAGGGCAAAGAAGGGCGGTTGCGAGCGGCTTTAATAGt                            <  1:4293/71‑1 (MQ=255)
                                      aTCCCCTTCATGCTGCGCCAGTAACTCTACTATCAGGGCAaagaa                                                      <  1:736753/45‑1 (MQ=255)
                                       tCGCCTTCATGCTGCGCCAGTAACTCTACTATCAGGGCAaag                                                        <  1:684828/42‑1 (MQ=255)
                                        cGCCTTCATGCTGCGCCAGTAACTCTACTATCAGGGCAaag                                                        <  1:208803/41‑1 (MQ=255)
                                        cGCCTTCATGCTGCGCCAGTAACTCTACTATCAGGGCAaag                                                        <  1:1514597/41‑1 (MQ=255)
                                        cGCCTTCATGCTGCGCCAGTAACTCTACTATCAGGGCAAAGAAGGGCGGTTGCGAGCGGCTTAAATAGTaa                          >  1:2216305/1‑71 (MQ=255)
                                         gCCTTCATGCTGCGCCAGTAACTCTACTGTCAGGGCAAAGAAGGGCGGTTGCGAGCGGCTTAAATAGTaa                          <  1:2075664/70‑1 (MQ=255)
                                         gCCTTCATGCTGCGCCAGTAACTCTACTATCAGGGCAAAGAAGGGCGGTTGCGAGCGGCTTAAATAGTAAc                         >  1:1998568/1‑71 (MQ=255)
                                           cTTCATGCTGCGCCAGTAACTCTACTATCAGGGCAAAGAAGGGCGGtttc                                            >  1:2488433/1‑48 (MQ=255)
                                              cATGCTGCGCCAGTAACTCTACTATCAGGGCAAAGAAGGGCGGt                                               >  1:319140/1‑44 (MQ=255)
                                              cATGCTGCGCCAGTAACTCTACTATCAGGGCAAAGAAGGGCGGt                                               >  1:2305478/1‑44 (MQ=255)
                                                  ctgcGCCAGTAACTCTACTATCAGGGCAAAGAAGGGCGGTTGCGAGCgg                                      >  1:1721021/1‑49 (MQ=255)
                                                  ctgcGCCAGTAACTCTACTATCAGGGCAAAGAAGGGCGGTTGCGAGCGGCTTAAATAGTaa                          >  1:1562831/1‑61 (MQ=255)
                                                     cgcCAGTAACTCTACTATCAGGGCAAAGAAGGGCGGTTGCGAGCGGCTTAAATAGTAACTGCGGTTGCCGt             <  1:2538757/71‑1 (MQ=255)
                                                       ccAGTAACTCTACTATCAGGGCAAAGAAGGGCGGTTGCGAGCGGCTTAAATa                              >  1:2602138/1‑52 (MQ=255)
                                                         aGTAACTCTACTATCAGGGCAAAGAAGGGCGGTTGCGAGCGGCTTAAATAGTAACTGCGGTTGCCGTTggg         <  1:570269/71‑1 (MQ=255)
                                                          gTAACTCTACTATCAGGGCAAAGAAGGGCGGTTGCGAGCGGCTTAAATAGTAACTGCGGTTGCCGTTggg         >  1:1363671/1‑70 (MQ=255)
                                                          gTAACTCTACTATCAGGGCAAAGAAGGGCGGTTGCGAGCGGCTTAAATAGTAACTGCGGTTGCCGTTGGGa        <  1:305344/71‑1 (MQ=255)
                                                          gTAACTCTACTATCAGGGCAAAGAAGGGCGGTTGCGAGCGGCTTAAATAGTAACTGCGGTTGCCGTTGGGa        <  1:613395/71‑1 (MQ=255)
                                                          gTAACTCTACTATCAGGGCAAAGAAGGGCGGTTGCGAGCGGCTTAAATAGTAACTGCGGTTGCCGTTGGGa        <  1:1904082/71‑1 (MQ=255)
                                                           tAACTCTACTATCAGGGCAAAGAAGGGCGGTTGCGAGCGGCTTAAATAGTAACTGCGGTTGc                <  1:2595292/62‑1 (MQ=255)
                                                           tAACTCTACTATCAGGGCAAAGAAGGGCGGTTGCGAGCGGCTTAAATAGTAACTGCGGTTGCCGTTGGGaa       >  1:1022819/1‑71 (MQ=255)
                                                           tAACTCTACTATCAGGGCAAAGAAGGGCGGTTGCGAGCGGCTTAAATAGTAACTGCGGTTGCCGTTGGGaa       >  1:895210/1‑71 (MQ=255)
                                                               tctACTATCAGGGCAAAGAAGGGCGGTTGCGAGCGGCTTAAATAGTAACTGCGGTTGCCGTTgg          >  1:1412777/1‑64 (MQ=255)
                                                                ctactaTCAGGGCAAAGAAGGGCGGTTGCGAGCGGCTTAAAAAGTAACTGCGGTTGCCGTTGGGAATATGa  >  1:918833/1‑71 (MQ=255)
                                                                    |                                                                  
TCTTTTTGCATTTGCGGCAGGTATTGCTTCAACGCGGCATCGCCTTCATGCTGCGCCAGTAACTCTACCATCAGGGCAAAGAAGGGCGGTTGCGAGCGGCTTAAATAGTAACTGCGGTTGCCGTTGGGAATATGA  >  minE/798881‑799015

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: