Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	422,067	C→G	6.5%	L627L (CTG→CTC) ‡	PP_0348 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	422,067	0	C	G	6.5%	195.9 / 4.0	62	L627L (CTG→CTC) ‡	PP_0348	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (29/29);  new base G (2/2);  total (31/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.35e-01

GATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAGGCAGCAATGGCAGGCGCTGGGCTACGTACGC  >  NC_002947/421920‑422213                                                                                                                                                    |                                                                                                                                                   gATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCa                                                                                                                                                   <  6:303901/149‑1 (MQ=255)   ttGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAgc                                                                                                                                                 <  6:243299/149‑1 (MQ=255)           cgtccAATCGACAGCACCAGCAACGCGCGGATGCGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCCTGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCggg                                                                                                                                          <  5:304601/145‑1 (MQ=255)           cGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCCGGGCGCTGATCGCCCAGGc                                                                                                                                         <  5:249469/149‑1 (MQ=255)            ggCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCCGGGCGCTGATCGCCCAGGcc                                                                                                                                        <  4:98897/149‑1 (MQ=255)            ggCCAATCGACAGCACCAGCAAAGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGcc                                                                                                                                        >  1:265560/1‑149 (MQ=255)             gCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCAGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCt                                                                                                                                       >  6:85795/1‑149 (MQ=255)               cAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCg                                                                                                                                     >  4:141670/1‑149 (MQ=255)                 aTCGACAGCACCAGCAACGCGCGGTTCGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGc                                                                                                                                   <  6:313203/149‑1 (MQ=255)                  tCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGc                                                                                                                                   >  2:161500/1‑148 (MQ=255)                       aGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATgta                                                                                                                             <  5:101621/149‑1 (MQ=255)                        gCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCACGGCCATgtag                                                                                                                            >  2:70497/1‑149 (MQ=255)                           ccAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGc                                                                                                                                   <  6:154679/139‑1 (MQ=255)                           ccAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCagatcg                                                                                                                            >  5:154679/1‑141 (MQ=255)                            cAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTagtagt                                                                                                                        >  5:169683/1‑149 (MQ=255)                                     gcgGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCgaagaa                                                                                                                 >  8:273644/1‑147 (MQ=255)                                       ggTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAgg                                                                                                             <  3:276037/149‑1 (MQ=255)                                         ttGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGgtc                                                                                                           <  6:169683/149‑1 (MQ=255)                                             gcagcagcaGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCt                                                                                                       >  4:87980/1‑149 (MQ=255)                                               agcagcagGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCttt                                                                                                     >  2:186403/1‑149 (MQ=255)                                                       gCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCgg                                                                                             >  1:5550/1‑149 (MQ=255)                                                        cTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCggg                                                                                            >  4:61731/1‑149 (MQ=255)                                                           atgaAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGc                                                                                           <  7:230929/147‑1 (MQ=255)                                                            tgaAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTg                                                                                        <  7:110473/149‑1 (MQ=255)                                                                       aaaCGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGt                                                                             <  7:273644/149‑1 (MQ=255)                                                                             cAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAg                                                                       <  3:87980/149‑1 (MQ=255)                                                                             cAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCACGCCCGAGGTCaaaga                                                                                >  4:129493/1‑138 (MQ=255)                                                                             cAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCACGCCCGAGGTCaaaga                                                                                <  3:129493/140‑3 (MQ=255)                                                                                        cATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGt                                                            >  8:260552/1‑149 (MQ=255)                                                                                         aTGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTc                                                           <  2:5550/149‑1 (MQ=255)                                                                                         aTGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCCCAGGCGCCCCTTGTc                                                           <  1:70497/149‑1 (MQ=255)                                                                                          tGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTcc                                                          <  2:265560/149‑1 (MQ=255)                                                                                          tGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTcc                                                          <  1:186403/149‑1 (MQ=255)                                                                                          tGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTc                                                           >  6:123812/1‑148 (MQ=255)                                                                                            aCCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCa                                                         <  3:61731/148‑1 (MQ=255)                                                                                             ccGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAg                                                        >  5:19903/1‑148 (MQ=255)                                                                                                    cGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCa                                                         <  3:219033/140‑1 (MQ=255)                                                                                                    cGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCa                                                         >  4:219033/1‑140 (MQ=255)                                                                                                           aCTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGt                                         >  7:137280/1‑149 (MQ=255)                                                                                                             ttCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTgg                                       >  3:32055/1‑149 (MQ=255)                                                                                                                  tGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGa                                    >  4:278121/1‑147 (MQ=255)                                                                                                                    gtcgtcATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACa                                  <  7:260552/147‑1 (MQ=255)                                                                                                                     tcgtcATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAgg                                >  5:68537/1‑148 (MQ=255)                                                                                                                       gtcATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAGgca                              <  3:141670/148‑1 (MQ=255)                                                                                                                          aTGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCACCCGCTGCACCGTGGCGACAGgcagca                           <  1:161500/148‑1 (MQ=255)                                                                                                                          aTGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAGGCAGCaa                          <  5:85795/149‑1 (MQ=255)                                                                                                                          aTGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAGGCAGCaa                          <  6:68537/149‑1 (MQ=255)                                                                                                                                   gcgtcggcCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAGGCAGCAATGGCAGgcg                 >  2:209500/1‑149 (MQ=255)                                                                                                                                   gcgtcggcCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAGGCAGCAATGGCAGgcg                 >  7:64869/1‑149 (MQ=255)                                                                                                                                    cgtcggcCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAGGCAGCAATGGCAGgcgc                <  1:209500/149‑1 (MQ=255)                                                                                                                                    ccgaggcCCGGGCGCTGATCGCCCAGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAGGCAGCAATGGCAGgcgc                >  3:265825/5‑149 (MQ=255)                                                                                                                                        ggcCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAGGCAGCAATGGCAGGCGCTgg             <  8:64869/148‑1 (MQ=255)                                                                                                                                         gcCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAGGCAGCAATGGCAGGCGCTGGGc           >  7:33004/1‑149 (MQ=255)                                                                                                                                         gcCCGGGCGCTGATCGCCCAGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAGGCAGCAATGGCAGGCGCTGGGc           >  4:284687/1‑149 (MQ=255)                                                                                                                                           cTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGa                                    <  4:43006/122‑1 (MQ=255)                                                                                                                                           cTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGa                                    >  3:43006/1‑122 (MQ=255)                                                                                                                                           cTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAGGCAGCAATGGCAGGCGCTGGGCta         >  1:15483/1‑149 (MQ=255)                                                                                                                                              ggCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAGGCAGCAATGGCAGGCGCTGGGCtacg       <  5:123812/148‑1 (MQ=255)                                                                                                                                                 gATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGcccc                                                                <  1:291436/88‑1 (MQ=255)                                                                                                                                                 gATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGcacc                                                                >  2:291436/1‑88 (MQ=255)                                                                                                                                                  aTCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAGGCAGCATTGGCAGGCGCTGGGCTAAGTACGc  >  8:115048/1‑149 (MQ=255)                                                                                                                                                  aTCAGCGACCGGGACTCGGCCATGTAGTAGTAGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAGGCAGCAATGGCAGGCGCTGGGCtacgtacg   >  4:292950/1‑148 (MQ=255)                                                                                                                                                    |                                                                                                                                                   GATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAGGCAGCAATGGCAGGCGCTGGGCTACGTACGC  >  NC_002947/421920‑422213

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