Predicted mutation
evidence seq id position mutation annotation gene description
RA minE 1,695,639 T→C E124G (GAG→GGG)  yfhQ ← predicted methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,695,6390TC100.0% 131.0 / NA 36E124G (GAG→GGG) yfhQpredicted methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (16/20);  total (16/20)

TCGCCAGGTTCAGCGAGCTGTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAGCGCCACCGGGGTATTTGCCGCTTCA  >  minE/1695571‑1695704
                                                                    |                                                                 
tCGCCAGGTTCAGCGAGCTGTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAAccct                                                                 >  1:57108/1‑71 (MQ=255)
     aGGTTCAGCGAGCTGTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAACCCTTCAtt                                                            >  1:1772447/1‑71 (MQ=255)
           aGCGAGCTGTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAACCCTTCATTGGTCaa                                                      <  1:2515871/71‑1 (MQ=255)
           aGCGAGCTGTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAACCCTTCATTGGTCaa                                                      <  1:189349/71‑1 (MQ=255)
           aGCGAGCTGTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAACCCTTCATTGGTCaa                                                      <  1:1753506/71‑1 (MQ=255)
             cGAGCTGTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAACCCTTCATTGGTCAAGc                                                    >  1:669515/1‑71 (MQ=255)
              gAGCTGTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAACCCTTCATTGGTCAAGcc                                                   <  1:72440/71‑1 (MQ=255)
                gCTGTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAACCCTTCATTGGTCAAGCCGa                                                 <  1:116321/71‑1 (MQ=255)
                 cTGTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAACCCTTCATTGGTCAAGCCGAc                                                <  1:2375290/71‑1 (MQ=255)
                  tGTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAACCCTTCATTGGTCAAGCCGAcg                                               >  1:1081751/1‑71 (MQ=255)
                            gTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAACCCTTCATTGGTCAAGCCGACGCGCTCGCGAc                                     >  1:2195811/1‑71 (MQ=255)
                               aGCCGCAATCGCGACATGATATTGGCATTTCTGCAACCCTTCATTGGTCAAGCCGACGCGCTCGCGACCa                                   <  1:2369862/70‑1 (MQ=255)
                               aGCCGCAATCGCGACATGATAATGGCATTTCTGCAACCCTTCATTGGTCAAGCCGACGCGCTCGCGACCa                                   <  1:2496433/70‑1 (MQ=255)
                                   gCAATCGCGACATGATAATGGCATTTCTGCAACCCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAacac                              <  1:1225690/71‑1 (MQ=255)
                                   gCAATCGCGACATGATAATGGCATTTCTGCAACCCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAacac                              <  1:1095582/71‑1 (MQ=255)
                                       tCGCGACATGATAATGGCATTTCTGCAACCCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAgc                          >  1:680865/1‑71 (MQ=255)
                                       tCGCGACATGATAATGGCATTTCTGCAACCCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAgc                          >  1:1341250/1‑71 (MQ=255)
                                       tCGCGACATGATAATGGCATTTCTGCAACCCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAgc                          >  1:2027396/1‑71 (MQ=255)
                                       tCGCGACATGATAATGGCATTTCTGCAACCCTTCATTGGTCAAGACGACGCGCTCGCGACCAAACACCAgc                          >  1:1701196/1‑71 (MQ=255)
                                          cgACATGATAATGGCATTTCTGCAACCCTTCATTGGTCAAGCCGACGCGCTc                                          >  1:782666/1‑52 (MQ=255)
                                          cgACATGATAATGGCATTTCTGCAACCCTTCATTGGTCAAGCCGACGCGCTc                                          >  1:912293/1‑52 (MQ=255)
                                            aCATGATAATGGCATTTCTGCAACCCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAGCGCCAc                     >  1:3073271/1‑71 (MQ=255)
                                            aCATGATAATGGCATTTCTGCAACCCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAGCGCCAc                     >  1:1902796/1‑71 (MQ=255)
                                            aCATGATAATGGCATTTCTGCAACCCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAGCGCCAc                     >  1:179659/1‑71 (MQ=255)
                                              aTGATAATGGCATTTCTGCAACCCTTCATTGGTCa                                                       <  1:1332336/35‑1 (MQ=255)
                                              aTGATAATGGCATTTCTGCAACCCTTCATTGGTCa                                                       <  1:1147683/35‑1 (MQ=255)
                                                gATAATGGCATTTCTGCAACCCTTCATTGGTCAAg                                                     >  1:2041835/1‑35 (MQ=255)
                                                gATAATGGCATTTCTGCAACCCTTCATTGGTCAAg                                                     >  1:3164191/1‑35 (MQ=255)
                                                gATAATGGCATTTCTGCAACCCTTCATTGGTCAAg                                                     >  1:863556/1‑35 (MQ=255)
                                                    aTGGCATTTCTGCAACCCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACAcc                             <  1:2116269/55‑1 (MQ=255)
                                                    aTGGCATTTCTGCAACCCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAGCGCCACCGGGGTAt             <  1:3118679/71‑1 (MQ=255)
                                                    aTGGCATTTCTGCAACCCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAGCGCCACCGGGGTAt             <  1:1346152/71‑1 (MQ=255)
                                                          tttCTGCAACCCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAGCGCCACCGGGGTATTTGc         <  1:2155700/69‑1 (MQ=255)
                                                                cAACCCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCa                            <  1:2739306/44‑1 (MQ=255)
                                                                cAACCCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCa                            <  1:8065/44‑1 (MQ=255)
                                                                cAACCCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAGCGCCACCGGGGTATTTGCCGCtt    <  1:2800117/68‑1 (MQ=255)
                                                                cAACCCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAGCGCCACCGGGGTATTTGCCGCTTCa  <  1:1276341/70‑1 (MQ=255)
                                                                    |                                                                 
TCGCCAGGTTCAGCGAGCTGTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAGCGCCACCGGGGTATTTGCCGCTTCA  >  minE/1695571‑1695704

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: