Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,179,175	T→G	D99E (GAT→GAG)	rplA →	50S ribosomal subunit protein L1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,179,175	0	T	G	100.0%	178.6 / NA	49	D99E (GAT→GAG)	rplA	50S ribosomal subunit protein L1
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (18/31);  total (18/31)

GACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCTGCAGGCGCAGAACTGGTAGGTATGGAAGATCTGGCTGACCAGATCAAGAAAGGCGAAATGAACTTTGACGTTGTTATTGCTTCTCCGGATGCAATGCGCGTTGTTGGCCAGCTGGGCCAGGTTCTGGGTCCGCGCGGCCTGATGCCAAACCCGAAAGTGGGTACTGTAACACCGAA  >  NC_000913/4179029‑4179321                                                                                                                                                   |                                                                                                                                                   gACGCTCGTCAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCTGCAGGCGCAGAACTGGTAGGTATGGAAGAGCTgg                                                                                                                                                <  2:84240/151‑1 (MQ=255) gACGCTCGTCAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCTGCAGGCGCAGAACTGGTAGGTATGGAAGAGCTgg                                                                                                                                                >  1:267263/1‑151 (MQ=255)  aCGCTCGTCAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCTGCAGGCGCAGAACTGGTAGGTATGGAAGAGCTGGc                                                                                                                                               <  1:274974/151‑1 (MQ=255)                     gAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCTGCAGGCGCAGAACTGGTAGGTATGGAAGAGCTGGCTGACCAGATCAAGAAAGGc                                                                                                                            <  1:422052/151‑1 (MQ=255)                                  gCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCTGCAGGCGCAGAACTGGTAGGTATGGAAGAGCTGGCTGACCAGATCAAGAAAGGCGAAATGAACTTTg                                                                                                               <  2:665431/151‑1 (MQ=255)                                  gCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCTGCAGGCGCAGAACTGGTAGGTATGGAAGAGCTGGCTGACCAGATCAAGAAAGGCGAAATGAACTTTg                                                                                                               <  1:456707/151‑1 (MQ=255)                                  gCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCTGCAGGCGCAGAACTGGTAGGTATGGAAGAGCTGGCTGACCAGATCAAGAAAGGCGAAATGAACTTTg                                                                                                               >  2:803080/1‑151 (MQ=255)                                    aaCTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCTGCAGGCGCAGAACTGGTAGGTATGGAAGAGCTGGCTGACCAGATCAAGAAAGGCGAAATGAACTTTGa                                                                                                              <  2:912644/150‑1 (MQ=255)                                     aCTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCTGCAGGCGCAGAACTGGTAGGTATGGAAGAGCTGGCTGACCAGATCAAGAAAGGCGAAATGAACTTTGACg                                                                                                            >  1:391016/1‑151 (MQ=255)                                     aCTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCTGCAGGCGCAGAACTGGTAGGTATGGAAGAGCTGGCTGACCAGATCAAGAAAGGCGAAATGAACTTTGACg                                                                                                            <  1:414753/151‑1 (MQ=255)                                        gTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCTGCAGGCGCAGAACTGGTAGGTATGGAAGAGCTGGCTGACCAGATCAAGAAAGGCGAAATGAACTTTGACgttg                                                                                                         <  2:260997/151‑1 (MQ=255)                                        gTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCTGCAGGCGCAGAACTGGTAGGTATGGAAGAGCTGGCTGACCAGATCAAGAAAGGCGAAATGAACTTTGACgt                                                                                                           <  2:71778/149‑1 (MQ=255)                                          aCTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCTGCAGGCGCAGAACTGGTAGGTATGGAAGAGCTGGCTGACCAGATCAAGAAAGGCGAAATGAACTTTGACgttgtt                                                                                                       >  1:571370/1‑151 (MQ=255)                                            tGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCTGCAGGCGCAGAACTGGTAGGTATGGAAGAGCTGGCTGACCAGATCAAGAAAGGCGAAATGAACTTTGACGTTGTTAt                                                                                                     >  2:726432/1‑151 (MQ=255)                                              ccGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCTGCAGGCGCAGAACTGGTAGGTATGGAAGAGCTGGCTGACCAGATCAAGAAAGGCGAAATGAACTTTGACGTTGTTATTg                                                                                                   <  2:412835/151‑1 (MQ=255)                                              ccGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCTGCAGGCGCAGAACTGGTAGGTATGGAAGAGCTGGCTGACCAGATCAAGAAAGGCGAAATGAACTTTGACGTTGTTATTg                                                                                                   >  2:513040/1‑151 (MQ=255)                                                gCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCTGCAGGCGCAGAACTGGTAGGTATGGAAGAGCTGGCTGACCAGATCAAGAAAGGCGAAATGAACTTTGACGTTGTTATTGCt                                                                                                 >  2:641133/1‑151 (MQ=255)                                                 cACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCTGCAGGCGCAGAACTAGTAGGTATGGAAGAGCTGGCTGACCAGATCAAGAAAGGCGAAATGAACTTTGACGTTGTTATTGCtt                                                                                                <  2:304387/151‑1 (MQ=255)                                                   cGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCTGCAGGCGCAGAACTGGTAGGTATGGAAGAGCTGGCTGACCAGATCAAGAAAGGCGAAATGAACTTTGACGTTGTTATTGCTtct                                                                                              >  1:353488/1‑151 (MQ=255)                                                     gTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCTGCAGGCGCAGAACTGGTAGGTATGGAAGAGCTGGCTGACCAGATCAAGAAAGGCGAAATGAACTTTGACGTTGTTATTGCTTCTcc                                                                                            >  1:191295/1‑151 (MQ=255)                                                       aCTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCTGCAGGCGCAGAACTGGTAGGTATGGAAGAGCTGGCTGACCAGATCAAGAAAGGCGAAATGAACTTTGACGTTGTTATTGCTTCTCCg                                                                                           <  2:174085/150‑1 (MQ=255)                                                        cTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCTGCAGGCGCAGAACTGGTAGGTATGGAAGAGCTGGCTGACCAGATCAAGAAAGGCGAAATGAACTTTGACGTTGTTATTGCTTCTCCgg                                                                                          <  1:785081/150‑1 (MQ=255)                                                          ggCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCTGCAGGCGCAGAACTGGTAGGTATGGAAGAGCTGGCTGACCAGATCAAGAAAGGCGAAATGAACTTTGACGTTGTTATTGCTTCTCCGGATg                                                                                       <  2:472118/151‑1 (MQ=255)                                                             cgttccgttcGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCTGCAGGCGCAGAACTGGTAGGTATGGAAGAGCTGGCTGACCAGATCAAGAAAGGCGAAATGAACTTTGACGTTGTTATTGCTTCTCCGGATGCaa                                                                                    <  1:570375/151‑1 (MQ=255)                                                             cgttccgttcGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCTGCAGGCGCAGAACTGGTAGGTATGGAAGAGCTGGCTGACCAGATCAAGAAAGGCGAAATGAACTTTGACGTTGTTATTGCTTCTCCGGATGCaa                                                                                    <  2:573144/151‑1 (MQ=255)                                                              gttccgttcGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCTGCAGGCGCAGAACTGGTAGGGATGGAAGAGCTGGCTGACCAGATCAAGAAAGGCGAAATGAACTTTGACGTTGTTATTGCTTCTCCGGATGCAAt                                                                                   >  1:799461/1‑151 (MQ=255)                                                                     tcGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCTGCAGGCGCAGAACTGGTAGGTATGGAAGAGCTGGCTGACCAGATCAAGAAAGGCGAAATGAACTTTGACGTTGTTATTGCTTCTCCGGATGCAATGCGCgtt                                                                            >  2:3379/1‑151 (MQ=255)                                                                     tcGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCTGCAGGCGCAGAACTGGTAGGTATGGAAGAGCTGGCTGACCAGATCAAGAAAGGCGAAATGAACTTTGACGTTGTTATTGCTTCTCCGGATGCAATGCGCgt                                                                             <  1:469120/150‑1 (MQ=255)                                                                        cgTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCTGCAGGCGCAGAACTGGTAGGTATGGAAGAGCTGGCTGACCAGATCAAGAAAGGCGAAATGAACTTTGACGTTGTTATTGCTTCTCCGGATGCAATGCGCgttgtt                                                                         <  1:380416/151‑1 (MQ=255)                                                                                 aTTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCTGCAGGCGCAGAACTGGTAGGTATGGAAGAGCTGGCTGACCAGATCAAGAAAGGCGAAATGAACTTTGACGTTGTTATTGCTTCTCCGGATGCAATGCGCGTTGTTGGCCAGCTg                                                                >  2:815760/1‑151 (MQ=255)                                                                                            gTGCAAACGCTGAAGCTGCTAAAGCTGCAGGCGCAGAACTGGTAGGTATGGAAGAGCTGGCTGACCAGATCAAGAAAGGCGAAATGAACTTTGACGTTGTTATTGCTTCTCCGGATGCAATGCGCGTTGTTGGCCAGCTGGGCCAGGTTCt                                                     <  1:514972/151‑1 (MQ=255)                                                                                               cAAACGCTGAAGCTGCTAAAGCTGCAGGCGCAGAACTGGTAGGTATGGAAGAGCTGGCTGACCAGATCAAGAAAGGCGAAATGAACTTTGACGTTGTTATTGCTTCTCCGGATGCAATGCGCGTTGTTGGCCAGCTGGGCCAGGTTCTggg                                                  <  1:91553/151‑1 (MQ=255)                                                                                                aaaCGCTGAAGCTGCTAAAGCTGCAGGCGCAGAACTGGTAGGTATGGAAGAGCTGGCTGACCAGATCAAGAAAGGCGAAATGAACTTTGACGTTGTTATTGCTTCTCCGGATGCAATGCGCGTTGTTGGCCAGCTGGGCCAGGTTCTGGGt                                                 <  2:347497/151‑1 (MQ=255)                                                                                                                  aGCTGCAGGCGCAGAACTGGTAGGTATGGAAGAGCTGGCTGACCAGATCAAGAAAGGCGAAATGAACTTTGACGTTGTTATTGCTTCTCCGGATGCAATGCGCGTTGTTGGCCAGCTGGGCCAGGTTCTGGGTCCGCGCGGCCTGATGcc                                >  1:169524/1‑150 (MQ=255)                                                                                                                     tGCAGGCGCAGAACTGGTAGGTATGGAAGAGCTGGCTGACCAGATCAAGAAAGGCGAAATGAACTTTGACGTTGTTATTGCTTCTCCGGATGCAATGCGCGTTGTTGGCCAGCTGGGCCAGGTTCTGGGTCCGCGCGGCCTGATGCCAAAc                            >  1:652442/1‑151 (MQ=255)                                                                                                                      gCAGGCGCAGAACTGGTAGGTATGGAAGAGCTGGCTGACCAGATCAAGAAAGGCGAAATGAACTTTGACGTTGTTATTGCTTCTCCGGATGCAATGCGCGTTGTTGGCCAGCTGGGCCAGGTTCTGGGTCCGCGCGGCCTGATGCCAAAc                            <  1:615274/150‑1 (MQ=255)                                                                                                                           cgcAGAACTGGTAGGTATGGAAGAGCTGGCTGACCAGATCAAGAAAGGCGAAATGAACTTTGACGTTGTTATTGCTTCTCCGGATGCAATGCGCGTTGTTGGCCAGCTGGGCCAGGTTCTGGGTCCGCGCGGCCTGATGCCAAACCCGaa                       >  1:140431/1‑150 (MQ=255)                                                                                                                            gcAGAACTGGTAGGTATGGAAGAGCTGGCTGACCAGATCAAGAAAGGCGAAATGAACTTTGACGTTGTTATTGCTTCTCCGGATGCAATGCGCGTTGTTGGCCAGCTGGGCCAGGTTCTGGGTCCGCGCGGCCTGATGCCAAACCCGAAAg                     <  1:620222/151‑1 (MQ=255)                                                                                                                            gcAGAACTGGTAGGTATGGAAGAGCTGGCTGACCAGATCAAGAAAGGCGAAATGAACTTTGACGTTGTTATTGCTTCTCCGGATGCAATGCGCGTTGTTGGCCAGCTGGGCCAGGTTCTGGGTCCGCGCGGCCTGATGCCAAACCCGAAAg                     <  1:329846/151‑1 (MQ=255)                                                                                                                             cAGAACTGGTAGGTATGGAAGAGCTGGCTGACCAGATCAAGAAAGGCGAAATGAACTTTGACGTTGTTATTGCTTCTCCGGATGCAATGCGCGTTGTTGGCCAGCTGGGCCAGGTTCTGGGTCCGCGCGGCCTGATGCCAAACCCGAAAg                     >  2:777305/1‑150 (MQ=255)                                                                                                                             cAGAACTGGTAGGTATGGAAGAGCTGGCTGACCAGATCAAGAAAGGCGAAATGAACTTTGACGTTGTTATTGCTTCTCCGGATGCAATGCGCGTTGTTGGCCAGCTGGGCCAGGTTCTGGGTCCGCGCGGCCTGATGCCAAACCCGAAAGt                    <  1:518987/151‑1 (MQ=255)                                                                                                                              aGAACTGGTAGGTATGGAAGAGCTGGCTGACCAGATCAAGAAAGGCGAAATGAACTTTGACGTTGTTATTGCTTCTCCGGATGCAATGCGCGTTGTTGGCCAGCTGGGCCAGGTTCTGGGTCCGCGCGGCCTGATGCCAAACCCGAAAGTg                   >  1:176491/1‑151 (MQ=255)                                                                                                                              aGAACTGGTAGGTATGGAAGAGCTGGCTGACCAGATCAAGAAAGGCGAAATGAACTTTGACGTTGTTATTGCTTCTCCGGATGCAATGCGCGTTGTTGGCCAGCTGGGCCAGGTTCTGGGTCCGCGCGGCCTGATGCCAAACCCGAAAGTg                   <  2:385807/151‑1 (MQ=255)                                                                                                                               gAACTGGTAGGTATGGAAGAGCTGGCTGACCAGATCAAGAAAGGCGAAATGAACTTTGACGTTGTTATTGCTTCTCCGGATGCAATGCGCGTTGTTGGCCAGCTGGGCCAGGTTCTGGGTCCGCGCGGCCTGATGCCAAACCCGAAAGTgg                  <  2:194371/151‑1 (MQ=255)                                                                                                                                aaCTGGTAGGTATGGAAGAGCTGGCTGACCAGATCAAGAAAGGCGAAATGAACTTTGACGTTGTTATTGCTTCTCCGGATGCAATGCGCGTTGTTGGCCAGCTGGGCCAGGTTCTGGGTCCGCGCGGCCTGATGCCAAACCCGAAAGTggg                 >  2:213299/1‑151 (MQ=255)                                                                                                                                 aCTGGTAGGTATGGAAGAGCTGGCTGACCAGATCAAGAAAGGCGAAATGAACTTTGACGTTGTTATTGCTTCTCCGGATGCAATGCGCGTTGTTGGCCAGCTGGGCCAGGTTCTGGGTCCGCGCGGCCTGATGCCAAACCCGAAAGTGGGt                <  2:755988/151‑1 (MQ=255)                                                                                                                                    ggtaggtaTGGAAGAGCTGGCTGACCAGATCAAGAAAGGCGAAATGAACTTTGACGTTGTTATTGCTTCTCCGGATGCAATGCGCGTTGTTGGCCAGCTGGGCCAGGTTCTGGGTCCGCGCGGCCTGATGCCAAACCCGAAAGTGGGTACt             <  2:607529/151‑1 (MQ=255)                                                                                                                                     gtaggtaTGGAAGAGCTGGCTGACCAGATCAAGAAAGGCGAAATGAACTTTGACGTTGTTATTGCTTCTCCGGATGCAATGCGCGTTGTTGGCCAGCTGGGCCAGGTTCTGGGTCCGCGCGGCCTGATGCCAAACCCGAAAGTGGGTACTg            <  2:755234/151‑1 (MQ=255)                                                                                                                                               aaGAGCTGGCTGACCAGATCAAGAAAGGCGAAATGAACTTTGACGTTGTTATTGCTTCTCCGGATGCAATGCGCGTTGTTGGCCAGCTGGGCCAGGTTCTGGGTCCGCGCGGCCTGATGCCAAACCCGAAAGTGGGTACTGTAACACCGaa  <  1:534685/151‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                   GACGCTCGTAAATCTGACCAGAACGTACGTGGTGCAACTGTACTGCCGCACGGTACTGGCCGTTCCGTTCGCGTAGCCGTATTTACCCAAGGTGCAAACGCTGAAGCTGCTAAAGCTGCAGGCGCAGAACTGGTAGGTATGGAAGATCTGGCTGACCAGATCAAGAAAGGCGAAATGAACTTTGACGTTGTTATTGCTTCTCCGGATGCAATGCGCGTTGTTGGCCAGCTGGGCCAGGTTCTGGGTCCGCGCGGCCTGATGCCAAACCCGAAAGTGGGTACTGTAACACCGAA  >  NC_000913/4179029‑4179321

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