Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,643,274	C→T	13.3%	intergenic (+59/‑64)	dbpA → / → yhiN	23S rRNA specific ATP‑dependent RNA helicase/FAD/NAD(P)‑binding domain‑containing oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,643,274	0	C	T	13.3%	112.8 / 7.8	45	intergenic (+59/‑64)	dbpA/yhiN	23S rRNA specific ATP‑dependent RNA helicase/FAD/NAD(P)‑binding domain‑containing oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (18/21);  new base T (4/2);  total (22/23)
Fisher's exact test for biased strand distribution p-value = 4.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.44e-01

CGTGGAGCGGGCGCTGGCCAAGCAGGCCATGCAGCGGCTGAACAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCAGCTCAGTGCCCGTAGCGGCCGCCG  >  NC_002947/5643128‑5643420                                                                                                                                                   |                                                                                                                                                   cGTGGAGCGGGCGCTGGCCAAGCAGGCCATGCAGCGGCTGAACAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCAt                                                                                                                                                  <  5:16781/149‑1 (MQ=255)   tGGAGCGGGCGCTGGCCAAGCAGGCCATGCAGCGGCTGAACAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCg                                                                                                                                                >  6:205899/1‑149 (MQ=255)            cgcTGGCCAAGCAGGCCATGCAGCGGCTGAACAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCgg                                                                                                                                       <  7:269152/149‑1 (MQ=255)            cgcTGGCCAAGCAGGCCATGCAGCGGCTGAACAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGATGGTGCATTCACTGCCGCATTCGCgg                                                                                                                                       <  5:58106/149‑1 (MQ=37)            cgagggCCAAGCAGGCCATGCAGCGGCTGAACAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCgg                                                                                                                                       <  6:166862/145‑1 (MQ=255)               tGGCCAAGCAGGCCATGCAGCGGCTGAACAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGt                                                                                                                                     >  5:290624/1‑148 (MQ=255)                   cAAGCAGGCCATGCAGCGGCTGAACAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAAcc                                                                                                                                <  7:253704/149‑1 (MQ=255)                   cAAGCAGGCCATGCAGCGGCTGAACAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAAcc                                                                                                                                >  3:60272/1‑149 (MQ=255)                          gCCATGCAGCGGCTGAACAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTcc                                                                                                                          >  3:217449/1‑148 (MQ=255)                                cAGCGGCTGAACAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAgg                                                                                                                    <  5:34025/148‑1 (MQ=255)                                    ggCTGAACAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCa                                                                                                               <  5:205899/149‑1 (MQ=255)                                          aCAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCAttt                                                                                                         <  6:290624/149‑1 (MQ=255)                                              cGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGtt                                                                                                                 <  1:49602/137‑1 (MQ=255)                                              cGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGtt                                                                                                                 >  2:49602/1‑137 (MQ=255)                                               ggCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGAGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAAt                                                                                                    <  6:309246/149‑1 (MQ=255)                                                gCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATc                                                                                                   <  8:85097/149‑1 (MQ=255)                                                 cAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTCGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCt                                                                                                  <  6:70058/149‑1 (MQ=255)                                                    gATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAAt                                                                                                    >  8:209184/1‑144 (MQ=255)                                                    gATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAAt                                                                                                    <  7:209184/144‑1 (MQ=255)                                                        aaGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGtt                                                                                           <  3:283392/149‑1 (MQ=255)                                                                     tGAAAGTCCGCATTATCTAAAGCCAGGCTAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGAGGCAGTGAATGCACCATCTCATTCGCGGGTGAACCAGCTCCCACAGGGTTCAGCATATCAAAACTGAAGGTTAATACAGTGCaa                                                                               >  1:176103/1‑147 (MQ=38)                                                                                 ttATCTAAAGCCCGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGa                                                                  <  4:60272/149‑1 (MQ=255)                                                                                          gCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACGGAc                                                                      >  7:236992/1‑136 (MQ=255)                                                                                          gCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACGGAc                                                                      >  7:232640/1‑136 (MQ=255)                                                                                          gCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGAc                                                                      <  8:236992/136‑1 (MQ=255)                                                                                            cAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTggc                                                       >  3:74870/1‑149 (MQ=255)                                                                                             aGGCAAGGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTggcg                                                      >  8:252606/1‑149 (MQ=255)                                                                                             aGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTc                                                                                                                <  4:13850/91‑1 (MQ=255)                                                                                             aGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTc                                                                                                                >  3:13850/1‑91 (MQ=255)                                                                                                        ccTGTGGGAGCGGGTTCACCCGCGAATGCGATGGTGCATTCACTGCCGCATTCGCGGGTAAACCCGCTCCTACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCgccgcc                                           >  5:224401/1‑149 (MQ=18)                                                                                                              ggAGCGGGTTCACCCGCGAATGCGATGGTGCATTCACTGCCGCATTCGCGGGTAAACCCGCTCCTACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCt                                      >  8:217043/1‑148 (MQ=18)                                                                                                               gAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGa                                    >  1:134487/1‑149 (MQ=255)                                                                                                                 gCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGAtg                                  >  1:66868/1‑149 (MQ=255)                                                                                                                       ttACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGAt                                   >  3:131970/1‑142 (MQ=255)                                                                                                                       ttACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGAt                                   <  4:131970/142‑1 (MQ=255)                                                                                                                          cccGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCAg                         <  2:134487/149‑1 (MQ=255)                                                                                                                           ccGCGAATGCGATGGTGCATTCACTGCCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGccc                           >  7:211521/1‑146 (MQ=37)                                                                                                                           ccGCGAATGCGATGGTGCATTCACTGCCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGccc                           <  8:211521/146‑1 (MQ=37)                                                                                                                             gcgAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCAGCt                       >  1:150424/1‑148 (MQ=255)                                                                                                                             gcgAATGCGATGGTGCATTCACTGCCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTAATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGTCCTGATTTGCGCCCAGCt                       >  7:201455/1‑148 (MQ=25)                                                                                                                              cgAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCAGCTCa                     <  7:252606/149‑1 (MQ=255)                                                                                                                              cgAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCAGCTCa                     <  8:25268/149‑1 (MQ=255)                                                                                                                                      gCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCAGCTCAGTGCCCg              >  5:7317/1‑148 (MQ=255)                                                                                                                                       cAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCAGCTCAGTGCCCGTa            <  4:217449/149‑1 (MQ=255)                                                                                                                                                 aCCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCAGCTCAGTGCCCGTAGCGGccgccg  >  2:221417/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                   CGTGGAGCGGGCGCTGGCCAAGCAGGCCATGCAGCGGCTGAACAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCAGCTCAGTGCCCGTAGCGGCCGCCG  >  NC_002947/5643128‑5643420

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