Predicted mutation
evidence seq id position mutation annotation gene description
RA minE 1,530,607 G→A T69I (ACT→ATT)  yfcA ← conserved inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,530,6070GA100.0% 115.7 / NA 36T69I (ACT→ATT) yfcAconserved inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (15/21);  total (15/21)

CATGGCGATATTGAGTTTCTGATCGCTTAAACTAACCACTTTGCGGCGAATAAAGTAGATAGTAGCGGAAATAGAGCCGCCGCAGGCTTGCAGTTTATTGGTTGCCAGCGCATT  >  minE/1530546‑1530659
                                                             |                                                    
cATGGCGATATTGAGTTTCTGATCGCTTAAACTAACCACTTTGCGGCGAATAAAGTAGATAATAGCGGaaa                                             >  1:2210337/1‑71 (MQ=255)
cATGGCGATATTGAGTTTCTGATCGCTTAAACTAACCACTTTGCGGCGAATAAAGTAGATAATAGCGGaaa                                             >  1:3148272/1‑71 (MQ=255)
cATGGCGATATTGAGTTTCTGATCGATTAAACTAACCACTTTGCGGCGAATAAAGTAGATAATAGCGGaaa                                             >  1:739745/1‑71 (MQ=255)
 aTGGCGATATTGAGTTTCTGATCGCTTAAACTAACCACTTTGCGGCTAATAAAGTAGATAATAGCGGaaa                                             <  1:1255595/70‑1 (MQ=255)
 aTGGCGATATTGAGTTTCTGATCGCTTAAACTAACCACTTTGCGGCGAATAAAGTAGATAATAGCGGaaa                                             <  1:1247641/70‑1 (MQ=255)
 aTGGCGATATTGAGTTTCTGATCGCTTAAACTAACCACTTTGCGGCGAATAAAGTAGATAATAGCGGaaa                                             <  1:821287/70‑1 (MQ=255)
 aTGGCGATATTGAGTTTCTGATCGCTTAAACTAACCACTTTGCGGCGAATAAAGTAGATAATAGCGGaaa                                             <  1:607608/70‑1 (MQ=255)
 aTGGCGATATTGAGTTTCTGATCGCTTAAACTAACCACTTTGCGGCGAATAAAGTAGATAATAGCGGaaa                                             <  1:521307/70‑1 (MQ=255)
 aTGGCGATATTGAGTTTCTGATCGCTTAAACTAACCACTTTGCGGCGAATAAAGTAGATAATAGCGGaaa                                             <  1:192800/70‑1 (MQ=255)
 aTGGCGATATTGAGTTTCTGATCGCTTAAACTAACCACTTTGCGGCGAATAAAGTAGATAATAGCGGaaa                                             <  1:465254/70‑1 (MQ=255)
 aTGGCGATATTGAGTTTCTGATCGCTTAAACTAACCACTTTGCGGCGAATAAAGTAGATAATAGCGGaaa                                             <  1:274549/70‑1 (MQ=255)
    gCGATATTGAGTTTCTTATCGCTTAAACTAACCACTTTGCGGCGAATAAAGTAGATAATAGCGGAAATaga                                         <  1:2004793/71‑1 (MQ=255)
     cGATATTGAGTTTCTGATCGCTTAAACTAACCACTTTGCGGCGAATAAAGTAGATAATAGCGGAAATagag                                        >  1:909012/1‑71 (MQ=255)
     cGATATTGAGTTTCTGATCGCTTAAACTAACCACTTTGCGGCGAATAAAGTAGATAATAGCGGAAATagag                                        >  1:4076495/1‑71 (MQ=255)
       atatTGAGTTTCTGATCGCTTAAACTAACCACTTTGCGGCGAATAAAGTAGATAATAGCGGAAATAGAgcc                                      >  1:1533550/1‑71 (MQ=255)
       atatTGAGTTTCTGATCGCTTAAACTAACCACTTTGCGGCGAATAAAGTAGATAATAGCGGAAATAGAgc                                       >  1:4001528/1‑70 (MQ=255)
             aGTTTCTGATCGCTTAAACTAACCACTTTGCGGCGAATAAAGTAGATAATAGCGGAAATAGAGCcgccgc                                 <  1:716815/70‑1 (MQ=255)
                 tCTGATCGCTTAAACTAACCACTTTGCGGCGAATAAAGTAGATAATAGCGGAAATAGAGCCGCCGCAGGc                             >  1:1446588/1‑70 (MQ=255)
                 tCTGATCGCTTAAACTAACCACTTTGCGGCGAATAAAGTAGATAATAGCGGAAATAGAGCCGCCGCAGGCt                            >  1:338239/1‑71 (MQ=255)
                 tCTGATCGCTTAAACTAACCACTTTGCGGCGAATAAAGTAGATAATAGCGGAAATAGAGCCGCCGCAGGCt                            >  1:3543942/1‑71 (MQ=255)
                 tCTGATCGCTTAAACTAACCACTTTGCGGCGAATAAAGTAGATAATAGCGGAAATAGAGCCGCCGCAGGCt                            >  1:3947933/1‑71 (MQ=255)
                 tCTGATCGCTTAAACTAACCACTTTGCGGCGAATAAAGTAGATAATAGCGGAAATAGAGCCGCCGCAGGCt                            >  1:4217846/1‑71 (MQ=255)
                 tCTGATCGCTTAAACTAACCACTTTGCGGCGAATAAAGTAGATAATAGCGGAAATAGAGCCGCCGCAGGCt                            >  1:4652523/1‑71 (MQ=255)
                 tCTGATCGCTTAAACTAACCACTTTGCGGCGAATAAAGTAGATAATAGCGGAAATAGAGCCGCCGCAGGCt                            >  1:4815413/1‑71 (MQ=255)
                 tCTGATCGCTTAAACTAACCACTTTGCGGCGAATAAAGTAGATAATAGCGGAAATAGAGCCGCCGCAGGCt                            >  1:1875137/1‑71 (MQ=255)
                     aTCGCTTAAACTAACCACTTTGCGGCGAATAAAGTAGATAATAGCGGAAATAGAGCCGCCGCAGGCTTg                          <  1:3326685/69‑1 (MQ=255)
                     aTCGCTTAAACTAACCACTTTGCGGCGAATAAAGTAGATAATAGCGGAAATAGAGCCGCCGCAGGCTTg                          <  1:114891/69‑1 (MQ=255)
                     aTCGCTTAAACTAACCACTTTGCGGCGAATAAAGTAGATAATAGCGGAAATAGAGCCGCCGCAGGCTTg                          <  1:4564089/69‑1 (MQ=255)
                     aTCGCTTAAACTAACCACTTTGCGGCGAATAAAGTAGATAATAGCGGAAATAGAGCCGCCGCAGGCTTg                          <  1:2914858/69‑1 (MQ=255)
                     aTCGCTTAAACTAACCACTTTGCGGCGAATAAAGTAGATAATAGCGGAAATAGAGCCGCCGCAGGCTTg                          <  1:2949241/69‑1 (MQ=255)
                      tCGCTTAAACTAACCACTTTGCGGCGAATAAAGTAGATAATAGCGGAAATAGAGCCGCCGCAGGCTTg                          <  1:1503639/68‑1 (MQ=255)
                             aaCTAACCACTTTGCGGCGAATAAAGTAGATAATAGCGGAAATAGGGCCGCCGCAGGCTTGCAGTTTATTg                <  1:3063858/71‑1 (MQ=255)
                             aaCTAACCACTTTGCGGCGAATAAAGTAGATAATAGCGGAAATAGAGCCGCCGCAGGCTTGCAGTTTATTg                <  1:2307538/71‑1 (MQ=255)
                             aaCTAACCACTTTGCGGCGAATAAAGTAGATAATAGCGGAAATAGAGCCGCCGCAGGCTTGCAGTTTATTg                <  1:4150617/71‑1 (MQ=255)
                                          gcggcgAATAAAGTAGATAATAGCGGAAATAGAGCCGCCGCAGGCTTGCAGTTTATTGGTTGCCAgcg      <  1:4545301/68‑1 (MQ=255)
                                           cggcgAATAAAGTAGATAATAGCGGAAATAGAGCCGCCGCAGGCTTGCAGTTTATTGGTTGCCAGCGCAtt  <  1:2337662/71‑1 (MQ=255)
                                                             |                                                    
CATGGCGATATTGAGTTTCTGATCGCTTAAACTAACCACTTTGCGGCGAATAAAGTAGATAGTAGCGGAAATAGAGCCGCCGCAGGCTTGCAGTTTATTGGTTGCCAGCGCATT  >  minE/1530546‑1530659

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: