BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,871,650	G→C	5.9%	L259V (CTG→GTG)	PP_3420 ←	sensor histidine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,871,650	0	G	C	5.9%	320.8 / 8.6	102	L259V (CTG→GTG)	PP_3420	sensor histidine kinase
Reads supporting (aligned to +/- strand): ref base G (45/51); new base C (3/3); total (48/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.38e-01

Only 100 of 102 total aligned reads displayed.

TCTGCGCAGCAGCGTTGGCGGTTTCCTGGGCAAGTTCGGCGCTGGAGCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCCGCGCTGGATGCGCCGGATGAGGTCTTGCGGCTCCTTGTAGATTTCCAC  >  NC_002947/3871503‑3871794

tctGCGCAGCAGCGTTGGCGGTTTCCTGGGCAAGTTCGGCGCTGGAGCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGt                                                                                                                                                <  1:267753/150‑1 (MQ=255)

    cgcAGCAGCGTTGGCGGTTTCCTGGGCAAGTTCGGCGCTGGAGCGTATGTTGGCCAGCGTGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTgcag                                                                                                                                            <  2:2330144/150‑1 (MQ=255)

     ccagcagCGTTGGCGGTTTCCTGGGCAAGTTCGGCGCTGGAGCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTgcagc                                                                                                                                           <  2:252444/149‑1 (MQ=255)

      cagcagCGTTGGCGGTTTCCTGGGCAAGTTCGGCGCTGGAGCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCACCTgcagca                                                                                                                                          <  1:1680603/150‑1 (MQ=255)

       agcagcGTTGGCGGTTTCCTTGGCAAGTTCGGCGCTGGAGCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCACCTgcagcag                                                                                                                                         <  2:880221/150‑2 (MQ=255)

         cagcGTTGGCGGTTTCCTGGGCAAGTTCGGCGCTGGAGCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCatat                                                                                                                                       <  1:564229/150‑1 (MQ=255)

           gcGTTGGCGGTTTCCTGGGCAAGTTCGGCGCTGGAGCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATAtgt                                                                                                                                     <  2:2123833/150‑1 (MQ=255)

             gTTGGCGGTTTCCTGGGCAAGTTCGGCGCTGGAGCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTgc                                                                                                                                   <  1:345056/150‑1 (MQ=255)

             gTTGGCGGTTTCCTGGGCAAGTTCGGCGCTGGAGCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTgc                                                                                                                                   <  2:802419/150‑1 (MQ=255)

                 gCGGTTTCCTGGGCAAGTTCGGCGCTGGAGCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTgccgcc                                                                                                                               <  2:2293338/150‑1 (MQ=255)

                  cGGTTTCCTGGGCAAGTTCGGCGCTGGAGCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGccc                                                                                                                              <  1:1257631/150‑1 (MQ=255)

                    gTTTCCTGGGCAAGTTCGGCGCTGGAGCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCAGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCgg                                                                                                                            >  1:1788083/1‑150 (MQ=255)

                     tttCCTGGGCAAGTTCGGCGCTGGAGCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGgc                                                                                                                           <  2:2170916/150‑1 (MQ=255)

                       tCCTGGGCAAGTTCGGCGCTGGAGCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTCGAATGCTCACAGCACGATCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGgcgc                                                                                                                         >  2:1880663/1‑150 (MQ=255)

                         cTGGGCAAGTTCGGCGCTGGAGCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTgca                                                                                                                                             <  2:1061289/128‑1 (MQ=255)

                         cTGGGCAAGTTCGGCGCTGGAGCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTgca                                                                                                                                             >  1:1061289/1‑128 (MQ=255)

                         cTGGGCAAGTTCGGCGCTGGAGCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCAc                                                                                                                       >  2:750485/1‑150 (MQ=255)

                         cTGGGCAAGTTCGGCGCTGGAGCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCAc                                                                                                                       >  2:105987/1‑150 (MQ=255)

                         cTGGGCAAGTTCGGCGCTGGAGCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCAc                                                                                                                       >  1:118420/1‑150 (MQ=255)

                          tGGGCAAGTTCGGCGCTGGAGCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACg                                                                                                                      >  1:536953/1‑150 (MQ=255)

                           gggCAAGTTCGGCGCTGGAGCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGa                                                                                                                     >  2:1444125/1‑150 (MQ=255)

                            ggCAAGTTCGGCGCTGGAGCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGAt                                                                                                                    >  2:2174644/1‑150 (MQ=255)

                            ggCAAGTTCGGCGCTGGAGCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGCTTCGCTTGCCACCAGCATGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGAt                                                                                                                    <  2:2203051/150‑1 (MQ=255)

                                 gTTCGGCGCTGGAGCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCAATGTAACGACCTTCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGa                                                                                                               >  2:664121/1‑150 (MQ=255)

                                  ttCGGCGCTGGAGCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGaa                                                                                                              <  2:1602593/150‑1 (MQ=255)

                                       cgcTGGAGCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGc                                                                                                         <  2:1135731/150‑1 (MQ=255)

                                        gcTGGAGCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGcc                                                                                                        <  1:2166616/150‑1 (MQ=255)

                                         cTGGAGCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCg                                                                                                       >  2:374390/1‑150 (MQ=255)

                                            gAGCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAg                                                                                                    <  1:1936261/150‑1 (MQ=255)

                                             aGCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGt                                                                                                   <  2:536953/150‑1 (MQ=255)

                                                 tATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGAt                                                                                               <  1:2303513/150‑1 (MQ=255)

                                                            gCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCcccac                                                                                    >  2:1943117/1‑150 (MQ=255)

                                                            gCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCcccac                                                                                    <  2:2156514/150‑1 (MQ=255)

                                                              gggTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCcccaccc                                                                                  >  1:794222/1‑150 (MQ=255)

                                                                  tACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCCGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACccacca                                                                              <  1:1991222/150‑1 (MQ=255)

                                                                     gCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAgt                                                                            <  2:506610/149‑1 (MQ=255)

                                                                     gCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAgt                                                                            >  1:506610/1‑149 (MQ=255)

                                                                       aTGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTc                                                                         >  2:1654076/1‑150 (MQ=255)

                                                                       aTGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTc                                                                         <  1:977622/150‑1 (MQ=255)

                                                                         gCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGa                                                                       <  2:794222/150‑1 (MQ=255)

                                                                         gCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGa                                                                       <  1:1654076/150‑1 (MQ=255)

                                                                          cTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAg                                                                      <  1:2346129/150‑1 (MQ=255)

                                                                           tgtgAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCatat                                                                                                                                       <  1:1214891/84‑1 (MQ=255)

                                                                           tgtgAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCatat                                                                                                                                       >  2:1214891/1‑84 (MQ=255)

                                                                           tgtgAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAgg                                                                     >  2:1193544/1‑150 (MQ=255)

                                                                           tgtgAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAgg                                                                     <  1:374390/150‑1 (MQ=255)

                                                                           tgtgAGCCACTGCCGACGACATTTCACCCAGCACAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAgg                                                                     <  1:12102/150‑1 (MQ=255)

                                                                               ggCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCa                                                                 <  2:1565930/149‑1 (MQ=255)

                                                                                gCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAg                                                                >  2:326245/1‑150 (MQ=255)

                                                                                 ccACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGa                                                               >  1:1014816/1‑150 (MQ=255)

                                                                                     tgcagacgacATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAg                                                           <  2:1788083/146‑1 (MQ=255)

                                                                                      gCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCTTGACCCACCAGTGTCGAGGCCCAGATCAGc                                                          <  1:664121/150‑1 (MQ=255)

                                                                                       ccGACGACATTTCACCCAGCGCGACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGca                                                         <  2:1114774/150‑1 (MQ=255)

                                                                                         gacgacATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGcaca                                                       >  2:2329514/1‑150 (MQ=255)

                                                                                         gacgacATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGcaca                                                       >  2:93205/1‑150 (MQ=255)

                                                                                         gacgacATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGcaca                                                       >  1:1655457/1‑150 (MQ=255)

                                                                                           cgacATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATa                                                     <  1:326245/150‑1 (MQ=255)

                                                                                            gacATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGCACCCGAGGAATTAGTTCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAg                                                    >  1:1743602/1‑150 (MQ=255)

                                                                                            gacATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAg                                                    >  1:34293/1‑150 (MQ=255)

                                                                                                tttCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCcgc                                                >  2:433285/1‑150 (MQ=255)

                                                                                                    aCCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGAt                                                                                               <  2:2077893/99‑1 (MQ=255)

                                                                                                    aCCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGAt                                                                                               >  1:2077893/1‑99 (MQ=255)

                                                                                                     cccAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCCGCGCTgg                                           >  2:717257/1‑150 (MQ=255)

                                                                                                      ccAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCCGATCAGCACATAGCCGCGCTGGa                                          >  1:493297/1‑150 (MQ=255)

                                                                                                      ccAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCa                                                            <  1:25970/132‑1 (MQ=255)

                                                                                                      ccAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCa                                                            >  2:25970/1‑132 (MQ=255)

                                                                                                       cAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCCGCGCTGGAt                                         <  2:1657228/150‑1 (MQ=255)

                                                                                                           gcAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACa                                                                                                            <  1:854914/79‑1 (MQ=255)

                                                                                                           gcAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACa                                                                                                            >  2:854914/1‑79 (MQ=255)

                                                                                                            cAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGTCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCCGCGCTGGATGCGcc                                    <  2:378839/150‑1 (MQ=255)

                                                                                                              aCGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCCGCGCTGGATGCGCCgg                                  >  1:978998/1‑150 (MQ=255)

                                                                                                                gTAGGTTTCGCTTGGCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGGGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCTCCCACCAGTGTCGAGGCCCAGATCAGCAGATAGCCGCGCTGGATGCGCCGGAt                                <  2:1625516/150‑1 (MQ=255)

                                                                                                                  aGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCCCGATCCAGGACAGGCCGGAGAAGCTCGCGGCGGGACGCGCCAGTGTCGAGGCCCAGATCAGCACATAGCCGCGCTGGGTGCGCCGGATGa                              >  1:168163/1‑150 (MQ=255)

                                                                                                                  aGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCCGCGCTGGATGCGCCGGATGa                              <  1:105987/150‑1 (MQ=255)

                                                                                                                  aGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCCGCGCTGGATGCGCCGGATGa                              <  2:1014816/150‑1 (MQ=255)

                                                                                                                   ggTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCCGCGCTGGATGCGCCGGATGAg                             <  1:750485/150‑1 (MQ=255)

                                                                                                                    gTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCCGCGCTGGATGCGCCGGATGAgg                            >  2:1312929/1‑150 (MQ=255)

                                                                                                                     tttCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGAACATAGACGCGCTGGAGGCGCCGGAAGAggg                           >  1:1698542/1‑149 (MQ=255)

                                                                                                                       tCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATCGACGAGCTGGTTGAGCCGGTTGAGATCt                         >  2:1095921/1‑150 (MQ=255)

                                                                                                                         gCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATGAGCACATAGCCGCGCTGGATGCGCCGGATGAGGTCTTg                       <  1:433285/150‑1 (MQ=255)

                                                                                                                         gCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCcgc                                                >  2:2452450/1‑125 (MQ=255)

                                                                                                                         gCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCcgc                                                <  1:2452450/125‑1 (MQ=255)

                                                                                                                         gCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCCGCGCTGGATGCGCCGGATGAGGTCTTg                       <  1:1943117/150‑1 (MQ=255)

                                                                                                                           ttGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAgtgt                                                                          >  1:1820929/1‑97 (MQ=255)

                                                                                                                           ttGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAgtgt                                                                          <  2:1820929/97‑1 (MQ=255)

                                                                                                                           ttGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCCGCGCTGGATGCGCCGGATGAGGTCTTGCg                     >  1:232142/1‑150 (MQ=255)

                                                                                                                              ccaccaGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCCGCGCTGGATGCGCCGGATGAGGTCTTGCGGCt                  >  1:410838/1‑150 (MQ=255)

                                                                                                                                accaGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCCGCGCTGGATGCGCCGGATGAGGTCTTGCGGCTcc                >  2:283053/1‑150 (MQ=255)

                                                                                                                                 ccaGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCCGCGCTGGATGCGCCGGATGAGGTCTTGCGGCTCCt               >  1:1221877/1‑150 (MQ=255)

                                                                                                                                     ccTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGATCAGGCCGAAGTAGATGAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCCGCGCTGGATGCGCCGGATGAGGTCTTGCGGCTCCTTGTa           <  2:1253463/150‑1 (MQ=255)

                                                                                                                                       tGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCCGCGCTGGATGCGCCGGATGAGGTCTTGCGGCTCCTTGTAGa         <  1:2174644/150‑1 (MQ=255)

                                                                                                                                       acatatgctgcACCTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGc                                                          >  1:2139190/6‑101 (MQ=255)

                                                                                                                                       acatatgctgcACCTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGc                                                          <  2:2139190/96‑1 (MQ=255)

                                                                                                                                        catatgctgcACCTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCCGCGCTGGATGCGCCGGATGAGGTCTTGCGGCTCCTTGTAGAt        >  1:1314783/5‑150 (MQ=255)

                                                                                                                                           ctgctgCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGACCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCCGCGCTGGATGCGCCGGATGAGGTCTTGCGGCTCCTTGTAGATTTc     >  1:1054978/1‑150 (MQ=255)

                                                                                                                                            tgctgcAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAGGTAGGTCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCCGCGCTGGCTGCGCCGGATGAGGTCTTGCGGGCCCTTCTAGATTTcc    >  2:603320/1‑150 (MQ=255)

                                                                                                                                            tgctgcAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGca                                                         <  1:96038/97‑1 (MQ=255)

                                                                                                                                            tgctgcAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGca                                                         >  2:96038/1‑97 (MQ=255)

                                                                                                                                            tgctgcACCTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCCGCGCTGGATGCGCCGGATGAGGTCTTGCGGCTCCTTGTAGATTTcc    >  2:1990060/1‑150 (MQ=255)

                                                                                                                                              ctgcAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCCGCGCTGGATGCGCCGGATGAGGTCTTGCGGCTCCTTGTAGATTTccac  <  1:1193544/150‑1 (MQ=255)

TCTGCGCAGCAGCGTTGGCGGTTTCCTGGGCAAGTTCGGCGCTGGAGCGTATGTTGGCCAGCGGGTTACGCATGCTGTGAGCCACTGCCGACGACATTTCACCCAGCGCAACGTAGGTTTCGCTTGCCACCAGCCTGTCCTGCTGCAGGTGCAGCATATGTGCCGCCCGGCGCACGATCCAGAACAGGCCGAAGTAGATCAACGCCCCACCCACCAGTGTCGAGGCCCAGATCAGCACATAGCCGCGCTGGATGCGCCGGATGAGGTCTTGCGGCTCCTTGTAGATTTCCAC  >  NC_002947/3871503‑3871794

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