Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,619,957 G→A *128* (TAG→TAA marA → multiple antibiotic resistance transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,619,9570GA100.0% 91.7 / NA 30*128* (TAG→TAAmarAmultiple antibiotic resistance transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (17/13);  total (17/13)

AGTCGCAACAAACTCTGACCCGAACCTTCAAAAATTACTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAGTTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCACGCAGCCAGTTGTTACTTCTTGTGCCAATG  >  NC_000913/1619833‑1620091
                                                                                                                            |                                                                                                                                      
aGTCGCAACAAACTCTGACCCGAACCTTCAAAAATTACTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACa                                                                                                                          >  1:260558/1‑139 (MQ=255)
aGTCGCAACAAACTCTGACCCGAACCTTCAAAAATTACTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACa                                                                                                                          >  2:721381/1‑139 (MQ=255)
       acaaACTCTGACCCGAACCTTCAAAAATTACTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCAc                                                                                                                   <  1:721381/139‑1 (MQ=255)
            ctctGACCCGAACCTTCAAAAATTACTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAgg                                                                                                              <  1:49600/139‑1 (MQ=255)
             tctGACCCGAACCTTCAAAAATTACTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGc                                                                                                             <  2:605502/139‑1 (MQ=255)
             tctGACCCGAACCTTCAAAAATTACTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGc                                                                                                             <  2:420388/139‑1 (MQ=255)
              ctGACCCGAACCTTCAAAAATTACTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCa                                                                                                            >  2:910598/1‑139 (MQ=255)
                  cccGAACCTTCAAAAATTACTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAACTCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCa                                                                                                        >  2:24264/1‑139 (MQ=255)
                  cccGAACCTTCAAAAATTACTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCa                                                                                                        <  1:910598/139‑1 (MQ=255)
                                                     cATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCtt                                                                     <  1:355494/139‑1 (MQ=255)
                                                       tAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTAt                                                                   >  1:324431/1‑139 (MQ=255)
                                                                    aCCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCcgcg                                                      >  2:425005/1‑139 (MQ=255)
                                                                     ccAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAAc                                                                                                                        <  2:858711/72‑1 (MQ=255)
                                                                     ccAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAAc                                                                                                                        >  1:858711/1‑72 (MQ=255)
                                                                      cAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCa                                                    <  2:123013/139‑1 (MQ=255)
                                                                               gggCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGc                                           >  2:9684/1‑139 (MQ=255)
                                                                                   gAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGaa                                       >  1:146773/1‑139 (MQ=255)
                                                                                    aaTCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCa                                                                                                            >  1:871422/1‑69 (MQ=255)
                                                                                    aaTCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCa                                                                                                            <  2:871422/69‑1 (MQ=255)
                                                                                               ttACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCACGCAg                           <  2:260558/139‑1 (MQ=255)
                                                                                                        ttAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCAt                                                                                                       >  2:447413/1‑54 (MQ=255)
                                                                                                        ttAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCAt                                                                                                       <  1:447413/54‑1 (MQ=255)
                                                                                                        ttAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCa                                                                                          >  1:7729/1‑67 (MQ=255)
                                                                                                        ttAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCa                                                                                          <  2:7729/67‑1 (MQ=255)
                                                                                                        ttAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCACGCAGCCAgttgtt                  >  2:382522/1‑139 (MQ=255)
                                                                                                        ttAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATCAAACCACTTTCATCCGCAATAGCAGCTg                                                                          >  2:247222/1‑83 (MQ=255)
                                                                                                        ttAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATCAAACCACTTTCATCCGCAATAGCAGCTg                                                                          <  1:247222/83‑1 (MQ=255)
                                                                                                                 tACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCACGCAGCCAGTTGTTACTTCTtgt         >  1:712454/1‑139 (MQ=255)
                                                                                                                 tACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCACGCAGCCAGTTGTTACTTCTtgt         >  1:469104/1‑139 (MQ=255)
                                                                                                                        gCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCACGCAGCCAGTTGTTACTTCTTGTGCCAATg  >  2:399577/1‑139 (MQ=255)
                                                                                                                            |                                                                                                                                      
AGTCGCAACAAACTCTGACCCGAACCTTCAAAAATTACTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAGTTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCACGCAGCCAGTTGTTACTTCTTGTGCCAATG  >  NC_000913/1619833‑1620091

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