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

TATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAG  >  NC_000913/2089915‑2090181                                                                                                                                    |                                                                                                                                        tataAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCa                                                                                                                                  >  1:355176/1‑139 (MQ=255) tataAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCa                                                                                                                                  >  1:326005/1‑139 (MQ=255)   taAAAAAGCCCATGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAgg                                                                                                                                <  2:1046585/139‑1 (MQ=255)      aaaaGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGa                                                                                                                             >  1:1221627/1‑139 (MQ=255)       aaaGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGAt                                                                                                                            <  2:1221627/139‑1 (MQ=255)       aaaGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGAt                                                                                                                            <  2:49952/139‑1 (MQ=255)       aaaGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGAt                                                                                                                            <  2:1527753/139‑1 (MQ=255)         aGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGAtgt                                                                                                                          >  2:167946/1‑139 (MQ=255)           cccTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGAtgtgt                                                                                                                        >  2:1324811/1‑139 (MQ=255)                 cTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGa                                                                                                                  <  2:326005/139‑1 (MQ=255)                   ttCTAAGGTGAAAGTGTTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAg                                                                                                                <  1:1426371/139‑1 (MQ=255)                   ttCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAg                                                                                                                <  1:167946/139‑1 (MQ=255)                     cTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGAc                                                                                                              >  1:175818/1‑139 (MQ=255)                           tGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAg                                                                                                        >  1:1212869/1‑139 (MQ=255)                                   gTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGaa                                                                                                                 <  2:719796/122‑1 (MQ=255)                                   gTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGaa                                                                                                                 >  1:719796/1‑122 (MQ=255)                                    tttAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAg                                                                                               <  1:1324811/139‑1 (MQ=255)                                       aGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTcc                                                                                                 <  1:1410502/134‑1 (MQ=255)                                       aGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTcc                                                                                                 >  2:1410502/1‑134 (MQ=255)                                            aaaaGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCAt                                                                                       <  2:175818/139‑1 (MQ=255)                                             aaaGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATg                                                                                      >  2:1406106/1‑139 (MQ=255)                                             aaaGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCGTCCAGTGGTGCATg                                                                                      <  1:1406106/139‑1 (MQ=255)                                               aGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGaa                                                                                    <  2:1423373/139‑1 (MQ=255)                                                   aTCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGAt                                                                                                                            >  2:1592064/1‑95 (MQ=255)                                                   aTCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGAt                                                                                                                            <  1:1592064/95‑1 (MQ=255)                                                        ttGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATgaga                                                                           <  2:355176/139‑1 (MQ=255)                                                          gAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGaaa                                                                         <  2:754330/139‑1 (MQ=255)                                                               aaCATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGcccc                                                                    <  2:973554/139‑1 (MQ=255)                                                                aCATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCAGGAGAAAGccccc                                                                   <  1:983844/139‑1 (MQ=255)                                                                  aTTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGAtgtgtg                                                                                                                       <  1:1469145/85‑1 (MQ=255)                                                                  aTTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGAtgtgtg                                                                                                                       >  2:1469145/1‑85 (MQ=255)                                                                  aTTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGAtgtg                                                                                                                         >  1:105809/1‑83 (MQ=255)                                                                  aTTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGAtgtg                                                                                                                         <  2:105809/83‑1 (MQ=255)                                                                          agaCTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGa                                                                                                       >  1:1415643/1‑93 (MQ=255)                                                                          agaCTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGa                                                                                                       <  2:1415643/93‑1 (MQ=255)                                                                          agaCTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGa                                                                >  1:1443746/1‑132 (MQ=255)                                                                          agaCTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGa                                                                <  2:1443746/132‑1 (MQ=255)                                                                           gaCTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGAtg                                                                                                                           >  1:1350137/1‑72 (MQ=255)                                                                           gaCTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGAtg                                                                                                                           <  2:1350137/72‑1 (MQ=255)                                                                                      cacGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGAt                                                                                                                            >  1:507678/1‑60 (MQ=255)                                                                                      cacGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGAt                                                                                                                            <  2:507678/60‑1 (MQ=255)                                                                                        cgcgTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGAtgtg                                                                                                                         >  1:1480771/1‑61 (MQ=255)                                                                                        cgcgTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGAtgtg                                                                                                                         <  2:1480771/61‑1 (MQ=255)                                                                                        cgcgTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGAt                                                                                                                            <  1:1513209/58‑1 (MQ=255)                                                                                        cgcgTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGAt                                                                                                                            >  2:1513209/1‑58 (MQ=255)                                                                                        cgcgTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAg                                                                                               >  1:236828/1‑87 (MQ=255)                                                                                        cgcgTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAg                                                                                               <  2:236828/87‑1 (MQ=255)                                                                                          cgTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAg                                                                                               <  2:1385052/85‑1 (MQ=255)                                                                                          cgTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAg                                                                                               >  1:1385052/1‑85 (MQ=255)                                                                                          cgTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCttttt                                         >  2:325367/1‑139 (MQ=255)                                                                                            ttCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCTAg                                                                                               >  2:1427047/1‑83 (MQ=255)                                                                                            ttCAATTTAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAg                                                                                               <  1:1427047/83‑1 (MQ=255)                                                                                                 tttAAACACCACCATCATCACCATCATCCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTgcg                                  <  1:325367/139‑1 (MQ=255)                                                                                                             catcatcaccatcatcCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCggg                                                 <  1:386149/112‑1 (MQ=255)                                                                                                             catcatcaccatcatcCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCggg                                                 >  2:386149/1‑112 (MQ=255)                                                                                                                        tcatcCTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGAt                                                            >  1:901501/1‑90 (MQ=255)                                                                                                                        tcatcCTGACTAATCTTTAAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGAt                                                            <  2:901501/90‑1 (MQ=255)                                                                                                                             cTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGa               <  1:1546681/130‑1 (MQ=255)                                                                                                                             cTGACTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGa               >  2:1546681/1‑130 (MQ=255)                                                                                                                                aCTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTgcg                                  <  1:314609/108‑1 (MQ=255)                                                                                                                                aCTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTgcg                                  >  2:314609/1‑108 (MQ=255)                                                                                                                                aCTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGa               <  2:917745/127‑1 (MQ=255)                                                                                                                                aCTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGa               >  1:917745/1‑127 (MQ=255)                                                                                                                                 cTAATCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAagag  <  2:1212869/139‑1 (MQ=255)                                                                                                                                    |                                                                                                                                        TATAAAAAAGCCCTTGCTTTCTAACGTGAAAGTGGTTTAGGTTAAAAGACATCAGTTGAATAAACATTCACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAG  >  NC_000913/2089915‑2090181

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