Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	523,741	(G)6→5	coding (481/4281 nt)	rhsD →	protein RhsD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	523,736	0	G	.	100.0%	162.8 / NA	39	coding (476/4281 nt)	rhsD	protein RhsD
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (14/25);  total (14/25)

ACAACGGCGGGCGGAGCATTCACTTTGAGCCGCTGCTGCCGGGGGAGGCGGTGTACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCTGTGGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCGTGGTGGATACTGGGGTGGTCTGAGCGGGTGCCGGGTGCTGAGGACGTACTGCCAGCGCCGCTGCC  >  NC_000913/523604‑523871                                                                                                                                     |                                                                                                                                        aCAACGGCGGGCGGAGCATTCACTTTGAGCCGCTGCTGCCGGGGGAGGCGGTGTACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCTGTgggggcgc                                                                                                                                  <  1:677811/140‑4 (MQ=255)       gcgggcggAGCATTCACTTTGAGCCGCTGCTGCCGGGGGAGGCGGTGTACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCTGT‑GGGGGCGCTgccgc                                                                                                                           <  2:989104/140‑1 (MQ=255)            cggAGCATTCACTTTGCGCCGCTGGTGCCGGGGGAGGCGGTGTACAGCCGCAGTGAGTCAATGTGGGTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCTGT‑GGGGGCGCTGCCGCCGGat                                                                                                                      <  1:348173/140‑1 (MQ=255)              gAGCATTCACTTTGAGCCGCTGCTGCCGGGGGAGGCGGTGTACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCTGT‑GGGGGCGCTGCCGCCGGatat                                                                                                                    <  1:1200782/140‑1 (MQ=255)                   ttCACTTTGAGCCGCTGCTGCCGGGGGAGGCGGTGTACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCTGT‑GGGGGCGCTGCCGCCGGATATCCGGt                                                                                                               <  1:1301083/140‑1 (MQ=255)                    tCACTTTGAGCCGCTGCTGCCGGGGGAGGCGGTGTACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCTGT‑GGGGGCGCTGCCGCCGGATATCCGGtt                                                                                                              <  1:628953/140‑1 (MQ=255)                         ttGAGCCGCTGCTGCCGGGGGAGGCGGTGTACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCTCGGCTGT‑GGGGGCGCTGCCGCCGGATATCCGGTTAAGcc                                                                                                         <  1:356404/140‑1 (MQ=255)                          tGAGCCGCTGCTGCCGGGGGAGGCGGTGTACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCTGT‑GGGGGCGCTGCCGCCGGATATCCGGTTAAGccc                                                                                                        <  1:519209/140‑1 (MQ=255)                                 ctgctgCCGGGGGAGGCGGTGTACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCTGT‑GGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCtt                                                                                                 >  2:868384/1‑140 (MQ=255)                                        cGGGGGAGGCGGTGTACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCTGT‑GGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTgg                                                                                          <  2:385881/140‑1 (MQ=255)                                            ggAGGCGGTGTACAGCCGCGGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCTGT‑GGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGAc                                                                                      >  2:291877/1‑140 (MQ=255)                                            ggAGGCGGTGTACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCTGT‑GGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGAc                                                                                      <  1:1322743/140‑1 (MQ=255)                                            ggAGGCGGTGTACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCTGT‑GGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCAACTTGACCTGGCGAc                                                                                      >  2:839277/1‑140 (MQ=255)                                               ggCGGTGTACAGCCGCGGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCTGT‑GGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCaa                                                                                   <  1:291877/140‑1 (MQ=255)                                                           ccGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCTGT‑GGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCCACAGCGCACAggg                                                                       >  1:339949/1‑140 (MQ=255)                                                              cAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCTGT‑GGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAgggg                                                                      <  1:86508/138‑1 (MQ=255)                                                              cAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCTGT‑GGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAgggg                                                                      >  2:86508/1‑138 (MQ=255)                                                               aGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCTGT‑GGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCg                                                                   <  2:341027/140‑1 (MQ=255)                                                                 tGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCTGT‑GGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCgtg                                                                 <  1:336559/140‑1 (MQ=255)                                                                   aGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCTGT‑GGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCgtggt                                                               >  2:1144455/1‑140 (MQ=255)                                                                   aGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCTGT‑GGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCgtggt                                                               <  1:1165603/140‑1 (MQ=255)                                                                   aGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCTGT‑GGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCgtggt                                                               <  1:868384/140‑1 (MQ=255)                                                                   aGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCTGT‑GGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCgtggt                                                               <  1:740071/140‑1 (MQ=255)                                                                         tgtgGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCTGT‑GGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGcc                                                                    <  1:755066/129‑1 (MQ=255)                                                                         tgtgGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCTGT‑GGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGcc                                                                    >  2:755066/1‑129 (MQ=255)                                                                           tgGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCTGT‑GGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCGTGGTGGATACTg                                                       >  2:56082/1‑140 (MQ=255)                                                                                 gTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCTGT‑GGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCGTGGTGGATACTGtggtgg                                                 >  2:579648/1‑140 (MQ=255)                                                                                  tGCGCGGTGGTAATGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCTGT‑GGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCGTGGTGGATACTGGggtggt                                                <  1:249454/140‑1 (MQ=255)                                                                                     gcgGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCTGT‑GGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCGTGGTGGATACTGGGGTGGTCTg                                             >  2:1240042/1‑140 (MQ=255)                                                                                      cgGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCTGT‑GGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCGTGGTGGATACTGGGGTGGTCTGa                                            >  2:548423/1‑140 (MQ=255)                                                                                                   gcaCAGCCGGACGGCCATACGCTGGCGCGGCTGT‑GGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCGTGGTGGATACTGGggtggt                                                <  1:202541/123‑1 (MQ=255)                                                                                                   gcaCAGCCGGACGGCCATACGCTGGCGCGGCTGT‑GGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCGTGGTGGATACTGGggtggt                                                >  2:202541/1‑123 (MQ=255)                                                                                                   gcaCAGCCGGACGGCCATACGCTGGCGCGGCTGT‑GGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCGTGGTGGATACTGGGGTGGTCTGAGcgggtgccgggg                               >  2:1031489/1‑139 (MQ=255)                                                                                                     acaGCCGGACGGCCATACGCTGGCGCGGCTGT‑GGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCGTGGTGGATACTGGGGTGGTCTGAGcgggtgcc                                   <  1:265403/134‑1 (MQ=255)                                                                                                     acaGCCGGACGGCCATACGCTGGCGCGGCTGT‑GGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCGTGGTGGATACTGGGGTGGTCTGAGcgggtgcc                                   >  2:265403/1‑134 (MQ=255)                                                                                                         ccGGACGGCCATACGCTGGCGCGGCTGT‑GGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCGTGGTGGATACTGGGGTGGTCTGAGCGGGTGCCGGGTGCTGAg                         <  2:339949/140‑1 (MQ=255)                                                                                                                   aTACGCTGGCGCGGCTGT‑GGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCGTGGTGGATACTGGGGTGGTCTGAGCGGGTGCCGGGTGCTGAGGACGTACTGc               <  1:548423/140‑1 (MQ=255)                                                                                                                            cgcgGCTGT‑GGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCGTGGTGGATACTGGGGTGGTCTGAGCGGGTGCCGGGTGCTGAGGACGTACTGCCAGcgccgc      <  1:56082/140‑1 (MQ=255)                                                                                                                            cgcgGCTGT‑GGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCGTGGTGGATACTGGGGTGGTCTGAGCGGGTGCCGGGTGCTGAGGACGTACTGCCAGcgccgc      <  1:839277/140‑1 (MQ=255)                                                                                                                                gCTGT‑GGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCGTGGTGGATACTGGGGTGGTCTGAGCGGGTGCCGGGTGCTGAGGACGTACTGCCAGCGCCGCTgcc  <  1:1240042/140‑1 (MQ=255)                                                                                                                                     |                                                                                                                                        ACAACGGCGGGCGGAGCATTCACTTTGAGCCGCTGCTGCCGGGGGAGGCGGTGTACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCTGTGGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCGTGGTGGATACTGGGGTGGTCTGAGCGGGTGCCGGGTGCTGAGGACGTACTGCCAGCGCCGCTGCC  >  NC_000913/523604‑523871

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