Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	422,059	T→C	6.0%	Q630R (CAG→CGG)	PP_0348 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	422,059	0	T	C	6.0%	299.6 / 5.7	100	Q630R (CAG→CGG)	PP_0348	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (49/45);  new base C (2/4);  total (51/49)
Fisher's exact test for biased strand distribution p-value = 4.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.99e-01

ATCGATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCA‑TGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAGGCAGCAATGGCAGGCGCTGGGCTA  >  NC_002947/421917‑422206                                                                                                                                                |                                                                                                                                                    atCGATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCA‑TGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCg                                                                                                                                                >  6:143102/1‑148 (MQ=255)    gATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCA‑TGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGAt                                                                                                                                              >  6:108953/1‑147 (MQ=255)    gATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCA‑TGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGAt                                                                                                                                              >  2:211558/1‑147 (MQ=255)          cATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCA‑TGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGccc                                                                                                                                      >  6:131239/1‑149 (MQ=255)           aTGCGGCCAATCGACAGCACCAGCAACGCGCGGATGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCA‑TGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCg                                                                                                                                     >  1:231028/1‑149 (MQ=255)              cGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCA‑TGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCCGGGCGCTGATCGCCCagg                                                                                                                                   <  1:62663/148‑1 (MQ=255)              cGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCA‑TGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCCGGGCGCTGATCGCCCAGGc                                                                                                                                  <  8:56648/149‑1 (MQ=255)               ggCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCA‑TGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGcc                                                                                                                                 >  8:253667/1‑149 (MQ=255)                gCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCA‑TGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCgctga                                                                                                                                            >  3:255397/1‑133 (MQ=255)                gCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCA‑TGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCgctga                                                                                                                                            <  4:255397/137‑5 (MQ=255)                gCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCA‑TGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCt                                                                                                                                >  2:297057/1‑149 (MQ=255)                  cAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCA‑TGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCg                                                                                                                              <  7:267759/149‑1 (MQ=255)                  cAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCA‑TGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCCGGGCGCTGATCGCCCAGGc                                                                                                                                  <  3:258352/145‑1 (MQ=255)                  cAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCA‑TGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCCGGGCGCTGATCGCCCAGGc                                                                                                                                  >  4:258352/1‑145 (MQ=255)                   aaTCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCA‑TGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCgg                                                                                                                             >  7:123842/1‑149 (MQ=255)                   aaTCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCA‑TGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCgg                                                                                                                             >  8:142634/1‑149 (MQ=255)                   aaTCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCA‑TGACCGACTCCGCATTCACTTCCATGGTAGTCATGCTGTCGGCGTCCGCCTGGGCGATCAGAGCCCGGGCCGCgg                                                                                                                             >  4:331857/1‑149 (MQ=255)                    aTCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCA‑TGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGc                                                                                                                            >  3:317502/1‑149 (MQ=255)                    aTCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCA‑TGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGc                                                                                                                            <  8:323940/149‑1 (MQ=255)                       gACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCA‑TGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCAt                                                                                                                         >  2:155235/1‑149 (MQ=255)                          aGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCA‑TGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATgta                                                                                                                      <  4:317502/149‑1 (MQ=255)                          aGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCA‑TGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATg                                                                                                                        >  7:74866/1‑147 (MQ=255)                           gCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCA‑TGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATgtag                                                                                                                     <  5:108953/149‑1 (MQ=255)                            cACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCA‑TGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATgtag                                                                                                                     >  1:166533/1‑148 (MQ=255)                               cAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCA‑TGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTagtagt                                                                                                                 <  5:133863/149‑1 (MQ=255)                                 gCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCA‑TGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTc                                                                                                                >  8:31684/1‑148 (MQ=255)                                 gCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCA‑TGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTc                                                                                                                >  6:350708/1‑148 (MQ=255)                                      gcgcgGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCA‑TGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCgaaga                                                                                                           >  1:333964/1‑148 (MQ=255)                                       cgcgGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCA‑TGCCCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCgaagaa                                                                                                          >  4:37375/1‑148 (MQ=255)                                       cgcgGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCA‑TGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCgaagaa                                                                                                          >  6:295817/1‑148 (MQ=255)                                           gTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCA‑TGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGgt                                                                                                     <  1:297057/149‑1 (MQ=255)                                            ttGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCA‑TGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGgtc                                                                                                    <  2:231028/149‑1 (MQ=255)                                              gggCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCA‑TGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGgtcg                                                                                                   >  8:265112/1‑148 (MQ=255)                                               ggCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCA‑TGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGgtcgtc                                                                                                 >  2:77159/1‑149 (MQ=255)                                               ggCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCA‑TGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGgtcgt                                                                                                  <  8:74866/148‑1 (MQ=255)                                                gctgctgcagGCTGATGAAATCCAGCAAACGCAAGCTCAGCGCCT‑TGACCTACTCCGAATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGTCCATGTAGTAGTCGAAGAAAAGGTCGTCt                                                                                                <  2:313511/143‑1 (MQ=255)                                                gcagcagcaGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCA‑TGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGgtcgt                                                                                                  >  3:325125/1‑147 (MQ=255)                                                gcagcagcaGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCA‑TGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGgtcgt                                                                                                  <  4:325125/147‑1 (MQ=255)                                                  agcagcagGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCA‑TGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCttt                                                                                              <  5:295817/149‑1 (MQ=255)                                                  agcagcagGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCA‑TGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCttt                                                                                              <  1:77159/149‑1 (MQ=255)                                                   gcagcagGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCA‑TGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTg                                                                                             <  7:31684/149‑1 (MQ=255)                                                    cagcagGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCA‑TGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGa                                                                                            <  2:311926/149‑1 (MQ=255)                                                          gCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCA‑TGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATgta                                                                                                                      <  8:6258/117‑1 (MQ=255)                                                          gCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCA‑TGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATgta                                                                                                                      >  7:6258/1‑117 (MQ=255)                                                          gCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCA‑TGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCgg                                                                                      <  5:350708/149‑1 (MQ=255)                                                             gatgaAATCCAGCAAACGCCAGCTCAGCGCCATTGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAgg                                                                                                      >  1:226691/1‑131 (MQ=255)                                                             gatgaAATCCAGCAAACGCCAGCTCAGCGCCA‑TGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAgg                                                                                                      <  2:226691/130‑1 (MQ=255)                                                              atgaAATCCAGCAAACGCCAGCTCAGCGCCA‑TGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGc                                                                                    <  3:331857/147‑1 (MQ=255)                                                               tgaAATCCAGCAAACGCCAGCTCAGCGCCA‑TGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTg                                                                                 <  8:123842/149‑1 (MQ=255)                                                               tgaAATCCAGCAAACGCCAGCTCAGCGCCA‑TGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTg                                                                                 <  7:253667/149‑1 (MQ=255)                                                                 aaaTCCAGCAAACGCCAGCTCAGCGCCA‑TGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTgg                                                                                <  1:211558/148‑1 (MQ=255)                                                                     ccAGCAAACGCCAGCTCAGCGCCA‑TGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATgtag                                                                                                                     <  2:97417/107‑1 (MQ=255)                                                                     ccAGCAAACGCCAGCTCAGCGCCA‑TGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCaaa                                                                           >  5:135012/1‑149 (MQ=255)                                                                      cAGCAAACGCCAGCTCAGCGCCA‑TGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCaaa                                                                           <  5:59005/148‑1 (MQ=255)                                                                      cAGCAAACGCCAGCTCAGCGCCA‑TGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCaaa                                                                           <  5:143102/148‑1 (MQ=255)                                                                      cAGCAAACGCCAGCTCAGCGCCA‑TGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCaaa                                                                           <  3:210994/148‑1 (MQ=255)                                                                      cAGCAAACGCCAGCTCAGCGCCA‑TGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCACGCCCGAGGTCaaa                                                                           <  5:105457/148‑1 (MQ=255)                                                                      cAGCAAACGCCAGCTCAGCGCCA‑TGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCACGCCCGAGGTCa                                                                             >  6:105457/1‑146 (MQ=255)                                                                        gCAAACGCCAGCTCAGCGCCA‑TGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGc                                                                        >  6:13299/1‑149 (MQ=255)                                                                            aCGCCAGCTCAGCGCCA‑TGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGc                                                                    <  7:142634/149‑1 (MQ=255)                                                                              cccAGCTCAGCGCCA‑TGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCCCCCGGCCCTCGGCCATGTACTATTAGAAGAAAAGGTCGTCATTGACCTCGTGCGTGGTCAAAtgcttgct                                                                   >  8:45813/2‑148 (MQ=255)                                                                                cAGCTCAGCGCCA‑TGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTGGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAg                                                                <  7:45813/149‑1 (MQ=255)                                                                                 aGCTCAGCGCCA‑TGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAgg                                                               >  3:358652/1‑149 (MQ=255)                                                                                    tCAGCGCCA‑TGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGgcgc                                                            >  1:132797/1‑149 (MQ=255)                                                                                      aGCGCCA‑TGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGcc                                                           <  7:265112/148‑1 (MQ=255)                                                                                      aGCGCCA‑TGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAACGCCCCGGCCTCGGGCAAGTAGTAGTCGAAGAAAAGGTCGTCTTTTACCTCCGGCGAGGTCAAATTCGTGCTAAGGCGccc                                                          >  2:300581/1‑149 (MQ=255)                                                                                         gcCA‑TGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCtt                                                       >  5:310708/1‑149 (MQ=255)                                                                                         gcCA‑TGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCtt                                                       <  5:13299/149‑1 (MQ=255)                                                                                                 ccGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGc                                                <  3:37375/149‑1 (MQ=255)                                                                                                    aCTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGgtc                                                                                                    <  6:110430/94‑1 (MQ=255)                                                                                                    aCTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGgtc                                                                                                    >  5:110430/1‑94 (MQ=255)                                                                                                    aCTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCGGGCGCCCCTTGTCCAGCAAc                                             <  2:333964/149‑1 (MQ=255)                                                                                                     cTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAAcc                                            >  3:252392/1‑149 (MQ=255)                                                                                                      tCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCg                                           >  7:138628/1‑149 (MQ=255)                                                                                                         gCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTg                                        <  4:358652/149‑1 (MQ=255)                                                                                                            ttCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAAc                                             <  2:332366/141‑1 (MQ=255)                                                                                                            ttCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAAc                                             >  1:332366/1‑141 (MQ=255)                                                                                                                     cTGGTCGTCATGCTGTCGGCGTCGGCCTTGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGAc                            <  2:166533/149‑1 (MQ=255)                                                                                                                      agtcggtcATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACa                           <  5:155043/144‑1 (MQ=255)                                                                                                                       ggTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTATTCGAAGAAAAGGTCGTCTTTTACCTCGGGCGTGGTAAAATGAGTGCTCAGGCGCCCCATGTACAGCAACCGCGGCACCGTGGCGACAg                          >  4:186292/1‑149 (MQ=255)                                                                                                                       ggTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAg                          >  4:279367/1‑149 (MQ=255)                                                                                                                        gtcgtcATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAg                          <  5:131239/148‑1 (MQ=255)                                                                                                                           gtcATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGc                                       <  1:288687/132‑1 (MQ=255)                                                                                                                           gtcATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGc                                       >  2:288687/1‑132 (MQ=255)                                                                                                                              aTGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGTTCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAGGCAGCaa                   <  1:300581/149‑1 (MQ=255)                                                                                                                              aTGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAGGCAGCaa                   >  2:282550/1‑149 (MQ=255)                                                                                                                              aTGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAGGCAGCaa                   >  5:232291/1‑149 (MQ=255)                                                                                                                               tGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAGGCAGCAAt                  <  1:155235/149‑1 (MQ=255)                                                                                                                                 cTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAGGCAGCAATg                 <  1:282550/148‑1 (MQ=255)                                                                                                                                    tcggcgtcggcCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAGGCAGCAATGGCAg             >  5:220324/1‑149 (MQ=255)                                                                                                                                         cgaggcCCGGGCGCTGATCGCCCAGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCGAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAGGCAGCAATGGCAGGCGCt        >  4:117409/4‑149 (MQ=255)                                                                                                                                          tcggcCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAGGCAGCAATGGCAGGCGCTg       >  6:232514/1‑149 (MQ=255)                                                                                                                                           cggcCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAGGCAGCaa                   <  5:333249/136‑1 (MQ=255)                                                                                                                                           cggcCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAGGCAGCaa                   <  3:268581/136‑1 (MQ=255)                                                                                                                                           cggcCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAGGCAGCaa                   >  4:268581/1‑136 (MQ=255)                                                                                                                                           cggcCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAGGCAGCaa                   >  6:333249/1‑136 (MQ=255)                                                                                                                                            ggcCCGGGCGCTGATCGCCCAGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCGAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAGGCAGCAATGGCAGGCGCTggg     <  3:117409/149‑1 (MQ=255)                                                                                                                                              ccTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAGGCAGCAATGGCAGGCGCTGGGc    >  2:350550/1‑148 (MQ=255)                                                                                                                                               cTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAGGCAGCAATGGCAGGCGCTGGGCta  >  7:7590/1‑149 (MQ=255)                                                                                                                                               cTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAGGCAGCAATGGCAGGCGCTGGGCta  >  4:311568/1‑149 (MQ=255)                                                                                                                                                |                                                                                                                                                    ATCGATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCA‑TGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCGCTGCACCGTGGCGACAGGCAGCAATGGCAGGCGCTGGGCTA  >  NC_002947/421917‑422206

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