Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	428,854	Δ1 bp	7.0%	coding (153/957 nt)	PP_0351 ←	transmembrane sensor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	428,854	0	G	.	7.0%	232.0 / 9.8	71	coding (153/957 nt)	PP_0351	transmembrane sensor
Reads supporting (aligned to +/- strand):  ref base G (35/31);  new base . (2/3);  total (37/34)
Fisher's exact test for biased strand distribution p-value = 6.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.19e-01

GACCCATGCCAGCAGCGCGCCGAGCAGAAACGTGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCGCTGCTCTCGCCAAGGCTGAGTTTCCAGGCAATCGCGGCTTCCAGGACCCGGGACGAAACCGGGAGG  >  NC_002947/428707‑428998                                                                                                                                                    |                                                                                                                                                  gACCCATGCCAGCAGCGCGCCGAGCAGAAACGTGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGa                                                                                                                                                  <  8:344334/149‑1 (MQ=255)    ccATGCCAGCAGCGCGCCGAGCAGAAACGTGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCt                                                                                                                                               >  5:126697/1‑149 (MQ=255)      aTGCCAGCAGCGCGCCGAGCAGAAACGTGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑gg                                                                                                                                            >  6:364602/1‑149 (MQ=255)        gCCAGCAGCGCGCCGAGCAGAAACGTGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑Ggcg                                                                                                                                          >  1:380057/1‑149 (MQ=255)           agcagcGCGCCGAGCAGAAACGTGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCa                                                                                                                                       >  8:295185/1‑149 (MQ=255)            gcagcGCGCCGAGCAGAAACGTGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCTCCAg                                                                                                                                      <  5:123797/149‑1 (MQ=255)                cgcgcCGAGCAGAAACGTGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCAc                                                                                                                                  >  1:272915/1‑149 (MQ=255)                cgcgcCGAGCAGAAACGTGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCAc                                                                                                                                  <  7:295185/149‑1 (MQ=255)                cgcgcCGAGCAGAAACGTGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCC‑AGCTCGGCGCCatggac                                                                                                                                  <  3:87675/149‑6 (MQ=255)                     cGAGCAGAAACGTGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑Ggcg                                                                                                                                          >  8:106535/1‑136 (MQ=255)                     cGAGCAGAAACGTGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑Ggcg                                                                                                                                          <  7:106535/136‑1 (MQ=255)                        gCAGAAACGTGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACgtgcgtgc                                                                                                                          >  3:4830/1‑149 (MQ=255)                            aaaCGTGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCtt                                                                                                                      >  6:89491/1‑149 (MQ=255)                               cGTGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGc                                                                                                                   <  2:272915/149‑1 (MQ=255)                                    cGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTggcg                                                                                                              >  5:126304/1‑149 (MQ=255)                                       ccAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTggcggct                                                                                                           <  5:89491/149‑2 (MQ=255)                                              gccgccAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAgc                                                                                                    <  5:364602/149‑1 (MQ=255)                                               ccgccAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAgcg                                                                                                   <  2:380057/149‑1 (MQ=255)                                                 gccAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCa                                                                                                 >  4:183017/1‑149 (MQ=255)                                                   cAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCAt                                                                                                          >  6:255467/1‑138 (MQ=255)                                                   cAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCAt                                                                                                          <  5:255467/138‑1 (MQ=255)                                                       cTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAAtt                                                                                           >  1:317470/1‑149 (MQ=255)                                                        ttGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTc                                                                                          >  1:201149/1‑149 (MQ=255)                                                          gCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCg                                                                                         <  4:377604/148‑1 (MQ=255)                                                           ccACTGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGt                                                                                        <  6:126697/148‑1 (MQ=255)                                                           ccAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGt                                                                                        <  3:133012/148‑1 (MQ=255)                                                             aaTGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCa                                                                                                 <  2:185364/137‑1 (MQ=255)                                                             aaTGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCa                                                                                                 >  1:185364/1‑137 (MQ=255)                                                             aaTGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGc                                                                                     <  6:126304/149‑1 (MQ=255)                                                             aaTGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGc                                                                                     >  5:157633/1‑149 (MQ=255)                                                                      ccgcAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAg                                                                            <  2:317470/149‑1 (MQ=255)                                                                      ccgcAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAg                                                                            >  3:262754/1‑149 (MQ=255)                                                                       cgcAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCC‑AGCTCGGCGCCat                                                                                                                                      >  8:336971/1‑89 (MQ=255)                                                                       cgcAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCC‑AGCTCGGCGCCat                                                                                                                                      <  7:336971/90‑2 (MQ=255)                                                                         cAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCg                                                                         >  3:304178/1‑149 (MQ=255)                                                                           ggCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTg                                                                       >  6:93027/1‑149 (MQ=255)                                                                                gcgcgcGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCgc                                                                   <  4:355334/148‑1 (MQ=255)                                                                                 cgcgcgATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCgctg                                                                 <  6:157633/149‑1 (MQ=255)                                                                                          gcagcaGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCGCTGCTCTCGcc                                                         <  5:93027/148‑1 (MQ=255)                                                                                           cagcagCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCGCTGCTCTCg                                                           <  8:168862/145‑1 (MQ=255)                                                                                           cagcagCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCGCTGCTCTCg                                                           >  7:168862/1‑145 (MQ=255)                                                                                           cagcagCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCGCTGCTCTCGcc                                                         >  5:43144/1‑147 (MQ=255)                                                                                             gcagcGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCGCTGCTCTCGCCAAgg                                                     >  5:321190/1‑149 (MQ=255)                                                                                             gcagcGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGGCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCGCTGCTCTCGCCAAgg                                                     <  3:29405/149‑1 (MQ=255)                                                                                                gcgcCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCGCTGCTCTCGCCAAGGCt                                                   >  1:313388/1‑148 (MQ=255)                                                                                                  gcCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCGCTGCTCTCGCCAAGGCTGAg                                                >  6:306664/1‑149 (MQ=255)                                                                                                   ccTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTc                                                                                          <  1:94893/106‑1 (MQ=255)                                                                                                   ccTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTc                                                                                          >  2:94893/1‑106 (MQ=255)                                                                                                    cTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGt                                                                        >  8:155181/1‑123 (MQ=255)                                                                                                    cTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGt                                                                        <  7:155181/123‑1 (MQ=255)                                                                                                       gCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCGCTGCTCTCGCCAAGGCTGAGTTTcc                                           <  3:183017/149‑1 (MQ=255)                                                                                                       gCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCGCTGCTCTCGCCAAGGCTGAGTTTc                                            >  7:48094/1‑148 (MQ=255)                                                                                                         gTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCGCTGCTCTCGCCAAGGCTGAGTTTCCAg                                         >  3:19582/1‑149 (MQ=255)                                                                                                         gTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCGCTGCTCTCGCCAAGGCTGAGTTTCCAg                                         <  8:48094/149‑1 (MQ=255)                                                                                                         gTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCGCTGCTCTCGCCAAGGCTGAGTTTCCAg                                         >  4:141097/1‑149 (MQ=255)                                                                                                         gTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCGCTGCTCTCGCCAAGGCTGAGTTTCCAg                                         <  4:262754/149‑1 (MQ=255)                                                                                                                  cAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCGCTGCTCTCGCCAAGGCTGAGTTTCCAGGCGATcgcg                                <  4:304178/149‑1 (MQ=255)                                                                                                                  cAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCGCTGCTCTCGCCAAGGCTGAGTTTCCAGGCAATcgcg                                <  5:306664/149‑1 (MQ=255)                                                                                                                   aGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑CGCGCCAGGCCCGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCGCTGCTCTCGCCAAGGCTGAGTTTCCAGGCAATCGCgg                               >  2:365678/1‑149 (MQ=255)                                                                                                                     cTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCGCTGCTCTCCCCAAGGCTGAGTTTCCAGGCAATCGCGGCt                             >  6:182944/1‑149 (MQ=255)                                                                                                                        ccGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCGCTGCTCTCGCCAAGGCTGAGTTTCCAGGCAATCGCGGCTTc                           <  2:201149/148‑1 (MQ=255)                                                                                                                        ccGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCGCTGCTCTCGCCAAGGCTGAGTTTCCAGGCAATCGCGGCTTc                           >  4:365973/1‑148 (MQ=255)                                                                                                                        ccGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCAAGAATTCGTTGCGTTCATCAGGCGTGCCGCTGCTCTCGCCAAGGCTGAGTTTCCAGGCAATCGCGGCTTc                           >  7:247560/1‑148 (MQ=255)                                                                                                                          gcgcTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCGCTGCTCTCGCCAAGGCTGAGTTTCCAGGCAATCGCGGCTTCCAg                        <  6:321190/149‑1 (MQ=255)                                                                                                                           cgcTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGtt                                                                                       <  4:127070/85‑1 (MQ=255)                                                                                                                           cgcTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGtt                                                                                       >  3:127070/1‑85 (MQ=255)                                                                                                                                    gTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCGCTGCTCTCGCCAAGGCTGAGTTT‑CAGGTAATCGCGGCTTCCAGGACCCGGGAc              >  1:281701/1‑148 (MQ=255)                                                                                                                                          gcctGGCGCC‑AGCTCGGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCGCTGCTCTCGCCAAGGCTGAGTTTCCAGGCAATCGCGGCTTCCa                         <  3:150526/129‑1 (MQ=255)                                                                                                                                          gcctGGCGCC‑AGCTCGGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCGCTGCTCTCGCCAAGGCTGAGTTTCCAGGCAATCGCGGCTTCCa                         >  4:150526/4‑132 (MQ=255)                                                                                                                                            aTGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCGCTGCTCTCGCCAAGGCTGAGTTTCCAGGCAATCGCGGCTTCCAGGACCCGGGACGAAACCgg      >  5:308686/1‑149 (MQ=255)                                                                                                                                                cgcCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCGCTGCTCTCGCCAAGGCTGAGTTTCCAGGCAATCGCGGCTTCCAGGACCCGGGACGAAACCGGGAgg  <  4:19582/149‑1 (MQ=255)                                                                                                                                                    gAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCGCTGCTCTCGCCaa                                                       >  3:34085/1‑92 (MQ=255)                                                                                                                                                    gAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCGCTGCTCTCGCCaa                                                       <  4:34085/92‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                  GACCCATGCCAGCAGCGCGCCGAGCAGAAACGTGCCGGCCAGCCCGCCGCCAACCTTGCCAATGCGTCGCCGCAGGCTGGCGCGCGATTGCAGCAGCGCCTGGCGTGCCGGGCCAGCTGCCGCGCTAACGCGTTGGTCCATGGCGCCGAGCT‑GGCGCCAGGCACGTGCGTGCTCTTCGCTGGCGGCATGCCAGCGCATGAATTCGTTGCGTTCATCAGGCGTGCCGCTGCTCTCGCCAAGGCTGAGTTTCCAGGCAATCGCGGCTTCCAGGACCCGGGACGAAACCGGGAGG  >  NC_002947/428707‑428998

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑