Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,968,718	T→C	D87G (GAC→GGC)	rppH ←	RNA pyrophosphohydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,968,718	0	T	C	100.0%	190.9 / NA	53	D87G (GAC→GGC)	rppH	RNA pyrophosphohydrolase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (28/25);  total (28/25)

CGGATACCAGTAACTTACCCATCGCCAGCCGTCAAACTCTGGTGTACTGCTGGTTTGCATATTGATTTCTGCATCGCCGCTCACCAGCTGCAAGAGAAACCATTTTTGTTTTTGGCCGATACAAACCGGCTTCGTGTCCCAACGCACCAAACGTTTCGGTAATTTGTAGCGCAACCAGTTACGCGTTGAAGCAAGGATTCGAACGTCTTTGCGGCTTAATCCTACTTCTTCAAACAATTCACGGTACATCGCCTGCTCTGCGGATTCTC  >  NC_000913/2968582‑2968850                                                                                                                                         |                                                                                                                                     cGGATACCAGTAACTTACCCATCGCCAGCCGTCAAACTCTGGTGTACTGCTGGTTTGCATATTGATTTCTGCATCGCCGCTCACCAGCTGCAAGAGAAACCATTTTTGTTTTTGGCCGATACAAACCGGCTTCGTGcccc                                                                                                                                   <  1:362789/140‑1 (MQ=255) cGGATACCAGTAACTTACCCATCGCCAGCCGTCAAACTCTGGTGTACTGCTGGTTTGCATATTGATTTCTGCATCGCCGCTCACCAGCTGCAAGAGAAACCATTTTTGTTTTTGGCCGATACAAACCGGCTTCGTGcccc                                                                                                                                   <  2:385646/140‑1 (MQ=255)   gATACCAGTAACTTACCCATCGCCAGCCGTCAAACTCTGGTGTACTGCTGGTTTGCATATTGATTTCTGCATCGCCGCTCACCAGCTGCAAGAGAAACCATTTTTGTTTTTGGCCGATACAAACCGGCTTCGTGCCCCaa                                                                                                                                 <  2:348301/140‑1 (MQ=255)   gATACCAGTAACTTACCCATCGCCAGCCGTCAAACTCTGGTGTACTGCTGGTTTGCATATTGATTTCTGCATCGCCGCTCACCAGCTGCAAGAGAAACCATTTTTGTTTTTGGCCGATACAAACCGGCTTCGTGCCCCaa                                                                                                                                 <  1:392824/140‑1 (MQ=255)   gATACCAGTAACTTACCCATCGCCAGCCGTCAAACTCTGGTGTACTGCTGGTTTGCATATTGATTTCTGCATCGCCGCTCACCAGCTGCAAGAGAAACCATTTTTGTTTTTGGCCGATACAAACCGGCTTCGTGCCCCaa                                                                                                                                 <  1:423553/140‑1 (MQ=255)     tACCAGTAACTTACCCATCGCCAGCCGTCAAACTCTGGTGTACTGCTGGTTTGCATATTGATTTCTGCATCGCCGCTCACCAGCTGCAAGAGAAACCATTTTTGTTTTTGGCCGATACAAACCGGCTTCGTGCCCCAACg                                                                                                                               >  1:573571/1‑140 (MQ=255)                       cGCCAGCCGTCAAACTCTGGTGTACTGCTGGTTTGCATATTGATTTCTGCATCGCCGCTCACCAGCTGCAAGAGAAACCATTTTTGTTTTTGGCCGATACAAACCGGCTTCGTGCCCCAACGCACCAAACGTTTCGGTaa                                                                                                             >  2:599853/1‑140 (MQ=255)                       cGCCAGCCGTCAAACTCTGGTGTACTGCTGGTTTGCATATTGATTTCTGCATCGCCGCTCACCAGCTGCAAGAGAAACCATTTTTGTTTTTGGCCGATACAAACCGGCTTCGTGCCCCAACGCACCAAACGTTTCGGTaa                                                                                                             >  2:96694/1‑140 (MQ=255)                       cGCCAGCCGTCAAACTCTGGTGTACTGCTGGTTTGCATATTGATTTCTGCATCGCCGCTCACCAGCTGCAAGAGAAACCATTTTTGTTTTTGGCCGATACAAACCGGCTTCGTGCCCCAACGCACCAAACGTTTCGGTaa                                                                                                             >  1:80410/1‑140 (MQ=255)                       cGCCAGCCGTCAAACTCTGGTGTACTGCTGGTTTGCATATTGATTTCTGCATCGCCGCTCACCAGCTGCAAGAGAAACCATTTTTGTTTTTGGCCGATACAAACCGGCTTCGTGCCCCAACGCACCAAACGTTTCGGTaa                                                                                                             >  1:685070/1‑140 (MQ=255)                         ccAGCCGTCAAACTCTGGTGTACTGCTGGTTTGCATATTGATTTCTGCATCGCCGCTCACCAGCTGCAAGAGAAACCATTTTTGTTTTTGGCCGATACAAACCGGCTTCGTGCCCCAACGCACCAAACGTTTCGGTAAtt                                                                                                           >  2:739123/1‑140 (MQ=255)                                tCAAACTCTGGTGTACTGCTGGTTTGCAAATTGATTTCTGCATCGCCGCTCACCAGCTGCAAGAGAAACCATTTTTGTTTTTGGCCGATACAAACCGGCTTCGTGCCCCAACGCACCAAACGTTTCGGTAATTTGTAgcg                                                                                                    <  1:96694/140‑1 (MQ=255)                                  aaaCTCTGGTGTACTGCTGGTTTGCATATTGATTTCTGCATCGCCGCTCACCAGCTGCAAGAGAAACCATTTTTGTTTTTGGCCGATACAAACCGGCTTCGTGCCCCAACGCACCAAACGTTTCGGTAATTTGTAGCGCa                                                                                                  <  2:211022/140‑1 (MQ=255)                                  aaaCTCTGGTGTACTGCTGGTTTGCATATTGATTTCTGCATCGCCGCTCACCAGCTGCAAGAGAAACCATTTTTGTTTTTGGCCGATACAAACCGGCTTCGTGCCCCAACGCACCAAACGTTTCGGTAATTTGTAGCGCa                                                                                                  >  2:252668/1‑140 (MQ=255)                                  aaaCTCTGGTGTACTGCTGGTTTGCATATTGATTTCTGCATCGCCGCTCACCAGCTGCAAGAGAAACCATTTTTGTTTTTGGCCGATACAAACCGGCTTCGTGCCCCAACGCACCAAACGTTTCGGTAATTTGTAGCGCa                                                                                                  <  1:375729/140‑1 (MQ=255)                                          gtgtACTGCTGGTTTGCATATTGATTTCTGCATCGCCGCTCACCAGCTGCAAGAGAAACCATTTTTGTTTTTGGCCGATACAAACCGGCTTCGTGCCCCAACGCACCAAACGTTTCGGTAATTTGTAGCGCAACCAGTTa                                                                                          <  2:74558/140‑1 (MQ=255)                                          gtgtACTGCTGGTTTGCATATTGATTTCTGCATCGCCGCTCACCAGCTGCAAGAGAAACCATTTTTGTTTTTGGCCGATACAAACCGGCTTCGTGCCCCAACGCACCAAACGTTTCGGTAATTTGTAGCGCAACCAGTTa                                                                                          >  2:673501/1‑140 (MQ=255)                                               ctgctgGTTTGCATATTGATTTCTGCATCGCCGCTCACCAGCTGCAAGAGAAACCATTTTTGTTTTTGGCCGATACAAACCGGCTTCGTGCCCCAACGCACCAAACGTTTCGGTAATTTGTAGCGCAACCAGTTACGCGt                                                                                     >  2:854965/1‑140 (MQ=255)                                                tgctgGTTTGCATATTGATTTCTGCATCGCCGCTCACCAGCTGCAAGAGAAACCATTTTTGTTTTTGGCCGATACAAACCGGCTTCGTGCCCCAACGCACCAAACGTTTCGGTAATTTGTAGCGCAACCAGTTACGCGtt                                                                                    <  2:80410/140‑1 (MQ=255)                                                    ggTTTGCATATTGATTTCTGCATCGCCGCTCACCAGCTGCAAGAGAAACCATTTTTGTTTTTGGCCGATACAAACCGGCTTCGTGCCCCAACGCACCAAACGTTTCGGTAATTTGTAGCGCAACCAGTTACGCGt                                                                                     <  1:613891/135‑1 (MQ=255)                                                    ggTTTGCATATTGATTTCTGCATCGCCGCTCACCAGCTGCAAGAGAAACCATTTTTGTTTTTGGCCGATACAAACCGGCTTCGTGCCCCAACGCACCAAACGTTTCGGTAATTTGTAGCGCAACCAGTTACGCGt                                                                                     >  2:613891/1‑135 (MQ=255)                                                          cATATTGATTTCTGCATCGCCGCTCACCAGCTGCAAGAGAAACCATTTTTGTTTTTGGCCGATACAAACCGGCTTCGTGCCCCAACGCACCAAACGTTTCGGTAATTTGTAGCGCAACCAGTTACGCGTTGAAGCa                                                                              <  1:704304/136‑1 (MQ=255)                                                          cATATTGATTTCTGCATCGCCGCTCACCAGCTGCAAGAGAAACCATTTTTGTTTTTGGCCGATACAAACCGGCTTCGTGCCCCAACGCACCAAACGTTTCGGTAATTTGTAGCGCAACCAGTTACGCGTTGAAGCa                                                                              >  2:704304/1‑136 (MQ=255)                                                            tatTGATTTCTGCATCGCCGCTCACCAGCTGCAAGAGAAACCATTTTTGTTTTTGGCCGATACAAACCGGCTTCGTGCCCCAACGCACCAAACGTTTCGGTAATTTGTAGCGCAACCAGTTACGCGTTGAAGCAAGGAtt                                                                        >  1:482917/1‑140 (MQ=255)                                                            tatTGATTTCTGCATCGCCGCTCACCAGCTGCAAGAGAAACCATTTTTGTTTTTGGCCGATACAAACCGGCTTCGTGCCCCAACGCACCAAACGTTTCGGTAATTTGTAGCGCAACCAGTTACGCGTTGAAGCAAGGAtt                                                                        >  1:413898/1‑140 (MQ=255)                                                               tGATTTCTGCATCGCCGCTCACCAGCTGCAAGAGAAACCATTTTTGTTTTTGGCCGATACAAACCGGCTTCGTGCCCCAACGCACCAAACGTTTCGGTAATTTGTAGCGCAACCAGTTACGCGTTGAAGCAAGGATTCGa                                                                     >  2:696067/1‑140 (MQ=255)                                                                    tCTGCATCGCCGCTCACCAGCTGCAAGAGAAACCATTTTTGTTTTTGGCCGATACAAACCGGCTTCGTGCCCCAACGCACCAAACGTTTCGGTAATTTGTAGCGCAACCAGTTGCGCGTTGAAGCAAGGATTCGGACGTc                                                                >  2:858656/1‑140 (MQ=255)                                                                        cATCGCCGCTCACCAGCTGCAAGAGAAACCATTTTTGTTTTTGGCCGATACAAACCGGCTTCGTGCCCCAACGCACCAAACGTTTCGGTAATTTGTAGCGCAACCAGTTACGCGTTGAAGCAAGGATTCGAACGTCTTTg                                                            >  1:741824/1‑140 (MQ=255)                                                                        cATCGCCGCTCACCAGCTGCAAGAGAAACCATTTTTGTTTTTGGCCGATACAAACCGGCTTCGTGCCCCAACGCACCAAACGTTTCGGTAATTTGTAGCGCAACCAGTTACGCGTTGAAGCAAGGATTCGAACGTCTTTg                                                            <  1:673501/140‑1 (MQ=255)                                                                        cATCGCCGCTCACCAGCTGCAAGAGAAACCATTTTTGTTTTTGGCCGATACAAACCGGCTTCGTGCCCCAACGCACCAAACGTTTCGGTAATTTGTAGCGCAACCAGTTACGCGTTGAAGCAAGGATTCGAACGTCTTTg                                                            <  1:599853/140‑1 (MQ=255)                                                                        cATCGCCGCTCACCAGCTGCAAGAGAAACCATTTTTGTTTTTGGCCGATACAAACCGGCTTCGTGCCCCAACGCACCAAACGTTTCGGTAATTTGTAGCGCAACCAGTTACGCGTTGAAGCAAGGATTCGAACGTCTTTg                                                            >  2:415855/1‑140 (MQ=255)                                                                         aTCGCCGCTCACCAGCTGCAAGAGAAACCATTTTTGTTTTTGGCCGATACAAACCGGCTTCGTGCCCCAACGCACCAAACGTTTCGGTAATTTGTAGCGCAACCAGTTACGCGTTGAAGCaa                                                                             <  1:769784/122‑1 (MQ=255)                                                                         aTCGCCGCTCACCAGCTGCAAGAGAAACCATTTTTGTTTTTGGCCGATACAAACCGGCTTCGTGCCCCAACGCACCAAACGTTTCGGTAATTTGTAGCGCAACCAGTTACGCGTTGAAGCaa                                                                             >  2:769784/1‑122 (MQ=255)                                                                             ccgcTCACCAGCTGCAAGAGAAACCATTTTTGTTTTTGGCCGATACAAACCGGCTTCGTGCCCCAATGCACCAAACGTTTCGGTAATTTGTAGCGCAACCAGTTACGCGTTGAAGCAAGGATTCGAACGTCTTTGCGGCt                                                       >  1:926053/1‑140 (MQ=255)                                                                             ccgcTCACCAGCTGCAAGAGAAACCATTTTTGTTTTTGGCCGATACAAACCGGCTTCGTGCCCCAACGCACCAAACGTTTCGGTAATTTGTAGCGCAACCAGTTACGCGTTGAAGCAAGGATTCGAACGTCttt                                                             <  1:25018/134‑1 (MQ=255)                                                                             ccgcTCACCAGCTGCAAGAGAAACCATTTTTGTTTTTGGCCGATACAAACCGGCTTCGTGCCCCAACGCACCAAACGTTTCGGTAATTTGTAGCGCAACCAGTTACGCGTTGAAGCAAGGATTCGAACGTCttt                                                             >  2:25018/1‑134 (MQ=255)                                                                             ccgcTCACCAGCTGCAAGAGAAACCATTTTTGTTTTTGGCCGATACAAACCGGCTTCGTGCCCCAACGCACCAAACGTTTCGGTAATTTGTAGCGCAACCAGTTACGCGTTGAAGCAAGGATTCGAACGTCTTTGCGGCt                                                       <  2:573571/140‑1 (MQ=255)                                                                                cTCACCAGCTGCAAGAGAAACCATTTTTGTTTTTGGCCGATACAAACCGGCTTCGTGCCCCAACGCACCAAACGTTTCGGTAATTTGTAGCGCAACCAGTTACGCGTTGAAGCAAGGATTCGAACGTCTTTGCGGCTTaa                                                    <  1:252668/140‑1 (MQ=255)                                                                                     cAGCTGCAAGAGAAACCATTTTTGTTTTTGGCCGATACAAACCGGCTTCGTGCCCCAACGCACCAAACGTTTCGGTAATTTGTAGCGCAACCAGTTACGCGTTGAAGCAAGGATTCGAACGTCTTTGCGGCTTAATCCTa                                               <  1:858656/140‑1 (MQ=255)                                                                                               agaAACCATTTTTGTTTTTGGCCGATACAAACCGGCTTCGTGCCCCAACGCACCAAACGTTTCGGTAATTTGTAGCGCAACCAGTTACGCGTTGAAGCAAGGATTCGAACGTCTTTGCGGCTTAATCCTACTTCTTCaaa                                     >  1:512376/1‑140 (MQ=255)                                                                                               agaAACCATTTTTGTTTTTGGCCGATACAAACCGGCTTCGTGCCCCAACGCACCAAACGTTTCGGTAATTTGTAGCGCAACCAGTTACGCGTTGAAGCAAGGATTCGAACGTCTTTGCGGCTTAATCCTACTTCTTCaaa                                     <  2:685070/140‑1 (MQ=255)                                                                                                  aaCCATTTTTGTTTTTGGCCGATACAAACCGGCTTCGTGCCCCAACGCACCAAACGTTTCGGTAATTTGTAGCGCAACCAGTTACGCGTTGAAGCAAGGATTCGAACGTCTTTGCGGCTTAATCCTACTTCTTCAAACaa                                  >  2:65012/1‑140 (MQ=255)                                                                                                   aCCATTTTTGTTTTTGGCCGATACAAACCGGCTTCGTGCCCCAACGCACCAAACGTTTCGGTAATTTGTAGCGCAACCAGTTACGCGTTGAAGCAAGGATTCGAACGTCTTTGCGGCTTAATCCTACTTCTTCAAACAAt                                 <  1:854965/140‑1 (MQ=255)                                                                                                    ccATTTTTGTTTTTGGCCGATACAAACCGGCTTCGTGCCCCAACGCACCAAACGTTTCGGTAATTTGTAGCGCAACCAGTTACGCGTTGAAGCAAGGATTCGAACGTCTTTGCGGCTTAATCCTACTTCTTCAAACAAtt                                <  1:696067/140‑1 (MQ=255)                                                                                                       tttttgtttttgGCCGATACAAACCGGCTTCGTGCCCCAACGCACCAAACGTTTCGGTAATTTGTAGCGCAACCAGTTACGCGTTGAAGCAAGGATTCGAACGTCTTTGCGGCTTAATCCTACTTCTTCAAACAATTCAc                             <  1:552907/140‑1 (MQ=255)                                                                                                       tttttgtttttgGCCGATACAAACCGGCTTCGTGCCCCAACGCACCAAACGTTTCGGTAATTTGTAGCGCAACCAGTTACGCGTTGAAGCAAGGATTCGAACGTCTTTGCGGCTTAATCCTACTTCTTCAAACAATTCAc                             <  1:65012/140‑1 (MQ=255)                                                                                                                  gtccGATACAAACCGGCTTCGTGCCCCAACGCACCAAACGTTTCGGTAATTTGTAGCGCAACCAGTTACGCGTTGAAGCAAGGATTCGAACGTCTTTGCGGCTTAATCCTACTTCTTCAAACAATTCACGGTACATCGcc                  >  1:673173/3‑140 (MQ=255)                                                                                                                             aCCGGCTTCGTGCCCCAACGCACCAAACGTTTCGGTAATTTGTAGCGCAACCAGTTACGCGTTGAAGCAAGGATTCGAACGTCTTTGCGGCTTAATCCTACTTCTTCAAACAATTCACGGTACATCGCCTGCTCTGCGGa       <  1:415855/140‑1 (MQ=255)                                                                                                                             aCCGGCTTCGTGCCCCAACGCACCAAACGTTTCGGTAATTTGTAGCGCAACCAGTTACGCGTTGAAGCAAGGATTCGAACGTCTTTGCGGCTTAATCCTACTTCTTCAAACAATTCACGGTACATCGCCTGCTCTGCGGa       >  2:555293/1‑140 (MQ=255)                                                                                                                             aCCGGCTTCGTGCCCCAACGCACCAAACGTTTCGGTAATTTGTAGCGCAACCAGTTACGCGTTGAAGCAAGGATTCGAACGTCTTTGCGGCTTAATCCTACTTCTTCAAACAATTCACGGTACATCGCCTGCTCTGCGGa       <  1:555293/140‑1 (MQ=255)                                                                                                                                  cTTCGTGCCCCAACGCACCAAACGTTTCGGTAATTTGTAGCGCAACCAGTTACGCGTTGAAGCAAGGATTCGAACGTCTTTGCGGCTTAATCCTACTTCTTCAAACAATTCACGGTACATCGCCTGc               <  2:877745/127‑1 (MQ=255)                                                                                                                                  cTTCGTGCCCCAACGCACCAAACGTTTCGGTAATTTGTAGCGCAACCAGTTACGCGTTGAAGCAAGGATTCGAACGTCTTTGCGGCTTAATCCTACTTCTTCAAACAATTCACGGTACATCGCCTGc               >  1:877745/1‑127 (MQ=255)                                                                                                                                  cTTCGTGCCCCAACGCACCAAACGTTTCGGTAATTTGTAGCGCAACCAGTTACGCGTTGAAGCAAGGATTCGAACGTCTTTGCGGCTTAATCCTACTTCTTCAAACAATTCACGGTACATCGCCTGCTCTGCGGATtctc  >  1:793797/1‑140 (MQ=255)                                                                                                                                  cTTCGTGCCCCAACGCACCAAACGTTTCGGTAATTTGTAGCGCAACCAGTTACGCGTTGAAGCAAGGATTCGAACGTCTTTGCGGCTTAATCCTACTTCTTCAAACAATTCACGGTACATCGCCTGCTCTGCGGATtctc  >  2:648830/1‑140 (MQ=255)                                                                                                                                  cTTCGTGCCCCAACGAACCAAACGTTTCGGTAATTTGTAGCGCAACCAGTTACGCGTTGAAGCAAGGATTCGAACGTCTTTGCGGCTTAATCCTACTTCTTCAAACAATTCACGGTACATCGCCTGCTCTGCGGATtctc  >  1:552592/1‑140 (MQ=255)                                                                                                                                         |                                                                                                                                     CGGATACCAGTAACTTACCCATCGCCAGCCGTCAAACTCTGGTGTACTGCTGGTTTGCATATTGATTTCTGCATCGCCGCTCACCAGCTGCAAGAGAAACCATTTTTGTTTTTGGCCGATACAAACCGGCTTCGTGTCCCAACGCACCAAACGTTTCGGTAATTTGTAGCGCAACCAGTTACGCGTTGAAGCAAGGATTCGAACGTCTTTGCGGCTTAATCCTACTTCTTCAAACAATTCACGGTACATCGCCTGCTCTGCGGATTCTC  >  NC_000913/2968582‑2968850

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: ‑