Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,054,198	Δ1 bp	7.2%	coding (772/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,198	0	A	.	7.2%	210.7 / 10.7	69	coding (772/1320 nt)	PP_3573	monooxygenase
Reads supporting (aligned to +/- strand):  ref base A (31/33);  new base . (3/2);  total (34/35)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.42e-01

CAACCGCGACTTCTGCTGGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGG‑ACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATCGAACGCAACGGCCTGGACCTGCC  >  NC_002947/4054050‑4054346                                                                                                                                                       |                                                                                                                                                     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                                                                                                                               >  2:42164/1‑149 (MQ=255)                        tGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGG‑ACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCg                                                                                                                               >  3:401589/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                                                                                                                        <  8:264271/146‑1 (MQ=255)                                  gggCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGG‑ACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGaaa                                                                                                                        >  7:264271/1‑146 (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                                                                                                                                  >  3:411443/1‑130 (MQ=255)                                        gTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGG‑ACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAg                                                                                                                                  <  4:411443/130‑1 (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                                                                                                        >  4:72450/1‑149 (MQ=255)                                               gcgGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGG‑ACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCtt                                                                                                        >  8:235699/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                          <  8:337253/149‑1 (MQ=255)                                                                                                                              cGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACAt                          <  6:294816/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)                                                                                                                                                      gACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATCGAACGCAACGGCCTGGACCTGc   >  6:14673/1‑149 (MQ=255)                                                                                                                                                       aCCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATCGAACGCAACGGCCTGGACCTGcc  <  1:275418/149‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                     CAACCGCGACTTCTGCTGGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGG‑ACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATCGAACGCAACGGCCTGGACCTGCC  >  NC_002947/4054050‑4054346

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