BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	422,065	A→C	5.9%	I628S (ATC→AGC)	PP_0348 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	422,065	0	A	C	5.9%	321.9 / 7.8	102	I628S (ATC→AGC)	PP_0348	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (49/47); new base C (2/4); total (51/51)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01

Only 100 of 102 total aligned reads displayed.

GATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCA‑TGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAA‑CCGCTGCACCGTGGCGACAGGCAGCAATGGCAGGCGCTGGGCTACGT  >  NC_002947/421920‑422209

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

                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)

                                            gctgctgcagGCTGATGAAATCCAGCAAACGCAAGCTCAGCGCCT‑TGACCTACTCCGAATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGTCCATGTAGTAGTCGAAGAAAAGGTCGTCt                                                                                                    <  2:313511/143‑1 (MQ=255)

                                            gcagcagcaGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCA‑TGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGgtcgt                                                                                                      <  4:325125/147‑1 (MQ=255)

                                            gcagcagcaGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCA‑TGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGgtcgt                                                                                                      >  3:325125/1‑147 (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                                                                                     <  7:253667/149‑1 (MQ=255)

                                                           tgaAATCCAGCAAACGCCAGCTCAGCGCCA‑TGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTg                                                                                     <  8:123842/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                                                                               <  3:210994/148‑1 (MQ=255)

                                                                  cAGCAAACGCCAGCTCAGCGCCA‑TGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCaaa                                                                               <  5:143102/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:13299/149‑1 (MQ=255)

                                                                                     gcCA‑TGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCtt                                                           >  5:310708/1‑149 (MQ=255)

                                                                                             ccGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGc                                                    <  3:37375/149‑1 (MQ=255)

                                                                                                aCTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGgtc                                                                                                        >  5:110430/1‑94 (MQ=255)

                                                                                                aCTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGgtc                                                                                                        <  6:110430/94‑1 (MQ=255)

                                                                                                aCTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCGGGCGCCCCTTGTCCAGCAA‑c                                                <  2:333964/149‑1 (MQ=255)

                                                                                                 cTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAA‑cc                                               >  3:252392/1‑149 (MQ=255)

                                                                                                  tCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAA‑CCg                                              >  7:138628/1‑149 (MQ=255)

                                                                                                     gCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAA‑CCGCTg                                           <  4:358652/149‑1 (MQ=255)

                                                                                                        ttCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAA‑c                                                >  1:332366/1‑141 (MQ=255)

                                                                                                        ttCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAA‑c                                                <  2:332366/141‑1 (MQ=255)

                                                                                                                 cTGGTCGTCATGCTGTCGGCGTCGGCCTTGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAA‑CCGCTGCACCGTGGCGAc                               <  2:166533/149‑1 (MQ=255)

                                                                                                                  agtcggtcATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAA‑CCGCTGCACCGTGGCGACa                              <  5:155043/144‑1 (MQ=255)

                                                                                                                   ggTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTATTCGAAGAAAAGGTCGTCTTTTACCTCGGGCGTGGTAAAATGAGTGCTCAGGCGCCCCATGTACAGCAA‑CCGCGGCACCGTGGCGACAg                             >  4:186292/1‑149 (MQ=255)

                                                                                                                   ggTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAA‑CCGCTGCACCGTGGCGACAg                             >  4:279367/1‑149 (MQ=255)

                                                                                                                    gtcgtcATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAA‑CCGCTGCACCGTGGCGACAg                             <  5:131239/148‑1 (MQ=255)

                                                                                                                       gtcATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAA‑CCGCTGc                                          <  1:288687/132‑1 (MQ=255)

                                                                                                                       gtcATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAA‑CCGCTGc                                          >  2:288687/1‑132 (MQ=255)

                                                                                                                          aTGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGTTCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAA‑CCGCTGCACCGTGGCGACAGGCAGCaa                      <  1:300581/149‑1 (MQ=255)

                                                                                                                          aTGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAA‑CCGCTGCACCGTGGCGACAGGCAGCaa                      >  5:232291/1‑149 (MQ=255)

                                                                                                                          aTGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAA‑CCGCTGCACCGTGGCGACAGGCAGCaa                      >  2:282550/1‑149 (MQ=255)

                                                                                                                           tGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAA‑CCGCTGCACCGTGGCGACAGGCAGCAAt                     <  1:155235/149‑1 (MQ=255)

                                                                                                                             cTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAA‑CCGCTGCACCGTGGCGACAGGCAGCAATg                    <  1:282550/148‑1 (MQ=255)

                                                                                                                                tcggcgtcggcCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAA‑CCGCTGCACCGTGGCGACAGGCAGCAATGGCAg                >  5:220324/1‑149 (MQ=255)

                                                                                                                                     cgaggcCCGGGCGCTGATCGCCCAGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCGAATGCGTGCTCAGGCGCCCCTTGTCCAGCAA‑CCGCTGCACCGTGGCGACAGGCAGCAATGGCAGGCGCt           >  4:117409/4‑149 (MQ=255)

                                                                                                                                      tcggcCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAA‑CCGCTGCACCGTGGCGACAGGCAGCAATGGCAGGCGCTg          >  6:232514/1‑149 (MQ=255)

                                                                                                                                       cggcCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAA‑CCGCTGCACCGTGGCGACAGGCAGCaa                      <  5:333249/136‑1 (MQ=255)

                                                                                                                                       cggcCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAA‑CCGCTGCACCGTGGCGACAGGCAGCaa                      >  4:268581/1‑136 (MQ=255)

                                                                                                                                       cggcCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAA‑CCGCTGCACCGTGGCGACAGGCAGCaa                      <  3:268581/136‑1 (MQ=255)

                                                                                                                                       cggcCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAA‑CCGCTGCACCGTGGCGACAGGCAGCaa                      >  6:333249/1‑136 (MQ=255)

                                                                                                                                        ggcCCGGGCGCTGATCGCCCAGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCGAATGCGTGCTCAGGCGCCCCTTGTCCAGCAA‑CCGCTGCACCGTGGCGACAGGCAGCAATGGCAGGCGCTggg        <  3:117409/149‑1 (MQ=255)

                                                                                                                                          ccTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAA‑CCGCTGCACCGTGGCGACAGGCAGCAATGGCAGGCGCTGGGc       >  2:350550/1‑148 (MQ=255)

                                                                                                                                           cTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAA‑CCGCTGCACCGTGGCGACAGGCAGCAATGGCAGGCGCTGGGCta     >  4:311568/1‑149 (MQ=255)

                                                                                                                                           cTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAA‑CCGCTGCACCGTGGCGACAGGCAGCAATGGCAGGCGCTGGGCta     >  7:7590/1‑149 (MQ=255)

                                                                                                                                              ggCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAA‑CCGCTGCACCGTGGCGACAGGCAGCAATGGCAGGCGCTGGGCtacgt  <  6:135012/149‑1 (MQ=255)

                                                                                                                                               gCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGcccc                                                             <  8:302131/90‑1 (MQ=255)

                                                                                                                                               gCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGcccc                                                             >  7:302131/1‑90 (MQ=255)

                                                                                                                                               gCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAACCCGCTGCACCGTGGCGACAGGCAGCAATGGCAGGCGCTGGGCtac    <  1:69186/147‑1 (MQ=255)

                                                                                                                                               gCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAA‑CCGCTGCACCGTGGCGACAGGCAGCAATGGCAGGCGCTGGGCtac    >  2:69186/1‑146 (MQ=255)

GATTGCCATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCA‑TGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAA‑CCGCTGCACCGTGGCGACAGGCAGCAATGGCAGGCGCTGGGCTACGT  >  NC_002947/421920‑422209

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