BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	422,069	G→T	6.4%	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.4%	321.4 / 8.1	101	L627M (CTG→ATG) ‡	PP_0348	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (47/47); new base T (2/4); total (50/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 = 5.70e-01

Only 100 of 101 total aligned reads displayed.

CATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCA‑TGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAA‑CCGCTGCACCGTGGCGACAGGCAGCAATGGCAGGCGCTGGGCTACGTACGCTC  >  NC_002947/421926‑422215

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)

         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                                                                                                                          >  6:350708/1‑148 (MQ=255)

                       gCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCA‑TGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTc                                                                                                                          >  8:31684/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                                                                                           <  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:13299/149‑1 (MQ=255)

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

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

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

                                                                                          aCTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGgtc                                                                                                              >  5:110430/1‑94 (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                                                >  2:288687/1‑132 (MQ=255)

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

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

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

                                                                                                                    aTGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAA‑CCGCTGCACCGTGGCGACAGGCAGCaa                            >  5:232291/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                            >  6:333249/1‑136 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

                                                                                                                                              cAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAA‑CCGCTGCACCGTGGCGACAGGCAGCAATGGCAGGCGCTGGGCTACGTACGctc  >  3:148363/1‑149 (MQ=255)

CATGCGGCCAATCGACAGCACCAGCAACGCGCGGTTGGGCAGCAGCAGGCTGATGAAATCCAGCAAACGCCAGCTCAGCGCCA‑TGACCGACTCCGCATTCACTTCCCTGGTCGTCATGCTGTCGGCGTCGGCCTGGGCGATCAGCGCCCGGGCCTCGGCCATGTAGTAGTCGAAGAAAAGGTCGTCTTTGACCTCGGGCGTGGTCAAATGCGTGCTCAGGCGCCCCTTGTCCAGCAA‑CCGCTGCACCGTGGCGACAGGCAGCAATGGCAGGCGCTGGGCTACGTACGCTC  >  NC_002947/421926‑422215

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