BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	961,713	G→C	5.0%	H235Q (CAC→CAG)	selD ←	selenide, water dikinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	961,713	0	G	C	5.0%	389.8 / 7.9	120	H235Q (CAC→CAG)	selD	selenide, water dikinase
Reads supporting (aligned to +/- strand): ref base G (57/56); new base C (2/4); total (59/61)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 121 total aligned reads displayed.

GATCTTGTGCCCGTAGCTGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGC‑GGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCCTTCTTTTCGGCAG  >  NC_002947/961566‑961861

gATCTTGTGCCCGTAGCTGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGC‑GGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTg                                                                                                                                                    <  2:2461824/150‑1 (MQ=255)

 aTCTTGTGCCCGTAGCTGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGC‑GGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTgg                                                                                                                                                   <  2:302510/150‑1 (MQ=255)

  tCTTGTGCCCGTAGCTGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGC‑GGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGc                                                                                                                                                  >  1:1879038/1‑150 (MQ=255)

   cTTGTGCCCGTAGCTGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGC‑GGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGcc                                                                                                                                                 <  2:2508699/150‑1 (MQ=255)

       tgCCCGTAGCTGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGC‑GGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGc                                                                                                                                             <  2:1879038/150‑1 (MQ=255)

       tgCCCGTAGCTGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGC‑GGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGc                                                                                                                                             >  2:1019194/1‑150 (MQ=255)

         cccGTAGCTGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGC‑GGTGAGGCCGCTGCCTTCGGCCTGCTCGGCCACCTGGTCGAGctg                                                                                                                                           >  2:2745593/1‑148 (MQ=255)

          ccGTAGCTGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGC‑GGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAgg                                                                                                                                          >  2:2678614/1‑150 (MQ=255)

           cggaGCTGTCGTAGTTGCGCAGGGTTCCAACGGGGATGCAGCCTTAGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGC‑GGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGc                                                                                                                                         <  2:1413789/147‑1 (MQ=255)

           cGTAGCTGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGC‑GGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGc                                                                                                                                         <  1:2723943/150‑1 (MQ=255)

               gCTGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGC‑GGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGaa                                                                                                                                     <  1:3043238/150‑1 (MQ=255)

                cTGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGC‑GGTGAGGCCGCTGCCTTCGGCCTGCTCGGCCACCTGGTCGAGCTGGCCGAAg                                                                                                                                    <  1:1179116/150‑1 (MQ=38)

                cTGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGCGACGGCACGGCGGCGTAGTTCAGGTGCGC‑GGTGAGGCCGCTGCCTTCGGCCTGCTCGGCCACCTGGTCGAGCTGGCCGAAg                                                                                                                                    <  2:247168/150‑1 (MQ=37)

                   tCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGC‑GGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGccagatcgg                                                                                                                                 >  1:661407/1‑143 (MQ=255)

                   tCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCATGGCACGGCGGCGTAGTTCAGGTGCGC‑GGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCg                                                                                                                                 >  2:590443/1‑150 (MQ=255)

                     gTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCATGGCACGGCGGCGTAGTTCAGGTGCGC‑GGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGt                                                                                                                               <  1:590443/150‑1 (MQ=255)

                      tAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGC‑GGTGAGGCCGCTGCCTTCGGCCTGCTCGGCCACCTGGTCGAGCTGGCCGAAggcagcc                                                                                                                              <  1:2745593/150‑7 (MQ=37)

                         ttGCGCAGGGTTCCACCGGGGATGCAGCAGCCGGCGAGGGAAGGGGCGACACGATGACGGAGTGGCACGGCGGAG‑ATCGGAAGAGCACACGTCTGACCTCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGACGGCTGAAGCCGGTGACg                                                                                                                           >  1:2996253/1‑150 (MQ=255)

                             gcAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGC‑GGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCgg                                                                                                                       <  2:1714287/150‑1 (MQ=255)

                               aGGGTTCCACCGGGGATGCAGCCTTCGGCCATGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGC‑GGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTc                                                                                                                     <  1:2777452/150‑1 (MQ=255)

                               aGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGC‑GGTGAGGCCGCTGCCTTCGGCCAGCTCGGCCACCTGGTCGAGCTGGCCGAAg                                                                                                                                    >  1:624948/1‑135 (MQ=38)

                               aGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGC‑GGTGAGGCCGCTGCCTTCGGCCAGCTCGGCCACCTGGTCGAGCTGGCCGAAg                                                                                                                                    <  2:624948/135‑1 (MQ=38)

                               aGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCAACACAAGGCAGGCGGGGCACGGCGGCGTAGTTAAGATGAAA‑GCGTTCGCCGCTGCCGTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTc                                                                                                                     >  1:700127/1‑150 (MQ=255)

                                gggTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGC‑GGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGt                                                                                                                               >  1:3030616/1‑139 (MQ=255)

                                gggTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGC‑GGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGt                                                                                                                               <  2:3030616/139‑1 (MQ=255)

                                 ggTTCGACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTCCGC‑GGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCAt                                                                                                                   >  1:504473/1‑150 (MQ=255)

                                 ggTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGC‑GGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCAt                                                                                                                   >  2:76606/1‑150 (MQ=255)

                                   ttCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGC‑GGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTc                                                                                                                         <  2:2722233/142‑1 (MQ=255)

                                   ttCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGC‑GGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATgg                                                                                                                 >  1:1666628/1‑150 (MQ=255)

                                     ccACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGC‑GGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGcc                                                                                                               >  1:3047959/1‑150 (MQ=255)

                                       aCCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGC‑GGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCtt                                                                                                             >  1:591286/1‑150 (MQ=255)

                                         cGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGC‑GGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGa                                                                                                           <  1:2856290/150‑1 (MQ=255)

                                          ggggATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGC‑GGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGAc                                                                                                          >  1:2164361/1‑150 (MQ=255)

                                                  aGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGC‑GGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTc                                                                                                                         >  1:2698342/1‑127 (MQ=255)

                                                  aGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGC‑GGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTc                                                                                                                         <  2:2698342/127‑1 (MQ=255)

                                                    ccTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGC‑GGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGt                                                                                                <  1:2702199/150‑1 (MQ=255)

                                                      ttCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGC‑GGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGttt                                                                                              >  1:2526801/1‑150 (MQ=255)

                                                       tCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGC‑GGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTg                                                                                             <  1:759780/150‑1 (MQ=255)

                                                          gCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGC‑GGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCg                                                                                          >  1:2153997/1‑150 (MQ=255)

                                                           ccAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGC‑GGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTGGAGTTTGCCGa                                                                                         >  2:795490/1‑150 (MQ=255)

                                                           ccAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGC‑GGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGa                                                                                         >  2:1146658/1‑150 (MQ=255)

                                                            cAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGC‑GGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGaa                                                                                        >  1:2550421/1‑150 (MQ=255)

                                                              ggTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGC‑GGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGc                                                                                      <  2:1145124/150‑1 (MQ=255)

                                                              ggTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGC‑GGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGc                                                                                      <  1:76606/150‑1 (MQ=255)

                                                               gTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGC‑GGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCg                                                                                     >  1:2325396/1‑150 (MQ=255)

                                                               gTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGC‑GGTGAGGCAGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCg                                                                                     >  1:2051804/1‑150 (MQ=255)

                                                                  aTGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGC‑GGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCt                                                                                  >  1:1241999/1‑150 (MQ=255)

                                                                  aTGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGC‑GGTGAGGCCGCTGCCTTCGGCCAGCGCGACCAGGTGGCCGAGCAGGCCGAAGGCGGTGACGTCGGTCATGGCCTTGCCCCCGTCCAGGTAGCCGAAGCGGCt                                                                                  <  2:1981329/150‑1 (MQ=255)

                                                                   tGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGC‑GGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGAAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTg                                                                                 >  1:32097/1‑150 (MQ=255)

                                                                     gTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGC‑GGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGcc                                                                               <  2:6814/150‑1 (MQ=255)

                                                                      tCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGC‑GGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGCCCAGTTTGCCGAAGCGGCTGCCg                                                                              >  1:755892/1‑150 (MQ=255)

                                                                       cGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGC‑GGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCgg                                                                             <  2:755892/150‑1 (MQ=255)

                                                                       cGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGC‑GGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCgg                                                                             <  1:1381829/150‑1 (MQ=255)

                                                                           acTAGGCGGGCGTGGCACGGCGGCGTAGTTCAGGTGCGC‑GGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCgt                                                                         <  2:1078785/150‑1 (MQ=255)

                                                                             tAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGC‑GGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGc                                                                                            >  2:254794/1‑129 (MQ=255)

                                                                             tAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGC‑GGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGc                                                                                            <  1:254794/129‑1 (MQ=255)

                                                                               ggcaggcGTGGCACGGCGGCGTAGTTCAGGTGCGC‑GGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGa                                                                                                           <  1:410808/112‑1 (MQ=255)

                                                                                gcaggcGTGGCACGGCGGCGTAGTTCAGGTGCGC‑GGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGa                                                                    >  2:505675/1‑150 (MQ=255)

                                                                                gcaggcGTGGCACGGCGGCGTAGTTCAGGTGCGC‑GGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCACGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGa                                                                    >  2:1453944/1‑150 (MQ=255)

                                                                                  aggcGTGGCACGGCGGCGTAGTTCAGGTGCGC‑GGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAgg                                                                  <  1:945186/150‑1 (MQ=255)

                                                                                  aggcGTGGCACGGCGGCGTAGTTCAGGTGCGC‑GGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAgg                                                                  <  1:505675/150‑1 (MQ=255)

                                                                                  aggcGTGGCACGGCGGCGTAGTTCAGGTGCGC‑GGTGAGGCAGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTAAgg                                                                  <  2:2051804/150‑1 (MQ=255)

                                                                                      gTGGCACGGCGGCGTAGTTCAGGTGCGC‑GGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGAc                                                                                                          <  1:2416414/106‑1 (MQ=255)

                                                                                      gTGGCACGGCGGCGTAGTTCAGGTGCGC‑GGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGc                                                                                            <  1:1969475/120‑1 (MQ=255)

                                                                                      gTGGCACGGCGGCGTAGTTCAGGTGCGC‑GGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGc                                                                                            >  2:1969475/1‑120 (MQ=255)

                                                                                       tGGCACGGCGGCGTAGTTCAGGTGCGC‑GGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGca                                                             <  2:2325396/150‑1 (MQ=255)

                                                                                           aCGGCGGCGTAGTTCAGGTGCGC‑GGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGc                                                                           >  1:2295300/1‑132 (MQ=255)

                                                                                           aCGGCGGCGTAGTTCAGGTGCGC‑GGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGc                                                                           <  2:2295300/132‑1 (MQ=255)

                                                                                           aCGGCGGCGTAGTTCAGGTGCGC‑GGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATc                                                         <  1:2471147/150‑1 (MQ=255)

                                                                                                ggcgTAGTTCAGGTGCGC‑GGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTc                                                    >  1:1051973/1‑150 (MQ=255)

                                                                                                   gTAGTTCAGGTGCGC‑GGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTcgcg                                                 >  2:375496/1‑150 (MQ=255)

                                                                                                      gTTCAGGTGCGC‑GGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCg                                                                                                                                 >  2:2501433/1‑67 (MQ=255)

                                                                                                         cAGGTGCGC‑GGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCTGCATGGCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAg                                           <  1:1753340/150‑1 (MQ=255)

                                                                                                         cAGGTGCGC‑GGTGAGGCCGCTGCATTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAg                                           >  1:921578/1‑150 (MQ=255)

                                                                                                         cAGGTGCGC‑GGTGAGGCCACTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAg                                           <  2:2996253/150‑1 (MQ=255)

                                                                                                              gcgc‑GGTGAGGCTGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCtt                                      <  1:2054217/150‑1 (MQ=255)

                                                                                                              gcgc‑GGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCtt                                      <  2:1051973/150‑1 (MQ=255)

                                                                                                                gc‑gGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTgg                                    >  2:2210494/1‑150 (MQ=255)

                                                                                                                cgatgtcaGGCCGCTGCCTTCGGCAAGCTCCTCGTGGTGCCTCAGAAGCGCGCAGCCGGTGACGTCGGTCATGGGCTTCACCCAGTCGAGGTTGACGAGGCGGCTGACGGGCGTGTTGTGGGTGCCCAGCTAGTCGCGAGCCAGGCCTTg                                     <  2:504473/143‑1 (MQ=255)

                                                                                                                     tGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCt                                <  1:1619701/150‑1 (MQ=255)

                                                                                                                         gCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCg                                                                                                                                 <  2:2030318/49‑1 (MQ=255)

                                                                                                                         gCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCg                                                                                                                                 >  1:2030318/1‑49 (MQ=255)

                                                                                                                         gCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCa                                            >  1:1860395/1‑134 (MQ=255)

                                                                                                                         gCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCa                                            <  2:1860395/134‑1 (MQ=255)

                                                                                                                          ccGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGcaca                           <  1:488843/150‑1 (MQ=255)

                                                                                                                           cGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCgg                                                                                                                       >  1:1098659/1‑57 (MQ=255)

                                                                                                                           cGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGACCTTGGTCCTGCacac                          <  1:805361/150‑1 (MQ=255)

                                                                                                                             cTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGc                        <  2:1247940/150‑1 (MQ=255)

                                                                                                                                  ttCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGTCCGTGACCCCGTCCATTGTGGCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCtt                   >  1:1912229/1‑150 (MQ=255)

                                                                                                                                  ttCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCtt                   >  2:1180574/1‑150 (MQ=255)

                                                                                                                                   tCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCTGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTg                  <  1:1146658/150‑1 (MQ=255)

                                                                                                                                      gCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGcc               >  2:2905808/1‑150 (MQ=255)

                                                                                                                                      gCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGcc               >  2:1862060/1‑150 (MQ=255)

                                                                                                                                       ccAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCct              <  2:1666628/150‑1 (MQ=255)

                                                                                                                                         tgCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACAt                                                          <  2:805747/103‑1 (MQ=255)

                                                                                                                                          gCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCcttct           >  1:549979/1‑150 (MQ=255)

                                                                                                                                          gCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCcttct           <  1:1862060/150‑1 (MQ=255)

                                                                                                                                          gCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCcttct           >  1:1482912/1‑150 (MQ=255)

                                                                                                                                            tCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCCTTCttt         >  1:2747586/1‑150 (MQ=255)

                                                                                                                                                   ggTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCCTTCTTTTCGGCAg  >  2:1564048/1‑150 (MQ=255)

GATCTTGTGCCCGTAGCTGTCGTAGTTGCGCAGGGTTCCACCGGGGATGCAGCCTTCGGCCAGGTAATGGTCGACACTAGGCAGGCGTGGCACGGCGGCGTAGTTCAGGTGCGC‑GGTGAGGCCGCTGCCTTCGGCCAGCTCGACCAGGTGGCCGAGCAGGCCGAAGCCGGTGACGTCGGTCATGGCCTTGACCCCGTCCAGTTTGCCGAAGCGGCTGCCGGGCGTGTTGAGGGTGCACATCCAGTCGCGAGCCAGGCCTTGGTCCTGCACACGCAGCTTGGCCTTCTTTTCGGCAG  >  NC_002947/961566‑961861

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 12 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`