Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,296,137	Δ1 bp	7.3%	coding (393/600 nt)	PP_2023 ←	glutathione S‑transferase family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,296,137	0	C	.	7.3%	169.4 / 8.2	55	coding (393/600 nt)	PP_2023	glutathione S‑transferase family protein
Reads supporting (aligned to +/- strand):  ref base C (27/24);  new base . (2/2);  total (29/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.08e-01

CGCACACGCGGATAAGCCGCCTCGGCGAACCACTGCGGTTCGACCCCGGCAAACTGGCGCATCAGGGGCAGCAACGCGGCATCGGCCATGCTCGGGTGGTCGG‑CCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCT‑GCTGGCGGTAATGCTCGCGGGAATGCTCCGGGTAACGTTCGGCATACTTGTACAGGTTCACGTGCGCCTTGAACGTGGTGTCATTGCGCGCAATCAGCGTTTCAGCCTGCTGCGCAGCCGCAGGGTTGGCCTGTA  >  NC_002947/2295988‑2296285                                                                                                                                                       |                                                                                                                                                      cGCACACGCGGATAAGCCGCCTCGGCGAACCACTGCGGTTCGACCCCGGCAAACTGGCGCATCAGGGGCAGCAACGCGGCATCGGCCATGCTCGGGTGGTCGG‑CCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCgccagc                                                                                                                                                       <  2:1178438/150‑1 (MQ=255)   cacaCGCGGATAAGCCGCCTCGGCGAACCACTGCGGTTCGACCACGGCAAACTGGCGCATCAGGGGCAGCAACGCGGCATCGGCCATGCTCGGGTGGTCGG‑CCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCgccagcca                                                                                                                                                     <  1:1036936/150‑1 (MQ=255)    acacGCGGATAAGCCGCCTCGGCGAACCACTGCGGTTCGACCCCGGCAAACTGGCGCATCAGGGGCAGCACAGCGGGAGCGGCCATGCTCCTTTGGTCGG‑CCAGCAGGTACGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGccagccag                                                                                                                                                    <  1:147518/150‑1 (MQ=255)    acacGCGGATAAGCCGCCTCGGCGAACCACTGCGGTTCGACCCCGGCAAACTGGCGCATCAGGGGCAGCAACGCGGCATCGGCCATGCTCGGGTGGTCGG‑CCAGCAGGTAAGCCCGCCCAGCCAGCAAGCCTTCGAGCTCCGccagccag                                                                                                                                                    <  1:933346/150‑1 (MQ=255)        gcgGATAAGCCGCCTCGGCGAACCACTGCGGTTCGACCCCGGCAAACTGGCGCATCAGGGGCAGCAACGCGGCATCGGCCATGCTCGGGTGGTCGG‑CCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCt                                                                                                                                                >  2:1005863/1‑150 (MQ=255)         cgGATAAGCCGCCTCGGCGAACCACTGCGGTTCGACCCCGGCAAACTGGCGCATCAGGGGCAGCAACGCGGCATCGGCCATGCTCGGGTGGTCGG‑CCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTc                                                                                                                                               <  1:1005863/150‑1 (MQ=255)            aTAAGCCGCCTCGGCGAACCACTGCGGTTCGACCCCGGCAAACTGGCGCATCAGGGGCAGCAACGCGGCATCGGCCATGCTCGGGTGGTCGG‑CCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGc                                                                                                                                            >  2:3942/1‑150 (MQ=255)                gccgccTCGGCGAACCACTGCGGTTCGACCCCGGCAAACTGGCGCATCAGGGGCAGCAACGCGGCATCGGCCATGCTCGGGTGGTCGG‑CCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCct‑gc                                                                                                                                       >  2:672670/1‑150 (MQ=255)                 ccgccTCGGCGAACCACTGCGGTTCGACCCCGGCAAACTGGCGCATCAGGGGCAGCAACGCGGCATCGGCCATGCTCGGGTGGTCGG‑CCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCct‑gct                                                                                                                                      >  1:64421/1‑150 (MQ=255)                 ccgccTCGGCGAACCACTGCGGTTCGACCCCGGCAAACTGGCGCATCAGGGGCAGCAACGCGGCATCGGCCATGCTCGGGTGGTCGG‑CCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCct‑gct                                                                                                                                      >  2:393264/1‑150 (MQ=255)                 ccgccTCGGCGAACCACTGCGGCTCGACCCCGGCAAACTGGCGCATCAGGGGCAGCAACGCGGCATCGGCCATGCTCGGGTGGTCGG‑CCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCct‑gct                                                                                                                                      >  1:567132/1‑150 (MQ=255)                   gccTCGGCGAACCACTGCGGTTCGACCCCGGCAAACTGGCGCATCAGGGGCAGCAACGCGGCATCGGCCATGCTCGGGTGGTCGG‑CCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCT‑GCTgg                                                                                                                                    <  2:862469/150‑1 (MQ=255)                       cGGCGAACCACTGCGGTTCGACCCCGGCAAACTGGCGCATCAGGGGCAGCAACGCGGCATCGGCCATGCTCGGGTGGTCGG‑CCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAG‑CAGGCCGAGGCCTGGCTGGCgga                                                                                                                                <  2:598076/150‑2 (MQ=255)                        ggCGAACCACTGCGGTTCGACCCCGGCAAACTGGCGCATCAGGGGCAGCAACGCGGCATCGGCCATGCTCGGGTGGTCGG‑CCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCT‑GCTGGCGGTa                                                                                                                               >  1:524889/1‑150 (MQ=255)                        ggCGAACCACTGCGGTTCGACCCCGGCAAACTGGCGCATCAGGGGCAGCAACGCGGCATCGGCCATGCTCGGGTGGTCGG‑CCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCT‑GCTGGCGGTa                                                                                                                               <  1:393264/150‑1 (MQ=255)                            aaCCACTGGGGTTCGACCCCGGCAAACTGGCGCATCAGGGGCAGCAACGCGGCATCGGCCATGCTCGGGTGGTCGG‑CCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCT‑GCTGGCGGTAATGc                                                                                                                           <  2:524889/150‑1 (MQ=255)                                  tGCGGTTCGACCCCGGCAAACTGGCGCATGAGGGGCAGCAACGCGGCATCGGCCATGCTCGGGTGGTCGG‑CCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCT‑GCTGGCGGTAATGCTCGCgg                                                                                                                     >  1:232174/1‑150 (MQ=255)                                  tGCGGTTCGACCCCGGCAAACTGGCGCATCAGGGGCAGCAACGCGGCATCGGCCATGCTCGGGTGGTCGG‑CCAGCAGGTCAAAACGCCCCGCCAGCAAGCCTTCGAGCTCCGACCGCCAGGCGGCGATGT‑GCTGGCGGTAAGGCTCGCgg                                                                                                                     >  1:327448/1‑150 (MQ=255)                                     ggctCGACCCCGGCAAACTGGCGCATCAGGGGCAGCAACGCGGCATCGGCCATGCTCGGGTGGTCGG‑CCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCT‑GCTGGCGGTAATGCTCGCGGGaa                                                                                                                  <  2:567132/147‑1 (MQ=255)                                     ggTTCGACCCCGGCAAACTGGCGCATCAGGGGCAGCAACGCGGCATCGGCCATGCTCGGGTGGTCGG‑CCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCT‑GCTGGCGGTAATGCTCGCGGGaa                                                                                                                  <  1:26519/150‑1 (MQ=255)                                          gaacccGGCAAACTGGCGCATCAGGGGCAGCAACGCGGCATCGGCCATGCTCGGGTGGTCGG‑CCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCT‑GCTGGCGGTAATGCTCGCGGGAATGCTc                                                                                                             >  2:1104145/4‑150 (MQ=255)                                            ccccGGCAAACTGGCGCATCAGGGGCAGCAACGCGGCATCGGCCATGCTCGGGTGGTCGG‑CCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCT‑GCTGGCGGTAATGCTCGCGGGAATGCTCCg                                                                                                           <  2:112084/150‑1 (MQ=255)                                             cccGGCAAACTGGCGCATCAGGGGCAGCAACGCGGCATCGGCCATGCTCGGGTGGTCGG‑CCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCT‑GCTGGCGGTAATGCTCGCGGGAATGCTCCgg                                                                                                          <  2:497975/150‑1 (MQ=255)                                              ccGGCAAACTGGCGCATCAGGGGCAGCAACGCGGCATCGGCCATGCTCGGGTGGTCGG‑CCAGCAGGTAAGCCCGCCCCGCTTGCTATCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCT‑GCTGGCGGTACTGCTCGCGGGAATGCTCCggg                                                                                                         >  1:380132/1‑150 (MQ=255)                                              ccGGCAAACTGGCGCATCAGGGGCAGCAACGAGGGATGGGACATGCTAGGGTGGTCGG‑CCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCT‑GCTGGCGGTAATGCTCGCGGGAATGCTCCggg                                                                                                         >  1:210718/1‑150 (MQ=255)                                                  cAAACTGGCGCATCAGGGGCAGCAACGCGGCATCGGCCATGCTCGGGTGGTCGG‑CCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCT‑GCTGGCGGTAATGCTCGCGGGAATGCTCCGGGTAAc                                                                                                     <  2:223143/150‑1 (MQ=255)                                                  cAAACTGGCGCATCAGGGGCAGCAACGCGGCATCGGCCATGCTCGGGTGGTCGG‑CCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCT‑GCTGGCGGTAATGCTCGCGGGAATGCTCCGGGTAAc                                                                                                     <  1:672670/150‑1 (MQ=255)                                                      cTGGCGCATCAGGGGCAGCAACGCGGCATCGGCCATGCTCGGGTGGTCGG‑CCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCT‑GCTGGCGGTAATGCTCGCGGGAATGCTCCGGGTAACGTTc                                                                                                 <  2:983974/150‑1 (MQ=255)                                                       tGGCGCATCAGGGGCAGCAACGCGGCATCGGCCATGCTTGGGCCGTTGGTCCCGCA‑GCAAGTCGTCCGCGTCATCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCT‑GCTGGCGGTAATGCTCGCGGGAATGCTCCGGGTAACGTTCg                                                                                                >  1:328007/1‑150 (MQ=255)                                                            cATCAGGGGCAGCAACGCGGCATCGGCCATGCTCGGGTGGTCGG‑CCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCT‑GCTGGCGGTAATGCTCGCGGGAATGCTCCGGGTAACGTTCGGCATa                                                                                           >  1:651503/1‑150 (MQ=255)                                                             aTCAGGGGCAGCAACGCGGCATCGGCCATGCTCGGGTGGTCGG‑CCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCT‑GCTGGCGGTAATGCTCGCGGGAATGCTCCGGGTAACGTTCGGCATAc                                                                                          <  2:713914/150‑1 (MQ=255)                                                                 ggggCAGCAACGCGGCATCGGCCATGCTCGGGTGGTCGG‑CCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCT‑GCTGGCGGTAATGCTCGCGGGAATGCTCCGGGTAACGTTCGGCATACTTGt                                                                                      >  1:460376/1‑150 (MQ=255)                                                                    gcagcaACGCGGCATCGGCCATGCTCGGGTGGTCGG‑CCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCT‑GCTGGCGGTAATGCTCGCGGGAATGCTCCGGGTAACGTTCGGCATACTTGTACa                                                                                   <  2:327448/150‑1 (MQ=255)                                                                       gcaACGCGGCATCGGCCATGCTCGGGTGGTCGG‑CCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCT‑GCTGGCGGTAATGCTCGCGGGAATGCTCCGGGTAACGTTCGGCATACTTGTACAGGt                                                                                <  1:467144/150‑1 (MQ=255)                                                                              ggCATCGGCCATGCTCGGGTGGTCGG‑CCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCT‑GCTGGc                                                                                                                                   <  2:436673/92‑1 (MQ=255)                                                                              ggCATCGGCCATGCTCGGGTGGTCGG‑CCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCT‑GCTGGc                                                                                                                                   >  1:436673/1‑92 (MQ=255)                                                                                      ccATGCTCGGGTGGTCGG‑CCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCT‑GCTGGCGGTAATGCTCGCGGGAATGCTCCGGGTAACGTTCGGCATACTTGTACAGGTTCACGTGCGCCTTGa                                                                 >  2:24256/1‑150 (MQ=255)                                                                                       cATGCTCGGGTGGTCGG‑CCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCT‑GCTGGCGGTAAt                                                                                                                             <  1:295841/89‑1 (MQ=255)                                                                                       cATGCTCGGGTGGTCGG‑CCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCT‑GCTGGCGGTAAt                                                                                                                             >  2:295841/1‑89 (MQ=255)                                                                                       cATGCTCGGGTGGTCGG‑CCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCT‑GCTGGCGGTAATGCTCGCGGGAATGCTCCGGGTAACGTTCGGCATACTTGTACAGGTTCACGTGCGCCTTGaa                                                                <  1:859191/150‑1 (MQ=255)                                                                                       cATGCTCGGGTGGTCGG‑CCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCT‑GCTGGCGGTAATGCTCGCGGGAATGCTCCGGGTAACGTTCGGCATACTTGTACAGGTTCACGTGCGCCTTGaa                                                                >  1:842041/1‑150 (MQ=255)                                                                                       cATGCTCGAGTGGTCCG‑ACAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCT‑GCTGGCGGTAATGCTCGCGGGAATGCTCCGGGTAACGTTCGGCATACTTGTACAGGTTCACGTGCGCCTTGaa                                                                >  2:596092/1‑150 (MQ=255)                                                                                            tCGGGTGGTCGG‑CCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAG‑CAGGCCGAGGCCtggctggc                                                                                                                                   <  2:655376/78‑1 (MQ=25)                                                                                            tCGGGTGGTCGG‑CCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAG‑CAGGCCGAGGCCtggctggc                                                                                                                                   >  1:655376/1‑78 (MQ=25)                                                                                             cGGGTGGTCGG‑CCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCT‑GCTGGCGGTAATGCTCGCGGGAATGCTCCGGGTAACGTTCGGCATACTTGTACAGGTTCACGTGCGCCTTGAACgtggt                                                          <  2:8624/150‑1 (MQ=255)                                                                                                gtggtCGG‑CCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCT‑GCTGGCGGTAATGCTCGCGGGAATGCTCCGGGTAACGTTCGGCATACTTGTACAgg                                                                                 >  2:178097/1‑124 (MQ=255)                                                                                                gtggtCGG‑CCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCT‑GCTGGCGGTAATGCTCGCGGGAATGCTCCGGGTAACGTTCGGCATACTTGTACAgg                                                                                 <  1:178097/124‑1 (MQ=255)                                                                                                      gg‑CCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCT‑GCTGGCGGTAATGCTCGCGGGAATGCTCCGGGTAACGTTCgg                                                                                               <  1:593855/104‑1 (MQ=255)                                                                                                      gg‑CCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCT‑GCTGGCGGTAATGCTCGCGGGAATGCTCCGGGTAACGTTCgg                                                                                               >  2:593855/1‑104 (MQ=255)                                                                                                       g‑CCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCT‑GCTGGCGGTAATGCTCGCGGGAATGCTCCGGGTAACGTTCGGCATACTTGTACAGGTTCACGTGCGCCTTGAACGTGGTGTCATTgcgc                                                <  2:380132/150‑1 (MQ=255)                                                                                                           agcagGTAAGCCCGCCCCGCCAGCAAGCCTTCGATCTCCGCCAGCCAGGCCTCGGCCT‑GCTGGCGGTAATGCTCGCGGGAATGCTCCGGGTAACGTTCGGCATACTTGTACAGGTTCACGTGCGCCTTGAACgtggt                                                          >  2:218704/1‑137 (MQ=255)                                                                                                           agcagGTAAGCCCGCCCCGCCAGCAAGCCTTCGATCTCCGCCAGCCAGGCCTCGGCCT‑GCTGGCGGTAATGCTCGCGGGAATGCTCCGGGTAACGTTCGGCATACTTGTACAGGTTCACGTGCGCCTTGAACgtggt                                                          <  1:218704/137‑1 (MQ=255)                                                                                                                gTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTTGGCCT‑GCTGGCGGTAATGCTCGCGGGAATGCTCCGGGTAACGTTCGGCATACTTGTACAGGTTCACGTGCGCCTTGAACGTGGTGTCATTGCGCGCAATCAg                                        >  1:468740/1‑150 (MQ=255)                                                                                                                       cgccccgccAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCT‑GCTGGCGGTAATGCTCGCGGGAATGCTCCGGGTAACGTTCGGCATACTTGTACAGGTTCACGTGCGCCTTGAACGTGGTGTCATTGCGCGCAATCAGCGTTTCa                                 <  2:210718/150‑1 (MQ=255)                                                                                                                        gccccgccAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCT‑GCTGGCGGTAATGCTCGCGGGAATGCTCCGGGTAACGTTCGGCATACTTGTACAGGTTCACGTGCGCCTTGAACGTGGTGTCATTGCGCGCAATCAGCGTTTCAg                                >  1:263567/1‑150 (MQ=255)                                                                                                                           ccgccAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCT‑GCTGGCGGTAATGCTCGCGGGAATGCTCCGGGTAACGTTCGGCATACTTGTACAGGTTCACGTGCGCCTTGAACGTGGTGTCATTGCGCGCAATCAGCGTTTCAGCct                             >  1:1120613/1‑150 (MQ=255)                                                                                                                                      ccTTCGAGCTCCGCCAGCCAGGCCTCGGCCT‑GCTGGCGGTAATGCTCGCGGGAATGCTCCGGGTAACGTTCGGCATACTTGTACAGGTTCACGTGCGCCTTGAACGTGGTGTCATTGCGCGCAATCAGCGTTTCAGCCTGCTGCGCAGcc                  <  2:571503/150‑1 (MQ=255)                                                                                                                                            attaccGCCAG‑CAGGCCGAGGCCTGGCTGGCGGTAATGCTCGCGGGAATGCTCCGGGTAACGTTCGGCATACTTGTACAGGTTCACGTGCGCCTTGAACGTGGTGTCATTGCGCGCAATCAGCGTTTCAGCCTGCTGCGCAGCCGCAggg            >  2:1155566/5‑150 (MQ=255)                                                                                                                                             gCTCCGCCAGCCAGGCCTCGGCCT‑GCTGGCGGTAATGCTCGCGGGAATGCTCCGGGTAACGTTCGGCATACTTGTACAGGTTCACGTGCGCCTTGAACGTGGTGTCATTGCGCGCAATCAGCGTTTCAGCCTGCTGCGCAGCCGCAGGGt           >  2:507644/1‑150 (MQ=255)                                                                                                                                                  gccagccaGGCCTCGGCCT‑GCTGGCGGTAATGCTCGCGGGAATGCTCCGGGTAACGTTCGGCATACTTGTACAGGTTCACGTGCGCCTTGAACGTGGTGTCATTGCGCGCAATCAGCGTTTCAGcc                              <  2:1029504/126‑1 (MQ=255)                                                                                                                                                  gccagccaGGCCTCGGCCT‑GCTGGCGGTAATGCTCGCGGGAATGCTCCGGGTAACGTTCGGCATACTTGTACAGGTTCACGTGCGCCTTGAACGTGGTGTCATTGCGCGCAATCAGCGTTTCAGcc                              >  1:1029504/1‑126 (MQ=255)                                                                                                                                                      gccagGCCTCGGCCT‑GCTGGCGGTAATGCTCGCGGGAATGCTCCGGGTAACGTTCGGCATACTTGTACAGGTTCACGTGCGCCTTGAACGTGGTGTCATTGCGCGCAATCAGCGTTTCAGCCTGCTGCGCAGCCGCAGGGTTGGCCTGTa  >  1:303446/1‑150 (MQ=255)                                                                                                                                                      gccagGCCTCGGCCT‑GCTGGCGGTAATGCTAGCGGGAATGCTCCGGGTAACGCTCGGCATACTTGTACAGGTTCACGTGCGCCTTGAAAGTGGTGTCATTGTGCGCCATCACCGTTTAAGCCTGCTGCGCAGCCGCAGGGTCGGCCTGTa  <  2:64421/150‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                      CGCACACGCGGATAAGCCGCCTCGGCGAACCACTGCGGTTCGACCCCGGCAAACTGGCGCATCAGGGGCAGCAACGCGGCATCGGCCATGCTCGGGTGGTCGG‑CCAGCAGGTAAGCCCGCCCCGCCAGCAAGCCTTCGAGCTCCGCCAGCCAGGCCTCGGCCT‑GCTGGCGGTAATGCTCGCGGGAATGCTCCGGGTAACGTTCGGCATACTTGTACAGGTTCACGTGCGCCTTGAACGTGGTGTCATTGCGCGCAATCAGCGTTTCAGCCTGCTGCGCAGCCGCAGGGTTGGCCTGTA  >  NC_002947/2295988‑2296285

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑