Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,810,336	C→A	10.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	10.1%	98.2 / 2.8	39	intergenic (+97/‑72)	eno/ftsB	enolase/cell division protein
Reads supporting (aligned to +/- strand):  ref base C (16/19);  new base A (2/2);  total (18/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 = 7.95e-01

ACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGCTGTTGGGTGGCCTGCAGTACCGCCTGT  >  NC_002947/1810191‑1810474                                                                                                                                                   |                                                                                                                                           aCTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGc                                                                                                                                         >  8:1863/1‑149 (MQ=255)        gcCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCGAGGCCAGGCtt                                                                                                                                  <  3:108987/149‑1 (MQ=255)          cAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCg                                                                                                                                >  4:88120/1‑149 (MQ=255)                     aCCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAg                                                                                                                     >  6:112703/1‑149 (MQ=255)                     aCCGCGGTCGTGCCGAGTTTAGCGTCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTATGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGa                                                                                                                      <  8:90229/148‑1 (MQ=255)                      ccGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAg                                                                                                                     <  2:141923/148‑1 (MQ=255)                      ccGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTAGAgg                                                                                                                    >  1:141923/1‑149 (MQ=255)                         cgGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGccc                                                                                                                 <  1:95683/149‑1 (MQ=255)                              gTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGaa                                                                                                             <  2:246343/148‑1 (MQ=255)                                 ccGAGTTTCGCGGCTAAGCAAGAGATTGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTaa                                                                                                         <  3:88120/149‑1 (MQ=255)                                 ccGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTaa                                                                                                         <  2:102691/149‑1 (MQ=255)                                 ccGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTaa                                                                                                         <  5:54370/149‑1 (MQ=255)                                   gAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAAc                                                                                                        >  5:132348/1‑148 (MQ=255)                                              cTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCgg                                                                                              >  8:186882/1‑147 (MQ=255)                                              cTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCa                                                                                            <  6:132348/149‑1 (MQ=255)                                              cTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCa                                                                                            <  5:112703/149‑1 (MQ=255)                                                      gagaTGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCttttt                                                                                    <  1:196064/149‑1 (MQ=255)                                                            gTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACCGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTg                                                                              <  4:14732/149‑1 (MQ=255)                                                               aaaaGGCAGCAGCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGa                                                                            <  7:1863/148‑1 (MQ=255)                                                                              gAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCAGGGCTCGGCCAGGCTTAGTGGGATTAGAGGCCCCCAATTAAAAATATAGCCGCAGGCGTTTTTAACAGGATAAAGTTATgccgagt                                                            >  4:1765/1‑144 (MQ=255)                                                                                   gTAGGAAAGTTCAGTAC‑GAACGTTCCTACGTACTAACGGGTTTCTGTCGACGAAAACAGTGCTAGGCCAAGGTTAGTCCTAATAGGCGCCCCGAAGAAAAAAGTACGGGGCAGACATTTTAAACAGGAATAAGTTATGAGGAGTccca                                                        >  2:115573/1‑147 (MQ=255)                                                                                                gTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTcc                                            <  7:186882/148‑1 (MQ=255)                                                                                                 tACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTc                                          >  7:126835/1‑149 (MQ=255)                                                                                                  aCTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCg                                         >  8:173483/1‑149 (MQ=255)                                                                                                    tGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGt                                        >  3:224696/1‑148 (MQ=255)                                                                                                                cGTACTATCGGGTTTCTGTCGTCGAAGCCTGGCCTCGGCCTGGCTTCGTGCTATTCGAGTCGCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTTATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTg                              <  7:97104/146‑1 (MQ=255)                                                                                                                    cTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGCTGTTgg                       <  1:210528/149‑1 (MQ=255)                                                                                                                         gggTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATa                                                                          >  4:85304/1‑93 (MQ=255)                                                                                                                         gggTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATa                                                                          <  3:85304/93‑1 (MQ=255)                                                                                                                         gggTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTACTCGTTCTGCTGTTGCTGTTGGGTTGc                  >  8:254302/1‑149 (MQ=255)                                                                                                                          ggTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGGTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCGTTTTAACCTGGATAAAGTGATGCGCAGTCGCTATTGGTTGTTCCTCGTTCCGCTCTTGCTGTTGGGTGGcc                 >  5:233921/1‑149 (MQ=255)                                                                                                                          ggTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGCTGTTGGGTGGcc                 <  7:173483/149‑1 (MQ=255)                                                                                                                                tGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATTACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGCTGTTGGGTGGCCTGCAg            >  2:156865/1‑148 (MQ=255)                                                                                                                                 tagcacgaAGCCTGGCCGAGGCCAGGCTTCATGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGgttgt                                               >  6:105722/5‑112 (MQ=39)                                                                                                                                 tagcacgaAGCCTGGCCGAGGCCAGGCTTCATGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGgttgt                                               <  5:105722/108‑1 (MQ=39)                                                                                                                                    gacgaAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGCTGTTGGGTGGCCTGCAGTAcc        <  4:224696/148‑1 (MQ=255)                                                                                                                                    gacgaAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGCTGTTGGGTGGCCTGCAGTACCg       >  3:78650/1‑149 (MQ=255)                                                                                                                                       gaAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCCGCGTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGCTGTTGGGTGGCCTGCAGTACCGCCt    >  1:260251/1‑149 (MQ=255)                                                                                                                                         aGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGCTGTTGGGTGGCCTGCAGTACCGCCtgt  <  3:1765/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                           ACTGGGCGCCAAAGCGGTTTACCGCGGTCGTGCCGAGTTTCGCGGCTAAGCAAGAGATGGTAAAAAGGCAGCAGCCGGAGTTGTAGGAACGTTC‑GTACTGAACTTTCCTACGTACTAACGGGTTTCTGTCGACGAAGCCTGGCCTCGGCCAGGCTTCGTGCTATTCGAGGCCCCGAAGTAACAAGTACGCGGCAGCCTTTTTAACCTGGATAACGTGATGCGCAGTCCCTATTGGTTGTTCCTCGTTCTGCTCTTGCTGTTGGGTGGCCTGCAGTACCGCCTGT  >  NC_002947/1810191‑1810474

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