Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 711,724 C→T A59T (GCG→ACG)  ybfE ← LexA‑regulated protein, CopB family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913711,7240CT100.0% 112.6 / NA 36A59T (GCG→ACG) ybfELexA‑regulated protein, CopB family
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (19/17);  total (19/17)

GATGCCCTGACTACGCAGCGCCGCCAGTTGCTGCATTAACATCTCTTCGATCAGTTCGCTACGGCTCATATTGCGCGACTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGCGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTTTATTAATACGCAGCTGTTCATCGCGCGAAAGCGGATTAGTTTTCGGTCGTCCCGGTCGACGCTCGTGCG  >  NC_000913/711611‑711855
                                                                                                                 |                                                                                                                                   
gATGCCCTGACTACGCAGCGCCGCCAGTTGCTGCATTAACATCTCTTCGATCAGTTCGCTACGGCTCATATTGCGCGACTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGACACGCTTa                                                                                                            <  2:30071/139‑1 (MQ=255)
  tGCCCTGACTACGCAGCGCCGCCAGTTGCTGCATTAACATCTCTTCGATCAGTTCGCTACGGCTCATATTGCGCGACTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGACACGCTTAAg                                                                                                          >  1:1551511/1‑139 (MQ=255)
     ccTGACTACGCAGCGCCGCCAGTTGCTGCATTAACATCTCTTCGATCAGTTCGCTACGGCTCATATTGCGCGACTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGACACGCTTAAGGcc                                                                                                       <  2:1551511/139‑1 (MQ=255)
            aCGCAGCGCCGCCAGTTGCTGCATTAACATCTCTTCGATCAGTTCGCTACGGCTCATATTGCGCGACTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACt                                                                                                >  2:731943/1‑139 (MQ=255)
                aGCGCCGCCAGTTGCTGCATTAACATCTCTTCGATCAGTACGCTACGGCTCATATTGCGCGACTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGt                                                                                            >  1:271697/1‑139 (MQ=255)
                   gccgctAGTTGCTGCATTAACATCTCTTCGATCAGTTCGCTACGGCTCATATTGCGCGACTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGaca                                                                                                                  <  2:858435/114‑1 (MQ=255)
                   gccgctAGTTGCTGCATTAACATCTCTTCGATCAGTTCGCTACGGCTCATATTGCGCGACTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGaca                                                                                                                  >  1:858435/1‑114 (MQ=255)
                   gccgccAGTTGCTGCATTAACATCTCTTCGATCAGTTCGCTACGGCTCATATTGCGCGACTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTcgc                                                                                         >  1:939555/1‑139 (MQ=255)
                                   ttAACATCTCTTCGATCAGTTCGCTACGGCTCATATTGCGCGACTCCGCCAGCTCGTTCAGTGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCg                                                                         <  1:1302074/139‑1 (MQ=255)
                                   ttAACATCTCTTCGATCAGTTCGCTACGGCTCATATTGCGCGACTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCg                                                                         <  2:939555/139‑1 (MQ=255)
                                    tAACATCTCTTCGATCAGTACGCTACGGCTCATATTGCGCGACTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGt                                                                        <  2:271697/139‑1 (MQ=255)
                                           tctTCGATCAGTTCGCTACGGCTCATATTGCGCGACTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTTTATTaa                                                                 >  2:1628116/1‑139 (MQ=255)
                                           tctTCGATCAGTTCGCTACGGCTCATATTGCGCGACTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTTTATTaa                                                                 >  2:1414131/1‑139 (MQ=255)
                                             ttCGATCAGTTCGCTACGGCTCATATTGCGCGACTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGt                                                                                                   >  2:1943254/1‑103 (MQ=255)
                                             ttCGATCAGTTCGCTACGGCTCATATTGCGCGACTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGt                                                                                                   <  1:1943254/103‑1 (MQ=255)
                                             ttCGATCAGTTCGCTACGGCTCATATTGCGCGACTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTTTATTAATa                                                               >  1:735093/1‑139 (MQ=255)
                                               cGATCAGTTCGCTACGGCTCATATTGCGCGACTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTTTATTAATACg                                                             <  1:731943/139‑1 (MQ=255)
                                                                cTCATATTGCGCGACTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTTTATTAATACGCAGCTGTTCATcgcgcg                                            >  1:1097075/1‑139 (MQ=255)
                                                                         cgcgACTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGa                                                                                                                    >  2:1105441/1‑58 (MQ=255)
                                                                         cgcgACTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGa                                                                                                                    <  1:1105441/58‑1 (MQ=255)
                                                                         cgcgACTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGACACGCTTa                                                                                                            <  1:1742210/66‑1 (MQ=255)
                                                                         cgcgACTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGACACGCTTa                                                                                                            >  2:1742210/1‑66 (MQ=255)
                                                                         cgcgACTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGc                                                                          >  1:2047348/1‑100 (MQ=255)
                                                                         cgcgACTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGc                                                                          <  2:2047348/100‑1 (MQ=255)
                                                                         cgcgACTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTTTATTAATACGCAGCTGTTCATCGCGCGAAAGCGGAt                                   <  2:224770/139‑1 (MQ=255)
                                                                          gcgACTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTTTATTAATa                                                               <  2:556825/110‑1 (MQ=255)
                                                                          gcgACTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTTTATTAATa                                                               >  1:556825/1‑110 (MQ=255)
                                                                              cTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTgg                                                                              <  2:998218/91‑1 (MQ=255)
                                                                              cTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTgg                                                                              >  1:998218/1‑91 (MQ=255)
                                                                              cTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTTTATTAATACGCAGCTGTTCATCGCGCGAAAGCGGATTAGtt                              >  1:970878/1‑139 (MQ=255)
                                                                                   ccAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTTTATTAATACGCAGCTGTTCATCGCGCGAAAGCGGATTAGTTTTCgg                         <  1:470069/139‑1 (MQ=255)
                                                                                    cAGCTCGTTCAGCGCCTCGACAGCTTCCGTGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTTTATTAATACGCAGCTGTTCATCGCGCGAAAGCGGATTAGTTTTCGgt                        <  1:208482/139‑1 (MQ=255)
                                                                                                     cGACAGCTTCCGTGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTTTATTAATACGCAGCTGTTCATCGCGCGAAAGCGGATTAGTTTTCGGTCGTCCCGGTCGACGCTc       >  2:37933/1‑139 (MQ=255)
                                                                                                          gCTTCCGTGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTTTATTAATACGCAGCTGTTCATCGCGCGAAAGCGGATTAGTTTTCGGTCGTCCCGGTCGACGCTCGTgcg  >  1:1076116/1‑139 (MQ=255)
                                                                                                             tCCGTGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGtt                                                                       >  2:447608/1‑67 (MQ=255)
                                                                                                             tCCGTGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGtt                                                                       <  1:447608/67‑1 (MQ=255)
                                                                                                                 |                                                                                                                                   
GATGCCCTGACTACGCAGCGCCGCCAGTTGCTGCATTAACATCTCTTCGATCAGTTCGCTACGGCTCATATTGCGCGACTCCGCCAGCTCGTTCAGCGCCTCGACAGCTTCCGCGTTCAGCTTCAGTTCGACACGCTTAAGGCCACGTACTTTGTCGCGTTTTAGCTGGTTGCGTTTATTAATACGCAGCTGTTCATCGCGCGAAAGCGGATTAGTTTTCGGTCGTCCCGGTCGACGCTCGTGCG  >  NC_000913/711611‑711855

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