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%	130.9 / NA	42	*17* (TAG→TAA)	hisL	his operon leader peptide
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (21/21);  total (21/21)

TATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTT  >  NC_000913/2089915‑2090175                                                                                                                                    |                                                                                                                                  tataAAAAAGCCCTTGCTTTATAAAGGGAAAGGGGAATAGGTTAAAAGAAATCAGTTGAATAAACATTCACAGAGACTTTAATGAAAAGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATATTTCa                                                                                                                            >  1:585193/1‑139 (MQ=255)           cccTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGAtgtgt                                                                                                                  >  1:198498/1‑139 (MQ=255)                   ttCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAg                                                                                                          <  2:198498/139‑1 (MQ=255)                           tGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAg                                                                                                  >  1:28382/1‑139 (MQ=255)                           tGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAg                                                                                                  >  2:641196/1‑139 (MQ=255)                                    tttAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAg                                                                                                                           >  2:4320/1‑105 (MQ=255)                                    tttAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAg                                                                                                                           <  1:4320/105‑1 (MQ=255)                                    tttAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAg                                                                                         <  2:491493/139‑1 (MQ=255)                                    tttAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAg                                                                                         <  1:641196/139‑1 (MQ=255)                                    tttAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAg                                                                                         <  1:327664/139‑1 (MQ=255)                                        ggTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAgtggt                                                                                     <  2:365326/139‑1 (MQ=255)                                         gTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCt                                                                                                               <  1:394568/112‑1 (MQ=255)                                         gTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCt                                                                                                               >  2:394568/1‑112 (MQ=255)                                            aaaaGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCAt                                                                                 <  2:479358/139‑1 (MQ=255)                                                gACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAAc                                                                             <  1:259301/139‑1 (MQ=255)                                                             tAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGcc                                                                >  2:86948/1‑139 (MQ=255)                                                               aaCATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGcccc                                                              <  2:212724/139‑1 (MQ=255)                                                                  aTTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCgg                                                           >  2:37295/1‑139 (MQ=255)                                                                          agaCTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGaa                                                                              >  1:338233/1‑112 (MQ=255)                                                                          agaCTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGaa                                                                              <  2:338233/112‑1 (MQ=255)                                                                           gaCTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCAcc                                                  >  2:74607/1‑139 (MQ=255)                                                                                  aTGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCggg                                           <  2:250127/139‑1 (MQ=255)                                                                                      cacGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCt                                       >  1:172973/1‑139 (MQ=255)                                                                                        cgcgTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAg                                                                                                          >  1:740310/1‑70 (MQ=255)                                                                                        cgcgTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAg                                                                                                          <  2:740310/70‑1 (MQ=255)                                                                                        cgcgTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATgaga                                                                     <  1:734619/107‑1 (MQ=255)                                                                                        cgcgTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATgaga                                                                     >  2:734619/1‑107 (MQ=255)                                                                                                         ccaccaTCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGAt                    >  2:336612/1‑139 (MQ=255)                                                                                                          caccatcatcaccatcatCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATa                   <  2:115281/139‑1 (MQ=255)                                                                                                          caccatcatcaccatcatCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATa                   <  2:464836/139‑1 (MQ=255)                                                                                                                  tcaccatcatcCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCa           <  1:336612/139‑1 (MQ=255)                                                                                                                    accatcatcCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGa         <  1:86948/139‑1 (MQ=255)                                                                                                                    accatcatcCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGa         <  1:74607/139‑1 (MQ=255)                                                                                                                    accatcatcCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGa         >  2:443638/1‑139 (MQ=255)                                                                                                                         catcCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGt    >  2:93384/1‑139 (MQ=255)                                                                                                                           tcCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGaa                                                                              <  1:42307/63‑1 (MQ=255)                                                                                                                           tcCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGaa                                                                              >  2:42307/1‑63 (MQ=255)                                                                                                                           tcCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGttt  >  2:689683/1‑139 (MQ=255)                                                                                                                           tcCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGttt  >  2:519437/1‑139 (MQ=255)                                                                                                                           tcCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGttt  >  1:293081/1‑139 (MQ=255)                                                                                                                            ccTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAAATCTTCCAg                                                                                         >  1:494680/1‑51 (MQ=39)                                                                                                                            ccTGACTAACCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAg                                                                                         <  2:494680/51‑1 (MQ=255)                                                                                                                                    |                                                                                                                                  TATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTT  >  NC_000913/2089915‑2090175

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