Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	955,121	C→T	5.4%	intergenic (+72/‑241)	cyoE → / → alaC	protoheme IX farnesyltransferase/aminotransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	955,121	0	C	T	5.4%	285.2 / 4.3	93	intergenic (+72/‑241)	cyoE/alaC	protoheme IX farnesyltransferase/aminotransferase
Reads supporting (aligned to +/- strand):  ref base C (26/61);  new base T (2/3);  total (28/65)
Fisher's exact test for biased strand distribution p-value = 6.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCTGTAGGA‑GCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTT  >  NC_002947/954972‑955263                                                                                                                                                         |                                                                                                                                                 aCCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGcc                                                                                                                                                  >  1:3381908/1‑150 (MQ=255)  ccGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCg                                                                                                                                                 <  2:1748906/150‑1 (MQ=255)       aCTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCTGGccgc                                                                                                                                            >  1:2415225/1‑148 (MQ=255)       aCTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑ACTGTGCTGGCctc                                                                                                                                            >  1:2415202/1‑150 (MQ=255)        cTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCctct                                                                                                                                           >  2:2702163/1‑150 (MQ=255)            ccTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTcgc                                                                                                                                       <  1:2422145/150‑1 (MQ=255)             cTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTcgc‑g                                                                                                                                     <  1:3686209/150‑1 (MQ=255)              tGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑gg                                                                                                                                    <  1:4408450/150‑1 (MQ=255)               gAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑ggg                                                                                                                                   <  2:2723763/150‑1 (MQ=255)               gAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑ggg                                                                                                                                   <  1:3411251/150‑1 (MQ=255)               gAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑ggg                                                                                                                                   <  2:131095/150‑1 (MQ=255)               gAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑ggg                                                                                                                                   <  2:131173/150‑1 (MQ=255)                      atgatgGGGGTGGACAGCCAGACCGCTGCGGATGTGCGGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAAccc                                                                                                                            <  2:1829358/150‑1 (MQ=255)                        gatgGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGc                                                                                                                          <  1:1758577/150‑1 (MQ=255)                        gatgGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGc                                                                                                                          <  2:24304/150‑1 (MQ=255)                         atgGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTGGC‑GGGTTAACCCGcg                                                                                                                         >  1:3589542/1‑149 (MQ=255)                          tgGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTc                                                                                                                        <  2:4101146/150‑1 (MQ=255)                          tgGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTc                                                                                                                        <  2:407320/150‑1 (MQ=255)                           gggTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTcc                                                                                                                       <  2:3792994/150‑1 (MQ=255)                            ggTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑‑TTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTccga                                                                                                                     >  2:3998760/1‑148 (MQ=255)                             gtgtGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCtc                                                                                                                     <  1:3352901/150‑2 (MQ=255)                             gtgtGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTa                                                                                                                     <  2:1903644/150‑1 (MQ=255)                              tgtgGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTAc                                                                                                                    <  2:211112/150‑1 (MQ=255)                                tgGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACaa                                                                                                                  <  2:2543490/150‑1 (MQ=255)                                 ggACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAg                                                                                                                 <  2:2268638/150‑1 (MQ=255)                                   aCAGCCAGACCGCTGCGGTTGTGCTGTTGAACTTCGCACGCTGAAACTGCTGCCTTTTTCACGAAAACCCGGGCAACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAAACCGCTCCTACAAGGt                                                                                                               <  1:155555/150‑1 (MQ=38)                                    cAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑ATTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACaggaac                                                                                                              <  2:1577800/150‑6 (MQ=255)                                     aGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAATCCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTc                                                                                                             <  2:3263366/150‑1 (MQ=255)                                     aGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTc                                                                                                             <  2:3263324/150‑1 (MQ=255)                                     aGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTc                                                                                                             <  2:2749605/150‑1 (MQ=255)                                     aGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCGCTTCGC‑GGGGGAACCCGCGCCGACCAgggtg                                                                                                             >  1:1735055/1‑147 (MQ=255)                                      gCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAAACCGCTCCTACAAGGttac                                                                                                            >  2:2790275/1‑148 (MQ=255)                                       ccAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCt                                                                                                           <  1:2918796/150‑1 (MQ=255)                                        cAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCTg                                                                                                          <  1:2852808/150‑1 (MQ=255)                                        cAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCTg                                                                                                          <  1:2561845/150‑1 (MQ=255)                                        cAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAAACCGCTCCAACAAGGTTccgg                                                                                                          >  2:3937322/1‑148 (MQ=255)                                         aGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAAGTGCTGCCTGTTTCACGAAAACCCCGGCCTCTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCGGGGCTCTTCGC‑GGGTAAACCCGCTCCTACCCGGTTCgtgt                                                                                                         <  2:494172/150‑1 (MQ=255)                                         aGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCTGt                                                                                                         <  2:2839473/150‑1 (MQ=255)                                         aGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCTGt                                                                                                         <  1:719996/150‑1 (MQ=255)                                           aCCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAAACCGCTCCGACAAGGTTGCTGTAg                                                                                                       >  2:4261689/1‑150 (MQ=255)                                            ccGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCTGTAgg                                                                                                      <  2:1485955/150‑1 (MQ=255)                                             cGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCTGTAGGa                                                                                                     <  2:4312527/150‑1 (MQ=255)                                             cGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCTGTAGGa                                                                                                     <  2:3010037/150‑1 (MQ=255)                                              gCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑‑CTGTGCTGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAGGAACCTTGTAGGA‑Gc                                                                                                  <  1:1886367/150‑1 (MQ=25)                                               cTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAGGAACCTTGTAGGA‑Gc                                                                                                  <  2:246154/150‑1 (MQ=38)                                               cTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAGGAACCTTGTAGGA‑Gc                                                                                                  <  2:2903532/150‑1 (MQ=38)                                               cTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCTGTAGGA‑Gc                                                                                                  <  1:475127/150‑1 (MQ=255)                                               cTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCTGTAGGA‑Gc                                                                                                  <  1:1158589/150‑1 (MQ=255)                                               cTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCTGTAGGA‑Gc                                                                                                  <  1:3093212/150‑1 (MQ=255)                                               cTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCTGTAGGA‑Gc                                                                                                  <  1:475235/150‑1 (MQ=255)                                                tGCGTATGAGCTGATGACCTACGCACGATGAAACTGCTGCCTGTTTCACTAAAACCCCGCCCACTGT‑GCCGGGG‑TATTTTATTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCTGTAGGA‑GCg                                                                                                 <  2:2569575/150‑1 (MQ=38)                                                tGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑‑CTGTGCTGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAGGAACCTTGTAGGA‑GCgg                                                                                                <  1:3794698/150‑1 (MQ=25)                                                 gCGGTTGTGCTGATGCCCTACGCACGCTGAAACTGCGGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGT‑TTTTTTTTTGGGGGTATTTTT‑TCTGTGCTGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAGGAACCTTGTAGGA‑GCgg                                                                                                <  1:4326311/150‑1 (MQ=25)                                                 gCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCTGTAGGA‑GCgg                                                                                                <  2:3590650/150‑1 (MQ=255)                                                 gCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCTGTAGGA‑GCgg                                                                                                <  1:3296377/150‑1 (MQ=255)                                                  cGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCTGTAGGA‑GCggg                                                                                               <  1:2237250/150‑1 (MQ=255)                                                  cGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCTGTAGGA‑GCggg                                                                                               <  1:4274752/150‑1 (MQ=255)                                                   ggATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCTGTAGGA‑GCGGGt                                                                                              <  2:613526/150‑1 (MQ=255)                                                    gATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCTGTAGGA‑GCGGGtt                                                                                             <  1:3050452/150‑1 (MQ=255)                                                         gCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGACAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCTGTAGGA‑GCGGGTTTAccc                                                                                        <  2:2231241/150‑1 (MQ=255)                                                           tgatgaCCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAGGAACCTTGTAGGA‑GCGGGTTTACCcgc                                                                                      <  2:272561/150‑1 (MQ=37)                                                           tgatgaCCGACGCACGCTGAAACTGCTGCCTGTATCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TATTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCTGTAGGA‑GCGGGTTTACCcgc                                                                                      <  2:3358492/150‑1 (MQ=255)                                                            gatgaCCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAAACCGCGCCGACAAGGTTCCGGTAGGA‑GCGGGGTTTACcgcg                                                                                     >  2:3832093/1‑150 (MQ=255)                                                                 ccTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAATCCCGCTCCTACAAGGTTCCTGGAGGT‑GCGGGTTTTCCCGCGAagag                                                                                >  2:3363419/1‑150 (MQ=255)                                                                 ccTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGGGAACCCGCGCCAACAAGGGTGCGGTAGGA‑GCGGGGTTTCCCCCGCagag                                                                                >  1:2253098/1‑150 (MQ=255)                                                                 ccTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTGTGGGGGTATTTTT‑CCAGTGCCGGCCTCTTCGC‑GGGGAAAACCCCAACAACGAAGTTACGGGAGGG‑GCGGGTGTTCCcgcgtaaag                                                                                >  1:2252846/1‑145 (MQ=255)                                                                  cTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTGAACCCCCGCCGACAAGGTTCCGGTAGGT‑GCGGGGTTTCCCGCGGAGAgg                                                                               >  1:714139/1‑150 (MQ=255)                                                                  cTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCGGTAGGG‑GCGGGTTTACCCGCGAAGAgg                                                                               >  2:713829/1‑150 (MQ=255)                                                                             gAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAAACCGCGCCTACAAGGTTCCTGGAGGA‑GCGGGTTTACCCGCGCAGAGGCCAGCGCCGaa                                                                    >  2:1746550/1‑150 (MQ=38)                                                                                aCTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCTGTAGGA‑GCGGGTTTACCCGCGAAGAGGCCAGCACAGaaaaa                                                                 <  2:2560407/150‑1 (MQ=255)                                                                                aCTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTGT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCTGTAGGA‑GCGGGTTTACCCGCGAAGAGGCCAGCACAGaaaaa                                                                 <  2:2635085/150‑1 (MQ=255)                                                                                aCTGCTGCCTGTCGCACGAAACCCCCGGCCAGTGT‑GCGGGGT‑TTTTTTTTGGGGGGTATTTTT‑CGGGGGCCGGCCTCTGGGC‑GGGTAAACCCGCTCCTCCCCGGTTCCTGTAGGA‑GCGGGTTTACCCGCGAAGAGGCCAGCACAGaaaaa                                                                 <  2:3810503/150‑1 (MQ=25)                                                                                     tgcCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGGTAACCCCCTCCGACAAGGTTACTGTAGGA‑GCGGGTTTTCCCGCGGAGAGGCCAGCCCCGacaaataaaa                                                            >  2:245371/1‑150 (MQ=255)                                                                                       ccTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTAGC‑GGGTTAACACCCGCCGACAAGGTTGCAGGTGGAGGGGGGTTTACCC‑CGAAGGGGGCCGCACACgcaaataaaaat                                                          >  1:1793927/1‑150 (MQ=255)                                                                                       ccTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTATTCGC‑GGGTAAACCCGCTCCTACCAGGTTCCGGTAGGA‑GCGGGTTTACCCGCGTAGTGGCCAGCAAAGaaaaataaaaat                                                          >  2:364766/1‑150 (MQ=38)                                                                                         tGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTAGCGGGGTATACCC‑CGCCCACCACGTTGGGGGTGGG‑GGGGGGGTTCCCGCGGAGAGGCCCGCCCAGcaaaaaaaaaaaaa                                                        >  1:3129213/1‑147 (MQ=255)                                                                                         tGTTTCACGAAAACCCCGGCCACGGT‑GCCGGGG‑TTTTTTTTGGGGGGTATTTTT‑CCGGTGCCGGCCTCTCCGC‑GGGTAAACCCGCCCCACCCCGGTTCCTGTAGGA‑GCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATAc                                                        <  2:856846/150‑1 (MQ=37)                                                                                           tttCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAAACCGCACCTACAAGGTTGCTGTATGA‑GCGGGTTTTCCCGCGAAGAGGCCAGCAAAGAAAAATAAAAATACCt                                                      >  2:4223411/1‑150 (MQ=38)                                                                                           ccaaaaaaaaaCCCCGG‑CACAGTGGCCGGGGTTTTTTTTTTGGTGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCGCCTTCAAGGTTCCGGTAGGG‑GCGGGGTTACCCGCGAAGAGGCCAGCACAGGAAAATAAAAATAcc                                                       >  2:2233506/8‑150 (MQ=21)                                                                                             tCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCACCTACAAAGTTCCGGTAGGA‑GCGGGTGTACCCCCCCACAGGCCCGCAACGAAAAATAAAAATAccagg                                                    >  2:789220/1‑147 (MQ=37)                                                                                               aCGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAAACCCCTCCGACCAGGTTGCGGGTGGG‑GCGGGGTTAACCGCGCAGCGGCCAGCACCGAAAAATAAAAAAACCTGCtt                                                  >  2:4159108/1‑150 (MQ=255)                                                                                                    aaCCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCGTCGC‑GGGGAAACCCGCACCCACAAGGGTCCGGGAGGC‑GCGGGGGTACCCGCGCAGAGGCCCCCCCCGCAAAAAAAAAAAACCGGATTTTCaa                                             >  2:1956709/1‑150 (MQ=255)                                                                                                      ccccGGCCCCTGT‑GCGGGTG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACACGGTTCCTTTTGGA‑GCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCaaaa                                           <  1:734189/150‑1 (MQ=37)                                                                                                              aCTGT‑GCCGGGG‑ATTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAAACCCCTCCTACAAGGGTGCTGTAGGG‑GCGGGTTTTCCCGCGACGAGGCCAGCACAGAAAAAAAAAAATACCTGATTTTGAAAAAATAGAta                                   >  2:2791069/1‑149 (MQ=37)                                                                                                               cTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCTGTAGGA‑GCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCt                                  <  2:3381908/150‑1 (MQ=255)                                                                                                                          gttttttttttGGGGTTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGTTAAACCCGCTCCTACACGGTTCCTGTAGGA‑GCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTc                         <  2:429023/148‑1 (MQ=255)                                                                                                                           ttttttttttGGGGTTATTTTT‑CCGGTGCCGGCCTCTTCGC‑GGTTAAACCCGCTCCTACAAGGTTCCTGTCGGA‑GCGGGTTAACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCt                        <  2:1877266/149‑1 (MQ=38)                                                                                                                               ttttttGGGGGTTTTTTTCCCGGTGCCGCCCTCTTCGC‑GGGAACCCCCGCTCCTACACGGTACGTGTAGGA‑GCGGGTTCACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCtaat                     <  1:2702163/150‑1 (MQ=37)                                                                                                                                    tGGGGGTATTTTC‑CCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCTGTAGGA‑GCGGTTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACa               <  2:1028763/150‑1 (MQ=255)                                                                                                                                        ggTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTTAACCCGCTCCTACAAGGTTGCTGGAGGA‑GCGGGTTTAACCGCGGAGGGGCCAGCACCGGAAAATAAAAATACCTGTTTTTAAAAAAATAGGTATTAAAAGTTGTAATAAAACAGGaa           >  2:3510506/1‑150 (MQ=25)                                                                                                                                        ggTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAAGCCGCTCCTACAATGTTGCTGTAGGT‑GCGGGTTTAGCCGCGTAGAGGCCAGCACAGAAAAATTAAAATACCTGATTTTGTAAAAATAGATCTGAAAAGTTGTAATAAAAGAGGaa           >  2:1441133/1‑150 (MQ=25)                                                                                                                                           aTTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCTGTAGGA‑GCGGTTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAAtt        <  2:3075426/150‑1 (MQ=255)                                                                                                                                           aTTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCTGTAGGA‑GCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAAtt        <  2:3075433/150‑1 (MQ=255)                                                                                                                                                  cgtgtgCCGCCCTCTGGGG‑GGGTAACCCCGCTCCTCCACGGTTCCTGTTGGC‑GCGGTTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCtt  <  1:3216641/148‑1 (MQ=255)                                                                                                                                                         |                                                                                                                                                 ACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTT‑CCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCTGTAGGA‑GCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTT  >  NC_002947/954972‑955263

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑