Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	523,512	G→A	K84K (AAG→AAA)	rhsD →	protein RhsD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	523,512	0	G	A	100.0%	76.9 / NA	25	K84K (AAG→AAA)	rhsD	protein RhsD
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (9/16);  total (9/16)

CGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGG‑CAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAGACGCCTGCACCGGTGGGCGTTTTCGGCCCCGGCTGGAAAGCGCCTTCTGATATCCGCTTACAGCTACGTGATGACGGACTGATACTCAACGACAACGGCGGGCGGAGCATTCACTTTGAGCCGCTGCTGCCGGGGGAGGCGGTGTACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCTGTGGGGGGCGCTGCCGCCGGAT  >  NC_000913/523282‑523755                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    cgTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGG‑CAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAAACGCCTGCACCGGTGGGCGt                                                                                                                                                                                                                                 <  1:347247/251‑1 (MQ=255)    cAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCACGGGTCGGCAGGAGTAAGAAATGGCGCGCCCACCGGCGTGGCGCGCTCGGGGGGTCCGGGGGGGGAGTCTTCGGGGCACCCCGG‑AAAACCGCTGCCGGGGGGGAAGGGGCTGCCCCGCGGGACGGGCCTTGCGCTGCCCGGCCCGCTGCCGTTCCTTCTCTCCCCCCCCCACAGCAGCTACCCGACGACAACGCCTGCCCCGGGGGGCGtttt                                                                                                                                                                                                                              >  1:344965/1‑251 (MQ=255)             atGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGG‑CAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAAACGCCTGCACCGGTGGGCGTTTTCGGCCCCg                                                                                                                                                                                                                      <  2:201190/250‑1 (MQ=255)                                         ggggTTCGGCAGGTGTTAGAATTGGAGCGCCCCCCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGG‑CACCCCTGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGAGAGACGGAGCTTGCGATGCCCGGCCCGCTGCCGGTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAAACGCCTGCACCGGTGGGCGTTTTCGGCCCCGGCTGGAAAGCGCCTTCTGATATCCGCTTa                                                                                                                                                                                         <  1:362536/250‑1 (MQ=255)                                                  cAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGG‑CAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAAACGCCTGCACCGGTGGGCGTTTTCGGCCCCGGCTGGAAAGCGCCTTCTGATATACGCTTACAGCTACGt                                                                                                                                                                                >  1:250031/1‑251 (MQ=255)                                                  cAGGAGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGG‑CAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCCGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCCACCGGACGCAAACGCCTGCACCGGGGGGGGTTTTCGGCCCCGGCTGGAAAGCGCCTTCTGATATCCGCTTACAGCtcagc                                                                                                                                                                                >  2:239201/1‑247 (MQ=38)                                                       gtAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTGCGGG‑CAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCGTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAAACGCCTGCACCGGTGGGCGTTTGCGGCCCCGGCTTGAAAGCGCCTTCTGATAGCCGCTTACAGCTACGtgatga                                                                                                                                                                           <  2:347776/251‑1 (MQ=255)                                                            aaTTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGG‑CAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAAACGCCTGCACCGGTGGGCGGTTTCGGCCCCGGCCGGAAAGCGCCTTCTGATATCCGCTTACAGCTACGTGACGACGGAc                                                                                                                                                                      >  2:273170/1‑251 (MQ=255)                                                                ggCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGG‑CAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAAACGCCTGCACCGGTGGGCGTTTTCGGCCCCGGCTGGAAAGCGCCTTCTGATATCCGCTTACAGCTACGTGATGACGGACTGAt                                                                                                                                                                  <  2:61217/251‑1 (MQ=255)                                                                  cgcgcCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGG‑CAACCCCGAAAATCCGCTGCTGGGGGCGAAGGTGCTGCACGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGACGTTCATTCTCTCCCGCACCTACAGCCGCTACCGGACGAAAACGCCTGCACCGGAGGGCGTTTTCGGCCCCGGCTTGAAAACGCCTTCTGCTATCCGCTTACAGCTACGGGATGACGGACTGATAc                                                                                                                                                                >  1:216189/1‑251 (MQ=255)                                                                           cGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGG‑CAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAAACGCCTGCACCGGTGGGCGTTTTCGGCCCCGGCTGGAAAGCGCCTTCTGATATCCGCTTACAGCTACGTGATGACGGACTGATACTCAACGAc                                                                                                                                                        >  1:141575/1‑250 (MQ=255)                                                                            ggcgtggcgtGCTCGGTGTGTCCGGGCGGGCTGAGTTTGGG‑CAACCCGGAAAATACGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACATTGCGCTGCCCAGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAAACGCCTGCACCGGAGGGCGTTTTCGGCCCCGGCTGGAAAGCGCCTTCTGATATCCGCTTACAGCTACGTGATGACGGACTGATACTCAACGACa                                                                                                                                                       >  1:342353/1‑250 (MQ=38)                                                                                              tgtCCGGGCGGGATGACTTCGGG‑CAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAAACGCCTGCACCGGTGGGCGTTTTCGGCCCCGGCTGGAAAGCGCCTTCTGATATCCGCTTACAGCTACGTGATGACGGACTGATACTCAACGACAACGgcgggcggggggctcc                                                                                                                                    >  2:46589/1‑243 (MQ=255)                                                                                               gtCCGGGCGGGATGACGTCGGG‑CAGCCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCGGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAAACGCCTGCACCGGTGGGCGTTTTCGGCCCCGGCTGGAAAGCGCCTTCTGATAGCCGCTTACAGCTACGTGATGACGGACTGATACTCAACGACAACGGCGGGCGGCGCATTCAc                                                                                                                                   <  2:200998/251‑1 (MQ=255)                                                                                                                tCGGG‑CAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAAACGCCTGCACCGGTGGGCGTTTTCGGCCCCGGCTGGAAAGCGCCTTCTGATATCCGCTTACAGCTACGTGATGACGGACTGATACTCAACGACAACGGCGGGCGGGGCAATCAccttttgcccccgctgcc                                                                                                                  >  1:340776/1‑234 (MQ=38)                                                                                                                   gg‑CAACCCGTGAAATCCGCTGCTGGGGGCGGAGGTGCTGCCCGGCGAGACGGACCTTGCGCGGCCCGGCCCGCGTCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAAACGCCTGCACCGGTGGGCGTTTTCGGCCCCGGCTGGAAAGCGCCTTCTGATATCCGCTTACAGCTACGTGATGACGGACTGATACTCAACGACAACGGCGGGCGGAGCATTCACTTTGAGCCGCTGCTGCCggg                                                                                                               <  2:321200/251‑1 (MQ=255)                                                                                                                     gcaccccGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAAACGCCTGCACCGGTGGGCGTTTTCGGCCCCGGCTGGAAAGCGCCTTCTGATATCCGCTTACAGCTACGTGATGACGGACTGATACTCAACGACAACGGCGGGCGGAGCATTCACTTTGAGCCGCTGCTGCCgggg                                                                                                              <  2:198999/247‑1 (MQ=255)                                                                                                                             gTAAAGCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAAACGCCTGCACCGGTGGGCGTTTTCGGCCCCGGCTGGAAAGCGCCTTCTGATATCCGCTTACAGCTACGTGATGACGGACTGATACTCAACGACAACGGCGGGCGGAGCATTCACTTTGAGCCGCTGCTGCCGGGGGAGGCGgt                                                                                                      <  1:147664/251‑1 (MQ=255)                                                                                                                                                ggAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAAACGCCTGCACCGGTGGGCGTTTTCGGCCCCGGCTGGAAAGCGCCTTCTGATATCCGCTTACAGCTACGTGATGACGGACTGATACTCAACGACAACGGCGGGCGGAGCATTCACTTTGAGCCGCTGCTGCCGGGGGAGGCGGTGTACAGCCGCAGTGAGTCa                                                                                   <  2:357411/250‑1 (MQ=255)                                                                                                                                                                      ggACCTGGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAAACGCCTGCACCGGTGGGCGTTTGCGGCCCCGGCTGGAAAGCGCCTTCTGATATCCGCTTACAGCTACGTGATGACGGACTGATACTCAACGACAACGGCGGGCGGAGCATTCACTTTGAGCCGCTGCTGCCGGGGGAGGCGGTGTACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTa                                                             <  2:250031/251‑1 (MQ=255)                                                                                                                                                                            tGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAAACGCCTGCACCGGTGGGCGTTTTCGGCCCCGGCTGGAAAGCGCCTTCTGATATCCGCTTACAGCTACGTGATGACGGACTGATACTCAACGACAACGGCGGGCGGAGCATTCACTTTGAGCCGCTGCTGCCGGGGGAGGCGGTGTACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGgca                                                        <  2:112732/250‑1 (MQ=255)                                                                                                                                                                                              gCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAAACGCCTGCACCGGTGGGCGTTTTCGGCCCCGGCTGGAAAGCGCCTTCTGATATCCGCTTACAGCTACGTGATGACGGACTGATACTCAACGACAACGGCGGGCGGAGCATTCACTTTGAGCCGCTGCTGCCGGGGGAGGCGGTGTACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATa                                     <  1:46589/251‑1 (MQ=255)                                                                                                                                                                                              gCCGGTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAAACGCCTGCACCGGTGGGCGTTTTCGGCGCCGTCGGGAAAGCGCCTGCTGATATCCGCGTACAGCTACGTGATGACGGACGGATGCTGAACGACAACGGCGGGCGGGGTATTCACGTTGGGCTGCTGCTGCGGGGGTAGGCGGGGTACGGCGGCAGTGAGCCGATGTTACTGGTGGGCGGTGGTAAGGCAGCACAGCCGGACGACCATa                                     <  1:174919/251‑1 (MQ=255)                                                                                                                                                                                                 gTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAAACGCCTGCACCGGTGGGCGTTTTCGGCCCCGGCTGGAAAGCGCCTTCTGATATCCGCTTACAGCTACGTGATGACGGACTGATACTCAACGACAACGGCGGGCGGAGCATTCACTTTGAGCCGCTGCTGCCGGGGGAGGCGGTGTACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGc                                  <  2:294323/251‑1 (MQ=255)                                                                                                                                                                                                                                 ggACGAAAACGCCTGCACCGGTGGGCGTTTTCGGCCCCGGCTGGAAAGCGCCTTCTGATATCCGCTTACAGCTACGTGATGACGGACTGATACTCAACGACAACGGCGGGCGGAGCATTCACTTTGAGCCGCTGCTGCCGGGGGAGGCGGTGTACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCTGTGGGGGGCGCTTCCGCCGGat  <  1:273170/251‑1 (MQ=255)                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    CGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGG‑CAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAGACGCCTGCACCGGTGGGCGTTTTCGGCCCCGGCTGGAAAGCGCCTTCTGATATCCGCTTACAGCTACGTGATGACGGACTGATACTCAACGACAACGGCGGGCGGAGCATTCACTTTGAGCCGCTGCTGCCGGGGGAGGCGGTGTACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCTGTGGGGGGCGCTGCCGCCGGAT  >  NC_000913/523282‑523755

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑