Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,387,095	G→A	A126T (GCA→ACA)	tyrR →	aromatic amino acid biosynthesis and transport regulon transcriptional regulator, autorepressor, ATPase, phosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,387,095	0	G	A	100.0%	116.7 / NA	37	A126T (GCA→ACA)	tyrR	aromatic amino acid biosynthesis and transport regulon transcriptional regulator, autorepressor, ATPase, phosphatase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (19/18);  total (19/18)

GCGTTGAGCGCGTTACTGGAGGCGTTGCCTGAACCTGTGCTCTCTGTCGATATGAAAAGCAAAGTGGATATGGCGAACCCGGCGAGCTGTCAGCTTTTTGGGCAAAAATTGGATCGCCTGCGCAACCATACCGCCGCACAATTGATTAACGGCTTTAATTTTTTACGTTGGCTGGAAAGCGAACCGCAAGATTCGCATAACGAGCATGTCGTTATTAATGGGCAGAATTTCCTGATGGAGATTACGCCTGTTTATCTT  >  NC_000913/1386960‑1387217                                                                                                                                        |                                                                                                                           gCGTTGAGCGCGTTACTGGAGGCGTTGCCTGAACCTGTGCTCTCTGTCGATATGAAAAGCAAAGTGGATATGGCGAACCCGGCGAGCTGTCAGCTTTTTGGGCAAAAATTGGATCGCCTGCGCAACCATACCGCCacac                                                                                                                         <  2:1249282/139‑1 (MQ=255)           cgTTACTGGAGGCGTTGCCTGAACCTGTGCTCTCTGTCGATATGAAAAGCAAAGTGGATATGGCGAACCCGGCGAGCTGTCAGCTTTTTGGGCAAAAATTGGATCGCCTGCGCAACCATACCGCCACACAATTGATTaa                                                                                                               >  1:563313/1‑139 (MQ=255)                          tGCCTGAACCTGTGCTCTCTGTCGATATGAAAAGCAAAGTGGATATGGCGAACCCGGCGAGCTGTCAGCTTTTTGGGCAAAAATTGGATCGCCTGCGCAACCATACCGCCACACAATTGATTAACGGCTTTAAtttttt                                                                                                >  1:1031687/1‑139 (MQ=255)                                      tgCTCTCTGTCGATATGAAAAGCAAAGTGGATATGGCGAACCCGGCGAGCTGTCAGCTTTTTGGGCAAAAATTGGATCGCCTGCGCAACCATACCGCCACAGAATTGATTAATCGCTTTAATTTTTTACGTTGGCTGGa                                                                                    >  1:1025919/1‑139 (MQ=255)                                      tgCTCTCTGTCGATATGAAAAGCAAAGTGGATATGGCGAACCCGGCGAGCTGTCAGCTTTTTGGGCAAAAATTGGATCGCCTGCGCAACCATACCGCCACACAATTGATTAACGGCTTTAATTTTTTACGTTGGCTGGa                                                                                    >  2:991174/1‑139 (MQ=255)                                      tgCTCTCTGTCGATATGAAAAGCAAAGTGGATATGGCGAACCCGGCGAGCTGTCAGCTTTTTGGGCAAAAATTGGATCGCCTGCGCAACCATACCGCCACACAATTGATTAACGGCTTTAATTTTTTACGTTGGCTGGa                                                                                    >  1:649338/1‑139 (MQ=255)                                                cGATATGAAAAGCAAAGTGGATATGGCGAACCCGGCGAGCTGTCAGCTTTTTGGGCAAAAATTGGATCGCCTGCGCAACCATACCGCCACACAATTGATTAACGGCTTTAATTTTTTACGTTGGCTGGAAAGCGAACCg                                                                          <  1:991174/139‑1 (MQ=255)                                                       aaaaGCAAAGTGGATATGGCGAACCCGGCGAGCTGTCAGCTTTTTGGGCAAAAATTGGATCGCCTGCGCAACCATACCGCCACACAATTGATTAACGGCTTTAATTTTTTACGTTGGCTGGAAAGCGAACCGCAAGAtt                                                                   >  1:342603/1‑139 (MQ=255)                                                              aaGTGGATATGGCGAACCCGGCGAGCTGTCAGCTTTTTGGGCAAAAATTGGATCGCCTGCGCAACCATACCGCCACACAATTGATTAACGGCTTTAATTTTTTACGTTGGCTGGAAAGCGAACCGCAAGATTCGCATaa                                                            >  2:960786/1‑139 (MQ=255)                                                              aaGTGGATATGGCGAACCCGGCGAGCTGTCAGCTTTTTGGGCAAAAATTGGATCGCCTGCGCAACCATACCGCCACACAATTGATTAACGGCTTTAATTTTTTACGTTGGCTGGAAAGCGAACCGCAAGATTCGCATaa                                                            <  2:1031005/139‑1 (MQ=255)                                                               aGTGGATATGGCGAACCCGGCGAGCTGTCAGCTTTTTGGGCAAAAATTGGATCGCCTGCGCAACCATACCGCCACACAATTGATTAACGGCTTTAATTTTTTACGTTGGCTGGAAAGCGAACCGCAAGATTCGCATAAc                                                           >  2:1302247/1‑139 (MQ=255)                                                                gTGGATATGGCGAACCCGGCGAGCTGTCAGCTTTTTGGGCAAAAATTGGATCGCCTGCGCAACCATACCGCCACACAATTGATTAACGGCTTTAATTTTTTACGTTGGCTGGAAAGCGAACCGCAAGATTCGCATAACg                                                          <  1:678781/139‑1 (MQ=255)                                                                         gCGAACCCGGCGAGCTGTCAGCTTTTTGGGCAAAAATTGGATCGCCTGCGCAACCATACCGCCACACaa                                                                                                                       <  1:1587800/69‑1 (MQ=255)                                                                         gCGAACCCGGCGAGCTGTCAGCTTTTTGGGCAAAAATTGGATCGCCTGCGCAACCATACCGCCACACaa                                                                                                                       >  2:1587800/1‑69 (MQ=255)                                                                          cGAACCCGGCGAGCTGTCAGCTTTTTGGGCAAAAATTGGATCGCCTGCGCAACCATACCGCCACACAATTGATTAACGGCTTTAATTTTTTACGTTGGCTGGAAAGCGAACCGCAAGATTCGCATAACGAGCATGTCGt                                                <  1:960786/139‑1 (MQ=255)                                                                           gAACCCGGCGAGCTGTCAGCTTTTTGGGCAAAAATTGGATCGCCTGCGCAACCATACCGCCACACAATTGATTAACGGCTTTAATTTTTTACGTTGGCTGGAAAGCGAACCGCAAGATTCGCATAACGAGCATGTCGtt                                               >  2:225144/1‑139 (MQ=255)                                                                           gAACCCGGCGAGCTGTCAGCTTTTTGGGCAAAAATTGGATCGCCTGCGCAACCATACCGCCACACAATTGATTAACGGCTTTAATTTTTTACGTTGGCTGGAAAGCGAACCGCAAGATTCGCATAACGAGCATGTCGtt                                               >  1:187067/1‑139 (MQ=255)                                                                             aCCCGGCGAGCTGTCAGCTTTTTGGGCAAAAATTGGATCGCCTGCGCAACCATACCGCCACACAATTGATTAGCGGCTTTAATTTTTTACGTTGGCTGGAAAGCGAACCGCAAGATTCGCATAACGAGCATGTCGttat                                             >  1:661229/1‑139 (MQ=255)                                                                              cccGGCGAGCTGTCAGCTTTTTGGGCGAAAATTGGATCGCCTGCGCAACCATACCGCCACACAATTGATTAACGGCTTTAATTTTTTACGTTGGCTGGAAAGCGAACCGCAAGATTCGCATAACGAGCATGTCGttatt                                            >  2:126771/1‑139 (MQ=255)                                                                                 ggCGAGCTGTCAGCTTTTTGGGCAAAAATTGGATCGCCTGCGCAACCATACCGCCACACAAt                                                                                                                      >  1:1247209/1‑62 (MQ=255)                                                                                 ggCGAGCTGTCAGCTTTTTGGGCAAAAATTGGATCGCCTGCGCAACCATACCGCCACACAAt                                                                                                                      <  2:1247209/62‑1 (MQ=255)                                                                                  gCGAGCTGTCAGCTTTTTGGGCAAAAATTGGATCGCCTGCGCAACCATACCGCCACACAATTGATTAACGGCTTTAATTTTTTACGTTGGCTGGAAAGCGAACCGCAAGATTCGCATAACGAGCATGTCGTTATTAATg                                        >  1:487192/1‑139 (MQ=255)                                                                                      gCTGTCAGCTTTTTGGGCAAAAATTGGATCGCCTGCGCAACCATACCGCCACACAATTGATTAGCGGCTTTAATTTTTTACGTTGGCTGGAAAGCGAACCGCAAGATTCGCATAACGAGCATGTCGTTATTAATGGGCa                                    <  2:661229/139‑1 (MQ=255)                                                                                            aGCTTTTTGGGCAAAAATTGGATCGCCTGCGCAACCATACCGCCACACAATTGATTAACGGCTTTAATTTTTTACGTTGGCTGGAAAGCg                                                                               >  1:885558/1‑90 (MQ=255)                                                                                            aGCTTTTTGGGCAAAAATTGGATCGCCTGCGCAACCATACCGCCACACAATTGATTAACGGCTTTAATTTTTTACGTTGGCTGGAAAGCg                                                                               <  2:885558/90‑1 (MQ=255)                                                                                             gCTTTTTGGGCAAAAATTGGATCGCCTGCGCAACCATACCGCCACACAATTGATTAACGGCTTTAATTTTTTACGTTGGCTGGAAAGCGAACCGCAAGATTCGCATAACGAGCATGTCGTTATTAATGGGCAGAATTTc                             <  2:1031687/139‑1 (MQ=255)                                                                                             gCTTTTTGGGCAAAAATTGGATCGCCTGCGCAACCATACCGCCACACAATTGATTAACGGCTTTAATTTTTTACGTTGGCTGGAAAGCGAACCGCAAGATTCGCATAACGAGCATGTCGTTATTAATGGGCAGAATTTc                             <  2:563313/139‑1 (MQ=255)                                                                                                 tttGGGCAAAAATTGGATCGCCTGCGCAACCATACCGCCACACAATTGATTAACGGCTTTAATTTTTTACGTTGGCTGGAAAGCGAACCGCAAGATTCGCATAACGAg                                                        <  1:816991/108‑1 (MQ=255)                                                                                                 tttGGGCAAAAATTGGATCGCCTGCGCAACCATACCGCCACACAATTGATTAACGGCTTTAATTTTTTACGTTGGCTGGAAAGCGAACCGCAAGATTCGCATAACGAg                                                        >  2:816991/1‑108 (MQ=255)                                                                                                   tGGGCAAAAATTGGATCGCCTGCGCAACCATACCGCCACACAATTGATTAACGGCTTTAATTTTTTACGTTGGCTGGAAAGCGAACCGCAAGATTCGCATAACGAGCATGTCGTTATTAATGGGCAGAATTTCCTGATg                       <  2:187067/139‑1 (MQ=255)                                                                                                   tGGGCAAAAATTGGATCGCCTGCGCAACCATACCGCCACACAATTGATTAACGGCTTTAATTTTTTACGTTGGCTGGAAAGCGAACCGCAACATTCGCATAACGAGCATGTCGTTATTAATGGGCAGAATTTCCTGATg                       <  1:83149/139‑1 (MQ=255)                                                                                                         aaaaTTGGATCGCCTGCGCAACCATACCGCCACACAATTGATTAACGGCTTTAATTTTTTACGTTGGCTGGAAAGCGAACCGCAAGATTCGCATAACGAGCATGTCGTTATTAATGGGCAGAATTTCCTGATGGAGAtt                 <  1:1302247/139‑1 (MQ=255)                                                                                                            aTTGGATCGCCTGCGCAACCATACCGCCACACAATTGATTAACGGCTTTAATTTTTTACGTTGGCTGGAAAGCGAACCGCAAGATTCGCATAACGAGCATGTCGTTATTAATGGGCAGAATTTCCTGATGGAGATTACg              <  1:225144/139‑1 (MQ=255)                                                                                                            aTTGGATCGCCTGCGCAACCATACCGCCACACAATTGATTAACGGCTTTAATTTTTTACGTTGGCTGGAAAGCGAACCGCAAGATTCGCATAACGAGCATGTCGTTATTAATGGGCAGAATTTCCTGATGGAGATTACg              <  2:487192/139‑1 (MQ=255)                                                                                                                     ccTGCGCAACCATACCGCCACACAATTGATTAACGGCTTTAATTTTTTACGTTGGCTGGAAAGCGAACCGCAAGATTCGCATAACGAGCATGTCGTTATTAATGGGTAGAATTTCCTGATGGAGATTACGCCTGTTTAt     >  1:828244/1‑139 (MQ=255)                                                                                                                      cTGCGCAACCATACCGCCACACAATTGATTAACGGCTTTAATTTTTTACGTTGGCTGGAAAGCGAACCGCAAGATTCGCATAACGAGCATGTCGTTATTAATGGGCAGAATTTCCTGATGGAGATTACGcc            <  2:44683/131‑1 (MQ=255)                                                                                                                      cTGCGCAACCATACCGCCACACAATTGATTAACGGCTTTAATTTTTTACGTTGGCTGGAAAGCGAACCGCAAGATTCGCATAACGAGCATGTCGTTATTAATGGGCAGAATTTCCTGATGGAGATTACGcc            >  1:44683/1‑131 (MQ=255)                                                                                                                        gcgcAACCATACCGCCACACAATTGATTAACGGCTTTAATTTTTTACGTTGGCTGGAAAGCGAACCGCAAGATTCGCATAACGAGCATGTCGTTATTAATGGGCAGAATTTCCTGATGGAGATTACGCCTGTTTATCtt  <  2:342603/139‑1 (MQ=255)                                                                                                                                        |                                                                                                                           GCGTTGAGCGCGTTACTGGAGGCGTTGCCTGAACCTGTGCTCTCTGTCGATATGAAAAGCAAAGTGGATATGGCGAACCCGGCGAGCTGTCAGCTTTTTGGGCAAAAATTGGATCGCCTGCGCAACCATACCGCCGCACAATTGATTAACGGCTTTAATTTTTTACGTTGGCTGGAAAGCGAACCGCAAGATTCGCATAACGAGCATGTCGTTATTAATGGGCAGAATTTCCTGATGGAGATTACGCCTGTTTATCTT  >  NC_000913/1386960‑1387217

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