Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,643,270	G→T	13.4%	intergenic (+55/‑68)	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,270	0	G	T	13.4%	115.3 / 10.2	45	intergenic (+55/‑68)	dbpA/yhiN	23S rRNA specific ATP‑dependent RNA helicase/FAD/NAD(P)‑binding domain‑containing oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (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 = 4.37e-01

GCCGTGGAGCGGGCGCTGGCCAAGCAGGCCATGCAGCGGCTGAACAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCAGCTCAGTGCCCGTA  >  NC_002947/5643126‑5643410                                                                                                                                                 |                                                                                                                                             gCCGTGGAGCGGGCGCTGGCCAAGCAGGCCATGCAGCGGCTGAACAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCAc                                                                                                                                           >  8:253704/1‑148 (MQ=255)   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)                                                                                                                                                 |                                                                                                                                             GCCGTGGAGCGGGCGCTGGCCAAGCAGGCCATGCAGCGGCTGAACAGCGGCAAGATCAAGGGTCGCGCCCTGAAAGTCCGCATTATCTAAAGCCAGGCAAAGCCCCCTGTAGGAGCGGGTTTACCCGCGAATGCGGCAGTGAATGCACCATCGCATTCGCGGGTGAACCCGCTCCCACAGGGTTCAGCATTTCAAATCTGAAGGTTCATACCGTGCACTCCACCGACGTGATCATCCTTGGCGCCGGCGCCGCCGGCCTGATGTGCGCCCAGCTCAGTGCCCGTA  >  NC_002947/5643126‑5643410

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