Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,938,662	T→A	11.5%	L408Q (CTG→CAG)	rbbA →	ribosome‑associated ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,938,662	0	T	A	11.5%	137.0 / 7.5	52	L408Q (CTG→CAG)	rbbA	ribosome‑associated ATPase
Reads supporting (aligned to +/- strand):  ref base T (22/24);  new base A (2/4);  total (24/28)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGCAGAACCTGGCCCTGCATGCCCGCCTGTTCGACTTGCCCAAAGCCGAGAGCGCCCAGC‑GCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGTCCCGCGAGCAGGGCGTGACCATTTTCCTCTCCA  >  NC_002947/5938514‑5938808                                                                                                                                                      |                                                                                                                                                   cGGCAGAACCTGGCCCTGCATGCCCGCCTGTTCGACTTGCCCAAAGCCGAGAGCGCCCAGC‑GCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCt                                                                                                                                                    <  7:259961/149‑1 (MQ=255)       aaCCTGGCCCTGCATGCCCGCCTGTTCGACTTGCCCAAAGCCGAGAGCGCCCAGC‑GCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCt                                                                                                                                              >  8:10090/1‑149 (MQ=255)         ccTGGCCCTGCATGCCCGCCTGTTCGACTTGCCCAAAGCCGAGAGCGCCCAGCGGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTg                                                                                                                                             >  1:262982/1‑149 (MQ=255)          cTGGCCCTGCATGCCCGCCTGTTCGACTTGCCCAAAGCCGAGAGCGCCCAGC‑GCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGACAGGCGCTGGc                                                                                                                                           <  3:135981/149‑1 (MQ=255)             gCCCTGCATGCCCGCCTGTTCGACTTGCCCAAAGCCGAGAGCGCCCAGC‑GCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTCTCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGACAGGCGCTGGcg                                                                                                                                          <  5:246304/147‑2 (MQ=255)               ccTGCATGCCCGCCTGTTCGACTTGCCCAAAGCCGAGAGCGCCCAGC‑GCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTgg                                                                                                                                      <  7:63740/149‑1 (MQ=255)                      gCCCGCCTGTTCGACTTGCCCAAAGCCGAGAGCGCCCAGC‑GCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCt                                                                                                                               >  5:338299/1‑149 (MQ=255)                                  gACTTGCCCAAAGCCGAGAGCGCCCAGC‑GCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGTAGTGGCGGTTCt                                                                                                                               <  5:181307/137‑1 (MQ=255)                                  gACTTGCCCAAAGCCGAGAGCGCCCAGC‑GCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGTAGTGGCGGTTCt                                                                                                                               >  6:181307/1‑137 (MQ=255)                                   aCTTGCCCAAAGCCGAGAGCGCCCAGC‑GCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGaa                                                                                                                  <  2:156025/149‑1 (MQ=255)                                   aCTTGCCCAAAGCCGAGAGCGCCCAGC‑GCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGaa                                                                                                                  <  7:52694/149‑1 (MQ=255)                                        cccAAAGCCGAGAGCGCCCAGC‑GCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCtt                                                                                                             >  1:97708/1‑149 (MQ=255)                                        cccAAAGCCGAGAGCGCCCAGC‑GCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCt                                                                                                              <  5:122618/148‑1 (MQ=255)                                          cAAAGCCGAGAGCGCCCAGC‑GCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGa                                                                                                           <  3:69639/149‑1 (MQ=255)                                                 gagaGCGCCCAGC‑GCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTcga                                                                                                    >  3:64558/1‑149 (MQ=255)                                                 gagaGCGCCCAGC‑GCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTcga                                                                                                    <  7:10090/149‑1 (MQ=255)                                                        cccAGC‑GCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTgg                                                                                                                                            <  2:945/102‑1 (MQ=255)                                                                 aTCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTc                                                                                                      >  4:193781/1‑132 (MQ=255)                                                                 aTCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTc                                                                                                      <  3:193781/132‑1 (MQ=255)                                                                 aTCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACAAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCgg                                                                                     >  8:244258/1‑149 (MQ=255)                                                                  tCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGACGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGc                                                                                    <  8:318011/149‑1 (MQ=255)                                                                     aTGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTc                                                                                 >  5:100793/1‑149 (MQ=255)                                                                       gAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCgg                                                                                                                    <  8:267933/112‑1 (MQ=255)                                                                       gAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCgg                                                                                                                    >  7:267933/1‑112 (MQ=255)                                                                        aGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGAc                                                                              <  2:262982/149‑1 (MQ=255)                                                                             aTCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCccgc                                                                         >  3:166674/1‑149 (MQ=255)                                                                                    gcTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGccgcccgc                                                                  <  2:280346/149‑1 (MQ=255)                                                                                     cTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTcc                                                                                       <  3:266913/127‑1 (MQ=255)                                                                                     cTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTCCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTcc                                                                                       >  4:266913/1‑127 (MQ=255)                                                                                      ttCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGAAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGcga                                                                >  4:285909/1‑149 (MQ=255)                                                                                        cGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGcga                                                                <  7:244258/147‑1 (MQ=255)                                                                                          aCCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCgacgac                                                            >  6:210309/1‑149 (MQ=255)                                                                                                         ccGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCctgctg                                             <  5:210309/149‑1 (MQ=255)                                                                                                         ccGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCctgctg                                             >  5:28709/1‑149 (MQ=255)                                                                                                         ccGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCctgctg                                             >  7:280425/1‑149 (MQ=255)                                                                                                          cGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTgg                                            <  4:64558/149‑1 (MQ=255)                                                                                                           gACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTgg                                            >  1:233491/1‑148 (MQ=255)                                                                                                            aCCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTgg                                            >  3:311540/1‑147 (MQ=255)                                                                                                                aCCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGaa                                                                                                <  5:295720/91‑1 (MQ=255)                                                                                                                aCCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGaa                                                                                                >  6:295720/1‑91 (MQ=255)                                                                                                                aCCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAAc                                      <  6:100793/149‑1 (MQ=255)                                                                                                                aCCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAAc                                      >  5:97490/1‑149 (MQ=255)                                                                                                                aCCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAAc                                      <  3:215908/149‑1 (MQ=255)                                                                                                                            ccTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGTCCCGCGAGc                          >  2:342931/1‑149 (MQ=255)                                                                                                                                  ccTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGTCCCGCGAGCAGGGCg                    <  8:280425/149‑1 (MQ=255)                                                                                                                                  ccTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGTCCCGCGAGCAGGGCg                    <  2:233491/149‑1 (MQ=255)                                                                                                                                   cTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGTCCCGCGAGCAGGGCGt                   >  1:14466/1‑149 (MQ=255)                                                                                                                                              ccAGCGACAGGCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGTCCcgc                              <  3:117102/127‑1 (MQ=255)                                                                                                                                              ccAGCGACAGGCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGTCCcgc                              >  4:117102/1‑127 (MQ=255)                                                                                                                                               cAGCGACAGGCGCTGGCGGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGTCCCGCGAGCAGGGCGTGa                 <  8:135983/139‑1 (MQ=255)                                                                                                                                               cAGCGACAGGCGCTGGCGGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGTCCCGCGAGCAGGGCGTGa                 >  7:135983/1‑139 (MQ=255)                                                                                                                                                  cgcCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGTCCCGCGAGCAGGGCGTGACCATTTTCCtctc    <  6:28709/149‑1 (MQ=255)                                                                                                                                                    ccTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGTCCCGCGAGCAGGGCGTGACCATTTTCCTCTcca  <  2:14466/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                   CGGCAGAACCTGGCCCTGCATGCCCGCCTGTTCGACTTGCCCAAAGCCGAGAGCGCCCAGC‑GCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGTCCCGCGAGCAGGGCGTGACCATTTTCCTCTCCA  >  NC_002947/5938514‑5938808

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