Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,810,336	C→A	9.1%	intergenic (+97/‑72)	eno → / → ftsB	enolase/cell division protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,810,336	0	C	A	9.1%	124.9 / 2.2	48	intergenic (+97/‑72)	eno/ftsB	enolase/cell division protein
Reads supporting (aligned to +/- strand):  ref base C (16/27);  new base A (3/2);  total (19/29)
Fisher's exact test for biased strand distribution p-value = 3.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.94e-01

AACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGC‑CTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGCTGTTGGGTGGCCTGCAGTACCGCCTGTGGGTGGGT  >  NC_002947/1810190‑1810482                                                                                                                                                    |                                                                                                                                                   aaCTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGC‑CTCgg                                                                                                                                                  >  6:66842/1‑149 (MQ=255) aaCCGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATTGTAAAAATGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGC‑CTCgg                                                                                                                                                  <  5:72711/149‑1 (MQ=255)          ccAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGC‑CTCGGCCAGGCTTc                                                                                                                                         <  5:66842/149‑1 (MQ=255)          ccAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGC‑CTCGGCCAGGCTTc                                                                                                                                         <  3:178873/149‑1 (MQ=255)           cAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGC‑CTCGGCCAGGCTTCg                                                                                                                                        <  8:266820/149‑1 (MQ=255)            aaaGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGC‑CGAGGCCAGGCTTCGt                                                                                                                                       <  3:217698/149‑1 (MQ=255)            aaaGCGGTTTACCGCGGTCGGGCCGAGTTTCGCGGCTAAGCACGAGATGGTAAAACGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGC‑CGAGGCCAGGCTTCGt                                                                                                                                       <  8:254139/149‑1 (MQ=255)             aaGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGC‑CTCGGCCAGGCTTCGTg                                                                                                                                      <  1:212225/149‑1 (MQ=255)              aGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCAGGC‑CTCGGCCAGGCATCGTgg                                                                                                                                     >  4:188999/1‑148 (MQ=255)                   tttACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGC‑CTCGGCCAGGCTTCGTGCTATTc                                                                                                                                <  1:229302/149‑1 (MQ=255)                        cgcgGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGC‑CTAGGCCAGCCTTCGTGCTATTGGAGGc                                                                                                                           >  4:238184/1‑149 (MQ=255)                          cgGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGC‑CTCGGCCAGGCTTCGTGCTATTCGAGGccc                                                                                                                         <  4:277984/149‑1 (MQ=255)                                tGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGC‑CTCGGCGAGGCTTCGTGCTATTCGAGGCCCCGAAg                                                                                                                    >  5:272090/1‑148 (MQ=255)                                  ccGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGTCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGC‑CTCTGCCAGGCTTCGAGCTATTCGAGGCCCCGAAGTaa                                                                                                                 <  3:238184/149‑1 (MQ=255)                                   cGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGC‑CTCGGCCAGGCGTCGTGCTATTGGAGGCCCCGAAGTaa                                                                                                                 >  3:147007/1‑148 (MQ=255)                                         tCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGC‑CTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTAc                                                                                                          >  8:186181/1‑149 (MQ=255)                                           gcgGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGC‑CTCGGGCAGGCTTCGTGCTATTCGa                                                                                                                              >  5:262292/1‑127 (MQ=255)                                           gcgGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGC‑CTCGGCCAGGCTTCGTGCTATTCGa                                                                                                                              <  6:262292/127‑1 (MQ=255)                                               cTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACTTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGC‑CTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCa                                                                                                    <  7:151565/149‑1 (MQ=255)                                                             gTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGC‑CTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCt                                                                                       <  7:235948/148‑1 (MQ=255)                                                                aaaaGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGC‑CTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGa                                                                                    <  3:188999/148‑1 (MQ=255)                                                                 taaGGCAGCAGCCGGAGTTGTAGGATCGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGC‑CTCGGCCTGGCTTCGTGCTATTCGAGGCCCCGAAGTAACaa                                                                                                              <  5:296692/120‑1 (MQ=255)                                                                    ggCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGC‑CTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACg                                                                               <  2:250677/149‑1 (MQ=255)                                                                     gcagcaGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGC‑CTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGt                                                                              >  7:180395/1‑149 (MQ=255)                                                                          agcCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCAGGC‑CTAGGCCAGGCTTAGTGCTATTAGAGGCCCCCAAATAACAAGTAAGCGGCAGCCTTTTTAACCTTGATAAACTGAttc                                                                         >  6:20360/1‑147 (MQ=255)                                                                               taGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGC‑CTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCttt                                                                                              <  6:231117/122‑1 (MQ=255)                                                                                  ttGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACCAAGCCTGGC‑CTCGGCCACGCTTCCTGCTATTATCGGCCCCGAACTAACAAGTACGCGGCAGCCTTTTTAAACTGGATAACATGATGCGCAgaac                                                                  >  6:237391/1‑145 (MQ=255)                                                                                       ggAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGC‑CTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTAtt                                                             <  5:20360/148‑1 (MQ=255)                                                                                       ggAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGC‑CTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTg                                                            >  8:215108/1‑149 (MQ=255)                                                                                        gAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGC‑CTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTgg                                                           <  7:215108/149‑1 (MQ=255)                                                                                                  aCTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGC‑CTCGGCCAGGCTTCGTGCTATTAGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCg                                                 >  3:113288/1‑149 (MQ=255)                                                                                                     gAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCTCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGt                                                >  1:63320/1‑148 (MQ=255)                                                                                                      aaCTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGC‑CTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCt                                             >  4:255407/1‑149 (MQ=255)                                                                                                            ccTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGC‑CTCGGCCAg                                                                                                                                              >  4:7982/1‑46 (MQ=255)                                                                                                            ccTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGC‑CTCGGCCAg                                                                                                                                              <  3:7982/46‑1 (MQ=255)                                                                                                               aCGTACTAACGGGTTTCTGTCGACGAAGCCTGGC‑CTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCa                                                                      >  1:125997/1‑115 (MQ=255)                                                                                                               aCGTACTAACGGGTTTCTGTCGACGAAGCCTGGC‑CTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCa                                                                      <  2:125997/115‑1 (MQ=255)                                                                                                                  tACTAACGGGTTTCTGTCGACGAAGCCTGGC‑CTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTg                                            <  4:59709/138‑1 (MQ=255)                                                                                                                  tACTAACGGGTTTCTGTCGACGAAGCCTGGC‑CTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTg                                            >  3:59709/1‑138 (MQ=255)                                                                                                                    cTAACGGGTTTCTGTCGACGAAGCCTGGC‑CTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGCTGTTgg                               <  3:255407/149‑1 (MQ=255)                                                                                                                          ggTTTCTGTCGACGAAGCCTGGC‑CTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGCTGTTGGGTGGcc                         <  8:180395/149‑1 (MQ=255)                                                                                                                                         aaccTGGC‑CGAGGCCAGGCATCGTGCTATTCGAGGCCACGAAGTAACAAGTACGCGGCAGCAATTTTAACCTGGATAACGTGATGCGAAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGCTGTTGGCTGGCCTGCAGTACCGCCtgt          >  1:295933/3‑149 (MQ=255)                                                                                                                                         aGCCTGGC‑CTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTTACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGCTGTTGGGTGGCCTGCAGTACCGCCtgt          <  1:191116/149‑1 (MQ=255)                                                                                                                                          gCCTGGC‑CTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGCTGTTGGGTGGCCTGCAGTACCGCCtgtg         >  3:251972/1‑149 (MQ=255)                                                                                                                                               gC‑CTCGTCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGTATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGCTGTTGGGTGGCCTGCAGTACCGCCtgtgttttg    <  4:147007/149‑6 (MQ=255)                                                                                                                                               gC‑CTCGGCCAGTCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGCTGTTGGGTGGCCTGCAGTACCGCCTgtgggtgg    <  5:237391/149‑1 (MQ=255)                                                                                                                                                c‑cTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGCTGTTGGGTGGCCTGCAGTACCGCCTGtgggtggg   <  4:251972/149‑1 (MQ=255)                                                                                                                                                  cTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGCTGTTGGGTGGCCTGCAGTACCGCCTGTgggtgggt  <  2:63320/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                   AACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTCGTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGC‑CTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGCTGTTGGGTGGCCTGCAGTACCGCCTGTGGGTGGGT  >  NC_002947/1810190‑1810482

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