Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 747,445 G→T P109T (CCT→ACT)  abrB ← putative regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913747,4450GT100.0% 111.0 / NA 31P109T (CCT→ACT) abrBputative regulator
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (15/16);  total (15/16)

CGCCGTAATCTTGCGCCATAGCGACCATTGCCGCCGCGCCGCCAGGAGAGGAGCCCCAGGCACCGGTATTTCCAGGCAATGAGCTATAGCGCACCAATAACCAACCCACAATGGCGCTGGAAAGCAGCGTTACCAGTAAAATTGCAAGCACG  >  NC_000913/747370‑747521
                                                                           |                                                                            
cgcCGTAATCTTGCGCCATAGCGACCATTGCCGCCGCGCCGCCAGGAGAGGAGCCCCAGGCACCGGTATTTCCAGTCAAt                                                                          >  1:153690/1‑80 (MQ=255)
cgcCGTAATCTTGCGCCATAGCGACCATTGCCGCCGCGCCGCCAGGAGAGGAGCCCCAGGCACCGGTATTTCCAGTCAAt                                                                          <  1:153647/80‑1 (MQ=255)
cgcCGTAATCTTGCGCCATAGCGACCATTGCCGCCGCGCCGCCAGGAGAGGAGCCCCAGGCACCGGTATTTCCAGTCAATg                                                                         <  2:153648/81‑1 (MQ=255)
    gTAATCTTGCGCCATAGCGACCATTGCCGCCGCGCCGCCAGGAGAGGAGCCCCAGGCACCGGTATTTCCAGTCAATGAGCt                                                                     <  2:153649/81‑1 (MQ=255)
         cTTGCGCCATAGCGACCATTGCCGCCGCGCCGCCAGGAGAGGAGCCCCAGGCACCGGTATTTCCAGTCAATGAGCTATAgc                                                                >  2:153695/1‑81 (MQ=255)
            gcgcCATAGCGACCATTGCCGCCGCGCCGCCAGGAGAGGAGCCCCAGGCACCGGTATTTCCAGTCAATGAGCTATAgcgc                                                              <  2:153651/80‑1 (MQ=255)
            gcgcCATAGCGACCATTGCCGCCGCGCCGCCAGGAGAGGAGCCCCAGGCACCGGTATTTCCAGTCAATGAGCTATAgc                                                                <  2:153650/78‑1 (MQ=255)
            gcgcCATAGCGACCATTGCCGCCGCGCCGCCAGGAGAGGAGCCCCAGGCACCGGTATTTCCAGTCAATGAGCTATAGCGCa                                                             >  2:153688/1‑81 (MQ=255)
             cgcCATAGCGACCATTGCCGCCGCGCCGCCAGGAGAGGAGCCCCAGGCACCGGTATTTCCAGTCAATGAGCTATAGCGCa                                                             <  1:153652/80‑1 (MQ=255)
             cgcCATAGCGACCATTGCCGCCGCGCCGCCAGGAGAGGAGCCCCAGGCACCGGTATTTCCAGTCAATGAGCTATAGCGCAc                                                            <  1:153653/81‑1 (MQ=255)
                     cGACCATTGCCGCCGCGCCGCCAGGAGAGGAGCCCCAGGCACCGGTATTTCCAGTCAATGAGCTATAGCGCACCAATaacc                                                    <  1:153654/81‑1 (MQ=255)
                       aCCATTGCCGCCGCGCCGCCAGGAGAGGAGCCCCAGGCACCGGTATTTCCAGTCAATGAGCTATAGCGCACCAATaaccaa                                                  >  1:153687/1‑81 (MQ=255)
                                    cgccgcCAGGAGAGGAGCCCCAGGCACCGGTATTTCCAGTCAATGAGCTATAGCGCACCAATAACCAACCCACAATGgcgc                                     >  1:153694/1‑81 (MQ=255)
                                      ccgccAGGAGAGGAGCCCCAGGCACCGGTATTTCCAGTCAATGAGCTATAGCGCACCAATAACCAACCCACAATGGCGCTg                                   <  1:153655/81‑1 (MQ=255)
                                      ccgccAGGAGAGGAGCCCCAGGCACCGGTATTTCCAGTCAATGAGCTATAGCGCACCAATAACCAACCCACAATGGCGCTg                                   >  1:153705/1‑81 (MQ=255)
                                       cgccAGGAGAGGAGCCCCAGGCACCGGTATTTCCAGTCAATGAGCTATAGCGCACCAATAACCAACCCACAATGGCGCTgg                                  <  2:153656/81‑1 (MQ=255)
                                           aggagaggagCCCCAGGCACCGGTATTTCCAGTCAATGAGCTATAGCGCACCAATAACCAACCCACAATGGCGCTGGAAag                              >  1:153711/1‑81 (MQ=255)
                                               gaggagCCCCAGGCACCGGTATTTCCAGTCAATGAGCTATAGCGCACCAATAACCAACCCACAATGGCGCTGGAAagcagc                          >  2:153697/1‑81 (MQ=255)
                                                 ggagCCCCAGGCACCGGTATTTCCAGTCAATGAGCTATAGCGCACCAATAACCAACCCACAATGGCGCTGGAAAGCAGCGt                        >  1:806754/1‑81 (MQ=255)
                                                   agCCCCAGGCACCGGTATTTCCAGTCAATGAGCTATAGCGCACCAATAACCAACCCACAATGGCGCTGGAAAGCAGCGTTa                      >  1:153699/1‑81 (MQ=255)
                                                      cccAGGCACCGGTATTTCCAGTCAATGAGCTATAGCGCACCAATAACCAACCCACAATGGCGCTGGAAAGCAGCGTTAc                     >  2:153725/1‑79 (MQ=255)
                                                              ccGGTATTTCCAGTCAATGAGCTATAGCGCACCAATAACCAACCCACAATGGCGCTGGAAAGCAGCGTTACCAGTAAAAtt           <  2:153657/81‑1 (MQ=255)
                                                                ggTATTTCCAGTCAATGAGCTATAGCGCACCAATAACCAACCCACAATGGCGCTGGAAAGCAGCGTTACCAGTAAAATTg          >  2:153708/1‑80 (MQ=255)
                                                                 gTATTTCCAGTCAATGAGCTATAGCGCACCAATAACCAACCCACAATGGCGCTGGAAAGCAGCGTTACCAGTAAAATTGCa        <  1:153658/81‑1 (MQ=255)
                                                                    tttCCAGTCAATGAGCTATAGCGCACCAATATCCAACCCACAATGGCGCTGGAAAGCAGCGTTACCAGTAAAATTGCAAg      >  2:153716/1‑80 (MQ=255)
                                                                    tttCCAGTCAATGAGCTATAGCGCACCAATAACCAACCCACAATGGCGCTGGAAAGCAGCGTTACCAGTAAAATTGCAAg      <  1:153659/80‑1 (MQ=255)
                                                                     ttCCAGTCAATGAGCTATAGCGCACCAATAACCAACCCACAATGGCGCTGGAAAGCAGCGTTACCAGTAAAATTGCAAGCa    >  1:153707/1‑81 (MQ=255)
                                                                      tCCAGTCAATGAGCTATAGCGCACCAATAACCAACCCACAATGGCGCTGGAAAGCAGCGTTACCAGTAAAATTGCAAGCAc   >  1:153703/1‑81 (MQ=255)
                                                                       ccAGTCAATGAGCTATAGCGCACCAATAACCAACCCACAATGGCGCTGGAAAGCAGCGTTACCAGTAAAATTGCAAGCAc   <  2:153660/80‑1 (MQ=255)
                                                                        cAGTCAATGAGCTATAGCGGACCAATAACCAACCCACAATGGCGCTGGAAAGCAGCGTTACCAGTAAAATTGCAAGCACg  <  1:153661/80‑1 (MQ=255)
                                                                        cAGTCAATGAGCTATAGCGCACCAATAACCAACCCACAATGGCGCTGGAAAGCAGCGTTACCAGTAAAATTGCAAGCACg  <  1:153662/80‑1 (MQ=255)
                                                                           |                                                                            
CGCCGTAATCTTGCGCCATAGCGACCATTGCCGCCGCGCCGCCAGGAGAGGAGCCCCAGGCACCGGTATTTCCAGGCAATGAGCTATAGCGCACCAATAACCAACCCACAATGGCGCTGGAAAGCAGCGTTACCAGTAAAATTGCAAGCACG  >  NC_000913/747370‑747521

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