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

CTTCGAGTCGCAACAAACTCTGACCCGAACCTTCAAAAATTACTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAGTTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCACGCAGCCAGTTGTTACTTCTTGT  >  NC_000913/1619828‑1620084                                                                                                                                  |                                                                                                                                ctacGAGTCGCAACAAACTCTGACCCGAACCTTCAAAAATTACTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACg                                                                                                                        >  2:265280/4‑139 (MQ=255)      aGTCGCAACAAACTCTGACCCGAACCTTCAAAAATTACTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACa                                                                                                                   >  2:915607/1‑139 (MQ=255)          gCAACAAACTCTGACCCGAACCTTCAAAAATTACTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGt                                                                                                               >  2:237732/1‑139 (MQ=255)          gCAACAAACTCTGACCCGAACCTTCAAAAATTACTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGt                                                                                                               >  1:909813/1‑139 (MQ=255)                  ctctGACCCGAACCTTCAAAAATTACTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAgg                                                                                                       <  2:909813/139‑1 (MQ=255)                  ctctGACCCGAACCTTCAAAAATTACTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAgg                                                                                                       <  1:915607/139‑1 (MQ=255)                        cccGAACCTTCAAAAATTACTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGaa                                                                                                                            >  1:1379452/1‑112 (MQ=255)                        cccGAACCTTCAAAAATTACTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGaa                                                                                                                            <  2:1379452/112‑1 (MQ=255)                            aaCCTTCAAAAATTACTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGaa                                                                                             >  1:564064/1‑139 (MQ=255)                                      aaTTACTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCa                                                                                   <  1:1422079/139‑1 (MQ=255)                                          aCTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGt                                                                                                               >  2:608416/1‑107 (MQ=255)                                          aCTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGt                                                                                                               <  1:608416/107‑1 (MQ=255)                                                           cATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATc                                                                                 <  2:520107/120‑1 (MQ=255)                                                           cATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATc                                                                                 >  1:520107/1‑120 (MQ=255)                                                             tAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTAt                                                            >  2:211096/1‑139 (MQ=255)                                                                       aTGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTcc                                                  >  1:303247/1‑139 (MQ=255)                                                                          aCCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCcgcg                                               <  1:265280/139‑1 (MQ=255)                                                                                 tGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCg                                        >  1:224510/1‑139 (MQ=255)                                                                                    aGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGc                                                                                                      >  2:1553674/1‑74 (MQ=255)                                                                                    aGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGc                                                                                                      <  1:1553674/74‑1 (MQ=255)                                                                                    aGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTg                                     <  2:303247/139‑1 (MQ=255)                                                                                       gCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGt                                                                                                                       >  2:317684/1‑54 (MQ=255)                                                                                       gCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGt                                                                                                                       <  1:317684/54‑1 (MQ=255)                                                                                        cGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGACACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCgaa                                 >  2:1636899/1‑137 (MQ=255)                                                                                        cGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTgc                                                                  <  2:242162/106‑1 (MQ=255)                                                                                        cGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTgc                                                                  >  1:242162/1‑106 (MQ=255)                                                                                           aTCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAAc                                                                                           >  1:963469/1‑78 (MQ=255)                                                                                           aTCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCtttt                                                    >  2:283415/1‑117 (MQ=255)                                                                                           aTCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCtttt                                                    <  1:283415/117‑1 (MQ=255)                                                                                           aTCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGa                                 <  2:1399365/136‑1 (MQ=255)                                                                                           aTCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGa                                 >  1:1399365/1‑136 (MQ=255)                                                                                           aTCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACa                              >  2:697034/1‑139 (MQ=255)                                                                                           aTCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATCAAAc                                                                                           <  2:963469/78‑1 (MQ=255)                                                                                                   ttttACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCACGc                      <  1:697034/139‑1 (MQ=255)                                                                                                     ttACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCACGCAg                    <  1:512366/139‑1 (MQ=255)                                                                                                         aTCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCACGCAGCCAg                <  2:564064/139‑1 (MQ=255)                                                                                                           ccATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTAACGGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCGCGCAGCCAgtt              <  2:857361/139‑1 (MQ=255)                                                                                                              ttAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCcgcg                                               >  1:1051555/1‑103 (MQ=255)                                                                                                              ttAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCcgcg                                               <  2:1051555/103‑1 (MQ=255)                                                                                                              ttAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCACGCAGCCAgttgtt           >  2:1213979/1‑139 (MQ=255)                                                                                                              ttAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCACGCAGCCAgttgtt           >  2:1488483/1‑139 (MQ=255)                                                                                                               tAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCACGCAGCCAGTTGTTa          <  1:1213979/139‑1 (MQ=255)                                                                                                                      ttACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGc                                                                              >  2:495803/1‑64 (MQ=255)                                                                                                                      ttACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGc                                                                              <  1:495803/64‑1 (MQ=255)                                                                                                                      ttACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCACGCAGCCAGTTGTTACTTCTtg   >  1:1355308/1‑139 (MQ=255)                                                                                                                       tACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCACGCAGCCAGTTGTTACTTCTtgt  >  1:1435569/1‑139 (MQ=255)                                                                                                                               cTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCtt                                                              <  2:362033/71‑1 (MQ=255)                                                                                                                               cTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCACGCAGcc                  >  2:1395599/1‑115 (MQ=255)                                                                                                                               cTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCACGCAGCCa                 >  2:1592819/1‑116 (MQ=255)                                                                                                                                  |                                                                                                                                CTTCGAGTCGCAACAAACTCTGACCCGAACCTTCAAAAATTACTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAGTTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCACGCAGCCAGTTGTTACTTCTTGT  >  NC_000913/1619828‑1620084

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