Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,627,990	G→A	Q517* (CAG→TAG)	rbbA ←	ribosome‑associated ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,627,990	0	G	A	100.0%	162.6 / NA	50	Q517* (CAG→TAG)	rbbA	ribosome‑associated ATPase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (26/24);  total (26/24)

AACGTACTGGATCGCGTCGCAGTTCCAGCGCTTCGCGGCGGCTGTAGCTAAACAGACGGCGCAGGCTAAATCCCTGACGCGGCGCGTGGGTGGTGTCGTGTACCACGGGCGGCGCTTCGGCTTCGTTGCTCTGCCCTGCCGCTTCCTGCAAATAGGCGATAAATGCCTCTTCCAGACTGGCGGCTCCGCGTTTCTCAACCAGTTCCTGCGGTGTACCGCTGGCAAGCACTTTTCCGGCGTGCATCAGTGAGATGCGGT  >  NC_000913/3627858‑3628115                                                                                                                                     |                                                                                                                              aaCGTACTGGATCGCGTCGCAGTTCCAGCGCTTCGCGGCGGCTGTAGCTAAACAGACGGTGCAGGCTAAATCCCTGACGCGGCGCGTGGGTGGTGTCGTGTACCACGGGCGGCGCTTCGGCTTCGTTGCTctaccctgcc                                                                                                                        <  2:528898/140‑1 (MQ=255) aaCGTACTGGATCGCGTCGCAGTTCCAGCGCTTCGCGGCGGCTGTAGCTAAACAGACGGCGCAGGCTAAATCCCTGACGCGGCGCGTGGGTGGTGTCGTGTACCACGGGCGGCGCTTCGGCTTCGTTGCTctaccctgcc                                                                                                                        <  2:722496/140‑1 (MQ=255) aaCGTACTGGATCGCGTCGCAGTTCCAGCGCTTCGCGGCGGCTGTAGCTAAACAGACGGCGCAGGCTAAATCCCTGACGCGGCGCGTGGGTGGTGTCGTGTACCACGGGCGGCGCTTCGGCTTCGTTGCTctaccctgcc                                                                                                                        >  1:528225/1‑140 (MQ=255)   cGTACTGGATCGCGTCGCAGTTCCAGCGCTTCGCGGCGGCTGTAGCTAAACAGACGGCGCAGGCTAAATCCCTGACGCGGCGCGTGGGTGGTGTCGTGTACCACGGGCGGCGCTTCGGCTTCGTTGCTCTACCCTGCCGc                                                                                                                      <  2:528225/140‑1 (MQ=255)         ggATCGCGTCGCAGTTCCAGCGCTTCGCGGCGGCTGTAGCTAAACAGACGGCGCAGGCTAAATCCCTGACGCGGCGCGTGGGTGGTGTCGTGTACCACGGGCGGCGCTTCGGCTTCGTTGCTCTACCCTGCCGCTTCCTg                                                                                                                >  1:917385/1‑140 (MQ=255)            tCGCGTCGCAGTTCCAGCGCTTCGCGGCGGCTGTAGCTAAACAGACGGCGCAGGCTAAATCCCTGACGCGGCGCGTGGGTGGTGTCGTGTACCACGGGCGGCGCTTCGGCTTCGTTGCTCTACCCTGCCGCTTCCTGCaa                                                                                                             <  2:810145/140‑1 (MQ=255)               cgTCGCAGTTCCAGCGCTTCGCGGCGGCTGTAGCTAAACAGACGGCGCAGGCTAAATCCCTGACGCGGCGCGTGGGTGGTGTCGTGTACCACGGGCGGCGCTTCGGCTTCGTTGCTCTACCCTGCCGCTTCCTGCAAATa                                                                                                          >  2:814760/1‑140 (MQ=255)               cgTCGCAGTTCCAGCGCTTCGCGGCGGCTGTAGCTAAACAGACGGCGCAGGCTAAATCCCTGACGCGGCGCGTGGGTGGTGTCGTGTACCACGGGCGGCGCTTCGGCTTCGTTGCTCTACCCTGCCGCTTCCTGCAAATa                                                                                                          >  1:868933/1‑140 (MQ=255)               cgTCGCAGTTCCAGCGCTTCGCGGCGGCTGTAGCTAAACAGACGGCGCAGGCTAAATCCCTGACGCGGCGCGTGGGTGGTGTCGTGTACCACGGGCGGCGCTTCGGCTTCGTTGCTCTACCCTGCCGCTTCCTGCAAATa                                                                                                          >  2:428519/1‑140 (MQ=255)               cgTCGCAGTTCCAGCGCTTCGCGGCGGCTGTAGCTAAACAGACGGCGCAGGCTAAATCCCTGACGCGGCGCGTGGGTGGTGTCGTGTACCACGGGCGGCGCTTCGGCTTCGTTGCTCTACCCTGCCGCTTCCTGCAAATa                                                                                                          >  1:40724/1‑140 (MQ=255)                gTCGCAGTTCCAGCGCTTCGCGGCGGCTGTAGCTAAACAGACGGCGCAGGCTAAATCCCTGACGCGGCGCGTGGGTGGTGTCGTGTACCACGGGCGGCGCTTCGGCTTCGTTGCTCTACCCTGCCGCTTCCTGCAAATAg                                                                                                         <  2:40724/140‑1 (MQ=255)                  cGCAGTTCCAGCGCTTCGCGGCGGCTGTAGCTAAACAGACGGCGCAGGCTAAATCCCTGACGCGGCGCGTGGGTGGTGTCGTGTACCACGGGCGGCGCTTCGGCTTCGTTGCTCTACCCTGCCGCTTCCTGCAAATAGGc                                                                                                       <  1:428519/140‑1 (MQ=255)                    cAGTTCCAGCGCTTCGCGGCGGCTGTAGCTAAACAGACGGCGCAGGCTAAATCCCTGACGCGGCGCGTGGGTGGTGTCGTGTACCACGGGCGGCGCTTCGGCTTCGTTGCTCTACCCTGCCGCTTCCTGCAAATAGGCGa                                                                                                     >  2:330612/1‑140 (MQ=255)                       ttCCAGCGCTTCGCGGCGGCTGTAGCTAAACAGACGGCGCAGGCTAAATCCCTGACGCGGCGCGTGGGTGGTGTCGTGTACCACGGGCGGCGCTTCGGCTTCGTTGCTCTACCCTGCCGCTTCCTGCAAATAGGCGATaa                                                                                                  <  2:868933/140‑1 (MQ=255)                           aGCGCTTCGCGGCGGCTGTAGCTAAACAGACGGCGCAGGCTAAATCCCTGACGCGGCGCGTGGGTGGTGTCGTGTACCACGGGCGGCGCTTCGGCTTCGTTGCTCTACCCTGCCGCTTCCTGCAAATAGGCGATAAATGc                                                                                              <  1:377068/140‑1 (MQ=255)                           aGCGCTTCGCGGCGGCTGTAGCTAAACAGACGGCGCAGGCTAAATCCCTGACGCGGCGCGTGGGTGGTGTCGTGTACCACGGGCGGCGCTTCGGCTTCGTTGCTCTACCCTGCCGCTTCCTGCAAATAGGCGATAAATGc                                                                                              <  2:917385/140‑1 (MQ=255)                               cTTCGCGGCGGCTGTAGCTAAACAGACGGCGCAGGCTAAATCCCTGACGCGGCGCGTGGGTGGTGTCGTGTACCACGGGCGGCGCTTCGGCTTCGTTGCTCTACCCTGCCGCTTCCTGCAAATAGGCGATAAATGActct                                                                                          >  1:650327/1‑140 (MQ=255)                                  cgcgGCGGCTGTAGCTAAACAGACGGCGCAGGCTAAATCCCTGACGCGGCGCGTGGGTGGTGTCGTGTACCACGGGCGGCGCTTCGGCTTCGTTGCTCTACCCTGCCGCTTCCTGCAAATAGGCGATAAATGCCTCTTcc                                                                                       >  2:758122/1‑140 (MQ=255)                                                 tAAACAGACGGCGCAGGCTAAATCCCTGACGCGGCGCGTGGGTGGTGTCGTGTACCACGGGCGGCGCTTCGGCTTCGTTGCTCTACCCTGCCGCTTCCTGCAAATAGGCGATAAATGCCTCTTCCAGACTGGCGGCTCcg                                                                        <  2:158007/140‑1 (MQ=255)                                                   aaCAGACGGCGCAGGCTAAATCCCTGACGCGGCGCGTGGGTGGTGTCGTGTACCACGGGCGGCGCTTCGGCTTCGTTGCTCTACCCTGCCGCTTCCTGCAAATAGGCGATAAATGCCTCTTCCAGACTGGCGGCTCcgcg                                                                      >  2:201615/1‑140 (MQ=255)                                                        aCGGCGCAGGCTAAATCCCTGACGCGGCGCGTGGGTGGTGTCGTGTACCACGGGCGGCGCTTCGGCTTCGTTGCTCTACCCTGCCGCTTCCTGCAAATAGGCGATAAATGCCTCTTCCAGACTGGCGGCTCCGCGTTtct                                                                 >  1:222046/1‑140 (MQ=255)                                                             gcAGGCTAAATCCCTGACGCGGCGCGTGGGTGGTGTCGTGTACCACGGGCGGCGCTTCGGCTTCGTTGCTCTACCCTGCCGCTTCCTGCAAATAGGCGATAAATGCCTCTTCCAGACTGGCGGCTCCGCGTTTCTCAAcc                                                            <  1:149060/140‑1 (MQ=255)                                                             gcAGGCTAAATCCCTGACGCGGCGCGTGGGTGGTGTCGTGTACCACGGGCGGCGCTTCGGCTTCGTTGCTCTACCCTGCCGCTTCCTGCAAATAGGCGATAAATGCCTCTTCCAGACTGGCGGCTCCGCGTTTCTCAAcc                                                            >  2:221780/1‑140 (MQ=255)                                                             gcAGGCTAAATCCCTGACGCGGCGCGTGGGTGGTGTCGTGTACCACGGGCGGCGCTTCGGCTTCGTTGCTCTACCCTGCCGCTTCCTGCAAATAGGCGATAAATGCCTCTTCCAGACTGGCGGCTCCGCGTTTCTCAAcc                                                            >  2:557242/1‑140 (MQ=255)                                                             gcAGGCTAAATCCCTGACGCGGCGCGTGGGTGGTGTCGTGTACCACGGGCGGCGCTTCGGCTTCGTTGCTCTACCCTGCCGCTTCCTGCAAATAGGCGATAAATGCCTCTTCCAGACTGGCGGCTCCGCGTTTCTCAAcc                                                            <  1:346776/140‑1 (MQ=255)                                                                     aaTCCCTGACGCGGCGCGTGGGTGGTGTCGTGTACCACGGGCGGCGCTTCGGCTTCGTTGCTCTACCCTGCCGCTTCCTGCAAATAGGCGATAAATGCCTCTTCCAGACTGGCGGCTCCGCGTTTCTCaa                                                              <  1:94118/130‑1 (MQ=255)                                                                     aaTCCCTGACGCGGCGCGTGGGTGGTGTCGTGTACCACGGGCGGCGCTTCGGCTTCGTTGCTCTACCCTGCCGCTTCCTGCAAATAGGCGATAAATGCCTCTTCCAGACTGGCGGCTCCGCGTTTCTCaa                                                              >  2:94118/1‑130 (MQ=255)                                                                     aaTCCCTGACGCGGCGCGTGGGTGGTGTCGTGTACCACGGGCGGCGCTTCGGCTTCGTTGCTCTACCCTGCCGCTTCCTGCAAATAGGCGATAAATGCCTCTTCCAGACTGGCGGCTCCGCGTTTCTCAACCAGTTCCTg                                                    <  2:954186/140‑1 (MQ=255)                                                                           tGACGCGGCGCGTGGGTGGTGTCGTGTACCACGGGCGGCGCTTCGGCTTCGTTGCTCTACCCTGCCGCTTCCTGCAAATAGGCGATAAATGCCTCTTCCAGACTGGCGGCTCCGCGTTTCTCAACCAGTTCCTGCGgtgt                                              >  1:739732/1‑140 (MQ=255)                                                                           tGACGCGGCGCGTGGGTGGTGTCGTGTACCACGGGCGGCGCTTCGGCTTCGTTGCTCTACCCTGCCGCTTCCTGCAAATAGGCGATAAATGCCTCTTCCAGACTGGCGGCTCCGCGTTTCTCAACCAGTTCCTGCGgtgt                                              >  2:219103/1‑140 (MQ=255)                                                                            gACGCGGCGCGTGGGTGGTGTCGTGTACCACGGGCGGCGCTTCGGCTTCGTTGCTCTACCCTGCCGCTTCCTGCAAATAGGCGATAAATGCCTCTTCCAGACTGGCGGCTCCGCGTTTCTCAACCAGTTCCTGCGGTGTa                                             <  2:222046/140‑1 (MQ=255)                                                                                 ggcgCGTGGGTGGTGTCGTGTACCACGGGCGGCGCTTCGGCTTCGTTGCTCTACCCTGCCGCTTCCTGCAAATAGGCGATAAATGCCTCTTCCAGACTGGCGGCTCCGCGTTTCTCAACCAGTTCCTGCGGTGTACCGCt                                        >  1:206813/1‑140 (MQ=255)                                                                                 ggcgCGTGGGTGGTGTCGTGTACCACGGGCGGCGCTTCGGCTTCGTTGCTCTACCCTGCCGCTTCCTGCAAATAGGCGATAAATGCCTCTTCCAGACTGGCGGCTCCGCGTTTCTCAACCAGTTCCTGCGGTGTACCGCt                                        >  1:924967/1‑140 (MQ=255)                                                                                       tgggtggTGTCGTGTACCACGGGCGGCGCTTCGGCTTCGTTGCTCTACCCTGCCGCTTCCTGCAAATAGGCGATAAATGCCTCTTCCAGACTGGCGGCTCCGCGTTTCTCAACCAGTTCCTGCGGTGTACCGCTGGCAAg                                  <  1:814760/140‑1 (MQ=255)                                                                                       tgggtggTGTCGTGTACCACGGGCGGCGCTTCGGCTTCGTTGCTCTACCCTGCCGCTTCCTGCAAATAGGCGATAAATGCCTCTTCCAGACTGGCGGCTCCGCGTTTCTCAACCAGTTCCTGCGGTGTACCGCTGGCAAg                                  <  1:330612/140‑1 (MQ=255)                                                                                                tCGTGTACCACGGGCGGCGCTTCGGCTTCGTTGCTCTACCCTGCCGCTTCCTGCAAATAGGCGATAAATGCCTCTTCCAGACTGGCGGCTCCGCGTTTCTCAACCAGTTCCTGCGGTGTACCGCTGGCAAGCACTTTTcc                         >  1:676310/1‑140 (MQ=255)                                                                                                tCGTGTACCACGGGCGGCGCTTCGGCTTCGTTGCTCTACCCTGCCGCTTCCTGCAAATAGGCGATAAATGCCTCTTCCAGACTGGCGGCTCCGCGTTTCTCAACCAGTTCCTGCGGTGTACCGCTGGCAAGCACTTTTcc                         >  2:855367/1‑140 (MQ=255)                                                                                                 cGTGTACCACGGGCGGCGCTTCGGCTTCGTTGCTCTACCCTGCCGCTTCCTGCAAATAGGCGATAAATGCCTCTTCCAGACTGGCGGCTCCGCGTTTCTCAACCAGTTCCTGCGGTGTACCGCTGGCAAGCACTTTTCCg                        <  1:221780/140‑1 (MQ=255)                                                                                                 cGTGTACCACGGGCGGCGCTTCGGCTTCGTTGCTCTACCCTGCCGCTTCCTGCAAATAGGCGATAAATGCCTCTTCCAGACTGGCGGCTCCGCGTTTCTCAACCAGTTCCTGCGGTGTACCGCTGGCAAGCACTTTTCCg                        <  1:557242/140‑1 (MQ=255)                                                                                                   tgtACCACGGGCGGCGCTTCGGCTTCGTTGCTCTACCCTGCCGCTTCCTGCAAATAGGCGATAAATGCCTCTTCCAGACTGGCGGCTCCGCGTTTCTCAACCAGTTCCTGCGGTGTACCGCTGGCAAGCACTTTTCCGGc                      >  2:455283/1‑140 (MQ=255)                                                                                                   tgtACCACGGGCGGCGCTTCGGCTTCGTTGCTCTACCCTGCCGCTTCCTGCAAATAGGCGATAAATGCCTCTTCCAGACTGGCGGCTCCGCGTTTCTCAACCAGTTCCTGCGGTGTACCGCTGGCAAGCACTTTTCCGGc                      >  1:26418/1‑140 (MQ=255)                                                                                                   tgtACCACGGGCGGCGCTTCGGCTTCGTTGCTCTACCCTGCCGCTTCCTGCAAATAGGCGATAAATGCCTCTTCCAGACTGGCGGCTCCGCGTTTCTCAACCAGTTCCTGCGGTGTACCACTGGCAAGCACTTTTCCGGc                      <  1:219103/140‑1 (MQ=255)                                                                                                      aCCACGGGCGGCGCTTCGGCTTCGTTGCTCTACCCTGCCGCTTCCTGCAAATAGGCGATAAATGCCTCTTCCAGACTGGCGGCTCCGCGTTTCTCAACCAGTTCCTGCGGTGTACCGCTGGCAAGCACTTTTCCGGCGTg                   <  2:676310/140‑1 (MQ=255)                                                                                                        cACGGGCGGCGCTTCGGCTTCGTTGCTCTACCCTGCCGCTTCCTGCAAATAGGCGATAAATGCCTCTTCCAGACTGGCGGCTCCGCGTTTCTCAACCAGTTCCTGCGGTGTACCGCTGGCAAGCACTTTTCCg                        >  1:524043/1‑133 (MQ=255)                                                                                                        cACGGGCGGCGCTTCGGCTTCGTTGCTCTACCCTGCCGCTTCCTGCAAATAGGCGATAAATGCCTCTTCCAGACTGGCGGCTCCGCGTTTCTCAACCAGTTCCTGCGGTGTACCGCTGGCAAGCACTTTTCCg                        <  2:524043/133‑1 (MQ=255)                                                                                                        cACGGGCGGCGCTTCGGCTTCGTTGCTCTACCCTGCCGCTTCCTGCAAATAGGCGATAAATGCCTCTTCCAGACTGGCGGCTCCGCGTTTCTCAACCAGTTCCTGCGGTGTACCGCTGGCAAGCACTTTTCCGGCGTg                   >  2:948298/1‑138 (MQ=255)                                                                                                        cACGGGCGGCGCTTCGGCTTCGTTGCTCTACCCTGCCGCTTCCTGCAAATAGGCGATAAATGCCTCTTCCAGACTGGCGGCTCCGCGTTTCTCAACCAGTTCCTGCGGTGTACCGCTGGCAAGCACTTTTCCGGCGTg                   <  1:948298/138‑1 (MQ=255)                                                                                                        cACGGGCGGCGCTTCGGCTTCGTTGCTCTACCCTGCCGCTTCCTGCAAATAGGCGATAAATGCCTCTTCCAGACTGGCGGCTCCGCGTTTCTCAACCAGTTCCTGCGGTGTACCGCTGGCAAGCACTTTTCCGGCGTGCa                 <  1:758122/140‑1 (MQ=255)                                                                                                           gggCGGCGCTTCGGCTTCGTTGCTCTACCCTGCCGCTTCCTGCAAATAGGCGATAAATGCCTCTTCCAGACTGGCGGCTCCGCGTTTCTCAACCAGTTCCTGCGGTGTACCGCTGGCAAGCACTTTTCCGGCGTGCATCa              >  1:489655/1‑140 (MQ=255)                                                                                                                 cgcTTCGGCTTCGTTGCTCTACCCTGCCGCTTCCTGCAAATAGGCGATAAATGCCTCTTCCAGACTGGCGGCTCCGCGTTTCTCAACCAGTTCCTGCGGTGTACCGCTGGCAAGCACTTTTCCGGCGTGCATCAGTgaga        <  1:455283/140‑1 (MQ=255)                                                                                                                   cttcggcttcgTTGCTCTACCCTGCCGCTTCCTGCAAATAGGCGATAAATGCCTCTTCCAGACTGGCGGCTCCGCGTTTCTCAACCAGTTCCTGCGGTGTACCGCTGGCAAGCACTTTTCCGGCGTGCATCAGTGAGATg      <  2:489655/140‑1 (MQ=255)                                                                                                                       ggcttcgTTGCTCTACCCTGCCGCTTCCTGCAAATAGGCGATAAATGCCTCTTCCAGACTGGCGGCTCCGCGTTTCTCAACCAGTTCCTGCGGTGTACCGCTGGCAAGCACTTTTCCGGCGTGCATCAGTGAGATGCGGt  >  2:146891/1‑140 (MQ=255)                                                                                                                       ggcttcgTTGCTCTACCCTGCCGCTTCCTGCAAATAGGCGATAAATGCCTCTTCCAGACTGGCGGCTCCGCGTTTCTCAACCAGTTCCTGCGGTGTACCGCTGGCAAGCACTTTTCCGGCGTGCATCAGTGAGATGCGGt  >  2:670973/1‑140 (MQ=255)                                                                                                                                     |                                                                                                                              AACGTACTGGATCGCGTCGCAGTTCCAGCGCTTCGCGGCGGCTGTAGCTAAACAGACGGCGCAGGCTAAATCCCTGACGCGGCGCGTGGGTGGTGTCGTGTACCACGGGCGGCGCTTCGGCTTCGTTGCTCTGCCCTGCCGCTTCCTGCAAATAGGCGATAAATGCCTCTTCCAGACTGGCGGCTCCGCGTTTCTCAACCAGTTCCTGCGGTGTACCGCTGGCAAGCACTTTTCCGGCGTGCATCAGTGAGATGCGGT  >  NC_000913/3627858‑3628115

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