Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,720,436	A→G	I130T (ATC→ACC)	slyA ←	DNA‑binding transcriptional dual regulator SlyA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,720,436	0	A	G	100.0%	110.1 / NA	31	I130T (ATC→ACC)	slyA	DNA‑binding transcriptional dual regulator SlyA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (16/15);  total (16/15)

TACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTTAATAACAGCTTCCATTTCGCTGATCAGCGGCTCTGCCTTTTCCGTCAGTTTAATACGTTTAGC  >  NC_000913/1720300‑1720566                                                                                                                                         |                                                                                                                                   tACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGGTGa                                                                                                                                 <  2:728687/140‑1 (MQ=255)      gCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGGTGAGCGTa                                                                                                                            >  2:381504/1‑140 (MQ=255)      gCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGGTGAGCGTa                                                                                                                            >  2:765370/1‑140 (MQ=255)         aGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGGTGAGTGTAATc                                                                                                                         <  1:914390/140‑1 (MQ=255)            aaTGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGGTGAGCGTAATCAGt                                                                                                                      <  2:781213/140‑1 (MQ=255)              tGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGGTGAGCGTAATCAGTTg                                                                                                                    >  1:502681/1‑140 (MQ=255)                    tCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGGTGAGCGTAATCAGTTGCTCCAg                                                                                                              <  1:490155/140‑1 (MQ=255)                       ggTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGGTGAGCGTAATCAGTTGCTCCAGTTc                                                                                                           <  2:874606/140‑1 (MQ=255)                                    ccGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGGTGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATg                                                                                              <  2:210996/140‑1 (MQ=255)                                       cgcTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGGTGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCa                                                                                           <  1:381504/140‑1 (MQ=255)                                                aGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGGTGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATAtt                                                                                  <  2:1061181/140‑1 (MQ=255)                                                      cGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGGTGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCcgcg                                                                            >  1:443285/1‑140 (MQ=255)                                                      cGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGGTGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCcgcg                                                                            <  1:120945/140‑1 (MQ=255)                                                         gtgGTCAGGTTACTGCCCCCACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGGTGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCgg                                                                         <  2:271294/140‑1 (MQ=255)                                                            gTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGGTGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGtt                                                                      <  2:443285/140‑1 (MQ=255)                                                            gTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGGTGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGtt                                                                      >  2:656162/1‑140 (MQ=255)                                                            gTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGGTGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGtt                                                                      <  1:936958/140‑1 (MQ=255)                                                                      tGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGGTGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTtaataa                                                            >  2:924990/1‑140 (MQ=255)                                                                       gACCACACGCCCCCTTCATTTCACCCTCTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGGTGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTTAATAAc                                                           >  1:470295/1‑140 (MQ=255)                                                                               gCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGGTGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTTAATAACAGCTTCCa                                                   <  1:715315/140‑1 (MQ=255)                                                                                  cccTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGGTGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTTAATAACAGCTTCCAttt                                                >  1:161258/1‑140 (MQ=255)                                                                                             aCCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGGTGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTTAATAACAGCTTCCATTTCGCTGATCAGc                                     >  2:995215/1‑140 (MQ=255)                                                                                                  ttGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGGTGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTTAATAACAGCTTCCATTTCGCTGATCAGCGGctc                                >  1:962733/1‑140 (MQ=255)                                                                                                    ggCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGGTGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTTAATAACAGCTTCCATTTCGCTGATCAGCGGCTCTg                              <  1:924990/140‑1 (MQ=255)                                                                                                    ggCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGGTGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTTAATAACAGCTTCCATTTCGCTGATCAGCGGCTCTg                              >  1:14990/1‑140 (MQ=255)                                                                                                      ccTGTAACTCAATGATATTATGCTCAAGTTTTGCGGTGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTTAATAACAGCTTCCATTTCGCTGATCAGCGGCTCTGcc                            >  1:381063/1‑140 (MQ=255)                                                                                                             cTCAATGATATTATGCTCAAGTTTTGCGGTGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTTAATAACAGCTTCCATTTCGCTGATCAGCGGCTCTGCCTTTTCCg                     >  2:364221/1‑140 (MQ=255)                                                                                                                        tatGCTCAAGTTTTGCGGTGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTTAATAACAGCTTCCATTTCGCTGATCAGCGGCTCTGCCTTTTCCGTCAGTTTAATa          >  1:240982/1‑140 (MQ=255)                                                                                                                            cTCAAGTTTTGCGGTGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTTAATAACAGCTTCCATTTCGCTGATCAGCGGCTCTGCCTTTTCCGTCAGTTTAATACGtt      <  1:656162/140‑1 (MQ=255)                                                                                                                            cTCAAGTTTTGCGGTGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTTAATAACAGCTTCCATTTCGCTGATCAGCGGCTCTGCCTTTTCCGTCAGTTTAATACGtt      >  1:477300/1‑140 (MQ=255)                                                                                                                                aGTTTTGCGGTGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTTAATAACAGCTTCCATTTCGCTGATCAGCGGCTCTGCCTTTTCCGTCAGTTTAATACGTTTAGc  >  1:665421/1‑140 (MQ=255)                                                                                                                                         |                                                                                                                                   TACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTTAATAACAGCTTCCATTTCGCTGATCAGCGGCTCTGCCTTTTCCGTCAGTTTAATACGTTTAGC  >  NC_000913/1720300‑1720566

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