Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,127,160	G→A	D50N (GAT→AAT)	rpmE →	50S ribosomal subunit protein L31

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,127,160	0	G	A	100.0%	127.7 / NA	35	D50N (GAT→AAT)	rpmE	50S ribosomal subunit protein L31
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (19/16);  total (19/16)

ATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGA  >  NC_000913/4127082‑4127232                                                                               |                                                                         aTCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTAAt                                                                        <  1:768614/81‑1 (MQ=255)  tCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTAATg                                                                       <  1:768615/81‑1 (MQ=255)   ccGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTAATGt                                                                      >  2:768676/1‑81 (MQ=255)    cGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTAATGtt                                                                     <  2:768615/81‑1 (MQ=255)    cGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTAATGtt                                                                     <  2:768616/81‑1 (MQ=255)    cGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTAATGtt                                                                     >  1:768658/1‑81 (MQ=255)      cTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTAATGTTGc                                                                   <  1:768618/81‑1 (MQ=255)        ccACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTAATGTTGCTa                                                                 >  1:768675/1‑81 (MQ=255)          aCCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTAATGTTGCTAcc                                                               <  1:768619/81‑1 (MQ=255)          aCCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTAATGTTGCTAcc                                                               >  1:768673/1‑81 (MQ=255)           ccGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTAATGTTGCTACCg                                                              >  1:768676/1‑81 (MQ=255)             gTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTAATGTTGCTACCGGt                                                            <  2:768619/81‑1 (MQ=255)                ggTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTAATGTTGCTACCGGTgg                                                          <  1:768621/80‑1 (MQ=255)                  tCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTAATGTTGCTACCGGTGGCCg                                                       <  1:768622/81‑1 (MQ=255)                  tCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTAATGTTGCTACCGGTGGCCg                                                       >  1:768664/1‑81 (MQ=255)                          tGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTAATGTTGCTACCGGTGGCCGTGTTGAcc                                               >  2:768673/1‑81 (MQ=255)                            aaCCTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTAATGTTGCTACCGGTGGCCGTGTTGACCGc                                             >  1:768659/1‑81 (MQ=255)                                    cGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTAATGTTGCTACCGGTGGCCGTGTTGACCGCTTcaacaa                                     <  1:768623/81‑1 (MQ=255)                                         gcagcaAGTGCCACCCGTTCTTCACTGGCAAACAGCGTAATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGtt                                >  2:768687/1‑81 (MQ=255)                                           agcaAGTGCCACCCGTTCTTCACTGGCAAACAGCGTAATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTc                              <  1:768624/81‑1 (MQ=255)                                            gcaAGTGCCACCCGTTCTTCACTGGCAAACAGCGTAATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCa                             <  2:768624/81‑1 (MQ=255)                                                  gCCACCCGTTCTTCACTGGCAAACAGCGTAATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATcc                       >  2:768666/1‑81 (MQ=255)                                                  gCCACCCGTTCTTCACTGGCAAACAGCGTAATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATcc                       >  1:768669/1‑81 (MQ=255)                                                   ccACCCGTTCTTCACTGGCAAACAGCGTAATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATccc                      >  2:768667/1‑81 (MQ=255)                                                   ccACCCGTTCTTCACTGGCAAACAGCGTAATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATccc                      >  2:768679/1‑81 (MQ=255)                                                    cACCCGTTCTTCACTGGCAAACAGCGTAATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCg                     >  2:768681/1‑81 (MQ=255)                                                     aCCCGTTCTTCACTGGCAAACAGCGTAATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCgg                    <  2:768625/81‑1 (MQ=255)                                                       ccGTTCTTCACTGGCAAACAGCGTAATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGgc                  <  2:768626/81‑1 (MQ=255)                                                         gTTCTTCACTGGCAAACAGCGTAATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGgcag                >  2:768670/1‑81 (MQ=255)                                                         gTTCTTCACTGGCAAACAGCGTAATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGgcag                <  1:768628/81‑1 (MQ=255)                                                           tcttcACTGGCAAACAGCGTAATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGgcagc               >  2:768669/1‑80 (MQ=255)                                                              tcACTGGCAAACAGCGTAATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCaaa            >  1:768678/1‑80 (MQ=255)                                                              tcACTGGCAAACAGCGTAATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCaaa            >  2:768628/1‑80 (MQ=255)                                                               cACTGGCAAACAGCGTAATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAAt           <  2:768629/80‑1 (MQ=255)                                                                       aaCAGCGTAATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGa  >  1:768672/1‑81 (MQ=255)                                                                               |                                                                         ATCCGCTCCACCGTTGGTCATGACCTGAACCTCGACGTGTGCAGCAAGTGCCACCCGTTCTTCACTGGCAAACAGCGTGATGTTGCTACCGGTGGCCGTGTTGACCGCTTCAACAAGCGTTTCAACATCCCGGGCAGCAAATAAGACCGGA  >  NC_000913/4127082‑4127232

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