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%	168.3 / NA	52	*128* (TAG→TAA)	marA	multiple antibiotic resistance transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (29/23);  total (29/23)

GGCTTCGAGTCGCAACAAACTCTGACCCGAACCTTCAAAAATTACTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAGTTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCACGCAGCCAGTTGTTACTTCTTGTGCCAA  >  NC_000913/1619826‑1620089                                                                                                                                    |                                                                                                                                     ggCTTCGAGTCGCAACAAACTCTGACCCGAACCTTCAAAAATTACTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGaaaa                                                                                                                               <  2:550007/139‑1 (MQ=255)    ttCGAGTCGCAACAAACTCTGACCCGAACCTTCAAAAATTACTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGt                                                                                                                            <  2:90481/139‑1 (MQ=255)            gCAACAAACTCTGACCCGAACCTTCAAAAATTACTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGt                                                                                                                    >  1:1666112/1‑139 (MQ=255)                 aaaCTCTGACCCGAACCTTCAAAAATTACTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTg                                                                                                               <  2:1421267/139‑1 (MQ=255)                               aCCTTCAAAAATTACTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGaaa                                                                                                 <  2:1666112/139‑1 (MQ=255)                                   tCAAAAATTACTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCAc                                                                                             >  1:85125/1‑139 (MQ=255)                                                  aTGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGc                                                                                                           >  2:587051/1‑110 (MQ=255)                                                  aTGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGc                                                                                                           <  1:587051/110‑1 (MQ=255)                                                      tCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATagcagc                                                                          <  1:1055305/139‑1 (MQ=255)                                                        cgccgcATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATagcagc                                                                          >  1:354139/1‑137 (MQ=255)                                                        cgccgcATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATagcagc                                                                          <  2:354139/137‑1 (MQ=255)                                                        cgccgcATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTg                                                                        <  2:1294068/139‑1 (MQ=255)                                                        cgccgcATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTg                                                                        >  1:419461/1‑139 (MQ=255)                                                          ccgcATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTgcg                                                                      >  2:442085/1‑139 (MQ=255)                                                                        gATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCAt                                                                                                     <  2:831202/94‑1 (MQ=255)                                                                        gATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCAt                                                                                                     >  1:831202/1‑94 (MQ=255)                                                                           gACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACaa                                                                                                                       <  2:356834/73‑1 (MQ=255)                                                                           gACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACaa                                                                                                                       >  1:356834/1‑73 (MQ=255)                                                                            aCCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGTAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCcgcg                                                    >  1:796827/1‑139 (MQ=255)                                                                            aCCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCcgcg                                                    >  2:638381/1‑139 (MQ=255)                                                                             ccAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACaa                                                                                                                       <  1:1209061/71‑1 (MQ=255)                                                                             ccAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACaa                                                                                                                       >  2:1209061/1‑71 (MQ=255)                                                                                 tatGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAggg                                               >  2:1240434/1‑139 (MQ=255)                                                                                        ggCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGGTTATTCTCTTTTCCGCGCAGGGCGTTGCg                                        >  1:1146382/1‑139 (MQ=255)                                                                                        ggCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCg                                        >  2:941048/1‑139 (MQ=255)                                                                                        ggCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCg                                        >  2:455916/1‑139 (MQ=255)                                                                                         gCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGt                                                                                                                            >  1:1686780/1‑54 (MQ=255)                                                                                         gCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGt                                                                                                                            <  2:1686780/54‑1 (MQ=255)                                                                                             aTCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAgg                                                <  1:551827/126‑1 (MQ=255)                                                                                             aTCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAgg                                                >  2:551827/1‑126 (MQ=255)                                                                                                  gcTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGaaa                                                                                                 >  2:260115/1‑72 (MQ=255)                                                                                                  gcTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGaaa                                                                                                 <  1:260115/72‑1 (MQ=255)                                                                                                    tttttACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCACg                            <  2:85125/139‑1 (MQ=255)                                                                                                              cATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCtt                                                                   >  2:713117/1‑90 (MQ=255)                                                                                                              cATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCtt                                                                   <  1:713117/90‑1 (MQ=255)                                                                                                                ttAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATagcagc                                                                          >  1:531384/1‑81 (MQ=255)                                                                                                                ttAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATagcagc                                                                          <  2:531384/81‑1 (MQ=255)                                                                                                                ttAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCcgcg                                                    >  2:407380/1‑103 (MQ=255)                                                                                                                ttAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCcgcg                                                    <  1:407380/103‑1 (MQ=255)                                                                                                                ttAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCACGCAGCCAgttgtt                >  1:1303049/1‑139 (MQ=255)                                                                                                                ttAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCACGCAGCCAgttgtt                >  2:907274/1‑139 (MQ=255)                                                                                                                ttAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCACGCAGCCAgttgtt                >  2:166282/1‑139 (MQ=255)                                                                                                                ttAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCACGCAGCCAgttgtt                >  2:1458341/1‑139 (MQ=255)                                                                                                                        ttACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCACGCAg                         <  2:1677604/122‑1 (MQ=255)                                                                                                                        ttACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCACGCAg                         >  1:1677604/1‑122 (MQ=255)                                                                                                                            aacaGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCACGCAGCCAGTTGTTACTTCTTGTGcc    <  2:419461/139‑1 (MQ=255)                                                                                                                             acaGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCACGCAGCCAGTTGTTACTTCTTGTGCCa   <  1:455916/139‑1 (MQ=255)                                                                                                                              caGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCACGCAGCCAGTTGTTACTTCTTGTGCCaa  >  1:36493/1‑139 (MQ=255)                                                                                                                                 cTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCa                                                                                        >  1:265834/1‑50 (MQ=255)                                                                                                                                 cTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCa                                                                                        <  2:265834/50‑1 (MQ=255)                                                                                                                                 cTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCtttt                                                         >  2:619080/1‑81 (MQ=255)                                                                                                                                 cTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCtttt                                                         <  1:619080/81‑1 (MQ=255)                                                                                                                                    |                                                                                                                                     GGCTTCGAGTCGCAACAAACTCTGACCCGAACCTTCAAAAATTACTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAGTTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCACGCAGCCAGTTGTTACTTCTTGTGCCAA  >  NC_000913/1619826‑1620089

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