Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,436,053	C→T	L283F (CTC→TTC)	rsmB →	16S rRNA m(5)C967 methyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,436,053	0	C	T	100.0%	189.4 / NA	54	L283F (CTC→TTC)	rsmB	16S rRNA m(5)C967 methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (29/25);  total (29/25)

TGCATGACCTGGCTTGCGCCACAAAACGGTGAACACATTTTGGATCTTTGTGCCGCCCCCGGCGGTAAAACAACGCATATCCTTGAGGTGGCACCAGAAGCGCAGGTTGTTGCGGTTGATATCGACGAACAGCGCCTCTCTCGGGTTTACGACAATTTAAAACGCCTTGGTATGAAGGCGACCGTGAAACAAGGTGATGGCCGTTACCCTTCTCAATGGTGTGGCGAGCAACAGTTTGATCGCATTTTATTAGATGCGCCTTGTT  >  NC_000913/3435918‑3436182                                                                                                                                        |                                                                                                                                  tGCATGACCTGGCTTGCGCCACAAAACGGTGAACACATTTTGGATCTTTGTGCCGCCCCCGGCGGTAAAACAACGCATATCCTTGAGGTGGCACCAGAAGCGCAGGTTGTTGCGGTTGATATCGACGAACAGCGCttctc                                                                                                                               <  1:130017/140‑1 (MQ=255)   cATGACCTGGCTTGCGCCACAAAACGGTGAACACATTTTGGATCTTTGTGCCGCCCCCGGCGGTAAAACAACGCATATCCTTGAGGTGGCACCAGAAGCGCAGGTTGTTGCGGTTGATATCGACGAACAGCGCTtctctc                                                                                                                             <  2:243506/140‑1 (MQ=255)                 cgcCACAAAACGGTGAACACATTTTGGATCTTTGTGCCGCCCCCGGCGGTAAAACAACGCATATCCTTGAGGTGGCACCAGAAGCGCAGGTTGTTGCGGTTGATATCGACGAACAGCGCTtctctc                                                                                                                             >  2:60047/1‑126 (MQ=255)                 cgcCACAAAACGGTGAACACATTTTGGATCTTTGTGCCGCCCCCGGCGGTAAAACAACGCATATCCTTGAGGTGGCACCAGAAGCGCAGGTTGTTGCGGTTGATATCGACGAACAGCGCTtctctc                                                                                                                             <  1:60047/126‑1 (MQ=255)                 cgcCACAAAACGGTGAACACATTTTGGATCTTTGTGCCGCCCCCGGCGGTAAAACAACGCATATCCTTGAGGTGGCACCAGAAGCGCAGGTTGTTGCGGTTGATATCGACGAACAGCGCTTCTCTCGGGTTTACGACAAt                                                                                                               >  2:357027/1‑140 (MQ=255)                  gcCACAAAACGGTGAACACATTTTGGATCTTTGTGCCGCCCCCGGCGGTAAAACAACGCATATCCTTGAGGTGGCACCAGAAGCGCAGGTTGTTGCGGTTGATATCGACGAACAGCGCTTCTCTCGGGTTTACGACAAtt                                                                                                              >  1:159857/1‑140 (MQ=255)                  gcCACAAAACGGTGAACACATTTTGGATCTTTGTGCCGCCCCCGGCGGTAAAACAACGCATATCCTTGAGGTGGCACCAGAAGCGCAGGTTGTTGCGGTTGATATCGACGAACAGCGCTTCTCTCGGGTTTACGACAAtt                                                                                                              >  2:310100/1‑140 (MQ=255)                       aaaaCGGTGAACACATTTTGGATCTTTGTGCCGCCCCCGGCGGTAAAACAACGCATATCCTTGAGGTGGCACCAGAAGCGCAGGTTGTTGCGGTTGATATCGACGAACAGCGCTTCTCTCGGGTTTACGACAATTTaaaa                                                                                                         <  1:357027/140‑1 (MQ=255)                                aaCACATTTTGGATCTTTGTGCCGCCCCCGGCGGTAAAACAACGCATATCCTTGAGGTGGCACCAGAAGCGCAGGTTGTTGCGGTTGATATCGACGAACAGCGCTTCTCTCGGGTTTACGACAATTTAAAACGCCTTGGt                                                                                                >  1:556422/1‑140 (MQ=255)                                      ttttGGATCTTTGTGCCGCCCCCGGCGGTAAAACAACGCATATCCTTGAGGTGGCACCAGAAGCGCAGGTTGTTGCGGTTGATATCGACGAACAGCGCTTCTCTCGGGTTTACGACAATTTAAAACGCCTTGGTATGAAg                                                                                          >  2:640917/1‑140 (MQ=255)                                        ttGGATCTTTGTGCCGCCCCCGGCGGTAAAACAACGCATATCCTTGAGGTTGCACCAGAAGCGCAGGTTGTTGCGGTTGATATCGACGAACAGCGCTTCTCTCGGGTTTACGACAATTTAAAACGCCTTGGTATGAAGGc                                                                                        <  1:745986/140‑1 (MQ=255)                                          ggATCTTTGTGCCGCCCCCGGCGGTAAAACAACGCATATCCTTGAGGTGGCACCAGAAGCGCAGTTTGTTGCGGTTGATATCGACGAACAGCGCTTCTCTCGGGTTTACGACAATTTAAAACGCCTTGGTATGAAGGCGa                                                                                      <  2:763733/140‑1 (MQ=255)                                          ggATCTTTGTGCCGCCCCCGGCGGTAAAACAACGCATATCCTTGAGGTGGCACCAGAAGCGCAGGTTGTTGCGGTTGATATCGACGAACAGCGCTTCTCTCGGGTTTACGACAATTTAAAACGCCTTGGTATGAAGGCGa                                                                                      <  2:507200/140‑1 (MQ=255)                                            aTCTTTGTGCCGCCCCCGGCGGTAAAACAACGCATATCCTTGAGGTGGCACCAGAAGCGCAGGTTGTTGCGGTTGATATCGACGAACAGCGCTTCTCTCGGGTTTACGACAATTTAAAACGCCTTGGTATGAAGGCGAcc                                                                                    <  1:485934/140‑1 (MQ=255)                                             tCTTTGTGCCGCCCCCGGCGGTAAAACAACGCATATCCTTGAGGTGGCACCAGAAGCGCAGGTTGTTGCGGTTGATATCGACGAACAGCGCTTCTCTCGGGTTTACGACAATTTAAAACGCCTTGGTATGAAGGCGACCg                                                                                   <  2:556422/140‑1 (MQ=255)                                                 tgtgCCGCCCCCGGCGGTAAAACAACGCATATCCTTGAGGTGGCACCAGAAGCGCAGGTTGTTGCGGTTGATATCGACGAACAGCGCTTCTCTCGGGTTTACGACAATTTAAAACGCCTTGGTATGAAGGCGACCGTGaa                                                                               >  1:154514/1‑140 (MQ=255)                                                     ccgccCCCGGCGGTAAAACAACGCATATCCTTGAGGTGGCACCAGAAGCGCAGGTTGTTGCGGTTGATATCGACGAACAGCGCTTCTCTCGGGTTTACGACAATTTAAAACGCCTTGGAATGAAGGCGACCGTGAAACaa                                                                           >  2:78469/1‑140 (MQ=255)                                                       gccCCCGGCGGTAAAACAACGCATATCCTTGAGGTGGCACCAGAAGCGCAGGTTGTTGCGGTTGATATCGACGAACAGGGCTTCTCTCGGGTTTACGACAATTTAAAACGCCTTGGTATGAAGGCGACCGTGAAACAAgg                                                                         >  2:715919/1‑140 (MQ=255)                                                        cccccGGCGGTAAAACAACGCATATCCTTGAGGTGGCACCAGAAGCGCAGGTTGTTGCGGTTGATATCGACGAACAGCGCTTCTCTCGGGTTTACGACAATTTAAAACGCCTTGGTATGAAGGCGACCGTGAAACaa                                                                           <  1:18058/137‑1 (MQ=255)                                                        cccccGGCGGTAAAACAACGCATATCCTTGAGGTGGCACCAGAAGCGCAGGTTGTTGCGGTTGATATCGACGAACAGCGCTTCTCTCGGGTTTACGACAATTTAAAACGCCTTGGTATGAAGGCGACCGTGAAACaa                                                                           >  2:18058/1‑137 (MQ=255)                                                        cccccGGCGGTAAAACAACGCATATCCTTGAGGTGGCACCAGAAGCGCAGGTTGTTGCGGTTGATATCGACGAACAGCGCTTCTCTCGGGTTTACGACAATTTAAAACGCCTTGGTATGAAGGCGACCGTGAAACAAGGt                                                                        >  2:311481/1‑140 (MQ=255)                                                        cccccGGCGGTAAAACAACGCATATCCTTGAGGTGGCACCAGAAGCGCAGGTTGTTGCGGTTGATATCGACGAACAGCGCTTCTCTCGGGTTTACGACAATTTAAAACGCCTTGGTATGAAGGCGACCGTGAAACAAGGt                                                                        >  2:227713/1‑140 (MQ=255)                                                         ccccGGCGGTAAAACAACGCATATCCTTGAGGTGGCACCAGAAGCGCAGGTTGTTGCGGTTGATATCGACGAACAGCGCTTCTCTCGGGTTTACGACAATTTAAAACGCCTTGGTATGAAGGCGACCGTGAAACAAGGTg                                                                       >  2:527259/1‑140 (MQ=255)                                                          cccGGCGGTAAAACAACGCATATCCTTGAGGTGGCACCAGAAGCGCAGGTTGTTGCGGTTGATATCGACGAACAGCGCTTCTCTCGGGTTTACGACAATTTAAAACGCCTTGGTATGAAGGCGACCGTGAAACAAGGTGa                                                                      <  1:640917/140‑1 (MQ=255)                                                                  tAAAACAACGCATATCCTTGAGGTGGCACCAGAAGCGCAGGTTGTTGCGGTTGATATCGACGAACAGCGCTTCTCTCGGGTTTACGACAATTTAAAACGCCTTGGTATGAAGGCGACCGTGAAACAAGGTGATGGCCGtt                                                              <  1:310100/140‑1 (MQ=255)                                                                   aaacCAACGCGTAGCCTTGAGGTGGCTCCTGAAGCGCAGGTTGTTGCGGTTGATATCTACGAACAGCGGTTCTCTCGGGTTTACGTCAGTTTAATGCGCCTTGGTATGAAGTCGACCGCTAAGCAAGGTGATGGCCgatc                                                             <  2:282958/140‑4 (MQ=255)                                                                           gCATATCCTTGAGGTGGCACCAGAAGCGCAGGTTGTTGCGGTTGATATCGACGAACAGCGCTTCTCTCGGGTTTACGACAATTTAAAACGCCTTGGTATGAAGGCGACCGTGAAACAAGGTGATg                                                                    >  1:2968/1‑125 (MQ=255)                                                                           gCATATCCTTGAGGTGGCACCAGAAGCGCAGGTTGTTGCGGTTGATATCGACGAACAGCGCTTCTCTCGGGTTTACGACAATTTAAAACGCCTTGGTATGAAGGCGACCGTGAAACAAGGTGATg                                                                    <  2:2968/125‑1 (MQ=255)                                                                            cATATCCTTGAGGTGGCACCAGAAGCGCAGGTTGTTGCGGTTGATATCGACGAACAGCGCTTCTCTCGGGTTTACGACAATTTAAAACGCCTTGGTATGAAGGCGACCGTGAAACAAGGTGATGGCCGTTACCCTTCTCa                                                    <  1:919234/140‑1 (MQ=255)                                                                                    tGAGGTGGCACCAGAAGCGCAGGTTGTTGCGGTTGATATCGACGAACAGCGCTTCTCTCGGGTTTACGACAATTTAAAACGCCTTGGTATGAAGGCGACCGTGAAACAAGGTGATGGCCGTTACCCTTCTCAATGGtgtg                                            <  1:715919/140‑1 (MQ=255)                                                                                     gAGGTGGCACCAGAAGCGCAGGTTGTTGCGGTTGATATCGACGAACAGCGCTTCTCTCGGGTTTACGACAATTTAAAACGCCTTGGTATGAAGGCGACCGTGAAACAAGGTGATGGCCGTTACCCTTCTCAATGGTGTgg                                           >  1:721465/1‑140 (MQ=255)                                                                                       ggTGGCACCAGAAGCGCAGGTTGTTGCGGTTGATATCGACGAACAGCGCTTCTCTCGGGTTTACGACAATTTAAAACGCCTTGGTATGAAGGCGACCGTGAAACAAGGTGATGGCCGTTACCCTTCTCAATGGTGTGGCg                                         >  2:883112/1‑140 (MQ=255)                                                                                         tGGCACCAGAAGCGCAGGTTGTTGCGGTTGATATCGACGAACAGCGCTTCTCTCGGGTTTACGACAATTTAAAACGCCTTGGTATGAAGGCGACCGTGAAACAAGGTGATGGCCGTTACCCTTCTCAATGGTGTGGCGAg                                       <  1:224742/140‑1 (MQ=255)                                                                                           gCACCAGAAGCGCAGGTTGTTGCGGTTGATATCGACGAACAGCGCTTCTCTCGGGTTTACGACAATTTAAAACGCCTTGGTATGAAGGCGACCGTGAAACAAGGTGATGGCCGTTACCCTTCTCAGTGGTGTGGCGAGCa                                     >  2:591724/1‑140 (MQ=255)                                                                                           gCACCAGAAGCGCAGGTTGTTGCGGTTGATATCGACGAACAGCGCTTCTCTCGGGTTTACGACAATTTAAAACGCCTTGGTATGAAGGCGACCGTGAAACAAGGTGATGGCCGTTACCCTTCTCAATGGTGTGGCGAGc                                      <  2:445399/139‑1 (MQ=255)                                                                                           gCACCAGAAGCGCAGGTTGTTGCGGTTGATATCGACGAACAGCGCTTCTCTCGGGTTTACGACAATTTAAAACGCCTTGGTATGAAGGCGACCGTGAAACAAGGTGATGGCCGTTACCCTTCTCAATGGTGTGGCGAGc                                      >  1:445399/1‑139 (MQ=255)                                                                                           gCACCAGAAGCGCAGGTTGTTGCGGTTGATATCGACGAACAGCGCTTCTCTCGGGTTTACGACAATTTAAAACGCCTTGGTATGAAGGCGACCGTGAAACAAGGTGATGGCCGTTACCCTTCTCAATGGTGTGGCGAGCa                                     >  2:251458/1‑140 (MQ=255)                                                                                           gCACCAGAAGCGCAGGTTGTTGCGGTTGATATCGACGAACAGCGCTTCTCTCGGGTTTACGACAATTTAAAACGCCTTGGTATGAAGGCGACCGTGAAACAAGGTGATGGCCGTTACCCTTCTCAATGGTGTGGCGAGCa                                     >  1:266053/1‑140 (MQ=255)                                                                                           gCACCAGAAGCGCAGGTTGTTGCGGTTGATATCGACGAACAGCGCTTCTCTCGGGTTTACGACAATTTAAAACGCCTTGGTATGAAGGCGACCGTGAAACAAGGTGATGGCCGTTACCCTTCTCAATGGTGTGGCGAGCa                                     >  2:333601/1‑140 (MQ=255)                                                                                           gCACCAGAAGCGCAGGTTGTTGCGGTTGATATCGACGAACAGCGCTTCTCTCGGGTTTACGACAATTTAAAACGCCTTGGTATGAAGGCGACCGTGAAACAAGGTGATGGCCGTTACCCTTCTCAATGGTGTGGCGAGCa                                     >  2:270204/1‑140 (MQ=255)                                                                                           gCACCAGAAGCGCAGGTTGTTGCGGTTGATATCGACGAACAGCGCTTCTCTCGGGTTTACGACAATTTAAAACGCCTTGGTATGAAGGCGACCGTGAAACAAGGTGATGGCCGTTACCCTTCTCAATGGTGTGGCGAGCa                                     >  1:416363/1‑140 (MQ=255)                                                                                             aCCAGAAGCGCAGGTTGTTGCGGTTGATATCGACGAACAGCGCTTCTCTCGGGTTTACGACAATTTAAAACGCCTTGGTATGAAGGCGACCGTGAAACAAGGTGATGGCCGTTACCCTTCTCAATGGTGTGGCGAGCAAc                                   <  1:568205/140‑1 (MQ=255)                                                                                                aGAAGCGCAGGTTGTTGCGGTTGATATCGACGAACAGCGCTTCTCTCGGGTTTACGACAATTTAAAACGCCTTGGTATGAAGGCGACCGTGAAACAAGGTGATGGCCGTTACCCTTCTCAATGGTGTGGCGAGCAACAGt                                <  2:721465/140‑1 (MQ=255)                                                                                                    gcgcAGGTTGTTGCGGTTGATATCGACGAACAGCGCTTCTCTCGGGTTTACGACAATTTAAAACGCCTTGGTATGAAGGCGACCGTGAAACAAGGTGATGGCCGTTACCCTTCTCAATGGTGTGGCGAGCAACAGTTTGa                            <  1:344783/140‑1 (MQ=255)                                                                                                       cAGGTTGTTGCGGTTGATATCGACGAACAGCGCTTCTCTCGGGTTTACGACAATTTAAAACGCCTTGGTATGAAGGCGACCGTGAAACAAGGTGATGGCCGTTACCCTTCTCAATGGTGTGGCGAGCAACAGTTTGATCg                         <  1:883112/140‑1 (MQ=255)                                                                                                       cAGGTTGTTGCGGTTGATATCGACGAACAGCGCTTCTCTCGGGTTTACGACAATTTAAAACGCCTTGGTATGAAGGCGACCGTGAAACAAGGTGATGGCCGTTACCCTTCTCAATGGTGTGGCGAGCAACAGTTTGATCg                         >  1:374729/1‑140 (MQ=255)                                                                                                        aGGTTGTTGCGGTTGATATCGACGAACAGCGCTTCTCTCGGGTTTACGACAATTTAAAACGCCTTGGTATGAAGGCGACCGTGAAACAAGGTGATGGCCGTTACCCTTCTCAATGGTGTGGCGAGCAACAGTTTGATCGc                        <  2:374729/140‑1 (MQ=255)                                                                                                          gttgttGCGGTTGATATCGACGAACAGCGCTTCTCTCGGGTTTACGACAATTTAAAACGCCTTGGTATGAAGGCGACCGTGAAACAAGGTGATGGCCGTTACCCTTCTCAATGGTGTGGCGAGCAACAGTTTGATCGCAt                      <  1:270204/140‑1 (MQ=255)                                                                                                          gttgttGCGGTTGATATCGACGAACAGCGCTTCTCTCGGGTTTACGACAATTTAAAACGCCTTGGTATGAAGGCGACCGTGAAACAAGGTGATGGCCGTTACCCTTCTCAATGGTGTGGCGAGCAACAGTTTGATCGCAt                      <  1:333601/140‑1 (MQ=255)                                                                                                                      gATATCGACGAACAGCGCTTCTCTCGGGTTTACGACAATTTAAAACGCCTTGGTATGAAGGCGACCGTGAAACAAGGTGATGGCCGTTACCCTTCTCAATGGTGTGGCGAGCAACAGTTTGATCGCATTTTATTAGATgc          >  2:884871/1‑140 (MQ=255)                                                                                                                             acgaACAGCGCTTCTCTCGGGTTTACGACAATTTAAAACGCCTTGGTATGAAGGCGACCGTGAAACAAGGTGATGGCCGTTACCCTTCTCAATGGTGTGGCGAGCAACAGTTTGATCGCATTTTATTAGATGCGCCTTGt   >  2:246393/1‑140 (MQ=255)                                                                                                                             acgaACAGCGCTTCTCTCGGGTTTACGACAATTTAAAACGCCTTGGTATGAAGGCGACCGTGAAACAAGGTGATGGCCGTTACCCTTCTCAATGGTGTGGCGAGCAACAGTTTGATCGCATTTTATTAGATGCGCCTTGt   >  1:938479/1‑140 (MQ=255)                                                                                                                              cgaACAGCGCTTCTCTCGGGTTTACGACAATTTAAAACGCCTTGGTATGAAGGCGACCGTGAAACAAGGTGATGGCCGTTACCCTTCTCAATGGTGTGGCGAGCAACAGTTTGATCGCATTTTATTAGATGCGCCTTGtt  >  1:392893/1‑140 (MQ=255)                                                                                                                                 aCAGCGCTTCTCTCGGGTTTACGACAATTTAAAACGCCTTGGTATGAAGGCGACCGTGAAACAAGGTGATGGCCGTTACCCTTCTCAATGGTGTGGCGAGCAACAGTTTGATCGCATTTTATTAGATgcg         >  2:497905/1‑130 (MQ=255)                                                                                                                                 aCAGCGCTTCTCTCGGGTTTACGACAATTTAAAACGCCTTGGTATGAAGGCGACCGTGAAACAAGGTGATGGCCGTTACCCTTCTCAATGGTGTGGCGAGCAACAGTTTGATCGCATTTTATTAGATgcg         <  1:497905/130‑1 (MQ=255)                                                                                                                                  cAGCGCTTCTCTCGGGTTTACGACAATTTAAAACGCCTTGGTATGAAGGCGACCGTGAAACAAGGTGATGGCCGTTACCCTTCTCAATGGTGTGGCGAGCAACAGTTTGATCGCATTTTATTAGATGCGcc       <  1:846829/131‑1 (MQ=255)                                                                                                                                  cAGCGCTTCTCTCGGGTTTACGACAATTTAAAACGCCTTGGTATGAAGGCGACCGTGAAACAAGGTGATGGCCGTTACCCTTCTCAATGGTGTGGCGAGCAACAGTTTGATCGCATTTTATTAGATGCGcc       >  2:846829/1‑131 (MQ=255)                                                                                                                                        |                                                                                                                                  TGCATGACCTGGCTTGCGCCACAAAACGGTGAACACATTTTGGATCTTTGTGCCGCCCCCGGCGGTAAAACAACGCATATCCTTGAGGTGGCACCAGAAGCGCAGGTTGTTGCGGTTGATATCGACGAACAGCGCCTCTCTCGGGTTTACGACAATTTAAAACGCCTTGGTATGAAGGCGACCGTGAAACAAGGTGATGGCCGTTACCCTTCTCAATGGTGTGGCGAGCAACAGTTTGATCGCATTTTATTAGATGCGCCTTGTT  >  NC_000913/3435918‑3436182

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