Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,761,201	A→G	6.1%	T85A (ACC→GCC)	gpmI →	2,3‑bisphosphoglycerate‑independent phosphoglycerate mutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,761,201	0	A	G	6.1%	187.0 / 2.4	66	T85A (ACC→GCC)	gpmI	2,3‑bisphosphoglycerate‑independent phosphoglycerate mutase
Reads supporting (aligned to +/- strand):  ref base A (31/30);  new base G (2/2);  total (33/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.35e-01

CGCCTGTGCGCCACCCAGCCGCATGGCCTCATTTCCGGCTCGGGCATGGATGTCGGCCTGCCGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTT‑CACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGG  >  NC_002947/5761057‑5761345                                                                                                                                                  |                                                                                                                                                 cGCCTGTGCGCCACCCAGCCGCATGGCCTCATTTCCGGCTCGGGCATGGATGTCGGCCTGCCGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTT‑CACGCGGGTGAcc                                                                                                                                                >  1:265614/1‑147 (MQ=255)   ccTGTGCGCCACCCAGCCGCATGGCCTCATTTCCGGCTCGGGCATGGATGTCGGCCTGCCGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTT‑CACGCGGGTGACCAAgg                                                                                                                                            <  1:176493/149‑1 (MQ=255)   ccTGTGCGCCACCCAGCCGCATGGCCTCATTTCCGGCTCGGGCATGGATGTCGGCCAGCCGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTT‑CACGCGGGTGACCAAgg                                                                                                                                            <  2:172405/149‑1 (MQ=255)    cTGTGCGCCACCCAGCCGCATGGCCTCATTTCCGGCTCGGGCATGGATGTCGGCCTGCCGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTT‑CACGCGGGTGACCAAGGc                                                                                                                                           >  5:194379/1‑149 (MQ=255)          gcCACCCAGCCGCATGGCCTCATTTCCGGCTCGGGCATGGATGTCGGCCTGCCGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTT‑CACGCGGGTGACCaa                                                                                                                                              >  2:59504/1‑140 (MQ=255)          gcCACCCAGCCGCATGGCCTCATTTCCGGCTCGGGCATGGATGTCGGCCTGCCGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTT‑CACGCGGGTGACCaa                                                                                                                                              <  1:59504/140‑1 (MQ=255)             aCCCAGCCGCATGGCCTCATTTCCGGCTCGGGCATGGATGTCGGCCTGCCGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTT‑CACGCGGGTGACCAAGGCCATCCGCGa                                                                                                                                  >  7:23020/1‑149 (MQ=255)             aCCCAGCCGCATGGCCTCATTTCCGGCTCGGGCATGGATGTCGGCCTGCCGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTT‑CACGCGGGTGACCAAGGCCATCCGCGa                                                                                                                                  <  2:265614/149‑1 (MQ=255)             aCCCAGCCGCATGGCCTCATTTCCGGCTCGGGCATGGATGTCGGCCTGCCGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTT‑CACGCGGATGGCCTTGGTCACCcgcgt                                                                                                                                  <  8:81127/149‑2 (MQ=255)             aCCCAGCCGCATGGCCTCATTTCCGGCTCGGGCATGGATGTCGGCCTGCCGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTT‑CACGCGGATGGCCTTGGTCACCcgcgt                                                                                                                                  <  7:130614/149‑2 (MQ=255)               ccAGCCGCATGGCCTCATTTCCGGCTCGGGCATGGATGTCGGCCTGCCGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTT‑CACGCGGGTGACCAAGGCCATCCGCGACg                                                                                                                                >  8:154178/1‑149 (MQ=255)               ccAGCCGCATGGCCTCATTTCCGGCTCGGGCATGGATGTCGGCCTGCCGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTT‑CACGCGGGTGACCAAGGCCATCCGCGACg                                                                                                                                >  1:338948/1‑149 (MQ=255)                cAGCCGCATGGCCTCATTTCCGGCTCGGGCATGGATGTCGGCCTGCCGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTT‑CACGCGGGTGACCAAGGCCATCCGCGACgg                                                                                                                               >  1:31126/1‑149 (MQ=255)                   ccGCATGGCCTCATTTCCGGCTCGGGCATGGATGTCGGCCTGCCGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTT‑CACGCGGGTGACCAAGGCCATCCGCTACGGCGa                                                                                                                            <  7:154178/149‑1 (MQ=255)                          gCCTCATTTCCGGCTCGGGCATGGATGTCGGCCTGCCGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTT‑CACGCGGGTGACCAAGGCCATCCGCGACGGCGAGttcttc                                                                                                                     <  4:153795/149‑1 (MQ=255)                          gCCTCATTTCCGGCTCGGGCATGGATGTCGGCCTGCCGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTT‑CACGCGGGTGACCAAGGCCATCCGCGACGGCGAGttcttc                                                                                                                     <  5:333513/149‑1 (MQ=255)                          gCCTCATTTCCGGCTCGGGCATGGATGTCGGCCTGCCGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTT‑CACGCGGGTGACCAAGGCCATCCGCGACGGCGAGttcttc                                                                                                                     <  6:194379/149‑1 (MQ=255)                           ccTCATTTCCGGCTCGGGCATGGATGTCGGCCTGCCGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTT‑CACGCGGGTGACCAAGGCCATCCGCGACGGCGAGttcttc                                                                                                                     <  3:147296/148‑1 (MQ=255)                           ccTCATTTCCGGCTCGGGCATGGATGTCGGCCTGCCGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTT‑CACGCGGGTGACCAAGGCCATCCGCGACGGCGAGttcttc                                                                                                                     >  4:147296/1‑148 (MQ=255)                               aTTTCCGGCTCGGGCATGGATGTCGGCCTGCCGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTT‑CACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGaaa                                                                                                                 <  1:296389/148‑1 (MQ=255)                                tttCCGGCTCGGGCATGGATGTCGGCCTGACGGAAGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTT‑CACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAAc                                                                                                               >  6:332084/1‑149 (MQ=255)                                      gCTCGGGCATGGATGTCGGCCTGCCGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTT‑CACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTa                                                                                                         >  3:169490/1‑149 (MQ=255)                                          gggCATGGATGTCGGCCTGCCGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTT‑CACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGa                                                                                                     <  8:23020/149‑1 (MQ=255)                                            gCATGGATGTCGGCCTGCCGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTT‑CACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGAcc                                                                                                   >  2:77570/1‑149 (MQ=255)                                            gCATGGATGTCGGCCTGCCGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTT‑CACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGAcc                                                                                                   >  6:125385/1‑149 (MQ=255)                                            gCATGGATGTCGGCCTGCCGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTT‑CACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGAcc                                                                                                   >  7:246526/1‑149 (MQ=255)                                                ggATGTCGGCCTGCCGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTT‑CACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGAc                                                                                                    >  8:215085/1‑144 (MQ=255)                                                ggATGTCGGCCTGCCGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTT‑CACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGAc                                                                                                    <  7:215085/144‑1 (MQ=255)                                                ggATGTCGGCCTGCCGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTT‑CACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCgg                                                                                                 <  3:203877/147‑1 (MQ=255)                                                ggATGTCGGCCTGCCGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTT‑CACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCgg                                                                                                 >  4:203877/1‑147 (MQ=255)                                                     tCGGCCTGCCGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTT‑CACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCtgt                                                                                           <  5:166059/148‑1 (MQ=255)                                                       ggCCTGCCGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTT‑CACGCGGGTGACCAAGGCCATCCGCGACAGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGa                                                                                        >  3:238509/1‑149 (MQ=255)                                                          cTGGCGGACGGGCAGATGGGCAACTCCGTAGTCGTTCATATGATAGTGGGCGCAGGTCGCGTTGTGTATGAGGACTT‑CACTCTGGTGTCCATGGCCATGCGTGACTGCGTGTTCTTGTATAAACCTGTAGTGTCCGGTGCTGTTGatta                                                                                     <  4:238509/149‑4 (MQ=255)                                                             ccGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTT‑CACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAgg                                                                                   <  2:338948/148‑1 (MQ=255)                                                             ccGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTT‑CACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAg                                                                                    >  1:131184/1‑147 (MQ=255)                                                             ccGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTT‑CACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAg                                                                                    <  2:131184/147‑1 (MQ=255)                                                              cGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTT‑CACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCa                                                                                 <  3:120010/149‑1 (MQ=255)                                                                  cGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTT‑CACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCg                                                                             <  4:169490/149‑1 (MQ=255)                                                                     gCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTT‑CACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGAcc                                                                                                   <  5:322479/124‑1 (MQ=255)                                                                     gCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTT‑CACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGAcc                                                                                                   >  6:322479/1‑124 (MQ=255)                                                                     gCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTT‑CACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCtgt                                                                                           <  6:117637/132‑1 (MQ=255)                                                                     gCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTT‑CACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCtgt                                                                                           >  5:117637/1‑132 (MQ=255)                                                                      cAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTT‑CACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGc                                                                         <  2:31126/149‑1 (MQ=255)                                                                       aGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTT‑CACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCt                                                                        >  2:188969/1‑149 (MQ=255)                                                                       aGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGAGTATCAGGACTT‑CACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCt                                                                        >  8:210021/1‑149 (MQ=255)                                                                          tGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTT‑CACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACaa                                                                                     >  4:198910/1‑133 (MQ=255)                                                                          tGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTT‑CACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACaa                                                                                     <  3:198910/133‑1 (MQ=255)                                                                           gggCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTT‑CACGCGGGTGACCAAGGCCATCCGCGACGGCGAGttcttc                                                                                                                     <  6:14475/100‑1 (MQ=255)                                                                           gggCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTT‑CACGCGGGTGACCAAGGCCATCCGCGACGGCGAGttcttc                                                                                                                     >  5:14475/1‑100 (MQ=255)                                                                            ggCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTT‑CACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCaa                                                                   <  5:332084/149‑1 (MQ=255)                                                                            ggCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTT‑CACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCaa                                                                   <  5:125385/149‑1 (MQ=255)                                                                             gCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTT‑CACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAg                                                                  >  5:282555/1‑149 (MQ=255)                                                                                                   gAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTTACACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGActg                                             >  1:56489/1‑149 (MQ=255)                                                                                                   gAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTT‑CACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAgg                                                                                   >  6:265739/1‑110 (MQ=255)                                                                                                   gAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTT‑CACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAgg                                                                                   <  5:265739/110‑1 (MQ=255)                                                                                                            cgcAGGTCGCGTCGTGTATCAGGACTT‑CACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGAtg                                   >  8:164326/1‑149 (MQ=255)                                                                                                               aGGTCGCGTCGTGTATCAGGACTT‑CACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGAtggtg                                <  1:77570/149‑1 (MQ=255)                                                                                                                          tgtATCAGGACTT‑CACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGc                     >  4:229458/1‑149 (MQ=255)                                                                                                                            tATCAGGACTT‑CACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGAtggtg                                >  7:286290/1‑136 (MQ=255)                                                                                                                            tATCAGGACTT‑CACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGAtggtg                                <  8:286290/136‑1 (MQ=255)                                                                                                                             aTCAGGACTT‑CACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGccac                  >  5:48969/1‑149 (MQ=255)                                                                                                                              tCAGGACTT‑CACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCCGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGccac                  <  8:246526/148‑1 (MQ=255)                                                                                                                                     tt‑CACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCa          >  3:189472/1‑149 (MQ=255)                                                                                                                                       cgtcgcgGATGGCCTTGGTCACCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCAc         >  7:19803/4‑149 (MQ=255)                                                                                                                                        gtcgcgGATGGCCTTGGTCACCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACt        >  1:292940/3‑149 (MQ=255)                                                                                                                                              ggTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTtggtgg  <  7:164326/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                 CGCCTGTGCGCCACCCAGCCGCATGGCCTCATTTCCGGCTCGGGCATGGATGTCGGCCTGCCGGACGGGCAGATGGGCAACTCCGAAGTCGGTCATATGAACCTCGGCGCAGGTCGCGTCGTGTATCAGGACTT‑CACGCGGGTGACCAAGGCCATCCGCGACGGCGAGTTCTTCGAAAACCCGGTACTGACCGGTGCTGTGGACAAGGCAGTCGGTGCTGGCAAGGCCGTGCACATCCTCGGACTGCTGTCCGATGGTGGTGTACACAGCCACCAGGACCACTTGGTGG  >  NC_002947/5761057‑5761345

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