Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,054,193:1	+T	7.0%	coding (767/1320 nt)	PP_3573 →	monooxygenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,054,193	1	.	T	7.0%	281.1 / 11.4	71	coding (767/1320 nt)	PP_3573	monooxygenase
Reads supporting (aligned to +/- strand):  ref base . (34/32);  new base T (3/2);  total (37/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 = 5.87e-01

CCGCAACCGCGACTTCTGCTGGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGG‑ACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATCGAACGCAACGGCCTGGAC  >  NC_002947/4054047‑4054341                                                                                                                                                     |                                                                                                                                                     ccGCAACCGCGACTTCTGCTGGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGG‑ACTTCCGTGCCCTCGCCCACCAGGG‑Ca                                                                                                                                                    >  5:321187/1‑149 (MQ=255)  cGCAACCGCGACTTCTGCTGGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGG‑ACTTCCGTGCCCTCGCCCACCAGGG‑CAt                                                                                                                                                   >  3:108494/1‑149 (MQ=255)  cGCAACCGCGACTTCTGCTGGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGG‑ACTTCCGTGCCCTCGCCCACCAGGG‑CAt                                                                                                                                                   >  7:187691/1‑149 (MQ=255)    cAACCGCGACTTCTGCTGGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGG‑ACTTCCGTGCCCTCGCCCACCAGGG‑CATGa                                                                                                                                                 >  5:76913/1‑149 (MQ=255)             cTTCTGCTGGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGG‑ACTTCCGTGCCCTCGCCCACCAGGGTCATG‑CCCTGGtgg                                                                                                                                        <  6:172126/149‑3 (MQ=255)                ctgctgGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGG‑ACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCgg                                                                                                                                       <  4:40603/147‑1 (MQ=255)                ctgctgGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGG‑ACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGc                                                                                                                                      <  6:76913/148‑1 (MQ=255)                  gctgGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGG‑ACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCt                                                                                                                                    >  6:351812/1‑148 (MQ=255)                   ctgGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGG‑ACTTCCGTGCCCTCGCCCACCAGGGTCATG‑CCCTGGTGGGcgag                                                                                                                                   <  1:239197/148‑4 (MQ=255)                   ctgGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGG‑ACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTg                                                                                                                                   <  5:163432/148‑1 (MQ=255)                    tggtggCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGG‑ACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGAc                                                                                                                                 >  3:299762/1‑149 (MQ=255)                      gtggCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGG‑ACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGAcc                                                                                                                                <  2:248835/148‑1 (MQ=255)                      gtggCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGG‑ACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGAc                                                                                                                                 <  6:221524/147‑1 (MQ=255)                           tGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGTCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGG‑ACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCg                                                                                                                          >  6:315307/1‑149 (MQ=255)                           tGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGG‑ACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCg                                                                                                                          >  3:401589/1‑149 (MQ=255)                           tGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGG‑ACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCg                                                                                                                          >  2:42164/1‑149 (MQ=255)                               cGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGTACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTc                                                                                                                       <  1:119041/149‑1 (MQ=255)                                    tGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGG‑ACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAAc                                                                                                                 >  5:358362/1‑149 (MQ=255)                                    tGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGG‑ACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAAc                                                                                                                 <  6:77542/149‑1 (MQ=255)                                     gggCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGG‑ACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGaaa                                                                                                                   >  7:264271/1‑146 (MQ=255)                                     gggCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGG‑ACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGaaa                                                                                                                   <  8:264271/146‑1 (MQ=255)                                     gggCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGG‑ACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACg                                                                                                                <  8:187691/149‑1 (MQ=255)                                      ggCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGG‑ACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACgg                                                                                                               >  2:253850/1‑149 (MQ=255)                                        cGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGG‑ACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCg                                                                                                             >  7:337253/1‑149 (MQ=255)                                           gTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGG‑ACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAg                                                                                                                             <  4:411443/130‑1 (MQ=255)                                           gTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGG‑ACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAg                                                                                                                             >  3:411443/1‑130 (MQ=255)                                              ggACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGTCACACCGTGG‑ACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTcgcgc                                                                                                       >  7:350486/1‑149 (MQ=255)                                              ggACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGG‑ACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTcgcgc                                                                                                       >  2:275418/1‑149 (MQ=255)                                               gACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGG‑ACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGaaaa                                                                                                                  <  6:277419/137‑1 (MQ=255)                                               gACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGG‑ACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGaaaa                                                                                                                  >  5:277419/1‑137 (MQ=255)                                                  gcgGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGG‑ACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCtt                                                                                                   >  8:235699/1‑149 (MQ=255)                                                  gcgGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGG‑ACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCtt                                                                                                   >  4:72450/1‑149 (MQ=255)                                                      aaaTTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGG‑ACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAg                                                                                               <  5:190713/149‑1 (MQ=255)                                                            ccAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGG‑ACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAAc                                                                                         >  3:82477/1‑149 (MQ=255)                                                              aaGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGG‑ACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGAc                                                                                            <  1:170914/144‑1 (MQ=255)                                                                        gcgAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGG‑ACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGaa                                                                              <  4:181987/148‑1 (MQ=255)                                                                          gAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGG‑ACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACAt                                                                           <  4:77565/149‑1 (MQ=255)                                                                                gtgACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGG‑ACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATcg                                                                     >  4:207291/1‑149 (MQ=255)                                                                                 tgACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGG‑ACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATcgc                                                                    <  5:351812/149‑1 (MQ=255)                                                                                   aCCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGG‑ACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATcgcg                                                                   <  5:315307/148‑1 (MQ=255)                                                                                   aCCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGG‑ACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCgcgg                                                                  <  1:42164/149‑1 (MQ=255)                                                                                   aCCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGG‑ACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCgcgg                                                                  >  6:362492/1‑149 (MQ=255)                                                                                         gcgGTCAGTGGCGCCCGTGGTGGCCACACCGTGG‑ACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGcgacga                                                            <  4:299762/149‑1 (MQ=255)                                                                                            ggcAGTGGCGCCCGTGGGGGCCACACCGGCG‑ACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACaa                                                         <  6:441702/147‑1 (MQ=255)                                                                                                gTGGCGCCCGTGGTGGCCACACCGTGG‑ACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTAc                                                     <  4:401589/149‑1 (MQ=255)                                                                                                 tGGCGCCCGTGGTGGCCACACCGTGG‑ACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTAcc                                                    >  4:289980/1‑149 (MQ=255)                                                                                                   gcgcCCGTGGTGGCCACACCGTGG‑ACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGA‑‑CAGTCGTTCGAAATCGGCGTCGCGCGCTTCCAGGACAACCTGGTCTTGTACATCTATCGCGGCGACGACAACTACCtcg                                                 <  4:362549/148‑3 (MQ=255)                                                                                                    cgcCCGTGGTGGTCACACCGTGG‑ACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGACGAGAACATCAATCGCGGCGACGACAACTACCTgg                                                 <  8:350486/149‑1 (MQ=255)                                                                                                     gcCCGTGGTGGCCACACCGTGG‑ACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGc                                                >  8:194207/1‑149 (MQ=255)                                                                                                     gcCCGTGGTGGCCACACCGTGG‑ACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGAGGCGACGACAACTACCTGGc                                                >  6:249686/1‑149 (MQ=255)                                                                                                               gCCACACCGTGG‑ACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCACTGCTGGAt                                      >  8:382274/1‑149 (MQ=255)                                                                                                                 cacaCCGTGG‑ACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGc                                    >  5:397035/1‑149 (MQ=255)                                                                                                                  acacCGTGG‑ACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCg                                   >  8:356623/1‑149 (MQ=255)                                                                                                                   cacCGTGG‑ACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGTCTACGACAACTACCTGGCCCTGCTGGTTGCgg                                  <  4:118542/149‑1 (MQ=255)                                                                                                                    ccggggg‑ACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGc                                 <  7:270186/144‑1 (MQ=255)                                                                                                                     ccGTGG‑ACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAAAGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGcc                                >  6:144209/1‑149 (MQ=255)                                                                                                                       gTGG‑ACTTCAGTGCCCTCGACCACGAGGG‑CATGACCCTGTTCGGCCTGACCCAGTCTTTCGAATACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGa                              <  6:358362/149‑1 (MQ=255)                                                                                                                        tGG‑ACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGAt                             <  5:362492/149‑1 (MQ=255)                                                                                                                          g‑ACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGc                           >  5:179140/1‑149 (MQ=255)                                                                                                                             cTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAACACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCt                         <  6:321187/149‑1 (MQ=255)                                                                                                                              ttCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTa                        >  4:314508/1‑149 (MQ=255)                                                                                                                                ccGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTAc                       <  5:187383/148‑1 (MQ=255)                                                                                                                                 cGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACAt                     <  6:294816/149‑1 (MQ=255)                                                                                                                                 cGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACAt                     <  8:337253/149‑1 (MQ=255)                                                                                                                                   tGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATCg                   <  3:289980/149‑1 (MQ=255)                                                                                                                                       cTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATCGAACg               >  7:37501/1‑149 (MQ=255)                                                                                                                                          gCCGACCAGGGTCATG‑CCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCTGCCGATGACTACATCGAACGCaa            >  8:139314/1‑149 (MQ=255)                                                                                                                                          gCCGACCAGGGTCATG‑CCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATCGAACGCaa            >  2:401174/1‑149 (MQ=255)                                                                                                                                           cccACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATCGAACGCAAc           >  4:426342/1‑149 (MQ=255)                                                                                                                                             gaccaGGGTCATG‑CCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATCGAACGCAACgg         >  1:318973/2‑149 (MQ=255)                                                                                                                                             caccaGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATCGAACGCAACgg         >  5:48656/1‑149 (MQ=255)                                                                                                                                                    g‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATCGAACGCAACGGCCTGGAc  <  6:48656/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                     CCGCAACCGCGACTTCTGCTGGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGG‑ACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATCGAACGCAACGGCCTGGAC  >  NC_002947/4054047‑4054341

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