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

TCGAGTCGCAACAAACTCTGACCCGAACCTTCAAAAATTACTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAGTTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCACGCAGCCAGTTGTTACTTCTTGTGCCAATG  >  NC_000913/1619830‑1620091                                                                                                                                |                                                                                                                                       tCGAGTCGCAACAAACTCTGACCCGAACCTTCAAAAATTACTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTg                                                                                                                             >  2:1344907/1‑139 (MQ=255)  cGAGTCGCAACAAACTCTGACCCGAACCTTCAAAAATTACTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGa                                                                                                                            <  1:1344907/139‑1 (MQ=255)   gAGTCGCAACAAACTCTGACCCGAACCTTCAAAATTTACTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGAc                                                                                                                           <  1:252895/139‑1 (MQ=255)    aGTCGCAACAAACTCTGACCCGAACCTTCAAAAATTACTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGAcc                                                                                                                          >  1:1177527/1‑138 (MQ=255)         gaacaaACTCTGACCCGAACCTTCAAAAATTACTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTc                                                                                                                     >  1:918396/2‑139 (MQ=255)                       ccGAACCTTCAAAAATTACTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCCCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAATGTCACTGAGGCAATCAt                                                                                                       <  2:918396/139‑1 (MQ=255)                                                    cgccgcATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTg                                                                          <  1:1075849/139‑1 (MQ=255)                                                      ccgcATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTgcg                                                                        >  1:1420993/1‑139 (MQ=255)                                                         taTAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCtt                                                                     <  1:57051/138‑1 (MQ=255)                                                         cATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCtt                                                                     <  2:859713/139‑1 (MQ=255)                                                                   ggATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATg                                                                                                      >  2:404618/1‑96 (MQ=255)                                                                   ggATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATg                                                                                                      <  1:404618/96‑1 (MQ=255)                                                                    gATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATcc                                                                                       <  1:1358941/110‑1 (MQ=255)                                                                    gATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATcc                                                                                       >  2:1358941/1‑110 (MQ=255)                                                                        aCCAATATGTAGGGCGAATAGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGa                                                                                                                                    <  1:906041/61‑1 (MQ=255)                                                                        aCCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGa                                                                                                                                    >  2:906041/1‑61 (MQ=255)                                                                         ccAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCa                                                                                                            >  1:672956/1‑84 (MQ=255)                                                                         ccAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCa                                                                                                            <  2:672956/84‑1 (MQ=255)                                                                            atatGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAgg                                                  >  1:481147/1‑139 (MQ=255)                                                                                  aGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTg                                            <  2:1420993/139‑1 (MQ=255)                                                                                  aGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTg                                            >  1:1068703/1‑139 (MQ=255)                                                                                   gggCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTg                                                                                                                 >  1:1460609/1‑69 (MQ=255)                                                                                   gggCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTg                                                                                                                 <  2:1460609/69‑1 (MQ=255)                                                                                   gggCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCa                                                                                          >  1:109712/1‑92 (MQ=255)                                                                                   gggCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCa                                                                                          <  2:109712/92‑1 (MQ=255)                                                                                    ggCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACaa                                                                                                                         >  2:1369187/1‑60 (MQ=255)                                                                                    ggCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACaa                                                                                                                         <  1:1369187/60‑1 (MQ=255)                                                                                    ggCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCg                                          >  2:558745/1‑139 (MQ=255)                                                                                    ggCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCg                                          >  1:1297517/1‑139 (MQ=255)                                                                                      cGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGa                                        >  1:503275/1‑139 (MQ=255)                                                                                         aTCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACa                                     >  1:137598/1‑139 (MQ=255)                                                                                         aTCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACa                                     <  2:1297517/139‑1 (MQ=255)                                                                                                            ttAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGc                                                                                                             <  2:245983/48‑1 (MQ=255)                                                                                                            ttAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGc                                                                                                             >  1:245983/1‑48 (MQ=255)                                                                                                            ttAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCtt                                                                     <  2:1016315/88‑1 (MQ=255)                                                                                                            ttAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCtt                                                                     >  1:1016315/1‑88 (MQ=255)                                                                                                            ttAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCACGCAGCCAgttgtt                  >  1:1427961/1‑139 (MQ=255)                                                                                                             tAAATCATTACAACAGCTAATTGAAAACGTGACaa                                                                                                                         >  2:823799/1‑35 (MQ=255)                                                                                                             tAAATCATTACAACAGCTAATTGAAAACGTGACaa                                                                                                                         <  1:823799/35‑1 (MQ=255)                                                                                                                    ttACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGaaa                                                                                                   <  2:175393/50‑1 (MQ=255)                                                                                                                    ttACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGaaa                                                                                                   >  1:175393/1‑50 (MQ=255)                                                                                                                     tACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCtt                                                                     >  2:44864/1‑79 (MQ=255)                                                                                                                     tACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCtt                                                                     <  1:44864/79‑1 (MQ=255)                                                                                                                     tACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCgcgc                                                     >  1:765945/1‑95 (MQ=255)                                                                                                                     tACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCgcgc                                                     <  2:765945/95‑1 (MQ=255)                                                                                                                     tACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCACGCAGCCAGTTGTTACTTCTtgt         >  1:871109/1‑139 (MQ=255)                                                                                                                     tACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCACGCAGCCAGTTGTTACTTCTtgt         >  2:801988/1‑139 (MQ=255)                                                                                                                          caGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCACGCAg                           >  2:104994/1‑116 (MQ=255)                                                                                                                          caGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCACGCAg                           <  1:104994/116‑1 (MQ=255)                                                                                                                            gCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCACGCAGCCAGTTGTTACTTCTTGTGCCAATg  >  1:133199/1‑139 (MQ=255)                                                                                                                                |                                                                                                                                       TCGAGTCGCAACAAACTCTGACCCGAACCTTCAAAAATTACTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAGTTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCACGCAGCCAGTTGTTACTTCTTGTGCCAATG  >  NC_000913/1619830‑1620091

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