Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	523,701	A→C	A147A (GCA→GCC)	rhsD →	protein RhsD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	523,701	0	A	C	100.0%	133.3 / NA	36	A147A (GCA→GCC)	rhsD	protein RhsD
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (15/21);  total (15/21)

TACAGCTACGTGATGACGGACTGATACTCAACGACAACGGCGGGCGGAGCATTCACTTTGAGCCGCTGCTGCCGGGGGAGGCGGTGTACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCTGTGGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCGTGGTGGATACTGGGGTGGTCTGAG  >  NC_000913/523571‑523830                                                                                                                                   |                                                                                                                                  tACAGCTACGTGATGACGGACTGATACTCAACGACAACGGCGGGCGGAGCATTCACTTTGAGCCGCTGCTGCCGGGGGAGGCGGTGTACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCCGCACAGCCg                                                                                                                          <  1:443226/140‑1 (MQ=255)   cAGCTACGTGATGACGGACTGATACTCAACGACAACGGCGGGCGGAGCATTCACTTTGAGCCGCTGCTGCCGGGGGAGGCGGTGTACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCCGCACAGCCGGa                                                                                                                        <  1:448882/140‑1 (MQ=255)   cAGCTACGTGATGACGGACTGATACTCAACGACAACGGCGGGCGGAGCATTCACTTTGAGCCGCTGCTGCCGGGGGAGGCGGTGTACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCCGCACAGCCGGa                                                                                                                        <  1:636372/140‑1 (MQ=255)                       gATACTCAACGACAACGGCGGGCGGAGCATTCACTTTGAGCCGCTGCTGCCGGGGGAGGCGGTGTACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCCGCACAGCCGGACGGCCATACGCTGGCGCGGc                                                                                                    <  1:741817/140‑1 (MQ=255)                                  aCAACGGCGGGCGGAGCATTCACTTTGAGCCGCTGCTGCCGGGGGAGGCGGTGTACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCCGCACAGCCGGACGGCCATACGCTGGCGCGGCTGTGGGGGgcg                                                                                         <  1:641997/140‑1 (MQ=255)                                         cgggcggAGCATTCACTTTGAGCCGCTGCTGCCGGGGGAGGCGGTGTACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCCGCACAGCCGGACGGCCATACGCTGGCGCGGCTGTGGGGGGCGCTgccgc                                                                                  <  1:900926/140‑1 (MQ=255)                                                aGCATTCACTTTGAGCCGCTGCTGCCGGGGGAGGCGGTGTACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCCGCACAGCCGGACGGCCATACGCTGGCGCGGCTGTGGGGGGCGCTGCCGCCGGatat                                                                           >  2:286428/1‑140 (MQ=255)                                                           tGAGCCGCTGCTGCCGGGGGAGGCGGTGTACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCCGCACAGCCGGACGGCCATACGCTGGCGCGGCTGTGGGGGGCGCTGCCGCCGGATATCCGGTTAAGcc                                                                <  1:286428/140‑1 (MQ=255)                                                           tGAGCCGCTGCTGCCGGGGGAGGCGGTGTACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCCGCACAGCCGGACGGCCATACGCTGGCGCGGCTGTGGGGGGCGCTGCCGCCGGATATCCGGTTAAGcc                                                                <  1:480032/140‑1 (MQ=255)                                                            gAGCCGCTGCTGCCGGGGGAGGCGGTGTACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCCGCACAGCCGGACGGCCATACGCTGGCGCGGCTGTGGGGGGCGCTGCCGCCGGATATCCGGTTAAGccc                                                               <  1:606504/140‑1 (MQ=255)                                                               ccGCTGCTGCCGGGGGAGGCGGTGTACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCCGCACAGCCGGACGGCCATACGCTGGCGCGGCTGTGGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCa                                                            >  2:851727/1‑140 (MQ=255)                                                                              gAGGCGGTGTACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCCGCACAGCCGGACGGCCATACGCTGGCGCGGCTGTGGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGAc                                             <  1:178398/140‑1 (MQ=255)                                                                                ggCGGTGTACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCCGCACAGCCGGACGGCCATACGCTGGCGCGGCTGTGGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCa                                           >  2:122253/1‑140 (MQ=255)                                                                                 gCGGTGTACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCCGCACAGCCGGACGGCCATACGCTGGCGCGGCTGTGGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCaa                                          <  2:498620/140‑1 (MQ=255)                                                                                 gCGGTGTACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCCGCACAGCCGGACGGCCATACGCTGGCGCGGCTGTGGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCaa                                          >  2:227120/1‑140 (MQ=255)                                                                                       tACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCCGCACAGCCGGACGGCCATACGCTGGCGCGGCTGTGGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTAcc                                                    <  1:488096/124‑1 (MQ=255)                                                                                       tACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCCGCACAGCCGGACGGCCATACGCTGGCGCGGCTGTGGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTAcc                                                    >  2:488096/1‑124 (MQ=255)                                                                                              gCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCCGCACAGCCGGACGGCCATACGCTGGCGCGGCTGTGGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAgggg                             <  1:122253/140‑1 (MQ=255)                                                                                               cAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCCGCACAGCCGGACGGCCATACGCTGGCGCGGCTGTGGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGc                            >  2:953349/1‑140 (MQ=255)                                                                                               cAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCCGCACAGCCGGACGGCCATACGCTGGCGCGGCTGTGGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGc                            >  2:156332/1‑140 (MQ=255)                                                                                               cAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCCGCACAGCCGGACGGCCATACGCTGGCGCGGCTGTGGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGc                            >  1:780186/1‑140 (MQ=255)                                                                                                aGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCCGCACAGCCGGACGGCCATACGCTGGCGCGGCTGTGGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGcc                           <  1:791547/140‑1 (MQ=255)                                                                                                 gTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCCGCACAGCCGGACGGCCATACGCTGGCGCGGCTGTGGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCg                          >  2:147009/1‑140 (MQ=255)                                                                                                 gTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCCGCACAGCCGGACGGCCATACGCTGGCGCGGCTGTGGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCg                          <  1:76316/140‑1 (MQ=255)                                                                                                 gTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCCGCACAGCCGGACGGCCATACGCTGGCGCGGCTGTGGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCCGGCGCCCCCCAGCGCACAGGGGccc                          >  2:600404/1‑139 (MQ=255)                                                                                                     gTCAATGTGGCTGGTGCGCGGTGGTAAGGCCGCACAGCCGGACGGCCATACGCTGGCGCGGCTGTGGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCgtggt                      <  1:573753/140‑1 (MQ=255)                                                                                                     gTCAATGTGGCTGGTGCGCGGTGGTAAGGCCGCACAGCCGGACGGCCATACGCTGGCGCGGCTGTGGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCgtggt                      <  1:227120/140‑1 (MQ=255)                                                                                                     gTCAATGTGGCTGGTGCGCGGTGGTAAGGCCGCACAGCCGGACGGCCATACGCTGGCGCGGCTGTGGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCgtggt                      <  1:662598/140‑1 (MQ=255)                                                                                                      tCAATGTGGCTGGTGCGCGGTGGTAAGGCCGCACAGCCGGACGGCCATACGCTGGCGCGGCTGTGGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCgtggtg                     >  1:440098/1‑140 (MQ=255)                                                                                                          tgtgGCTGGTGCGCGGTGGTAAGGCCGCACAGCCGGACGGCCATACGCTGGCGCGGCTGTGGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCGTGGTGGATa                 <  1:953349/140‑1 (MQ=255)                                                                                                              gCTGGTGCGCGGTGGTAAGGCCGCACAGCCGGACGGCCATACGCTGGCGCGGCTGTGGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCGTGGTGGATACTgg             >  1:799917/1‑140 (MQ=255)                                                                                                               cTGGTGCGCGGTGGTAAGGCCGCACAGCCGGACGGCCATACGCTGGCGCGGCTGTGGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCGTGGTGGATACTggg            <  2:296429/140‑1 (MQ=255)                                                                                                               cTGGTGCGCGGTGGTAAGGCCGCACAGCCGGACGGCCATACGCTGGCGCGGCTGTGGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCGTGGTGGATACTggg            >  2:169518/1‑140 (MQ=255)                                                                                                                tGGTGCGCGGTGGTAAGGCCGCACAGCCGGACGGCCATACGCTGGCGCGGCTGTGGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCGTGGTGGATACTgggg           >  1:1029/1‑140 (MQ=255)                                                                                                                    gcgcGGTGGTAAGGCCGCAAAGCCGGACGGCCATACGCTGGCGCGGCTGTGGGGGGCGCTGCCGCCCGATATCCGGTTAAGCCCGCATCTTTTCATGGCGGCCAACAGCGCCCATGGGCCGTGGGGGATACAGGggtggt       >  2:810162/1‑140 (MQ=255)                                                                                                                         gtggtAAGGCCGCACAGCCGGACGGCCATACGCTGGCGCGGCTGTGGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCGTGGTGGATACTGGGGTGGTCTGAg  <  2:440098/140‑1 (MQ=255)                                                                                                                         ggtggaaGGCCGCACAGCCGGACGGCCATACGCTGGCGCGGCTGTGGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCGTGGTGGATACTGGGGTGGTCTGAg  <  1:394153/135‑1 (MQ=255)                                                                                                                                   |                                                                                                                                  TACAGCTACGTGATGACGGACTGATACTCAACGACAACGGCGGGCGGAGCATTCACTTTGAGCCGCTGCTGCCGGGGGAGGCGGTGTACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCTGTGGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCGTGGTGGATACTGGGGTGGTCTGAG  >  NC_000913/523571‑523830

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