Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	5,674,751	+G	intergenic (+113/+27)	rhlE‑II → / ← yceI	ATP‑dependent RNA helicase/polyisoprene‑binding acidic stress response factor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,674,750	1	.	G	95.5%	135.0 / ‑2.1	44	intergenic (+112/+28)	rhlE‑II/yceI	ATP‑dependent RNA helicase/polyisoprene‑binding acidic stress response factor
Reads supporting (aligned to +/- strand):  ref base . (0/0);  major base G (20/22);  minor base A (2/0);  total (22/22)
Fisher's exact test for biased strand distribution p-value = 4.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.75e-01

GCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGC‑GCGCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAA  >  NC_002947/5674608‑5674888                                                                                                                                                |                                                                                                                                           gcgcACAGCGCTGCTGACCCGCAAGTGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAgc‑ggcc                                                                                                                                        <  6:114953/147‑4 (MQ=255)       aGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGCCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCgg                                                                                                                                  <  4:103812/148‑1 (MQ=255)         cgcTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGACTTTCTTGATCTCTCCAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGtt                                                                                                                                <  5:108245/148‑1 (MQ=37)            tgctgACCCGCAAGGGCTGATAGACGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACACGGCCGCTCCTGCAGAGAATTTCTTGCTCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTc                                                                                                                             <  1:114834/148‑1 (MQ=37)              ctgACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGCGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAg                                                                                                                           >  1:90576/1‑148 (MQ=255)               tgACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGc                                                                                                                          >  1:140811/1‑148 (MQ=255)                  cccGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCt                                                                                                                         <  3:78734/146‑1 (MQ=255)                         gggCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGaaa                                                                                                                 >  5:122837/1‑147 (MQ=255)                           gCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTa                                                                                                              >  3:60511/1‑148 (MQ=255)                               aTAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTc                                                                                                          >  5:109840/1‑148 (MQ=255)                                        ccGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTcaca                                                                                                 >  8:93627/1‑148 (MQ=255)                                         cGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCacac                                                                                                >  6:57946/1‑148 (MQ=255)                                                    ctTCGCAGCCCTATCACGACATAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTg                                                                                     >  6:140586/1‑148 (MQ=255)                                                     ttCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGCTCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCGCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCCCCCCTTCGACACTg                                                                                     >  7:12645/1‑147 (MQ=37)                                                        gCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGata                                                                                  <  3:165609/147‑1 (MQ=255)                                                                   aCGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATAtccaggtccagg                                                                      >  2:73253/1‑148 (MQ=255)                                                                   aCGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATAtccaggtccagg                                                                      >  2:68900/1‑148 (MQ=255)                                                                      acacAAGCCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTc                                                                   <  8:68106/148‑1 (MQ=255)                                                                        acaAGGCCGCTCCTGCAGAGAATTTCTTGATCTCAACAGGAGCGGACTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATAACTTCTGCTTCACACCTTCGACACTGATATCCAGGTACAGGGTCTg                                                                 >  6:49046/1‑148 (MQ=255)                                                                          aaGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACCGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGc                                                                >  3:68584/1‑147 (MQ=255)                                                                               cGCTGCTGCAGAGAATTTCTTGATCTCTACAGGTGAGGGCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTa                                                          <  6:105611/148‑1 (MQ=255)                                                                                gCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAg                                                         <  4:127708/148‑1 (MQ=255)                                                                                     tGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCAc                                                      <  8:95526/146‑1 (MQ=255)                                                                                      gCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCCCACCTTCGACACTGATATCCCGGTCAAGGGTCTGCGAAGTAGCACCCg                                                   >  6:167681/1‑148 (MQ=38)                                                                                         gagaATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGc                                                >  2:106113/1‑148 (MQ=255)                                                                                         gagaATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGc                                                >  2:116889/1‑148 (MQ=255)                                                                                         gagaATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCCTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGc                                                >  4:48190/1‑148 (MQ=255)                                                                                            aaTTTCTTGATCTCTACAGGAGCTGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCtt                                             <  7:82920/148‑1 (MQ=255)                                                                                                    gATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGcc                                      >  1:132083/1‑147 (MQ=255)                                                                                                           aCAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGGCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAg                                  <  3:27528/144‑1 (MQ=255)                                                                                                             aGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCg                               <  6:109840/145‑1 (MQ=255)                                                                                                              ggAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCgttca                           <  8:88081/148‑1 (MQ=255)                                                                                                                aGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCgttcagg                         <  3:87925/148‑1 (MQ=255)                                                                                                                          tgtCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCgt                 <  5:48713/146‑1 (MQ=255)                                                                                                                           gtCGTGTTAGGGCCGCAAAGCGGCCCCCTCGTTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCgtggt              <  2:113623/148‑1 (MQ=255)                                                                                                                           gtCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCgtggt              >  6:102252/1‑148 (MQ=255)                                                                                                                           gtCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCgtgg               <  2:132083/147‑1 (MQ=255)                                                                                                                               tGATAGGGCTGCGAAGCAGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCgtggtggtg          >  6:121850/1‑148 (MQ=255)                                                                                                                               tGATAGGGCTGCGAAGCAGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCgtggtggtg          >  6:119516/1‑148 (MQ=255)                                                                                                                               tGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCgtggtggt           <  2:90576/147‑1 (MQ=255)                                                                                                                                gATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGtggtggtgg         <  8:10071/148‑1 (MQ=255)                                                                                                                                  tAGGGCCGCAAAGCGGCCCCCTCGGTTTTCGGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCg       >  2:82506/1‑148 (MQ=255)                                                                                                                                  tAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGc        <  3:54943/147‑1 (MQ=255)                                                                                                                                      gCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGa   <  3:34900/148‑1 (MQ=255)                                                                                                                                       ccGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGaa  <  2:101830/148‑1 (MQ=255)                                                                                                                                                |                                                                                                                                           GCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGC‑GCGCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAA  >  NC_002947/5674608‑5674888

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