Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,743,335	T→C	12.3%	intergenic (‑168/+79)	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,335	0	T	C	12.3%	158.9 / 11.8	65	intergenic (‑168/+79)	fbp/glgP	class 1 fructose‑1,6‑bisphosphatase/glycogen phosphorylase
Reads supporting (aligned to +/- strand):  ref base T (29/28);  new base C (4/4);  total (33/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.37e-01

GGGGGATGAAAAACCCACGCAGTTTACAGGGAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATC  >  NC_002947/5743188‑5743483                                                                                                                                                    |                                                                                                                                                     gggggATGAAAAACCCACGCAGTTTACAGGGAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGcgctc                                                                                                                                                     >  2:14971/1‑149 (MQ=255)  ggggATGAAAAACCCACGCAGTTTACAGGGAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGcgctcg                                                                                                                                                    <  1:14971/149‑1 (MQ=255)   gggATGAAAAACCCACGCAGTTTACAGGGAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGcgctcgc                                                                                                                                                   >  2:129084/1‑149 (MQ=255)    ggATGAAAAACCCACGCAGTTTACAGGGAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTACATTCGCGGGCGcgctcgct                                                                                                                                                  >  6:55164/1‑149 (MQ=255)      aTGAAAAACCCACGCAGTTTACAGGGAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTcc                                                                                                                                                <  5:3319/149‑1 (MQ=255)        gAAAAACCCACGCAGTTTACAGGGAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCACGCCCGCTCCca                                                                                                                                              >  4:73465/1‑149 (MQ=255)         aaaaaCCCACGCAGTTTACAGGGAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCACGCCCGCTCCca                                                                                                                                              >  1:223990/1‑148 (MQ=255)          aaaaCCCACGCAGTTTACAGGGAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATATGTCTGTAGTTCAGCCTTAGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCcac                                                                                                                                             <  4:267858/148‑1 (MQ=255)            aaCCCACGCAGTTTACAGGGAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAgg                                                                                                                                          <  3:365590/149‑1 (MQ=255)               ccACGCAGTTTACAGGGAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGCCACAATAGGTCGGTAGTTCAGCCATCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTa                                                                                                                                        >  6:352154/1‑148 (MQ=255)                cACGCAGTTTACAGGGAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACCGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCGTTCGCGGGCACGCCCGCTCCCACAGGtgtc                                                                                                                                      <  2:143644/149‑4 (MQ=255)                          aCAGGGAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCt                                                                                                                            >  4:128162/1‑149 (MQ=255)                           cAGGGAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCtt                                                                                                                           <  5:303116/149‑1 (MQ=255)                               gagaGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAg                                                                                                                       <  5:55164/149‑1 (MQ=255)                                agagTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGc                                                                                                                      <  3:350805/149‑1 (MQ=255)                                   gTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTg                                                                                                                   >  1:277644/1‑149 (MQ=255)                                    tGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGc                                                                                                                  <  5:109763/149‑1 (MQ=255)                                    tGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGc                                                                                                                  <  3:352766/149‑1 (MQ=255)                                      aCGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCga                                                                                                                <  2:219359/149‑1 (MQ=255)                                        gggCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCgaga                                                                                                              >  8:101951/1‑149 (MQ=255)                                        gggCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCgaga                                                                                                              <  2:139135/149‑1 (MQ=255)                                         ggCTATCGAGTGATGAGACTCGCCAGGTCTCTCTGACGGCTCATTATTTCGGTAGTTCAGCCTTCGCCAGGCGGGTGGTCTGTTTGTGTGTCCATTGGCGTGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGac                                                                                                             <  7:34813/149‑1 (MQ=255)                                           cTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTAGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGacac                                                                                                           >  7:275094/1‑149 (MQ=255)                                            tAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTACCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGacac                                                                                                           >  2:323254/1‑148 (MQ=255)                                                 aGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTgg                                                                                                     <  8:275094/149‑1 (MQ=255)                                                 aGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTgg                                                                                                     <  3:128162/149‑1 (MQ=255)                                                     atgaGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGc                                                                                                 >  6:31104/1‑149 (MQ=255)                                                     atgaGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGc                                                                                                 >  4:221716/1‑149 (MQ=255)                                                      tgaGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGGCCGTTCGCGGGCACGCCCGCTCCCACAGGTGTCTCGCAAGCTTCAAGGCTTGCGCTGTACCTGTGGGAGCg                                                                                                >  8:250880/1‑149 (MQ=17)                                                         gaCTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGCTGTACCTGTGGgagcgagc                                                                                             <  8:247678/149‑1 (MQ=39)                                                         gaCTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGc                                                                                             <  7:101951/149‑1 (MQ=255)                                                           cTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGt                                                                                           <  3:221716/149‑1 (MQ=255)                                                           cTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGGCCGTTCGCGGGCACGCCCGCTCCCACAGGTGTCTCGCAAGCTTCAAGGCTTGCGCTGTACCTGTGGGagcgagcgc                                                                                           <  7:250880/149‑2 (MQ=14)                                                            tCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTg                                                                                          >  4:143484/1‑149 (MQ=255)                                                            tCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTg                                                                                          >  6:98519/1‑149 (MQ=255)                                                              gCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCACGCCCGCTCCCACAGGTGTCTCGCAAGCTTCAAGGCTTGCGCTGTACCTGTGGGAGCGAgcgcgc                                                                                         <  2:223990/148‑1 (MQ=17)                                                                cACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGt                                                                                                                                         >  3:15747/1‑98 (MQ=255)                                                                cACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGt                                                                                                                                         <  4:15747/98‑1 (MQ=255)                                                                     gtctaAGAGGGCAAAATAGGTCGGTAGTTCAGCCTTGGCCAGGCGGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACCGCGCAAGCCTTGACGCTTGGGAGACACGTGTGGGTGCGGGCGTGCCCGCGAAc                                                                                 <  8:129820/148‑1 (MQ=255)                                                                      tctcAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACg                                                                                >  7:198234/1‑149 (MQ=255)                                                                      tctcAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGGGGGAGCGGGCGTGCCCGCGAACg                                                                                >  3:231317/1‑149 (MQ=255)                                                                      tctcAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGCTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACg                                                                                >  7:337571/1‑149 (MQ=255)                                                                        tcAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGc                                                                              >  8:276594/1‑149 (MQ=255)                                                                                   aaTAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAAt                                                                   >  4:215751/1‑149 (MQ=255)                                                                                    aTAGGTCGGTCGTGCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGGGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGGGCCCGCGAACGGCCTCACCTCAATc                                                                  <  5:350386/149‑1 (MQ=255)                                                                                     tAGGTCGTTAGTTCAGCCTTCGCCTGGCGGGTGGTCTGGTTGTGTGTCCATTCGCGGGGGGGCTGGATCGCAGAGGTACTGCGCATGCCTTGATGGTTGGGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCaa                                                                    <  8:337571/146‑1 (MQ=255)                                                                                        gTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCcgc                                                                                     >  6:193970/1‑126 (MQ=255)                                                                                        gTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCcgc                                                                                     <  5:193970/126‑1 (MQ=255)                                                                                        gTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAAt                                                              >  2:335942/1‑149 (MQ=255)                                                                                                     ccTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCgg                                                                                               <  6:180569/103‑1 (MQ=255)                                                                                                     ccTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCgg                                                                                               >  5:180569/1‑103 (MQ=255)                                                                                                          gCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGcct                                                                            >  2:237872/1‑117 (MQ=255)                                                                                                          gCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGcct                                                                            <  1:237872/117‑1 (MQ=255)                                                                                                               gCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCata                                       >  6:315657/1‑149 (MQ=255)                                                                                                                ccGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCatat                                      >  6:377592/1‑149 (MQ=255)                                                                                                                              ttgtgcCCATGCGCGGGCGGGCTCGCTCCCACAGGTCCAGCCAATGCCTTTATGGTTGCGTGACACGTGTGGTTGGGGGGGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGc                         <  1:323254/147‑1 (MQ=255)                                                                                                                              gtgtgtCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCGTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCGCAATCCACTGCCTTCCAGATCTCGGTGGCATATTCGCGAAGGGTGCg                        >  5:262940/1‑149 (MQ=255)                                                                                                                                 tgtCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCg                        <  7:82903/146‑1 (MQ=255)                                                                                                                                 tgtCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCg                        >  8:82903/1‑146 (MQ=255)                                                                                                                                  gtCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGt                                                                                           >  4:236461/1‑78 (MQ=255)                                                                                                                                  gtCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGt                                                                                           <  3:236461/78‑1 (MQ=255)                                                                                                                                   tCCATTCGCCGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGCCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCAAATTAGCAAATGCTGAGGTcc                    >  6:203800/1‑148 (MQ=38)                                                                                                                                      aTTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCcgacg                >  1:174449/1‑149 (MQ=255)                                                                                                                                       ttCGCGGGCGCGCTCGCTCCCCCCCGTACAGGCCACGCCTTCAAGCTTGCGACAACCCTGTGGGAGGGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCcgacga               <  6:262940/149‑1 (MQ=255)                                                                                                                                       ttCGCGGGCACGCCCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCtgt                                                                                                       <  7:101436/61‑1 (MQ=255)                                                                                                                                       ttCGCGGGCACGCCCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCtgt                                                                                                       >  8:101436/1‑61 (MQ=255)                                                                                                                                          gcgGGCGCGCTCGCTCCCACAGGTACTGCGCAAGCCTTGATGCTTGCGAGACACCTGTGGGAGCGGGCGTGCGCGCGAACGGCCTCACCTCATTCCAATGCGTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACgaga             <  5:352154/148‑1 (MQ=255)                                                                                                                                          gcgGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGaa            >  4:155274/1‑149 (MQ=255)                                                                                                                                            gggCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAAc           >  8:103962/1‑148 (MQ=255)                                                                                                                                                    tcgctcCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGt                           <  4:159908/124‑1 (MQ=255)                                                                                                                                                    tcgctcCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGt                           >  3:159908/1‑124 (MQ=255)                                                                                                                                                    tcgctcCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATc  <  3:143484/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                     GGGGGATGAAAAACCCACGCAGTTTACAGGGAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATC  >  NC_002947/5743188‑5743483

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