BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,206,703	G→C	5.6%	H768Q (CAC→CAG)	PP_1058 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,206,703	0	G	C	5.6%	478.9 / 13.0	143	H768Q (CAC→CAG)	PP_1058	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (78/57); new base C (4/4); total (82/61)
Fisher's exact test for biased strand distribution p-value = 7.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 146 total aligned reads displayed.

CGGATGCGGCTGTCACGGGCGCCCTGGTGGCGGCTCTTGAACAGCCAGCCGTATACCTCCTGGCGTTCGAAACGACTGGCGTTGATCATCTCGGTCAGTGTCGAGCTCGGGTTCTTCAGCAGGTAGCCAAGCAGCCACAGGTCCAGCGGGTGGACCT‑TGGCGATGTAACCCTGGTCGGGCAGGTCGTACTTGCCGGGGCCGTAAGCTACGTACATCTCGGACAGACGGCCATCGGTCAGTTTGCCCAGCACGACCTTGTCACCCTTGAGGTGGGCACGCACGAACGCGTTGAAGGTTTC  >  NC_002947/1206554‑1206852

cGGATGCGGCTGTCACGGGCGCCCTGGTGGCGGCTCTTGAACAGCCAGCCGTATACCTCCTGGCGTTCGAAACGACTGGCGTTGATCATCTCGGTCAGTGTCGAGCTCGGGTTCTTCAGCAGGTAGCCAAGCAGCCACAGGTCCAGCggg                                                                                                                                                        >  1:504833/1‑150 (MQ=255)

 ggATGCGGCTGTCACGGGCGCCCTGGTGGCGGCTCTTGAACAGCCAGCCGTATACCTCCTGGCGTTCGAAACGACTGGCGTTGATCATCTCGGTCAGTGTCGAGCTCGGGTTCTTCAGCAGGTAGCCAAGCAGCCACAGGTCCAGCGGGt                                                                                                                                                       <  2:1861316/150‑1 (MQ=255)

     gCGGCTGTCACGGGCGCCCTGGTGGCGGCTCTTGAACAGCCAGCCGTATACCTCCTGGCGTTCGAAACGACTGGCGTTGATCATCTCGGTCAGTGTCGAGCTCGGGTTCTTCAGCAGGTAGCCAAGCAGCCACAGGTCCAGCGGGTGGAc                                                                                                                                                   <  2:5122332/150‑1 (MQ=255)

      cGGCTGTCACGGGCGCCCTGGTGGCGGCTCTTGAACAGCCAGCCGTATACCTCCTGGCGTTCGAAACGACTGGCGTTGATCATCTCGGTCAGTGTCGAGCTCGGGTTCTTCAGCAGGTAGCCAAGCAGCCACAGGTCCAGCGGGTGGAcc                                                                                                                                                  >  2:3721152/1‑150 (MQ=255)

       ggCTGTCACGGGCGCCCTGGTGGCGGCTCTTGAACAGCCAGCCGTATACCTCCTGGCGTTCGAAACGACTGGCGTTGATCATCTCGGGCAGTGTCGAGCTCGGGTTCTTCAGCAGGAAGCCAAGCAGCCACAGGTCCGGCGGGTGGACCt                                                                                                                                                 >  1:2346382/1‑150 (MQ=255)

         cTGTCACGGGCGCCCTGGTGGCGGCTCTTGAACAGCCAGCCGTATACCTCCTGGCGTTCGAAACGACTGGCGTTGATCATCTCGGTCAGTGTCGAGCTCGGGTTCTTCAGCAGGTAGCCAAGCAGCCACAGGTCCACCCGCTGGACCtgt                                                                                                                                               <  1:2579391/150‑3 (MQ=255)

         cTGTCACGGGCGCCCTGGTGGCGGCTCTTGAACAGCCAGCCGTATACCTCCTGGCGTTCGAAACGACTGGCGTTGATCATCTCGGTCAGTGTCGAGCTCGGGTTCTTCAGCAGGTAGCCAAGCAGCCACAGGTCCACCCGCTGGACCtgt                                                                                                                                               <  1:1766928/150‑3 (MQ=255)

          tGTCACGGGCGCCCTGGTGGCGGCTCTTGAACAGCCAGCCGTATACCTCCTGGCGTTCGAAACGACTGGCGTTGATCATCTCGGTCAGTGTCGAGCTCGGGTTCTTCAGCAGGTAGCCAAGCAGCCACAGGTCCAGCGGGTGGACCT‑Tgg                                                                                                                                             <  1:615713/150‑1 (MQ=255)

          tGTCACGGGCGCCCTGGTGGCGGCTCTTGAACAGCCAGCCGTATACCTCCTGGCGTTCGAAACGACTGGCGTTGATCATCTCGGTCAGTGTCGAGCTCGGGTTCTTCAGCAGGTAGCCAAGCAGCCACAGGTCCAGCGGGTGGACCT‑Tgg                                                                                                                                             <  2:4112610/150‑1 (MQ=255)

              aCGGGCGCCCTGGTGGCGGCTCTTGAACAGCCAGCCGTATACCTCCTGGCGTTCGAAACGACTGGCGTTGATCATCTCGGTCAGTGTCGAGCTCGGGTTCTTCAGCAGGTAGCCAAGCAGCCACAGGTCCAGCGGGTGGACCT‑TGGCGAt                                                                                                                                         >  1:4140854/1‑150 (MQ=255)

                 ggCGCCCTGGTGGCGGCTCTTGAACAGCCAGCCGTATACCTCCTGGCGTTCGAAACGACTGGCGTTGATCATCTCGGTCAGTGTCGAGCTCGGGTTCTTCAGCAGGTAGCCAAGCAGCCACAGGTCCAGCGGGTGGACCT‑TGGCGATGTa                                                                                                                                      >  2:2930846/1‑150 (MQ=255)

                      ccTGGTGGCGGCTCTTGAACAGCCAGCCGTATACCTCCTGGCGTTCGAAACGACTGGCGTTGATCATCTCGGTCAGTGTCGAGCTCGGGTTCTTCAGCAGGTAGCCAAGCAGCCACAGGTCCAGCGGGTGGACCT‑TGGCGATGTAACCCt                                                                                                                                 >  2:5086676/1‑150 (MQ=255)

                      ccTGGTGGCGGCTCTTGAACAGCCAGCCGTATACCTCCTGGCGTTCGAAACGACTGGCGTTGATCATCTCGGTCAGTGTCGAGCTCGGGTTCTTCAGCAGGTAGCCAAGCAGCCACAGGTCCAGCGGGTGGACCT‑TGGCGATGTAACCCt                                                                                                                                 <  2:225296/150‑1 (MQ=255)

                       cTGGTGGCGGCTCTTGAACAGCCAGCCGTATACCTCCTGGCGTTCGAAACGACTGGCGTTGATCATCTCGGTCAGTGTCGAGCTCGGGTTCTTCAGCAGGTAGCCAAGCAGCCACAGGTCCAGCGGGTGGACCT‑TGGCGATGTAACCCTg                                                                                                                                <  2:2803139/150‑1 (MQ=255)

                        tggtggCGGCTCTTGAACAGCCAGCCGTATACCTCCTGGCGTTCGAAACGACTGGCGTTGATCATCTCGGTCAGTGTCGAGCTCGGGTTCTTCAGCAGGTAGCCAAGCAGCCACAGGTCCAGCGGGTGGACCT‑TGGCGATGTAACCCTgg                                                                                                                               >  2:94972/1‑150 (MQ=255)

                        tggtggCGGCTCTTGAACAGCCAGCCGTATACCTCCTGGCGTTCGAAACGACTGGCGTTGATCATCTCGGTCAGTGTCGAGCTCGGGTTCTTCAGCAGGTAGCCAAGCAGCCACAGGTCCAGCGGGTGGACCT‑TGGCGATGTAACCCTgg                                                                                                                               >  2:63791/1‑150 (MQ=255)

                          gtggCGGCTCTTGAACAGCCAGCCGTATACCTCCTGGCGTTCGAAACGACTGGCGTTGATCATCTCGGTCAGTGTCGAGCTCGGGTTCTTCAGCAGGTAGCCAAGCAGCCACAGGTCCAGCGGGTGGACCT‑TGGCGATGTAACCCTGGTc                                                                                                                             <  2:504833/150‑1 (MQ=255)

                             gcggcTCTTGAACAGCCAGCCGTATACCTCCTGGCGTTCGAAACGACTGGCGTTGATCATCTCGGTCAGTGTCGAGCTCGGGTTCTTCAGCAGGTAGCCAAGCAGCCACAGGTCCAGCGGGTGGACCT‑TGGCGATGTAACCCTGGTCggg                                                                                                                          >  1:406858/1‑150 (MQ=255)

                              cggcTCTTGAACAGCCAGCCGTATACCTCCTGGCGTTCGAAACGACTGGCGTTGATCATCTCGGTCAGTGTCGAGCTCGGGTTCTTCAGCAGGTAGCCAAGCAGCCACAGGTCCAGCGGGTGGACCT‑TGGCGATGTAACCCTGGTCGGGc                                                                                                                         <  2:1501734/150‑1 (MQ=255)

                              cggcTCTTGAACAGCCAGCCGTAGACCTCCTGGCGTTCGAAACGACTGGCGTTGATCATCTCGGTCAGTGTCGAGCTCGGGTTCTTCAGCAGGTAGCCAAGCAGCCACAGGTCCAGCGGGTGGACCT‑TGGCGATGTAACCCTGGTCGGGc                                                                                                                         <  1:1248258/150‑1 (MQ=255)

                               ggcTCTTGAACAGCCAGCCGTATACCTCCTGGCGTTCGAAACGACTGGCGTTGATCATCTCGGTCAGTGTCGAGCTCGGGTTCTTCAGCAGGTAGCCAAGCAGCCACAGGTCCAGCGGGTGGACCT‑TGGCGATGTAACCCTGGTCGGGCa                                                                                                                        >  2:5002643/1‑150 (MQ=255)

                                  tctTGAACAGCCAGCCGTATACCTCCTGGCGTTCGAAACGACTGGCGTTGATCATCTCGGTCAGTGTCGAGCTCGGGTTCTTCAGCAGGTAGCCAAGCAGCTACAGGTCCAGCGGGTGGACCT‑TGGCGATGTAACCCTGGTCGGGCAGGt                                                                                                                     >  2:459614/1‑150 (MQ=255)

                                  tctTGAACAGCCAGCCGTATACCTCCTGGCGTTCGAAACGACTGGCGTTGATCATCTCGGTCAGTGTCGAGCTCGGGTTCTTCAGCAGGTAGCCAAGCAGCCACAGGTCCAGCGGGTGGACCT‑TGGCGATGTAACCCTGGTCGGGCAGGt                                                                                                                     >  2:343336/1‑150 (MQ=255)

                                      gAACAGCCAGCCGTATACCTCCTGGCGTTCGAAACGACTGGCGTTGATCATCTCGGTCAGTGTCGAGCTCGGGTTCTTCAGCAGGTAGCCAAGCAGCCACAGGTCCAGCGGGTGGACCT‑TGGCGATGTAACCCTGGTCGGGCAGGTCGTa                                                                                                                 <  2:3735984/150‑1 (MQ=255)

                                       aaCAGCCAGCCGTATACCTCCTGGCGTTCGAAACGACTGGCGTTGATCATCTCGGTCAGTGTCGAGCTCGGGTTCTTCAGCAGGTAGCCAAGCAGCCACAGGTCCAGCGGGTGGACCT‑TGGCGATGTAACCCTGGTCGGGCAGGTCGTAc                                                                                                                >  2:3355012/1‑150 (MQ=255)

                                        aCAGCCAGCCGTATACCTCCTGGCGTTCGAAACGACTGGCGTTGATCATCTCGGTCAGTGTCGAGCTCGGGTTCTTCAGCAGGTAGCCAAGCAGCCACAGGTCCAGCGGGTGGACCT‑TGGCGATGTAACCCTGGTCGGGCAGGTCGTACt                                                                                                               >  2:423483/1‑150 (MQ=255)

                                              agccGTATACCTCCTGGCGTTCGAAACGACTGGCGTTGATCATCTCGGTCAGTGTCGAGCTCGGGTTCTTCAGCAGGTAGCCAAGCAGCCACAGGTCCAGCGGGTGGACCT‑TGGCGATGTAACCCTGGTCGGGCAGGTCGTACTTGCCgg                                                                                                         <  2:3363463/150‑1 (MQ=255)

                                              agccGTATACCTCCTGGCGTTCGAAACGACTGGCGTTGATCATCTCGGTCAGTGTCGAGCTCGGGTTCTTCAGCAGGTAGCCAAGCAGCCACAGGTCCAGCGGGTGGACCT‑TGGCGATGTAACCCTGGTCGGGCAGGTCGTACTTGCCgg                                                                                                         <  1:145240/150‑1 (MQ=255)

                                                ccGTATACCTCCTGGCGTTCGAAACGACTGGCGTTGATCATCTCGGTCAGTGTCGAGCTCGGGTTCTTCAGCAGGTAGCCAAGCAGCCACAGGTCCAGCGGGTGGACCT‑TGGCGATGTAACCCTGGTCGGGCAGGTCGTACTTGCCgggg                                                                                                       <  1:1146793/150‑1 (MQ=255)

                                                 cGTATACCTCCTGGCGTTCGAAACGACTGGCGTTGATCATCTCGGTCAGTGTCGAGCTCGGGTTCTTCAGCAGGTAGCCAAGCAGCCACAGGTCCAGCGGGTGGACCT‑TGGCGATGTAACCCTGGTCGGGCAGGTCGTACTTGCCGGGGc                                                                                                      <  2:2306898/150‑1 (MQ=255)

                                                   tataCCTCCTGGCGTTCGAAACGACTGGCGTTGATCATCTCGGTCAGTGTCGAGCTCGGGTTCTTCAGCAGGTAGCCAAGCAGCCACAGGTCCAGCGGGTGGACCT‑TGGCGATGTAACCCTGGTCGGGCAGGTCGTACTTGCCGGGGCCg                                                                                                    >  1:4446061/1‑150 (MQ=255)

                                                    ataCCTCCTGGCGTTCGAAACGACTGGCGTTGATCATCTCGGTCAGTGTCGAGCTCGGGTTCTTCAGCAGGTAGCCAAGCAGCCACAGGTCCAGCGGGTGGACCT‑TGGCGATGTAACCCTGGTCGGGCAGGTCGTACTTGCCGGGGCCGt                                                                                                   >  2:4702401/1‑150 (MQ=255)

                                                         tcctGGCGTTCGAAACGACTGGCGTTGATCATCTCGGTCAGTGTCGAGCTCGGGTTCTTCAGCAGGTAGCCAAGCAGCCACAGGTCCAGCGGGTGGACCT‑TGGCGATGTAACCCTGGTCGGGCAGGTCGTACTTGCCGGGGCCGTAAGCt                                                                                              >  2:971836/1‑150 (MQ=255)

                                                           ctGGCGTTCGAAACGACTGGCGTTGATCATCTCGGTCAGTGTCGAGCTCGGGTTCTTCAGCAGGTAGCCAAGCAGCCACAGGTCCAGCGGGTGGACCT‑TGGCGATGTAACCCTGGTCGGGCAGGTCGTACTTGCCGGGGCCGTAAGCTAc                                                                                            <  2:1936790/150‑1 (MQ=255)

                                                              gCGTTCGAAACGACTGGCGTTGATCATCTCGGTCAGTGTCGAGCTCGGGTTCTTCAGCAGGTAGCCAAGCAGCCACAGGTCCAGCGGGTGGACCT‑TGGCGATGTAACCCTGGTCGGGCAGGTCGTACTTGCCGGGGCCGTAAGCTACGTa                                                                                         <  2:1639730/150‑1 (MQ=255)

                                                               cGTTCGAAACGACTGGCGTTGATCATCTCGGTCAGTGTCGAGCTCGGGTTCTTCAGCAGGTAGCCAAGCAGCCACAGGTCCAGCGGGTGGACCT‑TGGCGATGTAACCCTGGTCGGGCAGGTCGTACTTGCCGGGGCCGTAAGCTACGTAc                                                                                        <  1:3513729/150‑1 (MQ=255)

                                                                gTTCGAAACGACTGGCGTTGATCATCTCGGTCAGTGTCGAGCTCGGGTTCTTCAGCAGGTAGCCAAGCAGCCACAGGTCCAGCGGGTGGACCT‑TGGCGATGTAACCCTGGTCGGGCAGGTCGTACTTGCCGGGGCCGTAAGCTACGTACa                                                                                       >  1:3321900/1‑150 (MQ=255)

                                                                gTTCGAAACGACTGGCGTTGATCATCTCGGTCAGTGTCGAGCTCGGGTTCTTCAGCAGGTAGCCAAGCAGCCACAGGTCCAGCGGGTGGACCT‑TGGCGATGTAACCCTGGTCGGGCAGGTCGTACTTGCCGGGGCCGTAAGCTACGTACa                                                                                       >  2:1455035/1‑150 (MQ=255)

                                                                 ttCGAAACGACTGGCGTTGATCATCTCGGTCAGTGTCGAGCTCGGGTTCTTCAGCAGGTAGCCAAGCAGCCACAGGTCCAGCGGGTGGACCT‑TGGCGATGTAACCCTGGTCGGGCAGGTCGTACTTGCCGGGGCCGTAAGCTACGTACAt                                                                                      >  2:4996841/1‑150 (MQ=255)

                                                                    gAAACGACTGGCGTTGATCATCTCGGTCAGTGTCGAGCTCGGGTTCTTCAGCAGGTAGCCAAGCAGCCACAGGTCCAGCGGGTGGACCT‑TGGCGATGTAACCCTGGTCGGGCAGGTCGTACTTGCCGGGGCCGTAAGCTACGTACAtctc                                                                                   >  1:312209/1‑150 (MQ=255)

                                                                    gAAACGACTGGCGTTGATCATCTCGGTCAGTGTCGAGCTCGGGTTCTTCAGCAGGTAGCCAAGCAGCCACAGGTCCAGCGGGTGGACCT‑TGGCGATGTAACCCTGGTCGGGCAGGTCGTACTTGCCGGGGCCGTAAGCTACGTACAtctc                                                                                   <  1:5077572/150‑1 (MQ=255)

                                                                    gAAACGACTGGCGTTGATCATCTCGGTCAGTGTCGAGCTCGGGTTCTTCAGCAGGTAGCCAAGCAGCCACAGGTCCAGCGGGTGGACCT‑TGGCGATGTAACCCTGGTCGGGCAGGTCGTACTTGCCGGGGCCGTAAGCTACGTACAtctc                                                                                   >  1:312287/1‑150 (MQ=255)

                                                                    gAAACGACTGGCGTTGATCATCTCGGTCAGTGTCGAGCTCGGGTTCTTCAGCAGGTAGCCAAGCAGCCACAGGTCCAGCGGGTGGACCT‑TGGCGATGTAACCCTGGTCGGGCAGGTCGTACTTGCCGGGGCCGTAAGCTACGTACAtctc                                                                                   >  1:691993/1‑150 (MQ=255)

                                                                    gAAACGACTGGCGTTGATCATCTCGGTCAGTGTCGAGCTCGGGTTCTTCAGCAGGTAGCCAAGCAGCCACAGGTCCAGCGGGTGGACCT‑TGGCGATGTAACCCTGGTCGGGCAGGTCGTACTTGCCGGGGCCGTAAGCTACGTACAtctc                                                                                   >  1:1293739/1‑150 (MQ=255)

                                                                      aaCGACTGGCGTTGATCATCTCGGTCAGTGTCGAGCTCGGGTTCTTCAGCAGGTAGCCAAGCAGCCACAGGTCCAGCGGGTGGACCT‑TGGCGATGTAACCCTGGTCGGGCAGGTCGTACTTGCCGGGGCCGTAAGCTACGTACATCTCgg                                                                                 >  2:4996662/1‑150 (MQ=255)

                                                                      aaCGACTGGCGTTGATCATCTCGGTCAGTGTCGAGCTCGGGTTCTTCAGCAGGTAGCCAAGCAGCCACAGGTCCAGCGGGTGGACCT‑TGGCGATGTAACCCTGGTCGGGCAGGTCGTACTTGCCGGGGCCGTAAGCTACGTACATCTCgg                                                                                 >  2:164904/1‑150 (MQ=255)

                                                                          aCTGGCGTTGATCATCTCGGTCAGTGTCGAGCTCGGGTTCTTCAGCAGGTAGACAAGCAGCCACAGGTCCAGCGGGTGGACCT‑TGGCGATGTAACCCTGGTCGGGCAGGTCGTACTTGCCGGGGCCGTAAGCTACGGACATCTCGGACAg                                                                             >  2:1242803/1‑150 (MQ=255)

                                                                              gCGTTGACCATCTCGGTCAGTGTCGAGCTCGGGTTCTTCAGCAGGTAGCCAAGCAGCCACAGGTCCAGCGGGTGGACCT‑TGGCGATGTAACCCTGGTCGGGCAGGTCGTACTTGCCGGGGCCGTAAGCTACGTACATCTCGGACAGACgg                                                                         >  2:1691618/1‑150 (MQ=255)

                                                                               cGTTGATCATCTCGGTCAGTGTCGAGCTCGGGTTCTTCAGCAGGTAGCCAAGCAGCCACAGGTCCAGCGGGTGGACCT‑TGGCGATGTAACCCTGGTCGGGCAGGTCGTACTTGCCGGGGCCGTAAGCTACGTACATCTCGGACAGACGGc                                                                        >  1:4654377/1‑150 (MQ=255)

                                                                                    atcatcTCGGTCAGTGTCGAGCTCGGGTTCTTCAGCAGGTAGCCAAGCAGCCACAGGTCCAGCGGGTGGACCT‑TGGCGATGTAACCCTGGTCGGGCAGGTCGTACTTGCCGGGGCCGTAAGCTACGTACATCTCGGACAGACGGCCATCg                                                                   <  2:3321900/150‑1 (MQ=255)

                                                                                       atcTCGGTCAGTGTCGAGCTCGGGTTCTTCAGCAGGTAGCCAAGCAGCCACAGGTCCAGCGGGTGGACCT‑TGGCGATGTAACCCTGGTCGGGCAGGTCGTACTTGCCGGGGCCGTAAGCTACGTACATCTCGGACAGACGGCCATCGGTc                                                                <  1:1639772/150‑1 (MQ=255)

                                                                                        tctcGGTCAGTGTCGAGCTCGGGTTCTTCAGCAGGTAGCCAAGCAGCCACAGGTCCAGCGGGTGGACCT‑TGGCGATGTAACCCTGGTCGGGCAGGTCGTACTTGCCGGGGCCGTAAGCTACGTACATCTCGGACAGACGGCCATCGGTCa                                                               <  2:1389227/150‑1 (MQ=255)

                                                                                        tctcGGTCAGTGTCGAGCTCGGGTTCTTCAGCAGGTAGCCAAGCAGCCACAGGTCCAGCGGGTGGACCT‑TGGCGATGTAACCCTGGTCGGGCAGGTCGTACTTGCCGGGGCCGTAAGCTACGTACATCTCGGACAGACGGCCATCGGTCa                                                               >  1:2319854/1‑150 (MQ=255)

                                                                                          tcGGTCAGTGTCGAGCTCGGGTTCTTCAGCAGGTAGCCAAGCAGCCACAGGTCCAGCGGGTGGACCT‑TGGCGATGTAACCCTGGTCGGGCAGGTCGTACTTGCCGGGGCCGTAAGCTACGTACATCTCGGACAGACGGCCATCGGTCAGt                                                             <  2:617806/150‑1 (MQ=255)

                                                                                                 gtgtCGAGCTCGGGTTCTTCAGCAGGTAGCCAAGCAGCCACAGGTCCAGCGGGTGGACCT‑TGGCGATGTAACCCTGGTCGGGCAGGTCGTACTTGCCGGGGCCGTAAGCTACGTACATCTCGGACAGACGGCCATCGGTCAGTTTGCCCa                                                      >  2:4955286/1‑150 (MQ=255)

                                                                                                       aGCTCGGGTTCTTCAGCAGGTAGCCAAGCAGCCACAGGTCCAGCGGGTGGACCT‑TGGCGATGTAACCCTGGTCGGGCAGGTCGTACTTGCCGGGGCCGTAAGCTACGTACATCTCGGACAGACGGCCATCGGTCAGTTTGCCCAGCACGa                                                >  1:4937115/1‑150 (MQ=255)

                                                                                                        gCTCGGGTTCTTCAGCAGGTAGCCAAGCAGCCACAGGTCCAGCGGGTGGACCT‑TGGCGATGTAACCCTGGTCGGGCAGGTCGTACTTGCCGGGGCCGTAAGCTATGTACATCTCGGACAGACGGCCATCGGTCAGTTTGCCCAGCACGAc                                               >  1:3671147/1‑150 (MQ=255)

                                                                                                           cGGGTTCTTCAGCAGGTAGCCAAGCAGCCACAGGTCCAGCGGGTGGACCT‑TGGCGATGTAACCCTGGTCGGGCAGGTCGTACTTGCCGGGGCCGTAAGCTACGTACATCTCGGACAGACGGCCATCGGTCAGTTTGCCCAGCACGACCtt                                            <  2:4375929/150‑1 (MQ=255)

                                                                                                            gggTTCTTCAGCAGGTAGCCAAGCAGCCACAGGTCCAGCGGGTGGACCT‑TGGCGATGTAACCCTGGTCGGGCAGGTCGTACTTGCCGGGGCCGTAAGCTACGTACATCTCGGACGGACGGCCCTCGGTCAGTGTGCCCAGCACGACCTTg                                           >  2:3566069/1‑150 (MQ=255)

                                                                                                            gggTTCTTCAGCAGGTAGCCAAGCAGCCACAGGTCCAGCGGGTGGACCT‑TGGCGATGTAACCCTGGTCGGGCAGGTCGTACTTGCCGGGGCCGTAAGCTACGTACATCTCGGACAGACGGCCATCGGTCAGTTTGCCCAGCACGACCTTg                                           >  2:3566270/1‑150 (MQ=255)

                                                                                                            gggTTCTTCAGCAGGTAGCCAAGCAGCCACAGGTCCAGCGGGTGGACCT‑TGGCGATGTAACCCTGGTCGGGCAGGTCGTACTTGCCGGGGCCGTAAGCTACGTACATCTCGGACAGACGGCCATCGGTCAGTTTGCCCAGCACGACCTTg                                           >  1:473736/1‑150 (MQ=255)

                                                                                                             ggTTCTTCAGCAGGTAGCCAAGCAGCCACAGGTCCAGCGGGTGGACCT‑TGGCGATGTAACCCTGGTCGGGCAGGTCGTACTTGCCGGGGCCGTAAGCTACGTACATCTCGGACAGACGGCCATCGGTCAGTTTGCCCAGCACGACCTTGt                                          <  1:3963388/150‑1 (MQ=255)

                                                                                                                tcttcAGCAGGTAGCCAAGCAGCCACAGGTCCAGCGGGTGGACCT‑TGGCGATGTAACCCTGGTCGGGCAGGTCGTACTTGCCGGGGCCGTAAGCTACGTACATCTCGGACAGACGGCCATCGGTCAGTTTGCCCAGCACGACCTTGTCAc                                       >  1:1100950/1‑150 (MQ=255)

                                                                                                                tcttcAGCAGGTAGCCAAGCAGCCACAGGTCCAGCGGGTGGACCT‑TGGCGATGTAACCCTGGTCGGGCAGGTCGTACTTGCCGGGGCCGTAAGCTACGTACATCTCGGACAGACGGCCATCGGTCAGTTTGCCCAGCACGACCTTGTCAc                                       >  1:4762858/1‑150 (MQ=255)

                                                                                                                 cttcAGCAGGTAGCCAAGCAGCCACAGGTCCAGCGGGTGGACCT‑TGGCGATGTAACCCTGGTCGGGCAGGTCGTACTTGCCGGGGCCGTAAGCTACGTACATCTCGGACAGACGGCCATCGGTCAGTTTGCCCAGCACGACCTTGTCAcc                                      >  1:4857392/1‑150 (MQ=255)

                                                                                                                  ttcAGCAGGTAGCCAAGCAGCCACAGGTCCAGCGGGTGGACCT‑TGGCGATGTAACCCTGGTCGGGCAGGTCGTACTTGCCGGGGCCGTAAGCTACGTACATCTCGGACAGACGGCCATCGGTCAGTTTGcc                                                        <  2:216671/131‑1 (MQ=255)

                                                                                                                  ttcAGCAGGTAGCCAAGCAGCCACAGGTCCAGCGGGTGGACCT‑TGGCGATGTAACCCTGGTCGGGCAGGTCGTACTTGCCGGGGCCGTAAGCTACGTACATCTCGGACAGACGGCCATCGGTCAGTTTGcc                                                        >  1:216671/1‑131 (MQ=255)

                                                                                                                  ttcAGCAGGTAGCCAAGCAGCCACAGGTCCAGCGGGTGGACCT‑TGGCGATGTAACCCTGGTCGGGCAGGTCGTACTTGCCGGGGCCGTAAGCTACGTACATCTCGGACAGACGGCCATCGGTCAGTTTGCCCAGCACGACCTTGTCAccc                                     >  1:1174145/1‑150 (MQ=255)

                                                                                                                   tcAGCAGGTAGCCAAGCAGCCACAGGTCCAGCGGGTGGACCT‑TGGCGATGTAACCCTGGTCGGGCAGGTCGTACTTGCCGGGGCCGTAAGCTACGTACATCTCGGACAGACTGCCATCGGTCAGTTTGCCCAGCACGACCTTGTCACCCt                                    >  1:1731896/1‑150 (MQ=255)

                                                                                                                   tcAGCAGGTAGCCAAGCAGCCACAGGTCCAGCGGGTGGACCT‑TGGCGATGTAACCCTGGTCGGACAGGTCGTACTTGCCGGGGCCGTAAGCTACGTACATCTCGGACAGACGGCCATCGGTCAGTTTGCCCAGCACGACCTTGTCACCCt                                    >  2:4292230/1‑150 (MQ=255)

                                                                                                                       cagGTAGCCAAGCAGCCACAGGTCCAGCGGGTGGACCT‑TGGCGATGTAACCCTGGTCGGGCAGGTCGTACTTGCCGGGGCCGTAAGCTACGTACATCTCGGACAGACGGCCATCGGTCAGTTTGCCCAGCACGACCTTGTCACCCTTGAg                                <  1:1598764/150‑1 (MQ=255)

                                                                                                                         ggTAGCCAAGCAGCCACAGGTCCAGCGGGTGGACCT‑TGGCGATGTAACCCTGGTCGGGCAGGTCGTACTTGCCGGGGCCGTAAGCTACGTACATCTCGGACAGACGGCCATCGGTCAGTTTGCCCAGCACGACCTTGTCACCCTTGAGGt                              >  1:898862/1‑150 (MQ=255)

                                                                                                                            aGCCAAGCAGCCACAGGTCCAGCGGGTGGACCT‑TGGCGATGTAACCCTGGTCGGGCAGGTCGTACTTGCCGGGGCCGTAAGCTACGTACATCTCGGACAGACGGCCATCGGTCAGTTTGCCCAGCACGACCTTGTCACCCTTGAGGTggg                           >  2:454406/1‑150 (MQ=255)

                                                                                                                             gCCAAGCAGCCACAGGTCCAGCGGGTGGACCT‑TGGCGATGTAACCCTGGTCGGGCAGGTCGTACTTGCCGGGGCCGTAAGCTACGTACATCTCGGACAGACGGCCATCGGTCAGTTTGCCCAGCACGACCTTGTCACCCTTGAGGTGGgc                          <  1:708591/150‑1 (MQ=255)

                                                                                                                               cAAGCAGCCACAGGTCCAGCGGGTGGACCT‑TGGCGATGTAACCCTGGTCGGGCAGGTCGTACTTGCCGGGGCCGTAAGCTCCGTACATCTCGGACAGACGGCCATCGGTCAGTTTGCCCAGCACGACCTTGTCACCCTGGAGGTGGgcac                        >  2:1433142/1‑150 (MQ=255)

                                                                                                                               cAAGCAGCCACAGGTCCAGCGGGTGGACCT‑TGGCGATGTAACCCTGGTCGGGCAGGTCGTACTTGCCGGGGCCGTAAGCTACGTACATCTCGGACAGACGGCCATCGGTCAGTTTGCCCAGCACGACCTTGTCACCCTTGAGGTGGgcac                        <  1:3355012/150‑1 (MQ=255)

                                                                                                                               cAAGCAGCCACAGGTCCAGCGGGTGGACCT‑TGGCGATGTAACCCTGGTCGGGCAGGTCGTACTTGCCGGGGCCGTAAGCTACGTACATCTCGGACAGACGGCCATCGGTCAGTTTGCCCAGCACGACCTTGTCACCCTTGAGGTGGgcac                        >  2:3557737/1‑150 (MQ=255)

                                                                                                                                   atcgccaaGGTCCACCCGCTGGACCTGTGGCGATGTAACCCTGGTCGGGCAGGTCGTACTTGCCGGGGCCGTAAGCTACGTACATCTCGGACAGACGGCCATCGGTCAGTTTGCCCAGCACGACCTTGTCACCCTTGAGGTGGgcacgca                     >  1:4335130/8‑150 (MQ=255)

                                                                                                                                      ccACAGGTCCAGCGGGTGGACCT‑TGGCGATGTAACCCTGGTCGGGCAGGTCGTACTTGCCGGGGCCGTAAGCTACGTACATCTCGGACAGACGGCCATCGGTCAGTTTGCCCAGCACGACCTTGTCACCCTTGAGGTGGGCACGCACGaa                 >  1:4030390/1‑150 (MQ=255)

                                                                                                                                       cacaGGTCCAGCGGGTGGACCT‑TGGCGATGTAACCCTGGTCGGGCAGGTCGTACTTGCCGGGGCCGTAAGCTACGTACATCTCGGACAGACGGCCATCGGTCAGTTTGCCCAGCACGACCTTGTCACCCTTGAGGTGGGCACGCACGAAc                >  2:2490101/1‑150 (MQ=255)

                                                                                                                                       cacaGGTCCAGCGGGTGGACCT‑TGGCGATGTAACCCTGGTCGGGCAGGTCGTACTTGCCGGGGCCGTAAGCTACGTACATCTCGGACAGACGGCCATCGGTCAGTTTGCCCAGCACGACCTTGTCACCCTTGAGGTGGGCACGCACGAAc                <  1:1479729/150‑1 (MQ=255)

                                                                                                                                        caaGGTCCACCCGCTGGACCTGTGGCGATGTAACCCTGGTCGGGCAGGTCGTACTTGCCGGGGCCGTAAGCTACGTACATCTCGGACAGACGGCCATCGGTCAGTTTGCCCAGCACGACCTTGTCACCCTTGAGGTGGGCACGCACGAAc                >  1:3703291/3‑150 (MQ=255)

                                                                                                                                        caaGGTCCACCCGCTGGACCTGTGGCGATGTAACCCTGGTCGGGCAGGTCGTACTTGCCGGGGCCGTAAGCTACGTACATCTCGGACAGACGGCCATCGGTCAGTTTGCCCAGCACGACCTTGTCACCCTTGAGGTGGGCACGCACGAAc                >  1:3703268/3‑150 (MQ=255)

                                                                                                                                        acaGGTCCAGCGGGTGGACCT‑TGGCGATGTAACCCTGGTCGGGCAGGTCGTACTTGCCGGGGCCGTAAGCTACGTACATCTCGGACAGACGGCCATCGGTCAGTTTGCCCAGCACGACCTTGTCACCCTTGAGGTGGGCACGCACGAAcg               >  1:836259/1‑150 (MQ=255)

                                                                                                                                        acaGGTCCAGCGGGTGGACCT‑TGGCGATGTAACCCTGGTCGGGCAGGTCGTACTTGCCGGGGCCGTAAGCTACGTACATCTCGGACAGACGGCCATCGGTCAGTTTGCCCAGCACGACCTTGTCACCCTTGAGGTGGGCACGCACGAAcg               >  1:1186476/1‑150 (MQ=255)

                                                                                                                                        acaGGTCCAGCGGGTGGACCT‑TGGCGATGTAACCCTGGTCGGGCAGGTCGTACTTGCCGGGGCCGTAAGCTACGTACATCTCGGACAGACGGCCATCGGTCAGTTTGCCCAGCACGACCTTGTCACCCTTGAGGTGGGCACGCACGAAcg               >  1:131608/1‑150 (MQ=255)

                                                                                                                                         caGGTCCAGCGGGTGGACCT‑TGGCGATGTAACCCTGGTCGGGCAGGTCGTACTTGCCGGGGCCGTAAGCTACGTACATCTCGGACAGACGGCCCTCGGTCCGTGTGCCCAGCACGACCGTGTCACCCTGGAGGTGGGCACGCACGAAcgc              >  2:1748411/1‑150 (MQ=255)

                                                                                                                                         caGGTCCAGCGGGTGGACCT‑TGGCGATGTAACCCTGGTCGGGCAGGTCGTACTTGCCGGGGCCGTAAGCTACGTACATCTCGGACAGACGGCCATCGGTCAGTTTGCCCAGCACGACCTTGTCACCCTTGAGGTGGGCACGCACGAAcgc              >  1:957538/1‑150 (MQ=255)

                                                                                                                                            gTCCAGCGGGTGGACCT‑TGGCGATGTAACCCTGGTCGGGCAGGTCGTACTTGCCGGGGCCGTAAGCTACGTACATCTCGGACAGACGGCCATCGGTCAGTTTGCCCAGCACGACCTTGTCACCCTTGAGGTGGGCACGCACGAACGCGtt           <  1:2930846/150‑1 (MQ=255)

                                                                                                                                              ccAGCGGGTGGACCT‑TGGCGATGTAACCCTGGTCGGGCAGGTCGTACTTGCCGGGGCCGTAAGCTACGTACATCTCGGACAGACGGCCATCGGTCAGTTTGCCCAGCACGACCTTGTCACCCTTGAGGTGGGCACGCACGAACGCGTTGa         <  2:2591284/150‑1 (MQ=255)

                                                                                                                                               cAGCGGGTGGACCT‑TGGCGATGTAACCCTGGTCGGGCAGGTCGTACTTGCCGGGGCCGTAAGCTACGTACATCTCGGACAGACGGCCATCGGTCAGTTTGCCCAGCACGACCTTGTCACCCTTGAGGTGGGCACGCACGAACGCGTTGaa        <  2:1969860/150‑1 (MQ=255)

                                                                                                                                               cAGCGGGTGGACCT‑TGGCGATGTAACCCTGGTCGGGCAGGTCGTACTTGCCGGGGCCGTAAGCTACGTACATCTCGGACAGACGGCCATCGGTCAGTTTGCCCAGCACGACCTTGTCACCCTTGAGGTGGGCACGCACGAACGCGTTGaa        <  2:836259/150‑1 (MQ=255)

                                                                                                                                               cAGCGGGTGGACCT‑TGGCGATGTAACCCTGGTCGGGCAGGTCGTACTTGCCGGGGCCGTAAGCTACGTACATCTCGGACAGACGGCCATCGGTCAGTTTGCCCAGCACGACCTTGTCACCCTTGAGGTGGGCACGCACGAACGCGTTGaa        <  1:1958294/150‑1 (MQ=255)

                                                                                                                                               cAGCGGGTGGACCT‑TGGCGATGTAACCCTGGTCGGGCAGGTCGTACTTGCCGGGGCCGTAAGCTACGTACATCTCGGACAGACGGCCATCGGTCAGTTTGCCCAGCACGACCTTGTCACCCTTGAGGTGGGCACGCACGAACGCGTTGaa        <  2:4654377/150‑1 (MQ=255)

                                                                                                                                               cAGCGGGTGGACCT‑TGGCGATGTAACCCTGGTCGGGCAGGTCGTACTTGCCGGGGCCGTAAGCTACGTACATCTCGGACAGACGGCCATCGGTCAGTTTGCCCAGCACGACCTTGTCACCCTTGAGGTGGGCACGCACGAACGCGTTGaa        <  1:1797248/150‑1 (MQ=255)

                                                                                                                                                aGCGGGTGGACCT‑TGGCGATGTAACCCTGGTCGGGCAGGTCGTACTTGCCGGGGCCGTAAGCTACGTACATCTCGGACAGACGGCCATCGGTCAGTTTGCCCAGCACGACCTTGTCACCCTTGAGGTGGGCACGCACGAACGCGTTGAAg       <  2:1975573/150‑1 (MQ=255)

                                                                                                                                                aGCGGGTGGACCT‑TGGCGATGTAACCCTGGTCGGGCAGGTCGTACTTGCCGGGGCCGTAAGCTACGTACATCTCGGACAGACGGCCATCGGTCAGTTTGCCCAGCACGACCTTGTCACCCTTGAGGTGGGCACGCACGAACGCGTTGAAg       >  2:2026441/1‑150 (MQ=255)

                                                                                                                                                 gCGGGTGGACCT‑TGGCGATGTAACCCTGGTCGGGCAGGTCGTACTTGCCGGGGCCGTAAGCTACGTACATCTCGGACAGACGGCCATCGGTCAGTTTGCCCAGCACGACCTTGTCACCCTTGAGGTGGGCACGCACGAACGCGTTGAAgg      <  2:2160202/150‑1 (MQ=255)

                                                                                                                                                 gCGGGTGGACCT‑TGGCGATGTAACCCTGGTCGGGCAGGTCGTACTTGCCGGGGCCGTAAGCTACGTACATCTCGGACAGACGGCCATCGGTCAGTTTGCCCAGCACGACCTTGTCACCCTTGAGGTGGGCACGCACGAACGCGTTGAAgg      >  2:743363/1‑150 (MQ=255)

                                                                                                                                                     gTGGACCT‑TGGCGATGTAACCCTGGTCGGGCAGGTCGTACTTGCCGGGGCCGTAAGCTACGTACATCTCGGACAGACGGCCATCGGTCAGTTTGCCCAGCACGACCTTGTCACCCTTGAGGTGGGCACGCACGAACGCGTTGAAGGTTTc  <  1:343336/150‑1 (MQ=255)

CGGATGCGGCTGTCACGGGCGCCCTGGTGGCGGCTCTTGAACAGCCAGCCGTATACCTCCTGGCGTTCGAAACGACTGGCGTTGATCATCTCGGTCAGTGTCGAGCTCGGGTTCTTCAGCAGGTAGCCAAGCAGCCACAGGTCCAGCGGGTGGACCT‑TGGCGATGTAACCCTGGTCGGGCAGGTCGTACTTGCCGGGGCCGTAAGCTACGTACATCTCGGACAGACGGCCATCGGTCAGTTTGCCCAGCACGACCTTGTCACCCTTGAGGTGGGCACGCACGAACGCGTTGAAGGTTTC  >  NC_002947/1206554‑1206852

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