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_000913	1,619,957	0	G	A	100.0%	164.8 / NA	52	*128* (TAG→TAA)	marA	multiple antibiotic resistance transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (27/25);  total (27/25)

TCGAGTCGCAACAAACTCTGACCCGAACCTTCAAAAATTACTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAGTTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCACGCAGCCAGTTGTTACTTCTTGT  >  NC_000913/1619830‑1620084                                                                                                                                |                                                                                                                                tCGAGTCGCAACAAACTCTGACCCGAACCTTCAAAAATTACTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTg                                                                                                                      >  1:1292380/1‑139 (MQ=255)    aGTCGCAACAAACTCTGACCCGAACCTTCAAAAATTACTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACa                                                                                                                   >  2:646034/1‑139 (MQ=255)     gTCGCAACAAACTCTGACCCGAACCTTCAAAAATTACTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACaa                                                                                                                  <  2:1194064/139‑1 (MQ=255)                  ctGACCCGAACCTTCAAAAATTACTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCa                                                                                                     <  2:500317/139‑1 (MQ=255)                    gACCCGAACCTTCAAAAATTACTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACaa                                                                                                                  <  1:773715/124‑1 (MQ=255)                    gACCCGAACCTTCAAAAATTACTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACaa                                                                                                                  >  2:773715/1‑124 (MQ=255)                       ccGAACCTTCAAAAATTACTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCAt                                                                                                <  2:655784/139‑1 (MQ=255)                          aaCCTTCAAAAATTACTTTGATGTTCCGTCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAAc                                                                                                                 <  2:1232251/119‑1 (MQ=255)                          aaCCTTCAAAAATTACTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAAc                                                                                                                 >  1:1232251/1‑119 (MQ=255)                               tCAAAAATTACTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCAc                                                                                        >  2:123645/1‑139 (MQ=255)                                                   ccgccgCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCt                                                                    <  1:151437/139‑1 (MQ=255)                                                        gcATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAgg                                                                                                       <  1:766751/99‑1 (MQ=255)                                                        gcATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAgg                                                                                                       >  2:766751/1‑99 (MQ=255)                                                         cATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCtt                                                              >  2:913963/1‑139 (MQ=255)                                                         cATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCtt                                                              <  2:1292380/139‑1 (MQ=255)                                                                   ggATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGaa                                                                                                                            >  2:1089213/1‑67 (MQ=255)                                                                   ggATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGaa                                                                                                                            <  1:1089213/67‑1 (MQ=255)                                                                    gATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAAt                                                                                                   >  1:1235116/1‑91 (MQ=255)                                                                    gATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAAt                                                                                                   <  2:1235116/91‑1 (MQ=255)                                                                    gATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTc                                                   >  2:452508/1‑139 (MQ=255)                                                                      tGACCAATATGCAGGGCGCATGGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCcg                                                 >  1:661572/1‑139 (MQ=255)                                                                        aCCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCcgcg                                               >  1:133360/1‑139 (MQ=255)                                                                         ccAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCgcgc                                              >  1:1255140/1‑139 (MQ=255)                                                                                  aGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTg                                     <  2:1255140/139‑1 (MQ=255)                                                                                    ggCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAAtt                                                                                                                               <  1:1059160/47‑1 (MQ=255)                                                                                    ggCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCa                                                                                                     >  2:1175028/1‑73 (MQ=255)                                                                                    ggCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCa                                                                                                     <  1:1175028/73‑1 (MQ=255)                                                                                    ggCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATagcagc                                                                     <  1:1355868/105‑1 (MQ=255)                                                                                    ggCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATagcagc                                                                     >  2:1355868/1‑105 (MQ=255)                                                                                    ggCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCg                                   >  1:1019697/1‑139 (MQ=255)                                                                                    ggCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCg                                   >  1:731112/1‑139 (MQ=255)                                                                                          tCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACaa                             <  1:913963/139‑1 (MQ=255)                                                                                                tttttACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCACg                       <  1:452508/139‑1 (MQ=255)                                                                                                    tACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCACGCAGc                   >  2:1082735/1‑139 (MQ=255)                                                                                                      cATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCACGCAGCCa                 <  2:133360/139‑1 (MQ=255)                                                                                                         ccATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCACGCAGCCAgtt              <  2:1019697/139‑1 (MQ=255)                                                                                                            ttAAATCATTACAACAGCTAATTGAAAACGTGACaa                                                                                                                  <  2:789984/36‑1 (MQ=255)                                                                                                            ttAAATCATTACAACAGCTAATTGAAAACGTGACaa                                                                                                                  >  1:789984/1‑36 (MQ=255)                                                                                                            ttAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATcc                                                                                >  1:665200/1‑70 (MQ=255)                                                                                                            ttAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATcc                                                                                <  2:665200/70‑1 (MQ=255)                                                                                                            ttAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTg                                                                   >  1:633671/1‑83 (MQ=255)                                                                                                            ttAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTg                                                                   <  2:633671/83‑1 (MQ=255)                                                                                                            ttAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTg                                     >  2:1316851/1‑113 (MQ=255)                                                                                                            ttAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTg                                     <  1:1316851/113‑1 (MQ=255)                                                                                                                    ttACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCAc                        >  1:541901/1‑118 (MQ=255)                                                                                                                    ttACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCAc                        <  2:541901/118‑1 (MQ=255)                                                                                                                    ttACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCACGCAGCCAGTTGTTACTTCTtg   >  2:763790/1‑139 (MQ=255)                                                                                                                     tACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTg                                     <  1:472574/104‑1 (MQ=255)                                                                                                                     tACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTg                                     >  2:472574/1‑104 (MQ=255)                                                                                                                     tACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCACGCAGCCAGTTGTTACTTCTtgt  >  2:527040/1‑139 (MQ=255)                                                                                                                     tACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCACGCAGCCAGTTGTTACTTCTtgt  >  2:1259081/1‑139 (MQ=255)                                                                                                                             cTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCAGTTTCATc                                                                                 <  2:936181/52‑1 (MQ=255)                                                                                                                                |                                                                                                                                TCGAGTCGCAACAAACTCTGACCCGAACCTTCAAAAATTACTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAGTTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCACGCAGCCAGTTGTTACTTCTTGT  >  NC_000913/1619830‑1620084

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