Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,426,086	G→C	9.7%	G208R (GGC→CGC)	PP_4765 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,426,086	0	G	C	9.7%	213.8 / 12.4	72	G208R (GGC→CGC)	PP_4765	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (36/29);  new base C (3/4);  total (39/33)
Fisher's exact test for biased strand distribution p-value = 6.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.63e-01

AGGTGAAAGCCTCGGTCGTGGTGCTGGCGCTGGGTGGTGCCAGCTGGGCGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAGGCGTGGAAGGCAGCCTGGTGTATGCCTGGTCGGCAC  >  NC_002947/5425940‑5426233                                                                                                                                                   |                                                                                                                                                    aGGTGAAAGCCTCGGTCGTGGTGCTGGCGCTGGGTGGTGCCAGCTGGGCGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCggg                                                                                                                                                    <  4:291547/148‑1 (MQ=255)   gTGAAAGCCTCGGTCGTGGTGCTGGCGCTGGGTGGTGCCAGCTGGGCGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTg                                                                                                                                                 <  3:186019/149‑1 (MQ=255)    tGAAAGCCTCGGTCGTGGTGCTGGCGCTGGGTGGTGCCAGCTGGGCGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTg                                                                                                                                                 >  1:162825/1‑148 (MQ=255)            tCGGTCGTGGTGCTGGCGCTGGGTGGTGCCAGCTGGGCGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAgcgcgc                                                                                                                                         >  3:255372/1‑148 (MQ=255)             cGGTCGTGGTGCTGGCGCTGGGTGGTGCCAGCTGGGCGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTa                                                                                                                                       <  4:255372/149‑1 (MQ=255)               gTCGGGGGGCTGGCGCGGGGTGGTGCCAGCTGTGCGCGCCTGGGGTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCGGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACt                                                                                                                                     <  3:229745/149‑1 (MQ=255)                   tggtgCTGGCGCTGGGTGGTGCCAGCTGGGCGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGCGCTCCAGCCCGCTACttcaa                                                                                                                                 <  6:71403/149‑5 (MQ=255)                     gtgCTGGCGCTGGGTGGTGCCAGCTGGGCGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCggg                                                                                                                                                    <  4:144466/128‑1 (MQ=255)                     gtgCTGGCGCTGGGTGGTGCCAGCTGGGCGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCggg                                                                                                                                                    >  3:144466/1‑128 (MQ=255)                     gtgCTGGCGCTGGGTGGTGCCAGCTGGGCGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAgg                                                                                                                                <  7:90269/148‑1 (MQ=255)                      tgCTGGCGCTGGGTGGTGCCAGCTGGGCGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGa                                                                                                                               >  4:78163/1‑148 (MQ=255)                      tgCTGGCGCTGGGTGGTGCCAGCTGGGCGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGa                                                                                                                               >  7:2396/1‑148 (MQ=255)                         tGGCGCTGGGTGGTGCCAGCTGGGCGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAg                                                                                                                           >  7:95588/1‑149 (MQ=255)                             gcTGGGTGGTGCCAGCTGGGCGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCgg                                                                                                                                                     >  2:362091/1‑119 (MQ=255)                             gcTGGGTGGTGCCAGCTGGGCGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTcg                                                                                                                       >  5:19149/1‑149 (MQ=255)                               tGGGTGGTGCCAGCTGGGCGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTcgcg                                                                                                                     >  8:62727/1‑149 (MQ=255)                               tGGGTGGTGCCAGCTGGGCGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGCGCTCCAGCCCGCTACttcaaaa                                                                                                                               <  2:161791/139‑7 (MQ=255)                               tGGGTGGTGCCAGCTGGGCGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGCGCTCCAGCCCGCTACttcaaaa                                                                                                                               >  1:161791/1‑133 (MQ=255)                               tGGGTGGTGCCAGCTGGGAGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGAAACTGCGGTTTTGAAGTAGCGGGCTGGCGCGCGCTACTCAAGGACAAGTTcgcg                                                                                                                     >  8:200820/1‑149 (MQ=255)                                      tgCCAGCTGGGCGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCAc                                                                                                              >  2:184933/1‑149 (MQ=255)                                      tgCCAGCTGGGCGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCAc                                                                                                              <  6:84819/149‑1 (MQ=255)                                               ggCGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTgaaga                                                                                                     <  2:154823/149‑1 (MQ=255)                                                gcgcgcCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTgaagaa                                                                                                    <  1:184933/149‑1 (MQ=255)                                                gcgcgcCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTgaagaa                                                                                                    <  6:131037/149‑1 (MQ=255)                                                          gTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCt                                                                                           >  3:183475/1‑148 (MQ=255)                                                               gACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGt                                                                                     >  5:36851/1‑149 (MQ=255)                                                                aCGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTc                                                                                    >  8:37893/1‑149 (MQ=255)                                                                      gCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCggg                                                                               >  2:283543/1‑148 (MQ=255)                                                                       cTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCa                                                                             <  4:183475/149‑1 (MQ=255)                                                                          gCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAgcg                                                                          >  2:363728/1‑149 (MQ=255)                                                                          gCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAg                                                                            <  7:62727/147‑1 (MQ=255)                                                                           cAACCGCTGCTGGCAGAACGGGGTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGAGGTTTTGAAGTAGCGGGCTGGATAGCGATACTCAATGACAAGTTCGCGTGCTCACAGCTGAAGAAAATCTCCCTAAGCTTCACCTGCAgcg                                                                          >  7:135939/1‑148 (MQ=255)                                                                             aCCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCAc                                                                       >  5:3600/1‑149 (MQ=255)                                                                                  tgctgGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACcgcg                                                                   >  5:167242/1‑148 (MQ=255)                                                                                    ctgGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCgcgc                                                                  >  8:300424/1‑147 (MQ=255)                                                                                     tgGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAg                                                               >  6:98315/1‑149 (MQ=255)                                                                                       gCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAggg                                                             >  4:211878/1‑149 (MQ=255)                                                                                        cAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGc                                                            <  1:363728/149‑1 (MQ=255)                                                                                          gAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGa                                                          <  8:95588/149‑1 (MQ=255)                                                                                          gAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGa                                                          <  7:24469/149‑1 (MQ=255)                                                                                           aaCGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAg                                                         >  6:120265/1‑149 (MQ=255)                                                                                                  tgtgGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAg                                                         <  8:279405/142‑1 (MQ=255)                                                                                                  tgtgGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAg                                                         >  7:279405/1‑142 (MQ=255)                                                                                                   gtgGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATtct                                                 >  2:312843/1‑149 (MQ=255)                                                                                                   gtgGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATtc                                                  >  6:18232/1‑148 (MQ=255)                                                                                                    tgGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATtct                                                 >  1:249500/1‑148 (MQ=255)                                                                                                        tatCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCAcc                                             >  6:281392/1‑148 (MQ=255)                                                                                                         atCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCAtt                                                   <  8:121150/141‑1 (MQ=255)                                                                                                         atCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCAtt                                                   >  7:121150/1‑141 (MQ=255)                                                                                                            tcGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACTTCGCCCTCAGCGTCCCggg                                                                               >  5:183568/1‑110 (MQ=255)                                                                                                            tcGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACTTCGCCCTCAGCGTCCCggg                                                                               <  6:183568/110‑1 (MQ=255)                                                                                                              gCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAg                                      >  8:148646/1‑149 (MQ=255)                                                                                                               ccTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCAtt                                                   >  5:3210/1‑135 (MQ=255)                                                                                                               ccTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCAtt                                                   <  6:3210/135‑1 (MQ=255)                                                                                                                cTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGc                                                            <  6:340197/125‑1 (MQ=255)                                                                                                                cTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGc                                                            >  5:340197/1‑125 (MQ=255)                                                                                                                   tGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAGGCGt                                  <  5:98315/148‑1 (MQ=255)                                                                                                                   tGCAGCCGAGCAACTGAGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGGGAGTTCATTCTCACCGCTCAAGGCGt                                  <  8:135939/148‑1 (MQ=255)                                                                                                                   caaaggcGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAGGCGt                                  <  6:88220/142‑1 (MQ=255)                                                                                                                    gCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAGGCGTg                                 >  4:268709/1‑148 (MQ=255)                                                                                                                     cAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAGGCGTGGa                               <  6:36851/149‑1 (MQ=255)                                                                                                                          gAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAGGCGTGGAAGGCa                          <  1:312843/149‑1 (MQ=255)                                                                                                                                 tGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCCCTCAAGGCGTGGAAGGCAGCCTGgt                   >  2:347235/1‑149 (MQ=255)                                                                                                                                   cGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAGGCGTGGAAGGCAGCCTGgt                   <  6:167242/147‑1 (MQ=255)                                                                                                                                    ggTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAGGCGTGGAAGGCAGCCTGGTGTa                >  8:315671/1‑149 (MQ=255)                                                                                                                                    ggTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAGGCGTGGAAGGCAGCCTGGTGTa                <  5:281392/149‑1 (MQ=255)                                                                                                                                     gTTTTGAAGTAGCGGTCTGGAGCGCGGTACTCAAGCAGAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTGACCGCTCAAGGCGTGGAAGGCAGGCTGGTGTAt               <  8:51045/149‑1 (MQ=255)                                                                                                                                      ttttGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAGGCGTGGAAGGCAGCCTGGTGTATg              >  4:94604/1‑149 (MQ=255)                                                                                                                                        cttgaGTAGCGCGCTCCAGCCCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAGGCGTGGaa                              >  7:294371/5‑131 (MQ=255)                                                                                                                                        cttgaGTAGCGCGCTCCAGCCCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAGGCGTGGaa                              <  8:294371/127‑1 (MQ=255)                                                                                                                                         ttgaGTAGCGCGCTCCAGCCCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAggg                                                             <  4:289244/96‑1 (MQ=38)                                                                                                                                         ttgaGTAGCGCGCTCCAGCCCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAggg                                                             >  3:289244/4‑99 (MQ=38)                                                                                                                                         tGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAGGCGTGGAAGGCAGCCTGGTGTATGcc            <  2:162825/148‑1 (MQ=255)                                                                                                                                           aaGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAGGCGTGGAAGGCAGCCTGGTGTATGCCTg          <  1:340101/148‑1 (MQ=255)                                                                                                                                             gTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAGGCGTGGAAGGCAGCCTGGTGTATGCCTGGTc       <  1:283543/149‑1 (MQ=255)                                                                                                                                              tAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCCGCCAAGGGCGAGTTCATTCTCCCCGCTCAAGCCGTGGAAGGCAGCCTGGTGTAt               >  7:195928/1‑140 (MQ=255)                                                                                                                                                  gggCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAGGCGTGGAAGGCAGCCTGGTGTATGCCTGGTCGGCAc  <  7:37893/149‑1 (MQ=255)                                                                                                                                                  gggCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAGGCGTGGAAGGCAGCCTGGTGTATGCCTGGTCGGCAc  <  1:275487/149‑1 (MQ=255)                                                                                                                                                  gggCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACAGCTGAAGAAAATCGACCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAGGCGTGGAAGGCAGCCTGGTGTATGCCTGGTCGGCAc  >  8:115495/1‑149 (MQ=255)                                                                                                                                                   ggCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAGGCGTGGAAGGCAGc                        >  3:238271/1‑126 (MQ=255)                                                                                                                                                   ggCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAGGCGTGGAAGGCAGc                        <  4:238271/126‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                    AGGTGAAAGCCTCGGTCGTGGTGCTGGCGCTGGGTGGTGCCAGCTGGGCGCGCCTGGGTTCTGACGGCAGCTGGCAACCGCTGCTGGCAGAACGGGCTGTGGATATCTCGCCTTTGCAGCCGAGCAACTGCGGTTTTGAAGTAGCGGGCTGGAGCGCGCTACTCAAGGACAAGTTCGCGGGCGCACCGCTGAAGAACATCGCCCTCAGCGTCCCGGGCAGCGCACCGCGCAAGGGCGAGTTCATTCTCACCGCTCAAGGCGTGGAAGGCAGCCTGGTGTATGCCTGGTCGGCAC  >  NC_002947/5425940‑5426233

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