Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,312,151	G→C	Q281E (CAA→GAA)	alsR ←	DNA‑binding transcriptional repressor AlsR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,312,151	0	G	C	100.0%	113.6 / NA	31	Q281E (CAA→GAA)	alsR	DNA‑binding transcriptional repressor AlsR
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (17/14);  total (17/14)

GCCTTTTTTCATATCTACCCAAAATGGGTTGGAGAGGGTTTTCAATACGACAGCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCG‑CCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCA‑GCAAAGTTAGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGG‑CGTTTCCGGGGCTGGTGAGCAAATAATATAATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGG  >  NC_000913/4311908‑4312390                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                  cgctccagttaatcaaatcacgcatacgcgctcgcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCG‑CCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCatattaa                                                                                                                                                                                                                                           <  2:101179‑M2/102‑1 (MQ=255)                   cacgcatacgcgctctcgactacagtatgcatcttttatgccacattttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgtGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCG‑CCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGt                                                                                                                                                                                                                         >  2:152189‑M2/132‑251 (MQ=255)                        atacgcgctctcgactacagtatgcatcttttatgccacattttatgtggagtcattccctgatattactggctctatttagtcaaaactctgacgaaaaacaggcggtcgtcagtttttacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACTTTTTTCTTATTGCAGCGCG‑GGTGGTGAGAAGAAATCAACAATTGCGCTGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGc                                                                                                                                                                                                                    <  2:516995‑M2/125‑1 (MQ=255)                                                               attttatgtggggtcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCG‑CCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCA‑Gcaaggttggt                                                                                                                                                                            >  2:537251‑M2/88‑248 (MQ=255)                                                                ttttatgtggggtccttccctgatattacgggcactatttattcacaactctgccgaaaaacaggctgtcgtcagttttgacgtgaGGTAATGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGAG‑CCTGGTGAGAAGAAATCAACCATTGCGCCGGTTTTTTCCATATTAAGATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCA‑GCAAggttggtt                                                                                                                                                                           >  2:166302‑M2/87‑251 (MQ=255)                                                                            tcattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCG‑CCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCA‑GCAAGGTTGGTTGTAATATTCTTg                                                                                                                                                               >  2:171240‑M2/75‑251 (MQ=255)                                                                             cattccctgatattacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCG‑CCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCA‑GCAAGGTTGGTTGTAATATTCTTGc                                                                                                                                                              <  1:89801‑M2/178‑1 (MQ=255)                                                                                         ttacgggcactatttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCG‑CCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCA‑GCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACg                                                                                                                                                  >  2:111305‑M2/62‑251 (MQ=255)                                                                                                    atttattcaaaactctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCG‑CCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCA‑GCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACg                                                                                                                                       <  2:295170‑M2/201‑1 (MQ=255)                                                                                                                ctctgacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCG‑CCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCA‑GCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGG‑CGTTTCCgggg                                                                                                                          >  1:285417‑M2/39‑251 (MQ=255)                                                                                                                    gacgaaaaacaggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTGCTTATTTCAGCGCG‑CCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCA‑GCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGG‑CGTTTCCGGGGCTg                                                                                                                       <  1:430851‑M2/216‑1 (MQ=255)                                                                                                                           aacacgctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCG‑CCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATCTTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCA‑GCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGG‑CGTTTCCGGGGCTGGTGAgct                                                                                                                >  2:544250‑M2/28‑249 (MQ=255)                                                                                                                               ggctgtcgtcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCG‑CCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCA‑GCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGG‑CGTTTCCGGGGCTGGTGAGCAaataa                                                                                                           <  2:6429‑M2/228‑1 (MQ=255)                                                                                                                                   gtcgtcagttttgacgtgaGGTAAGGGTACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCG‑CCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCA‑GCAAGGTTGGTTGTAATATTCTTGCCGATGCATTACGACCTAATAACGG‑CGTTTCCGGTGCTGGTGAGCAaataatatg                                                                                                       <  2:549969‑M2/232‑2 (MQ=255)                                                                                                                                       tcagttttgacgtgaGGTAAGGGGACA‑GTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCG‑CCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCA‑GCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGG‑CGTTTCCGGGGCTGGTGAGCAAATAATATGATCg                                                                                                   >  1:185317‑M2/16‑251 (MQ=255)                                                                                                                                                           ggggACAG‑TAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCG‑CCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCA‑GCAAGGTTGGTTGTAatat                                                                                                                                                                    >  2:90275/1‑167 (MQ=255)                                                                                                                                                           ggggACAG‑TAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCG‑CCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCA‑GCAAGGTTGGTTGTAatat                                                                                                                                                                    <  1:90275/167‑1 (MQ=255)                                                                                                                                                            gggACAG‑TAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCG‑CCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCA‑GCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGG‑CGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGa                                                                               >  1:381647/1‑250 (MQ=255)                                                                                                                                                                   gTAACTCCGCGATTGCGGAATTACCTTTTTCTTATTTCAGCGCG‑CCTGGTGAGAAGAAATCAACAATTGCGACGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCA‑GCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGG‑CGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAg                                                                        <  2:438450/251‑1 (MQ=255)                                                                                                                                                                   gTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCG‑CCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCA‑GCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGG‑CGTTTCCGGGGCTGGt                                                                                                                     >  2:311484/1‑206 (MQ=255)                                                                                                                                                                   gTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCG‑CCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCA‑GCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGG‑CGTTTCCGGGGCTGGt                                                                                                                     <  1:311484/206‑1 (MQ=255)                                                                                                                                                                           gcgATGGCGGAATTTCCTTTTTCTTATTGCAGCGCG‑CCAGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCAGGGCGACAGAGACGAAAAAAGCGTCCA‑GCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGG‑CGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTGCGCTGTCGGTGAATGGTAGCTATGGGt                                                                <  2:204691/251‑1 (MQ=255)                                                                                                                                                                                          ccTTTTTCTTATTTTAGCGCGGCCTG‑TGAGGAGGAATCCACAATTTCGCCGGTTTTTTCCATATTAATAATAGCTTGTTCAATGTTGAGCTGGGGGACAGGAACAAAAAAAGCGTCCA‑GCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGG‑CGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTAGGGGTTATACAAATAATCt                                                  <  2:518494/250‑1 (MQ=255)                                                                                                                                                                                                   tATTTCAGCGCG‑CCTGGTGAGAAGAAATCAACAATTGCGCTGGTTTTTTCCATATTACTATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAACAAAGCGTCCA‑GCAAGGTTGGTTGTAATATTCTTGACGAGGCATTAAGACCTAATAACGG‑CGTTTCCGGGGCTGGTGAGCAAATAATATGCTCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGtt                                        >  2:83762/1‑251 (MQ=255)                                                                                                                                                                                                         aGCGCG‑CCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCA‑GCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGG‑CGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCttttt                                  <  2:12206/251‑1 (MQ=255)                                                                                                                                                                                                           cgcg‑cCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCA‑GCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGG‑CGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTAAGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTg                                >  2:52043/1‑251 (MQ=255)                                                                                                                                                                                                                           gaaATCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCA‑GCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGG‑CGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGAAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCg                 >  1:259513/1‑251 (MQ=255)                                                                                                                                                                                                                               tCAACAATTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCA‑GCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGG‑CGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCttt             >  2:388671/1‑251 (MQ=255)                                                                                                                                                                                                                                 aacaaTTGCGCCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCA‑GCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGG‑CGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCtttt            >  1:217614/1‑250 (MQ=255)                                                                                                                                                                                                                                    aatggggcCGGTTTTTTCCATATAAATATTAGCTGGTTCAAGGTTGAGCTGGGCGAAAGAGACAAAAAAAGC‑TCCAAGAAAGGTTGGTTGTAATATCCTGGCGGAGGAATTCCGCCCTAA‑AACGGGCGTTTCCGGGGCTGGTGGGCAAAAAATATGACGGGCCGGTTCGGCTACCGTTGAATGGTAGCAATGGGT‑ATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATc       <  2:185317/245‑1 (MQ=255)                                                                                                                                                                                                                                          gcCGGTTTTTTCCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCA‑GCAAGGTTGGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGG‑CGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTgg  >  1:293748/1‑251 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                  GCCTTTTTTCATATCTACCCAAAATGGGTTGGAGAGGGTTTTCAATACGACAGCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCG‑CCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCA‑GCAAAGTTAGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGG‑CGTTTCCGGGGCTGGTGAGCAAATAATATAATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGG  >  NC_000913/4311908‑4312390

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