Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,685,610	G→A	7.1%	T390I (ACC→ATC)	patD ←	medium chain aldehyde dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,685,610	0	G	A	7.1%	206.0 / 3.9	70	T390I (ACC→ATC)	patD	medium chain aldehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (22/43);  new base A (2/3);  total (24/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.38e-01

ACCAGCATGAAGTGGGTATTCACCCAGGTGCAGCCATACTGCAAGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCG‑CACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCTCGAAGAAAAAGCCATCGCCCGCTATTGGCTTGCCACCGGTTACCAGACGAATGTGCGGCTG  >  NC_002947/1685463‑1685758                                                                                                                                                    |                                                                                                                                                      gccAGCATGAAGTGGGTATTCACCCAGGTGCAGCCATACTGCAAGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCgg                                                                                                                                                       <  1:270728/147‑1 (MQ=255)    aGCATGAAGTGGGTATTCACCCAGGTGCAGCCATACTGCAAGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGaa                                                                                                                                                   <  2:10666/149‑1 (MQ=255)        tgatctGGGTAATCACCCAGGTGCAGACATACTGCAAGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCgg                                                                                                                                               <  8:90425/144‑1 (MQ=255)        tGAAGTGGGTATTCACCGAGGTGCAGGCATACTGCAAGCGCGCCGCCAGGGGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCGGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCACCGATGAAGCgc                                                                                                                                               <  4:355469/149‑2 (MQ=255)          aaGTGGGTATTCACCCAGTTGCAGCCATACTGCAAGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCggg                                                                                                                                              <  5:353160/148‑1 (MQ=255)          aaGTGGGTATTCACCCAGGTGCAGCCATACTGCAAGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCggg                                                                                                                                              <  8:149250/148‑1 (MQ=255)          aaGTGGGTATTCACCCAGGTGCAGCCATACTGCAAGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGt                                                                                                                                             <  7:277325/149‑1 (MQ=255)           aGTGGGTATTCACCCAGGTGCAGCCATACTGCAAGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGt                                                                                                                                             <  2:320890/148‑1 (MQ=255)           aGTGGGTATTCACCCAGGTGCAGCCATACTGCAAGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCACCGATGAAGCgcagg                                                                                                                                            <  3:134360/149‑5 (MQ=255)           aGTGGGTATTCACCCAGGTGCAGCCATACTGCAAGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCACCGATGAAGCgcagg                                                                                                                                            <  2:38967/149‑5 (MQ=255)            gTGGGTATTCACCCAGGTGCAGCCATACTGCAAGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGTTGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGa                                                                                                                                           >  6:363517/1‑149 (MQ=255)             tGGGTATTCACCCAGGTGCAGCCATACTGCAAGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGAc                                                                                                                                          >  2:181633/1‑149 (MQ=255)                     cACCCAGGTGCAGCCATACTGCAAGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAacca                                                                                                                                  >  8:184934/1‑149 (MQ=255)                        ccAGGTGCAGCCATACTGCAAGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCg                                                                                                                               <  4:358653/149‑1 (MQ=255)                        ccAGGTGCAGCCATACTGCAAGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCg                                                                                                                               >  8:107839/1‑149 (MQ=255)                         cAGGTGCAGCCATACTGCCAGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGGCGCGGGGCCATACCGATGCGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCgg                                                                                                                              <  1:363523/149‑1 (MQ=255)                                aGCCATACTGCAAGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTTCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAac                                                                                                                       <  1:73511/149‑1 (MQ=255)                                   cATACTGCAAGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACt                                                                                                                    <  6:121620/149‑1 (MQ=255)                                      aCTGCAAGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCg                                                                                                                 >  2:334500/1‑149 (MQ=255)                                         gCAAGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTg                                                                                                              >  7:144178/1‑149 (MQ=255)                                          cAAGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGtt                                                                                                               <  5:360239/147‑1 (MQ=255)                                          cAAGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGtt                                                                                                               <  7:184934/147‑1 (MQ=255)                                          cAAGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGtt                                                                                                               <  7:326328/147‑1 (MQ=255)                                          cAAGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGtt                                                                                                               <  1:181633/147‑1 (MQ=255)                                          cAAGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTgc                                                                                                             <  6:6639/149‑1 (MQ=255)                                           aaGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTgcg                                                                                                            >  4:177827/1‑149 (MQ=255)                                                cgccgcCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCATGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACATCGTTGCG‑CACGa                                                                                                      >  3:257853/1‑149 (MQ=255)                                                cgccgcCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCG‑CACGa                                                                                                      <  3:61519/149‑1 (MQ=255)                                                  ccgccAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGGTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCG‑CACGAtc                                                                                                    <  4:157869/149‑1 (MQ=255)                                                   cgccAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCG‑CACGAtct                                                                                                   >  1:236122/1‑149 (MQ=255)                                                           ggTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCG‑CACGATCTCGTCGTcc                                                                                           >  3:354774/1‑149 (MQ=255)                                                                   cgGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCACATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCG‑CACGATCTCGTCGTCCTGCAGGGc                                                                                   <  6:71041/149‑1 (MQ=255)                                                                     gCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCG‑CACGATCTCGTCGTCCTGCAGGGCGt                                                                                 <  3:213527/149‑1 (MQ=255)                                                                      ccGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCGCCACGATCTCGTCGTCCTGCAGGGCGt                                                                                 <  1:196279/149‑1 (MQ=255)                                                                        ggTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCG‑CACGATCTCGTCGTCCTGCAGGGCGTCgg                                                                              <  2:80350/149‑1 (MQ=255)                                                                            tCGCGGGTCCATACCGATGCGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCG‑CACGATCTCGTCGTCCTGCAGGGCGTCGGCCa                                                                           <  3:200814/148‑1 (MQ=255)                                                                                gggTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCG‑TGCG‑CACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACgg                                                                     <  1:136507/149‑1 (MQ=255)                                                                                     cATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCG‑CACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAg                                                                  <  1:334500/148‑1 (MQ=255)                                                                                       tACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCG‑CACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCt                                                               <  5:210805/149‑1 (MQ=255)                                                                                        aCCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCG‑CACGATCTAGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCTc                                                              >  6:25025/1‑149 (MQ=255)                                                                                           gATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCG‑CACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCTCgaa                                                           <  5:176284/149‑1 (MQ=255)                                                                                                      gccGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCG‑CACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCTCGAAGAAAAAGCCAt                                                <  4:354774/149‑1 (MQ=255)                                                                                                       ccGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCG‑CACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCTCGAAGAAAAAGCCATc                                               <  8:144178/149‑1 (MQ=255)                                                                                                        cGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCG‑CACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCTCGAAGAAAAAGCCATCg                                              >  5:130660/1‑149 (MQ=255)                                                                                                          tAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCG‑CACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCTCgaa                                                           >  5:66658/1‑134 (MQ=255)                                                                                                          tAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCG‑CACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCTCgaa                                                           <  6:66658/134‑1 (MQ=255)                                                                                                            gTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCG‑CACGAtctc                                                                                                  <  2:74753/93‑1 (MQ=255)                                                                                                            gTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCG‑CACGAtctc                                                                                                  >  1:74753/1‑93 (MQ=255)                                                                                                            gTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCG‑CACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCTCGAAGAAAAAGCCATCGCCCg                                          <  7:107839/149‑1 (MQ=255)                                                                                                            gTTGGAGGCGTTGGCCCATTCCAGTGCCTGCTCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCG‑GACGATCTGGGCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCTCGAAGAAAAAGCCATCGCCCg                                          <  5:377190/149‑1 (MQ=255)                                                                                                             ttGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCG‑CACGATCTCGTCGTCCTGCAGGGCGTCGGc                                                                             <  8:217889/113‑1 (MQ=255)                                                                                                             ttGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCG‑CACGATCTCGTCGTCCTGCAGGGCGTCGGc                                                                             >  7:217889/1‑113 (MQ=255)                                                                                                               ggAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCG‑CACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCTCGAAGAAAAAGCCATCGCCCGCTa                                       >  2:20402/1‑149 (MQ=255)                                                                                                                gAGTCGTTGGCCCATTCCTGTGCGTGGTCTTCATCGGTGAAGCGGGTGACCTAAACCTCGGGGCCAAACACTTCGTTGCG‑AACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTATGCTCGAAGAAAAAGCCATCGCCCGCTa                                       <  3:39770/148‑1 (MQ=255)                                                                                                                 aGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCG‑CACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCTCGAAGAAAAAGCCATCGCCCGCTAtt                                     <  5:363517/149‑1 (MQ=255)                                                                                                                  gTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCG‑CACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCTCGAAGAAAAAGCCATCGCCCGCTATTg                                    <  8:29534/149‑1 (MQ=255)                                                                                                                      ttGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCG‑CACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCTCGAAGAAAAAGCCATCGCCCGCTATTGGCtt                                <  2:256249/149‑1 (MQ=255)                                                                                                                        ggCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCG‑CACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCTCGAAGAAAAAGCCATCGCCCGCTATTGGCTTg                               <  7:34298/148‑1 (MQ=255)                                                                                                                         gCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCG‑CACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCTCGAAGAAAAAGCCATCGCCCGCTATTGGCTTGc                              <  1:20402/148‑1 (MQ=255)                                                                                                                             aTTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCG‑CACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCTCGAAGAAAAAGCCATCGCCCGCTATTGGCTTGCCACCg                         <  2:236122/149‑1 (MQ=255)                                                                                                                                 cAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCG‑CACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCTCGAAGAAAAAGCCATCGCCCGCTATTGGCTTGCCACCGGTTa                     >  8:261795/1‑149 (MQ=255)                                                                                                                                  aGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCG‑CACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCTCGAAGAAAAAGCCATCGCCCGCTATTGGCTTGCCACCGGTTAc                    <  3:177827/149‑1 (MQ=255)                                                                                                                                   gTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCG‑CACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCTCGAAGAAAAAGCCATCGCCCGCTATTGGCTTGCCACCGGTTAc                    >  3:81771/1‑148 (MQ=255)                                                                                                                                    cgCCTGCTCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCG‑CACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCTCGAAGAAAAAGCCATCGCCCGCTATTGGCTTGCCACCGGTTACCa                  <  4:106091/148‑1 (MQ=255)                                                                                                                                        tGCGCTTCATCGGTGAAGCGGGTGACCGAAACCCCCGGGCCAAACACTTCGTTGCG‑CACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCTCGAAGAAAAAGCCATCGCCCGCTATTGGCTTGCCACCGGTTACCAGAc               >  5:119183/1‑148 (MQ=255)                                                                                                                                        tGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTACG‑CACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCTCGAAGAAAAAGCCATCGCCCGCTATTGGCTTGCCACCGGTTACCAGAc               >  2:253605/1‑148 (MQ=255)                                                                                                                                         ccgcTTCACCGATGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCG‑CACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCTCGAAGAAAAAGCCATCGCCCGCTATTGGCTTGCCACCGGTTACCAGACg              >  1:337324/2‑148 (MQ=255)                                                                                                                                         ccgcTTCACCGATGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCG‑CACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCTCGAAGAAAAAGCCATCGCCCGCTATTGGCTTGCCACCGGTTACCAGACg              >  2:41543/2‑148 (MQ=255)                                                                                                                                                aTCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCG‑CACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCTCGAAGAAAAAGCCATCGCCCGCTATTGGCTTGCCACCGGTTACCAGACGAATGTGc       <  4:113074/148‑1 (MQ=255)                                                                                                                                                  cGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCG‑CACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCTCGAAGAAAAAGCCATCGCCCGCTATTGGCTTGCCACCGGTTACCAGACGAATGTGCGGc    >  1:161666/1‑149 (MQ=255)                                                                                                                                                  cGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCG‑CACGATCTCGTAGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCTCGAAGAAAAAGCCATCGCCCGCTATTGGCTTGCCACCGGTTACCAGACGAATGTGCGGc    >  7:24752/1‑149 (MQ=255)                                                                                                                                                    gTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCG‑CACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCTCGAAGAAAAAGCCATCGCCCGCTATTGGCTTGCCACCGGTTACCAGACGAATGTGCGGCTg  >  3:226780/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                      ACCAGCATGAAGTGGGTATTCACCCAGGTGCAGCCATACTGCAAGCGCGCCGCCAGGCGGTGGGCGCGGCCGGTGTCGCGGGTCCATACCGATGAGGCCAGGCCGTAGTTGGAGTCGTTGGCCCATTCCAGTGCCTGCGCTTCATCGGTGAAGCGGGTGACCGAAACCACCGGGCCAAACACTTCGTTGCG‑CACGATCTCGTCGTCCTGCAGGGCGTCGGCCAGCACGGTAGGCTCGAAGAAAAAGCCATCGCCCGCTATTGGCTTGCCACCGGTTACCAGACGAATGTGCGGCTG  >  NC_002947/1685463‑1685758

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