Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,629,478	G→A	H21Y (CAT→TAT)	rbbA ←	ribosome‑associated ATPase: ATP‑binding protein/ATP‑binding membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,629,478	0	G	A	100.0%	176.6 / NA	55	H21Y (CAT→TAT)	rbbA	ribosome‑associated ATPase: ATP‑binding protein/ATP‑binding membrane protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (26/29);  total (26/29)

CGGGCACCGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAA  >  NC_000913/3629354‑3629613                                                                                                                             |                                                                                                                                        cgGGCACCGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATACTGGCTCAcgcccg                                                                                                                           >  1:872859/1‑139 (MQ=255)   ggCACCGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCGATAATACTGGCTCAcgaccgcc                                                                                                                         <  1:1653959/139‑1 (MQ=255)                   gACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATACTGGCTCACGCCCGCCAGTTGCGCGACAGGcg                                                                                                         <  2:811639/139‑1 (MQ=255)                           gCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATACTGGCTCACGCCCGCCAGTTGCGCGACAGGcgggacggga                                                                                                 <  1:907208/139‑1 (MQ=255)                                aCTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATACTGGCTCACGCCCGCCAGTTGCGCGACAGGcgggacgg                                                                                                   <  2:143256/132‑1 (MQ=255)                                aCTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATACTGGCTCACGCCCGCCAGTTGCGCGACAGGcgggacgg                                                                                                   >  1:143256/1‑132 (MQ=255)                                           ccGTCCGGGCCAATCAGCCCGACCATACAGCGGGCTGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATACTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCg                                                                                 <  2:1010957/139‑1 (MQ=255)                                           ccGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATACTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCg                                                                                 >  2:148047/1‑139 (MQ=255)                                            cGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATACTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAAcc                                                                                              >  1:787794/1‑125 (MQ=255)                                            cGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATACTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAAcc                                                                                              <  2:787794/125‑1 (MQ=255)                                              tCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATACTGGCTCa                                                                                                                                 >  2:636360/1‑88 (MQ=255)                                              tCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATACTGGCTCa                                                                                                                                 <  1:636360/88‑1 (MQ=255)                                              tCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATACTGGCTCACGCCCGCCAGTTGCGCGACAGGcggg                                                                                                       >  1:1300782/1‑114 (MQ=255)                                              tCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATACTGGCTCACGCCCGCCAGTTGCGCGACAGGcggg                                                                                                       <  2:1300782/114‑1 (MQ=255)                                              tCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATACTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCa                                                                              >  2:1228768/1‑139 (MQ=255)                                                    ccAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATACTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCgg                                                                        <  2:872859/139‑1 (MQ=255)                                                                aCCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATACTGGCTc                                                                                                                                  >  2:1473382/1‑69 (MQ=255)                                                                aCCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATACTGGCTc                                                                                                                                  <  1:1473382/69‑1 (MQ=255)                                                                   aTACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATACTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAg                                                         <  1:148047/139‑1 (MQ=255)                                                                    tACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATACTGGCTCACGCCCGCCAGTTGCGCGa                                                                                                               >  1:1519132/1‑84 (MQ=255)                                                                    tACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATACTGGCTCACGCCCGCCAGTTGCGCGa                                                                                                               <  2:1519132/84‑1 (MQ=255)                                                                     aCAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATACTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGa                                                          >  1:1592736/1‑136 (MQ=255)                                                                     aCAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATACTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGa                                                          <  2:1592736/136‑1 (MQ=255)                                                                        gCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATACTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGgtcgt                                                    <  1:1228768/139‑1 (MQ=255)                                                                           ggCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATACTGGCTCACGCCCGCCAGt                                                                                                                      <  2:1030866/70‑1 (MQ=255)                                                                           ggCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATACTGGCTCACGCCCGCCAGt                                                                                                                      >  1:1030866/1‑70 (MQ=255)                                                                           ggCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATACTGGCTCACGCCCGCCAGTTGCGCGAc                                                                                                              >  1:454843/1‑78 (MQ=255)                                                                           ggCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATACTGGCTCACGCCCGCCAGTTGCGCGAc                                                                                                              <  2:454843/78‑1 (MQ=255)                                                                           ggCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATACTGGCTCACGCCCGCCAGTTGCGCGACAGGcgggacgg                                                                                                   <  1:921770/89‑1 (MQ=255)                                                                           ggCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATACTGGCTCACGCCCGCCAGTTGCGCGACAGGcgggacgg                                                                                                   >  2:921770/1‑89 (MQ=255)                                                                           ggCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATACTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTc                                                                                         <  1:1104397/99‑1 (MQ=255)                                                                           ggCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATACTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTc                                                                                         >  2:1104397/1‑99 (MQ=255)                                                                           ggCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATACTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGgtcgt                                                    <  1:1372651/136‑1 (MQ=255)                                                                           ggCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATACTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGgtcgt                                                    >  2:1372651/1‑136 (MQ=255)                                                                           ggCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATACTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAg                                                 >  2:1425818/1‑139 (MQ=255)                                                                              cGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATACTGGCTCa                                                                                                                                 <  2:1265827/56‑1 (MQ=255)                                                                              cGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATACTGGCTCa                                                                                                                                 >  1:1265827/1‑56 (MQ=255)                                                                               ggAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATACTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCa                                             <  1:238222/139‑1 (MQ=255)                                                                                  atatCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATACTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACgg                                          <  2:1422772/139‑1 (MQ=255)                                                                                             gATATTGTTCAGCGCAACGGTTTTTCCATAATACTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCttcat                               >  2:133714/1‑139 (MQ=255)                                                                                                  tGTTCAGCGCAACGGTTTTTCCATAATACTGGCTCACGCCCGCCAGTTGCGCGACAGGcggga                                                                                                      >  2:1454993/1‑63 (MQ=255)                                                                                                  tGTTCAGCGCAACGGTTTTTCCATAATACTGGCTCACGCCCGCCAGTTGCGCGACAGGcggga                                                                                                      <  1:1454993/63‑1 (MQ=255)                                                                                                       aGCGCAACGGTTTTTCCATAATACTGGCTCACGCCCGCCAGTTGCGCGACAGGcg                                                                                                         <  1:637487/55‑1 (MQ=255)                                                                                                       aGCGCAACGGTTTTTCCATAATACTGGCTCACGCCCGCCAGTTGCGCGACAGGcg                                                                                                         >  2:637487/1‑55 (MQ=255)                                                                                                         cgcAACGGTTTTTCCATAATACTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGc                                                                       >  1:1281849/1‑87 (MQ=255)                                                                                                         cgcAACGGTTTTTCCATAATACTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGc                                                                       <  2:1281849/87‑1 (MQ=255)                                                                                                         cgcAACGGTTTTTCCATAATACTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGgtcgt                                                    <  1:503909/106‑1 (MQ=255)                                                                                                         cgcAACGGTTTTTCCATAATACTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGgtcgt                                                    >  2:503909/1‑106 (MQ=255)                                                                                                          gcAACGGTTTTTCCATAATACTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCa                                             <  1:374191/112‑1 (MQ=255)                                                                                                          gcAACGGTTTTTCCATAATACTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCa                                             >  2:374191/1‑112 (MQ=255)                                                                                                           cAACGGTTTTTCCATAATACTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCAc                 >  2:1010235/1‑139 (MQ=255)                                                                                                                         tAATACTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCaa   >  1:522900/1‑139 (MQ=255)                                                                                                                         tAATACTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGC‑ACGCCCGGCaaa  <  1:133714/139‑1 (MQ=255)                                                                                                                          aaTACTGGCTCACGCCCGCCAGTTGCGCGACAGGcgggacgg                                                                                                   >  1:521494/1‑42 (MQ=255)                                                                                                                          aaTACTGGCTCACGCCCGCCAGTTGCGCGACAGGcgggacgg                                                                                                   <  2:521494/42‑1 (MQ=255)                                                                                                                             |                                                                                                                                        CGGGCACCGGAAATCAACGACAACAAGCTCGACTTCCCGACGCCGTCCGGGCCAATCAGCCCGACCATACAGCGGGCCGGAATATCGAGAGTGATATTGTTCAGCGCAACGGTTTTTCCATAATGCTGGCTCACGCCCGCCAGTTGCGCGACAGGCGGGACGGGAACCAGTTCCAGATGCGTCATTGCGGCAACCTCACCACGAGGTCGTCAGGCCACGGAAGTTCTTCATTCACCCGCACCCACGCTACGCCCGGCAAA  >  NC_000913/3629354‑3629613

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