Predicted mutation
evidence seq id position mutation annotation gene description
RA minE 2,279,749 Δ1 bp coding (633/1101 nt) trmA → tRNA (uracil‑5‑)‑methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,279,7460A.100.0% 168.4 / NA 41coding (630/1101 nt)trmAtRNA (uracil‑5‑)‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (22/19);  total (22/19)

CTTTACCCAGCCGAACGCGGCGATGAATATTCAGATGCTGGAATGGGCGCTGGACGTAACCAAAGGCTCAAAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACTTTTCATTAGCGCTGGCGCGTAATT  >  minE/2279677‑2279810
                                                                     |                                                                
cTTTACCCAGCCGAACGCGGCGATGAATATTCAGATGCTGGAATGGGCGCTGGACGTAACCAAAGGCTCaa                                                                 <  1:138024/71‑1 (MQ=255)
 tttACCCAGCCGAACGCGGCGATGAATATTCAGATGCTGGAATGGGCGCTGGACGTAACCAAAGGCTCa                                                                  >  1:902724/1‑69 (MQ=255)
  ttACCCAGCCGAACGCGGCGATGAATATTCAGATGCTGGAATGGGCGCTGGACGTAACCAAAGGCTCaaag                                                               >  1:2065599/1‑70 (MQ=255)
     cccAGCCGAACGCGGCGATGAATATTCAGATGCTGGAATGGGCGCTGGACGTAACCAAAGGCTCaa                                                                 <  1:1137642/66‑1 (MQ=255)
      ccAGCCGAACTCGGCGATGAATATTCAGATGCTGGAATGGGCGCTGGACGTAACCAAAGGCTCAAAGGCGa                                                          <  1:49140/71‑1 (MQ=255)
      ccAGCCGAACGCGGCGATGAATATTCAGATGCTGGAATGGGCGCTGGACGTAACCAAAGGCTCAAAGGCGa                                                          <  1:486080/71‑1 (MQ=255)
         gCCGAACGCGGCGATGAATATTCAGATGCTGGAATGGGCGCTGGAGGTAACCAAAGGCTCaaa                                                                >  1:2079147/1‑63 (MQ=255)
         gCCGAACGCGGCGATGAATATTCAGATGCTGGAATGGGCGCTGGACGTAACCAAAGGCTCaaa                                                                >  1:773499/1‑63 (MQ=255)
         gCCGAACGCGGCGATGAATATTCAGATGCTGGAATGGGCGCTGGACGTAACCAAAGGCTCaaa                                                                >  1:2870589/1‑63 (MQ=255)
         gCCGAACGCGGCGATGAATATTCAGATGCTGGAATGGGCGCTGGACGTAACCAAAGGCTCAAAGGCGAtt                                                        >  1:526759/1‑70 (MQ=255)
            gAACGCGGCGATGAATATTCAGATGCTGGAATGGGCGCTGGACGTAACCAAAGGCTCaa                                                                 >  1:627544/1‑59 (MQ=255)
            gAACGCGGCGATGAATATTCAGATGCTGGAATGGGCGCTGGACGTAACCAAAGGCTCaa                                                                 >  1:623925/1‑59 (MQ=255)
            gAACGCGGCGATGAATATTCAGATGCTGGAATGGGCGCTGGACGTAACCAAAGGCTCAAAGGCGATTTACt                                                    <  1:504121/71‑1 (MQ=255)
             aaCGCGGCGATGAATATTCAGATGCTGTAATGGGCGCTGGACGTAACCAAAGGCTCaaagg                                                              <  1:1305348/61‑3 (MQ=255)
             aaCGCGGCGATGAATATTCAGATGCTGGAATGGGCGCTGGACGTAACCAAAGGCTCaaagg                                                              >  1:326786/1‑59 (MQ=255)
             aaCGCGGCGATGAATATTCAGATGCTGGAATGGGCGCTGGACGTAACCAAAGGCTCAAAGGCGATTTACTg                                                   >  1:1908356/1‑71 (MQ=255)
                    cgATGAATATTCAGATGCTGGAATGGGCGCTGGACGTAACCAAAGGCTCAAAGGCGATTTACTGGAGCTg                                             >  1:2782639/1‑70 (MQ=255)
                      aTGAATATTCAGATGCTGGAATGGGCGCTGGACGTAACCAAAGGCTCAAAGGCGATTTACTGGAGCTGTa                                           >  1:2434595/1‑70 (MQ=255)
                       tGAATATTCAGATGCTGGAATGGGCGCTGGACGTAACCAAAGGCTCAAAGGCGATTTACTGGAGCTGTACt                                         >  1:2292245/1‑71 (MQ=255)
                           tatTCAGATGCTGGAATGGGCGCTGGACGTAACCAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCgg                                     >  1:1832283/1‑71 (MQ=255)
                           tatTCAGATGCTGGAATGGGCGCTGGACGTAACCAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCgg                                     <  1:605004/71‑1 (MQ=255)
                            atTCAGATGCTGGAATGGGCGCTGGACGTAACCAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCgg                                     >  1:2413812/1‑70 (MQ=255)
                             ttCAGATGCTGGAATGGGCGCTGGACGTAACCAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGc                                    <  1:486195/70‑1 (MQ=255)
                             ttCAGATGCTGGAATGGGCGCTGGACGTAACCAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGc                                    >  1:1614837/1‑70 (MQ=255)
                              tCAGATGCTGGAATGGGCGCTGGACGTAACCAAAGGCTCAAAGGCGATTTACTGGa                                                 >  1:1085051/1‑56 (MQ=255)
                              tCAGATGCTGGAATGGGCGCTGGACGTAACCAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGCaa                                  >  1:1426746/1‑71 (MQ=255)
                                aGATGCTGGAATGGGCGCTGGACGTAACCAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGCAACg                                >  1:2810234/1‑71 (MQ=255)
                                  aTGCTGGAATGGGCGCTGGACGTAACCAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGt                              <  1:1823034/71‑1 (MQ=255)
                                   tGCTGGAATGGGCGCTGGACGTAACCAAAGGCTCaaa                                                                >  1:1285997/1‑37 (MQ=255)
                                   tGCTGGAATGGGCGCTGGACGTAACCAAAGGCTCaaa                                                                >  1:1524275/1‑37 (MQ=255)
                                   tGCTGGAATGGGCGCTGGACGTAACCAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGt                              >  1:2846077/1‑70 (MQ=255)
                                    gCTGGAATGGGCGCTGGACGTAACCAAAGGCTCAAAGGCGATTTAc                                                     <  1:663136/46‑1 (MQ=39)
                                        gAATGGGCGCTGGACGTAACCAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACttt                        <  1:2091749/71‑1 (MQ=255)
                                        gAATGGGCGCTGGACGTAACCAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACt                          <  1:2824682/69‑1 (MQ=255)
                                          aTGGGCGCTGGACGTAACCAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACtttt                       >  1:2288780/1‑70 (MQ=255)
                                           tGGGCGCTGGACGTAACCAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGCaa                                  <  1:2604592/58‑1 (MQ=255)
                                            gggCGCTGGACGTAACCAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACTTTTCAt                    >  1:2696620/1‑71 (MQ=255)
                                              gcgcTGGACGTAACCAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACTTTTCATTa                  <  1:805630/71‑1 (MQ=255)
                                               cgcTGGACGTAACCAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACTTTTCATTAg                 <  1:672176/71‑1 (MQ=255)
                                                 cTGGACGTAACCAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACTTTTCATTAgcg               <  1:368073/71‑1 (MQ=255)
                                                  tGGACGTAACCAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACtt                         <  1:2137988/60‑1 (MQ=255)
                                                    gACGTAACCAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGc                                    <  1:2319796/47‑1 (MQ=38)
                                                       gTAACCAAAGGCTCAAAGGCGATTTACTGGAGCTg                                             <  1:1862145/35‑1 (MQ=38)
                                                        tAACCAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGCaa                                  <  1:128745/45‑1 (MQ=38)
                                                        tAACCAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTa                             <  1:509650/50‑1 (MQ=39)
                                                         aaCCAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACTTTTCATTAGCGCTGGcgcg       >  1:2839954/1‑71 (MQ=255)
                                                         aaCCAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACTTTTCATTAGCGCTGGcgcg       >  1:1551248/1‑71 (MQ=255)
                                                          aCCAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAAc                           >  1:2528020/1‑50 (MQ=39)
                                                           ccAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACTTTTCATTAGCGCTGGCGCGt      >  1:361146/1‑70 (MQ=255)
                                                           ccAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACTTTTCATTAGCGCTGGCGCGt      >  1:1313903/1‑70 (MQ=255)
                                                            cAAAGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACTTTTCATTAGCGCTGGCGCGTaa    <  1:2849682/71‑1 (MQ=255)
                                                              aaGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACTTTTCATTAGCGCTGGCGCGTAAtt  >  1:1558746/1‑71 (MQ=255)
                                                              aaGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACTTTTCATTAGCGCTGGCGCGTAAtt  >  1:899044/1‑71 (MQ=255)
                                                              aaGGCTCAAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACTTTTCATTAGCGCTGGCGCGTAAtt  >  1:902787/1‑71 (MQ=255)
                                                                     |                                                                
CTTTACCCAGCCGAACGCGGCGATGAATATTCAGATGCTGGAATGGGCGCTGGACGTAACCAAAGGCTCAAAAGGCGATTTACTGGAGCTGTACTGCGGCAACGGTAACTTTTCATTAGCGCTGGCGCGTAATT  >  minE/2279677‑2279810

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: