Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,409,596	C→A	A176E (GCG→GAG)	prmA →	methyltransferase for 50S ribosomal subunit protein L11

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,409,596	0	C	A	100.0%	182.1 / NA	54	A176E (GCG→GAG)	prmA	methyltransferase for 50S ribosomal subunit protein L11
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (29/25);  total (29/25)

CCAGGGCTGGCGTTTGGTACGGGTACCCATCCAACCACCTCTCTGTGCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGGCGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGTAATGGCGTTTCTGACCGTCTGGAACTCTACTTACCGAAAGATCAGCCAGAAG  >  NC_000913/3409466‑3409730                                                                                                                                   |                                                                                                                                       ccAGGGCTGGCGTTTGGTACGGGTACCCATCCAACCACCTCTCTGTGCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGGAGCTGAAACt                                                                                                                               >  1:519591/1‑140 (MQ=255)   aGGGCTGGCGTTTGGTACGGGTACCCATCCAACCACCTCTCTGTGCCTGCAATGGCTCTACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGGAGCTGAAACTgg                                                                                                                             >  2:313762/1‑140 (MQ=255)   aGGGCTGGCGTTTGGTACGGGTACCCATCCAACCACCTCTCTGTGCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGGAGCTGAAACt                                                                                                                               >  1:211099/1‑138 (MQ=255)   aGGGCTGGCGTTTGGTACGGGTACCCATCCAACCACCTCTCTGTGCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGGAGCTGAAACt                                                                                                                               <  2:211099/138‑1 (MQ=255)   aGGGCTGGCGTTTGGTACGGGTACCCATCCAACCACCTCTCTGTGCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGGAGCTGAAACTgg                                                                                                                             >  1:731381/1‑140 (MQ=255)   aGGGCTGGCGTTTGGTACGGGTACCCATCCAACCACCTCTCTGTGCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGGAGCTGAAACTgg                                                                                                                             >  2:813599/1‑140 (MQ=255)     ggCTGGCGTTTGGTACGGGTACCCATCCAACCACCTCTCTGTGCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGGAGCTGAAACTGGGt                                                                                                                           >  2:375750/1‑140 (MQ=255)          gCGTTTGGTACGGGTACCCATCCAACCACCTCTCTGTGCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGGAGCTGAAACTGGGTgcagc                                                                                                                      >  1:838766/1‑140 (MQ=255)                  tACGGGTACCCATCCACCCACCTCTCTGTGCCTGCAATGGCTCGACAGCCTGGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGGAGCTGAAACTGGGTGCAGCAAAAGCCa                                                                                                              <  1:180933/140‑1 (MQ=255)                  tACGGGTACCCATCCAACCACCTCTCTGTGCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGGAGCTGAAACTGGGTGCAGCAAAAGCCa                                                                                                              <  1:375750/140‑1 (MQ=255)                   tcGGGTACCCATCCAACCACCTCTCTGTGCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGGAGCTGAAACTGGGTGCAGCAAAAGCCAt                                                                                                             <  1:2608/139‑1 (MQ=255)                   aCGGGTACCCATCCAACCACCTCTCTGTGCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGGAGCTGAAACTGGGTGCAGCAAAAGCCAt                                                                                                             <  1:452384/140‑1 (MQ=255)                    cGGGTACCCATCCAACCACCTCTCTGTGCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGGAGCTGAAACTGGGTGCAGCAAAAGCCAtt                                                                                                            <  2:729153/140‑1 (MQ=255)                     gggTACCCATCCAACCACCTCTCTGTGCCTGCAATGGCTCGACAGCCTCGATTTAACCGGAAAAACCGTCATCGAATTTAGCTCTGGTTCCCGCATTCCGGCGATCGCGGACCTTCAACTGGCTGCAGCAAAACCCATTg                                                                                                           >  2:180933/1‑140 (MQ=255)                          cccATCCAACCACCTCTCTGTGCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGGAGCTGAAACTGGGTGCAGCAAAAGCCATTGGTAtt                                                                                                      >  2:455526/1‑140 (MQ=255)                           ccATCCAACCACCTCTCTGTGCCTGCAATGGCTCTACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGGAGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTg                                                                                                     <  1:313762/140‑1 (MQ=255)                                cAACCACCTCTCTGTGCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGGAGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATc                                                                                                >  1:225203/1‑140 (MQ=255)                                  accaccTCTCTGTGCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGGAGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGa                                                                                              >  2:514330/1‑140 (MQ=255)                                         ctctGTGCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGGAGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAg                                                                                       >  2:120949/1‑140 (MQ=255)                                           ctGTGCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGGAGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGc                                                                                     >  1:703960/1‑140 (MQ=255)                                           ctGTGCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGGAGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGc                                                                                     >  1:255663/1‑140 (MQ=255)                                              tgCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGGAGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGAt                                                                                  >  1:633145/1‑140 (MQ=255)                                                   gCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGGAGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAgg                                                                             <  2:703960/140‑1 (MQ=255)                                                   gCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGGAGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAgg                                                                             >  2:564009/1‑140 (MQ=255)                                                      aTGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGGAGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCa                                                                          >  1:166286/1‑140 (MQ=255)                                                        ggCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGGAGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGc                                                                        <  2:703271/140‑1 (MQ=255)                                                          cTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGGAGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCcg                                                                      >  2:353078/1‑140 (MQ=255)                                                              aCAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGGAGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGAt                                                                  <  1:762701/140‑1 (MQ=255)                                                                  ccTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGGAGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACg                                                              <  1:564009/140‑1 (MQ=255)                                                                  ccTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGGAGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACg                                                              >  2:33665/1‑140 (MQ=255)                                                                   cTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGGAGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGc                                                             >  2:433941/1‑140 (MQ=255)                                                                   cTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGGAGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGc                                                             >  1:488966/1‑140 (MQ=255)                                                                    tCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGGAGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGc                                                                        <  2:779791/128‑1 (MQ=255)                                                                    tCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGGAGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGc                                                                        >  1:779791/1‑128 (MQ=255)                                                                    tCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGGAGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGcc                                                            <  1:514330/140‑1 (MQ=255)                                                                    tCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGGAGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGcc                                                            <  2:519591/140‑1 (MQ=255)                                                                      gATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGGAGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGa                                                          <  1:33665/140‑1 (MQ=255)                                                                        tttAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGGAGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAAc                                                        <  2:731381/140‑1 (MQ=255)                                                                              cGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGGAGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGTAATg                                                  <  1:353078/140‑1 (MQ=255)                                                                                 tAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGGAGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGTAATGGCg                                               >  1:584296/1‑140 (MQ=255)                                                                                  aaaaCAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGGAGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGTAATGGCGt                                              <  1:813599/140‑1 (MQ=255)                                                                                        gTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGGAGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGTAATGGCGTTTCTGa                                        <  1:433941/140‑1 (MQ=255)                                                                                         tCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGGAGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGTAATGGCGTTTCTGAc                                       <  2:166286/140‑1 (MQ=255)                                                                                          cATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGGAGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGTAATGGCGTTTCTGAcc                                      <  2:633145/140‑1 (MQ=255)                                                                                           aTCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGGAGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGTAATGGCGTTTCTGACCg                                     >  1:819404/1‑140 (MQ=255)                                                                                             cGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGGAGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGTAATGGCGTTTCTGACCGTc                                   <  2:626494/140‑1 (MQ=255)                                                                                             cGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGGAGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGTAATGGCGTTTCTGACCGTc                                   <  1:445780/140‑1 (MQ=255)                                                                                              gACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGGAGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGTAATGGCGTTTCTGACCGTCt                                  >  1:81805/1‑140 (MQ=255)                                                                                                   tGGCTGTGGTTCCGGCATTCTGGCGATCGCGGAGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGTAATGGCGTTTCTGACCGTCTGGAAc                             >  1:835108/1‑140 (MQ=255)                                                                                                         tgGTTCCGGCATTCTGGCGATCGCGGAGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGTAATGGCGTTTCTGACCGTCTGGAACTCTACt                       <  2:584296/140‑1 (MQ=255)                                                                                                         tgGTTCCGGCATTCTGGCGATCGCGGAGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGTAATGGCGTTTCTGACCGTCTGGAACTCTACt                       >  1:675381/1‑140 (MQ=255)                                                                                                         tgGTTCCGGCATTCTGGCGATCGCGGAGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGTAATGGCGTTTCTGACCGTCTGGAACTCTACt                       >  1:496271/1‑140 (MQ=255)                                                                                                                    ttCTGGCGATCGCGGAGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGTAATGGCGTTTCTGACCGTCTGGAACTCTACTTACCGAAAGAt            <  2:488966/140‑1 (MQ=255)                                                                                                                    ttCTGGCGAGCGAGGAGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCTATAACGCCGAACGTAATGGCGGTTCTGACCGTCTGGAACTCTACTTACCGAAAGAt            <  1:501535/140‑1 (MQ=255)                                                                                                                     tCTGGCGATCGCGGAGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGTAATGGCGTTTCTGACCGTCTGGAACTCTACTTACCGAAAGATc           >  2:375825/1‑140 (MQ=255)                                                                                                                              cgcgGAGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGGAATGGCGTTTCTGACCGTCTGGAACTCTACTTACCGAAAGATCAGCCagaag  >  2:868398/1‑140 (MQ=255)                                                                                                                                   |                                                                                                                                       CCAGGGCTGGCGTTTGGTACGGGTACCCATCCAACCACCTCTCTGTGCCTGCAATGGCTCGACAGCCTCGATTTAACCGGTAAAACAGTCATCGACTTTGGCTGTGGTTCCGGCATTCTGGCGATCGCGGCGCTGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGTAATGGCGTTTCTGACCGTCTGGAACTCTACTTACCGAAAGATCAGCCAGAAG  >  NC_000913/3409466‑3409730

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