Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	5,674,753	G→C	intergenic (+115/+25)	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,753	0	G	C	100.0%	99.1 / NA	30	intergenic (+115/+25)	rhlE‑II/yceI	ATP‑dependent RNA helicase/polyisoprene‑binding acidic stress response factor
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (17/13);  total (17/13)

CCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGC‑G‑CGCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGC  >  NC_002947/5674626‑5674890                                                                                                                                  |                                                                                                                                          ccGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGCATTTCTTGATCTCTACTGGAGCGGCCTTGTGTAGTGATAGGGCCGCAAAGCGG‑CCCCCTAGGTTTTCAGCTGGc                                                                                                                        >  3:99675/1‑148 (MQ=37) ccGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CCCCCTCGGTTTTCAGCTg                                                                                                                          >  6:133486/1‑146 (MQ=255)      aGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CCCCCTCGGTTTTCAGCTGGCAGaaa                                                                                                                   <  3:168682/148‑1 (MQ=255)            gATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACtt                                                                                                             <  6:30457/148‑1 (MQ=255)                    tACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTca                                                                                                     <  5:53172/148‑1 (MQ=255)                                                      acaAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTg                                                                   >  6:153843/1‑148 (MQ=255)                                                          ggCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGaa                                                               <  4:142455/148‑1 (MQ=255)                                                             cGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTa                                                            >  5:193359/1‑148 (MQ=255)                                                                               cTTGATCTCTACAGGCGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGGCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGAt                                           <  2:150294/147‑1 (MQ=39)                                                                                ttGATCTCTACAGGAGCGGCCTTGTGTAGTGATAGGGCCGCAAAGCGG‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGc                                         >  6:167366/1‑148 (MQ=255)                                                                                   aTCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCg                                       <  7:100663/147‑1 (MQ=255)                                                                                      tctACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTc                                                                     >  7:176569/1‑114 (MQ=255)                                                                                      tctACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGt                                   >  2:199778/1‑148 (MQ=255)                                                                                       ctACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGC‑GTCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCt                                                                    <  8:176569/114‑1 (MQ=255)                                                                                         aCAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGaa                                     <  7:121435/143‑1 (MQ=255)                                                                                         aCAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGaa                                     <  2:50378/143‑1 (MQ=255)                                                                                         aCAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGCCCTTTGATGCCGAAGTcc                                 >  1:169524/1‑146 (MQ=255)                                                                                          cAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCg                                                     >  7:66384/1‑126 (MQ=255)                                                                                           aGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCATACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCgg                                                    <  8:66384/126‑1 (MQ=255)                                                                                                            gTGATAGGGCCGCAAAGCGG‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCgtggtggt             >  5:169495/1‑148 (MQ=255)                                                                                                              gATAGGGCCGCAAAGCGG‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGtggtggtgg           >  2:189313/1‑148 (MQ=255)                                                                                                               aTAGGGCCGCAAAGCGG‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGc          >  2:204235/1‑148 (MQ=255)                                                                                                                 aGGGCCGCAAAGCGG‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGt        <  2:169524/148‑1 (MQ=255)                                                                                                                  gggCCGCAAAGCGG‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGt        >  2:163061/1‑147 (MQ=255)                                                                                                                  gggCCGCAAAGCGG‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGt        >  2:126051/1‑147 (MQ=255)                                                                                                                   ggCCGCAAAGCGG‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGtt       >  2:56733/1‑147 (MQ=255)                                                                                                                      cGCAAAGCGG‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTATGCTTCACACCTTCGACACTGATATCCCGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGACGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAg   >  5:106841/1‑148 (MQ=255)                                                                                                                       gCAAAGCGG‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCATCGTGGTGGTGGCGTTGAAGc  >  8:154031/1‑148 (MQ=255)                                                                                                                       gCAAAGCGG‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGc  <  5:133486/148‑1 (MQ=255)                                                                                                                       gCAAAGCGG‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGc  <  6:193359/148‑1 (MQ=255)                                                                                                                                  |                                                                                                                                          CCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGC‑G‑CGCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGC  >  NC_002947/5674626‑5674890

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