Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,810,335	T→G	13.4%	intergenic (+96/‑73)	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,335	0	T	G	13.4%	97.3 / 10.8	45	intergenic (+96/‑73)	eno/ftsB	enolase/cell division protein
Reads supporting (aligned to +/- strand):  ref base T (21/18);  new base G (3/3);  total (24/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.05e-01

AGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGG‑CAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGCTGTTGGGTGGCCTGCAGTACCGCCTGTGGGTGGGT  >  NC_002947/1810187‑1810482                                                                                                                                                      |                                                                                                                                                     agCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCt                                                                                                                                                      <  3:221019/149‑1 (MQ=255) agCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCt                                                                                                                                                      >  4:221019/1‑149 (MQ=255)         ggCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCa                                                                                                                                                <  2:135254/147‑1 (MQ=255)              cAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCg                                                                                                                                         <  3:39864/149‑1 (MQ=255)                 aGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGc                                                                                                                                      <  7:248014/149‑1 (MQ=255)                     gTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGt                                                                                                                                        <  4:158057/143‑1 (MQ=255)                     gTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGt                                                                                                                                        >  3:158057/1‑143 (MQ=255)                          ccGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCGAGGCCAGGCt                                                                                                                                            >  4:148689/1‑134 (MQ=255)                          ccGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCGAGGCCAGGCt                                                                                                                                            <  3:148689/134‑1 (MQ=255)                                     ccGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTa                                                                                                                   >  5:306110/1‑148 (MQ=255)                                        aGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACaa                                                                                                               <  1:259174/149‑1 (MQ=255)                                                  cTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGG‑Ca                                                                                                    <  7:95134/149‑1 (MQ=255)                                                     aGCAAGAGATGGTAAAAAGGCAGCATCCGTAGTTGTAGGAACGTTC‑GTTCTGAACTTTCCTACGTACTAACGGGTTTCTGTCGAAGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGG‑CAGcc                                                                                                 <  4:163588/149‑1 (MQ=255)                                                      gCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGG‑CAGCCt                                                                                                >  2:252423/1‑149 (MQ=255)                                                      gCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCGAGGCCAGGCTTCGt                                                                                                                                        <  4:242810/110‑1 (MQ=255)                                                      gCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCGAGGCCAGGCTTCGt                                                                                                                                        >  3:242810/1‑110 (MQ=255)                                                      gCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCAGGGCTACGCCAGGCGAAGTGCTAATGGAGGCCCCGAAGTAAACAGTAAGCGG‑GAGCCt                                                                                                >  1:295860/1‑149 (MQ=255)                                                          gagaTGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGt                                                                                                                    <  5:21400/126‑1 (MQ=255)                                                                 tAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGG‑CAGCCTTTTTAACCTgg                                                                                     <  7:114701/149‑1 (MQ=255)                                                                 tAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGG‑CAGCCTTTTTAACCTg                                                                                      >  1:119538/1‑148 (MQ=255)                                                                    aaaGGCAGCAGCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCAGGCCTAGGCCAGGCTTCGTGGTTTTAGAGGCCCCGAACTAACAAGTAAGCGG‑CAGCCTTTTTaaacgtgat                                                                                   >  6:177250/1‑141 (MQ=255)                                                                        gcagcaGCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGG‑CAGCCTTTTTAACCTGGATAACGt                                                                              <  2:119538/149‑1 (MQ=255)                                                                          agcagcCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGG‑CAGCCTTTTTAACCTGGATAACGTGa                                                                            >  7:240829/1‑149 (MQ=255)                                                                                cGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGG‑CAGCCTTTTTAACCTGGATAACGTGATGCGCa                                                                      <  7:139307/149‑1 (MQ=255)                                                                                       gTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGGCCCGAAGTAACAAGTACGCGG‑CAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTa                                                               >  3:285754/1‑149 (MQ=255)                                                                                        tAGGAAAGTTCAGTAC‑GAACGTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGG‑CAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTAt                                                              >  8:218142/1‑149 (MQ=255)                                                                                          ggAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGG‑CAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTg                                                            <  8:213155/149‑1 (MQ=255)                                                                                          ggAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTAAGCGG‑CCGCCTTTTTAACCTGGATAACGTGATGCGCAGTcccgaatg                                                            >  6:294310/1‑144 (MQ=255)                                                                                                      aCTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCAGGCCTCGGCCAGGCTTCGTGGTATTCGAGGCCCCGAAGTAACAAATAAGCGG‑GAGCCGTTTTTACCTGGATAACGTGATTCGCAGTCCCGATTGgttgttgagcg                                                 >  4:54612/1‑144 (MQ=255)                                                                                                          aaCTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGACAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGtt                                               <  1:252423/148‑1 (MQ=255)                                                                                                          aaCTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGG‑CAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCt                                             >  7:15836/1‑149 (MQ=255)                                                                                                                     gTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGG‑CAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGCTGt                                  >  8:101863/1‑149 (MQ=255)                                                                                                                      tACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGG‑CAGCCTTTTTAACCTGGATAACGTGATGCACAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGCTGtt                                 >  8:86592/1‑149 (MQ=255)                                                                                                                           aCGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGGGCTATTAGAGGCCCCGAAGTAACAAGTACGCGG‑CAGCCTTTTTAACCGGGATAACGTGATGCGCAGTCCCGATTGGTTGTTCCTCGTTCTGCTCTTGCTGTTGGGTg                            >  3:136677/1‑149 (MQ=255)                                                                                                                                  tCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGG‑CAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGCTGTTGGGTGGCCTGc                      <  5:294310/148‑1 (MQ=255)                                                                                                                                   cTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGG‑CAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGCTGTTGGGTGGCCTGCa                     <  6:306110/148‑1 (MQ=255)                                                                                                                                    tGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGGTATTCGAGGCCCCGAAGTAACAAGTACGCGG‑GAGCCGTTTTAACCGGGATAACGTGCTGCGCCGTCCCTATTGGTTGTTCCTGGTTGTGCTCTTGCGGTTGGGTGGCCTGGAGt                   >  4:101572/1‑149 (MQ=255)                                                                                                                                     gTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGG‑CAGCCTTTTTAAACTGGATAACGTGATGCGCAGTCCCTATTGGGTGTTCCTCGTTCTGCTCTTGCTGTTGGGTGGCCTGCAGt                   >  3:184471/1‑148 (MQ=255)                                                                                                                                         acgaAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGG‑CAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGCTGTTGGGTGGCCTGCAGTACCGc              <  5:177250/149‑1 (MQ=255)                                                                                                                                           gaAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGG‑CAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGCTGTTGGGTGGCCTGCAGTACCGCCt            >  8:114824/1‑149 (MQ=255)                                                                                                                                           gaAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGG‑CAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGCTGTTGGGTGGCCTGCAGTACCGCCt            >  3:108841/1‑149 (MQ=255)                                                                                                                                           gaAGCCTGGCCGAGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGG‑CAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGCTGTTGGGTGGCCTGCAGTACCGcc             >  7:110382/1‑148 (MQ=255)                                                                                                                                           gaAGCCTGGCCGAGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGG‑CAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGCTGTTGGGTGGCCTGCAGTACCGcc             <  8:110382/148‑1 (MQ=255)                                                                                                                                              gCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGG‑CAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGCTGTTGGGTGGCCTGCAGTACCGCCtgtg         <  8:15836/149‑1 (MQ=255)                                                                                                                                               ccTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGG‑CAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGCTGTTGGGTGGCCTGCAGTACCGCCTgtgg        <  7:218142/149‑1 (MQ=255)                                                                                                                                               ccTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGG‑CAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGCTGTTGGGTGGCCTGCAGTACCGCCTgtgg        >  4:219246/1‑149 (MQ=255)                                                                                                                                                    ccTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGG‑CAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGCTGTTGGGTGGCCTGCAGTACCGCCTGtgggtggg   >  4:185989/1‑149 (MQ=255)                                                                                                                                                      tCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGG‑CAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGCTGTTGGGTGGCCTGCAGTACCGCCTGTgggtgggt  >  1:139830/1‑148 (MQ=255)                                                                                                                                                      |                                                                                                                                                     AGCAACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGG‑CAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGCTGTTGGGTGGCCTGCAGTACCGCCTGTGGGTGGGT  >  NC_002947/1810187‑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: ‑