Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	422,069	G→T	6.7%	L627M (CTG→ATG) ‡	PP_0348 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	422,069	0	G	T	6.7%	190.2 / 3.7	60	L627M (CTG→ATG) ‡	PP_0348	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (29/27);  new base T (2/2);  total (31/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.44e-01

TTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAGGCAGCAATGGCAGGCGCTGGGCTACGTACGC  >  NC_002947/421922‑422213                                                                                                                                                    |                                                                                                                                                 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                                                          <  1:186403/149‑1 (MQ=255)                                                                                        tGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTcc                                                          <  2:265560/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                                                         >  4:219033/1‑140 (MQ=255)                                                                                                  cGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCa                                                         <  3:219033/140‑1 (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                 >  7:64869/1‑149 (MQ=255)                                                                                                                                 gcgtcggcCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAGGCAGCAATGGCAGgcg                 >  2:209500/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)                                                                                                                                                    |                                                                                                                                                 TTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCATGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAGGCAGCAATGGCAGGCGCTGGGCTACGTACGC  >  NC_002947/421922‑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: ‑