Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,743,350	C→T	8.1%	intergenic (‑183/+64)	fbp ← / ← glgP	class 1 fructose‑1,6‑bisphosphatase/glycogen phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,743,350	0	C	T	8.1%	221.5 / 6.0	75	intergenic (‑183/+64)	fbp/glgP	class 1 fructose‑1,6‑bisphosphatase/glycogen phosphorylase
Reads supporting (aligned to +/- strand):  ref base C (38/30);  new base T (3/3);  total (41/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

CACGCAGTTTACAGGGAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCC‑TTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCA  >  NC_002947/5743203‑5743496                                                                                                                                                     |                                                                                                                                                   cACGCAGTTTACAGGGAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCC‑TTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTCCa                                                                                                                                                   >  7:99991/1‑149 (MQ=255)  aCGCAGTTTACAGGGAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCC‑TTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAg                                                                                                                                                  >  7:218237/1‑149 (MQ=255)           aCAGGGAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCC‑TTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGcc                                                                                                                                          <  2:3062/148‑1 (MQ=255)             aGGGAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCC‑TTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCtt                                                                                                                                        >  3:202238/1‑148 (MQ=255)              gggAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCC‑TTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGa                                                                                                                                      <  7:199141/149‑1 (MQ=255)               ggAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCC‑TTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTg                                                                                                                                       >  2:375856/1‑147 (MQ=255)               ggAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCC‑TTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTg                                                                                                                                       >  7:147660/1‑147 (MQ=255)               ggAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCC‑TTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTg                                                                                                                                       >  1:258396/1‑147 (MQ=255)                gagaGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCC‑TTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGaa                                                                                                                                     <  7:286797/148‑1 (MQ=255)                gagaGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCC‑TTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAg                                                                                                                                    >  4:313315/1‑149 (MQ=255)                      gACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCC‑TTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGc                                                                                                                               <  8:40846/148‑1 (MQ=255)                        cGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCC‑TTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTAGAGCGCAAGCCTTGAAGCTTGCga                                                                                                                             <  6:323901/148‑1 (MQ=255)                        cGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCC‑TTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCga                                                                                                                             <  8:147660/148‑1 (MQ=255)                            cTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCC‑TTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCGTTCGCGGGCACGCCCGCTCCCACAGGTGTCTCGCAAGCTTCAAGGCTTGCgctgtac                                                                                                                        <  5:97857/149‑7 (MQ=21)                             tAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCC‑TTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACAcc                                                                                                                       >  8:171899/1‑149 (MQ=255)                                  aGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCC‑TTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGTCCGCTCCCACAGGTGTCTCGCAAGCTTCAAGGCTTGCGCTGTACCTGTgg                                                                                                                  <  6:320911/149‑1 (MQ=21)                                    tgatgaGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCC‑TTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCCCGCTCCCACAGGTGTCTCGCAAGCTTCAAGGCTTGCGCTGTACCTGTGGGa                                                                                                                <  2:168572/149‑1 (MQ=21)                                      atgaGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCAATCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGa                                                                                                                <  1:324073/148‑1 (MQ=255)                                      atgaGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCC‑TTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGc                                                                                                              >  1:370381/1‑149 (MQ=255)                                      atgaGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCC‑TACGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTAGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGc                                                                                                              >  1:96476/1‑149 (MQ=255)                                         agaCTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCC‑TTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCgg                                                                                                            <  5:6296/148‑1 (MQ=255)                                         agaCTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCC‑TTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCgg                                                                                                            <  1:375856/148‑1 (MQ=255)                                         agaCTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCC‑TTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCACTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCgg                                                                                                            <  5:42945/148‑1 (MQ=255)                                          gaCTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCC‑TTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGTGAGACACCTGTGGGAGCGGGc                                                                                                          <  8:317617/149‑1 (MQ=255)                                          gaCTCGCCACGGCTCTCAGACGGCACAATAGGGCGGTAGTTCAGCC‑TTCGCCAGGCCGGTCGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGc                                                                                                          <  2:153006/149‑1 (MQ=255)                                           aCTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGAC‑TTCGCCCGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGGACATCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGTCGCGGGCg                                                                                                         >  5:204470/1‑149 (MQ=255)                                               gCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCC‑TTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGc                                                                                                      <  6:28416/148‑1 (MQ=255)                                               gCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCC‑TTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGc                                                                                                      <  1:30000/148‑1 (MQ=255)                                                    ggCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCC‑TTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCcgcg                                                                                                 >  8:289833/1‑148 (MQ=255)                                                         tcAGACGGCACAATAGGTCGGTAGTTCAGCC‑TTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCGCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGc                                                                                           >  3:46627/1‑149 (MQ=255)                                                           aGACGGCACAATAGGTCGGTAGTTCAGCC‑TTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCGCCCACCGTTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGc                                                                                           <  3:54336/147‑1 (MQ=255)                                                            gACGGCACAATAGGTCGGTAGTTCAGCC‑TTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCg                                                                                                                              <  3:54597/111‑1 (MQ=255)                                                            gACGGCACAATAGGTCGGTAGTTCAGCC‑TTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCg                                                                                                                              >  4:54597/1‑111 (MQ=255)                                                             aCGGCACAATAGGTCGGTAGTTCAGCC‑TTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGcctca                                                                                       >  2:37896/1‑149 (MQ=255)                                                               ggCACAATAGGTCGGTAGTTCAGCC‑TTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGcctcacc                                                                                     >  4:336332/1‑149 (MQ=255)                                                                    aaTAGGTCGGTAGTTCAGCC‑TTCACCAGGCTGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAAt                                                                                <  2:66869/149‑1 (MQ=255)                                                                        ggTCGGTAGTTCAGCC‑TTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCaa                                                                            <  2:96476/149‑1 (MQ=255)                                                                         gTCGGTAGTTCAGCC‑TTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAAt                                                                           <  7:171899/149‑1 (MQ=255)                                                                         gTCGGTAGTTCAGCC‑TTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAAt                                                                           <  7:289833/149‑1 (MQ=255)                                                                          tCGGTAGTTCAGCC‑TTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATg                                                                          <  3:313315/149‑1 (MQ=255)                                                                           cGGTAGTTCAGCC‑TGCGCCAGGCCGGTGGTCTGTTGGTGTGTCCATTCGCGGGCGCGGTCGCCCCCACTGGTACTGCGCAAGCGTTGAAGCTTGCGAGACACCTGTGGGAGGGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGc                                                                         <  7:155470/149‑1 (MQ=255)                                                                                 ttCAGCC‑TTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTcc                                                                    >  8:206184/1‑148 (MQ=255)                                                                                      cc‑TTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACAGCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGAtct                                                              <  2:258396/149‑1 (MQ=255)                                                                                            gCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTTAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGt                                                          <  4:164785/148‑1 (MQ=255)                                                                                             ccAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTgg                                                        >  4:129706/1‑149 (MQ=255)                                                                                             ccAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTgg                                                        <  8:218237/149‑1 (MQ=255)                                                                                                     gtggtCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTcg                                                >  3:70388/1‑149 (MQ=255)                                                                                                      tggtCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTcgc                                               <  3:129706/149‑1 (MQ=255)                                                                                                             gTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAGCGGCCGCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGt                                        >  1:197251/1‑149 (MQ=255)                                                                                                             gTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCGCCTCACTGCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATggg                                        >  5:165103/1‑148 (MQ=255)                                                                                                             gTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGt                                        >  6:44865/1‑149 (MQ=255)                                                                                                             gTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGt                                        >  4:10983/1‑149 (MQ=255)                                                                                                                gtgtgtCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCg                                     >  2:65143/1‑149 (MQ=255)                                                                                                                    gtCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTcc                                 >  7:357448/1‑149 (MQ=255)                                                                                                                    gtCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTcc                                 >  2:51443/1‑149 (MQ=255)                                                                                                                        aTTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCcgac                              >  6:126978/1‑148 (MQ=255)                                                                                                                        aTTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCcgac                              >  2:355598/1‑148 (MQ=255)                                                                                                                        aTTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTGCCAGATCTCGGTGGCATATTCGAGACTGGTGCGGTAcgac                              >  6:220390/1‑148 (MQ=255)                                                                                                                        aTTCGCGGGCGCGCTCGCT‑CCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTCGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCcgac                              >  8:210593/1‑147 (MQ=255)                                                                                                                         ttCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCg                                     <  5:148228/140‑1 (MQ=255)                                                                                                                         ttCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCg                                     >  6:148228/1‑140 (MQ=255)                                                                                                                              gggCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTg                                       >  1:208923/1‑133 (MQ=255)                                                                                                                              gggCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTg                                       <  2:208923/133‑1 (MQ=255)                                                                                                                                gcgcgcTCGCTCCCACAGGTACAGCCCAATCCTTCAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCTCGAACGGCCTCACCTCAATCCAATGCCTTCCAGTTCTCGTTGGCATATTCTCCTTTGTTGCGGTCCGACGAGAACCAg                     <  6:165103/149‑1 (MQ=255)                                                                                                                                  gcgcTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACgag                           >  7:15732/1‑141 (MQ=255)                                                                                                                                  gcgcTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACgag                           <  8:15732/141‑1 (MQ=255)                                                                                                                                    gcccgctcCCACAGGTGTCTCGCAAGCTTCAAGGCTTGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCACTCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGccc                  >  7:209922/4‑148 (MQ=18)                                                                                                                                     ctcgctcCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCAt                >  8:119084/1‑149 (MQ=255)                                                                                                                                      tcgctcCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATc               >  5:239103/1‑149 (MQ=255)                                                                                                                                        gctcCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGAGTTCCCGCGAACGGACTAACCTCAATCCAATGCCGTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCAGACGAGAACCAGCACATccg             >  5:276760/1‑149 (MQ=255)                                                                                                                                          tcCCACAGGTGTCTCGCAAGCTTCAAGGCTTGCGCTGTACCTGTGGGAGCGAGCGCGCCCGCGAATGGACACACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATccgtg           >  5:111876/1‑149 (MQ=17)                                                                                                                                          tcCCACAGGTGTCTCGCAAGCTTCAAGGCTTGCGCTGTACCTGTGGGAGCGAGCGCGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATccgtg           >  5:256017/1‑149 (MQ=18)                                                                                                                                              acaGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATccgtgccgt       >  2:376175/1‑149 (MQ=255)                                                                                                                                                   tACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTc   >  3:155752/1‑148 (MQ=255)                                                                                                                                                   tACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCa  <  8:99991/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                   CACGCAGTTTACAGGGAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCC‑TTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCA  >  NC_002947/5743203‑5743496

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