Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,090,046	G→A	*17* (TAG→TAA)	hisL →	his operon leader peptide

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,090,046	0	G	A	100.0%	141.7 / NA	44	*17* (TAG→TAA)	hisL	his operon leader peptide
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (19/25);  total (19/25)

ATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAG  >  NC_000913/2089916‑2090179                                                                                                                                   |                                                                                                                                      ataAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAg                                                                                                                               >  1:436678/1‑139 (MQ=255)     aaaaGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGa                                                                                                                           >  1:408253/1‑139 (MQ=255)        aGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGAtgt                                                                                                                        >  2:289575/1‑139 (MQ=255)         gCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGAtgtg                                                                                                                       <  2:436678/139‑1 (MQ=255)             ttGCTTTCTAACGTGGAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCt                                                                                                                   >  2:759604/1‑139 (MQ=255)                 tttCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGaa                                                                                                               <  2:691607/139‑1 (MQ=255)                       aCGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCg                                                                                                                            >  1:15451/1‑120 (MQ=255)                       aCGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCg                                                                                                                            <  2:15451/120‑1 (MQ=255)                           gAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGa                                                                                                     <  1:270401/139‑1 (MQ=255)                                   tttAGGTTAAAAGACATCAGTTGAATATATATTCACAGAGACTTTTATGACACGTGTTCAATTTAAACACCACCATCATCACCATCATCCTGATTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAg                                                                                             <  1:243641/139‑1 (MQ=255)                                   tttAGGTTAAAAGACATCAGTTGAATAAACATTCACAGATACTTTTATGACACGCGTTCAATTTAACCACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACCTTCAGATCTTCCAg                                                                                             <  2:486157/139‑1 (MQ=255)                                   tttAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAg                                                                                             <  2:195339/139‑1 (MQ=255)                                   tttAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAg                                                                                             <  2:623849/139‑1 (MQ=255)                                            aaaGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAg                                                                                             <  2:614884/130‑1 (MQ=255)                                            aaaGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAg                                                                                             >  1:614884/1‑130 (MQ=255)                                            aaaGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATg                                                                                    >  2:616572/1‑139 (MQ=255)                                              aGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGaa                                                                                  <  2:494434/139‑1 (MQ=255)                                                  aTCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGAtgt                                                                                                                        <  2:527967/97‑1 (MQ=255)                                                  aTCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGAtgt                                                                                                                        >  1:527967/1‑97 (MQ=255)                                                  aTCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGaa                                                                                  <  2:265651/135‑1 (MQ=255)                                                  aTCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGaa                                                                                  >  1:265651/1‑135 (MQ=255)                                                  aTCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCa                                                                              >  1:597244/1‑139 (MQ=255)                                                  aTCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCa                                                                              >  1:417658/1‑139 (MQ=255)                                                       ttGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATgaga                                                                         <  1:289575/139‑1 (MQ=255)                                                                 aTTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCgg                                                               >  1:47734/1‑139 (MQ=255)                                                                    cacaGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACAt                                                                                                          <  1:293809/93‑1 (MQ=255)                                                                    cacaGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACAt                                                                                                          >  2:293809/1‑93 (MQ=255)                                                                         agaCTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGa                                                                                                     >  2:800059/1‑93 (MQ=255)                                                                         agaCTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGa                                                                                                     <  1:800059/93‑1 (MQ=255)                                                                           aCTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCt                                                     >  2:457041/1‑139 (MQ=255)                                                                             ttttATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGAt                                                                                                                          <  2:609527/68‑1 (MQ=255)                                                                             ttttATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGAt                                                                                                                          >  1:609527/1‑68 (MQ=255)                                                                                   gACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCggggg                                             <  2:417658/139‑1 (MQ=255)                                                                                                  tAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGcgcg                              <  1:616572/139‑1 (MQ=255)                                                                                                          accatcatcaccatcatcCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATaa                      <  1:23500/139‑1 (MQ=255)                                                                                                                   accatcatcCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGa             <  2:597244/139‑1 (MQ=255)                                                                                                                   accatcatcCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCAGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGa             <  1:75496/139‑1 (MQ=255)                                                                                                                          tcCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGa                                                              >  2:368112/1‑83 (MQ=255)                                                                                                                          tcCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGa                                                              <  1:368112/83‑1 (MQ=255)                                                                                                                              gACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAag  <  2:408253/139‑1 (MQ=255)                                                                                                                               aCTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTc                                                                                                <  1:763934/44‑1 (MQ=255)                                                                                                                               aCTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTc                                                                                                >  2:763934/1‑44 (MQ=255)                                                                                                                               aCTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATa                       <  2:59913/117‑1 (MQ=255)                                                                                                                               aCTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATa                       >  1:59913/1‑117 (MQ=255)                                                                                                                                   |                                                                                                                                      ATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAG  >  NC_000913/2089916‑2090179

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: ‑