BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,113,166	G→C	5.6%	P108A (CCA→GCA)	PP_5361 ←	protein P47K

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,113,166	0	G	C	5.6%	540.4 / 16.2	161	P108A (CCA→GCA)	PP_5361	protein P47K
Reads supporting (aligned to +/- strand): ref base G (67/85); new base C (4/5); total (71/90)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

Only 100 of 162 total aligned reads displayed.

GTTTCGCCCCGGCTGGCCAGGCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATGCAGCCGTTGCTCATCTCGACCAGTTTCTCTTCAGCGCGGT  >  NC_002947/6113018‑6113315

gTTTCGCCCCGGCTGGCCAGGCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGc                                                                                                                                                      >  2:2421937/1‑150 (MQ=255)

 tttCGCCCCGGCTGGCCAGGCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCa                                                                                                                                                     <  1:2055659/150‑1 (MQ=255)

     gCCCCGGCTGGCCAGGCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTgg                                                                                                                                                 <  2:69256/150‑1 (MQ=255)

     gCCCCGGCTGGCCAGGCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTgg                                                                                                                                                 >  1:2869457/1‑150 (MQ=255)

     gCCCCGGCTGGCCAGGCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTgg                                                                                                                                                 >  1:677214/1‑150 (MQ=255)

     gCCCCGGCTGGCCAGGCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTgg                                                                                                                                                 <  2:69249/150‑1 (MQ=255)

      ccccGGCTGGCCAGGCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGc                                                                                                                                                >  2:4815039/1‑150 (MQ=255)

      ccccGGCTGGCCAGGCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGc                                                                                                                                                >  1:1803131/1‑150 (MQ=255)

      ccccGGCTGGCCAGGCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGc                                                                                                                                                <  1:4727794/150‑1 (MQ=255)

      ccccGGCTGGCCAGGCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGc                                                                                                                                                <  1:4727783/150‑1 (MQ=255)

       cccGGCTGGCCAGGCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCt                                                                                                                                               >  1:4525728/1‑150 (MQ=255)

       cccGGCTGGCCAGGCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCt                                                                                                                                               <  1:2480293/150‑1 (MQ=255)

        ccGGCTGGCCAGGCTGTCCGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGCCAGTGGCtg                                                                                                                                              <  1:1809242/150‑2 (MQ=255)

         cGGCTGGCCAGGCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCa                                                                                                                                             <  1:4553220/150‑1 (MQ=255)

         cGGCTGGCCAGGCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGCCAGTGGCtga                                                                                                                                             <  1:3025516/150‑3 (MQ=255)

           gCTGGCCAGGCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGTGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCaga                                                                                                                                           <  2:2946456/150‑1 (MQ=255)

           gCTGGCCAGGCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCaga                                                                                                                                           <  1:2371594/150‑1 (MQ=255)

           gCTGGCCAGGCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCaga                                                                                                                                           <  1:1202652/150‑1 (MQ=255)

              ggCCAGGCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGAt                                                                                                                                        <  2:1672795/150‑1 (MQ=255)

                ccAGGCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATAc                                                                                                                                      >  1:1251628/1‑150 (MQ=255)

                  aGGCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCg                                                                                                                                    >  2:1403511/1‑150 (MQ=255)

                   ggCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCgg                                                                                                                                   <  1:4215942/150‑1 (MQ=255)

                    gCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGt                                                                                                                                  <  1:2286451/150‑1 (MQ=255)

                    gCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGt                                                                                                                                  <  2:677214/150‑1 (MQ=255)

                    gCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGt                                                                                                                                  <  2:395105/150‑1 (MQ=255)

                    gCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGt                                                                                                                                  <  2:3515697/150‑1 (MQ=255)

                    gCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGt                                                                                                                                  <  1:3646931/150‑1 (MQ=255)

                       gTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGa                                                                                                                               <  2:1371458/150‑1 (MQ=255)

                          tgctgcCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTc                                                                                                                            >  2:4967694/1‑150 (MQ=255)

                          tgctgcCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTc                                                                                                                            <  1:138548/150‑1 (MQ=255)

                          tgctgcCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTc                                                                                                                            >  1:1453331/1‑150 (MQ=255)

                          tgctgcCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTc                                                                                                                            >  1:2218249/1‑150 (MQ=255)

                              gcCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATc                                                                                                                        <  1:3338480/150‑1 (MQ=255)

                              gcCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATc                                                                                                                        <  1:2128498/150‑1 (MQ=255)

                                 tGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGc                                                                                                                     >  1:3444331/1‑150 (MQ=255)

                                 tGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGc                                                                                                                     >  2:3904362/1‑150 (MQ=255)

                                       tCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCCaataa                                                                                                               >  2:3319862/1‑150 (MQ=255)

                                         gcgcAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATc                                                                                                             <  1:99535/150‑1 (MQ=255)

                                         gcgcAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATc                                                                                                             <  1:1472853/150‑1 (MQ=255)

                                             aGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGaaa                                                                                                         <  2:3067768/150‑1 (MQ=255)

                                             aGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGaaa                                                                                                         <  1:4643521/150‑1 (MQ=255)

                                             aGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGaaa                                                                                                         >  1:2217419/1‑150 (MQ=255)

                                              ggAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTTAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAAc                                                                                                        <  1:2997320/150‑1 (MQ=255)

                                                      aaGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCctct                                                                                                >  2:1552519/1‑150 (MQ=255)

                                                       aGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCa                                                                                                                       >  2:5108813/1‑126 (MQ=255)

                                                       aGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCa                                                                                                                       <  1:5108813/126‑1 (MQ=255)

                                                        gCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTc                                                                                              >  2:1797728/1‑150 (MQ=255)

                                                         ccATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCg                                                                                             <  1:3818447/150‑1 (MQ=255)

                                                          cATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCgg                                                                                            >  1:1956851/1‑150 (MQ=255)

                                                          cATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCgg                                                                                            <  2:4224014/150‑1 (MQ=255)

                                                           aTCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGc                                                                                           <  2:3679931/150‑1 (MQ=255)

                                                                ccACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTc                                                                                      <  2:1134222/150‑1 (MQ=255)

                                                                          cATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAAcctc                                                                            >  1:2982514/1‑150 (MQ=255)

                                                                               tgtCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGa                                                                       <  2:4394188/150‑1 (MQ=255)

                                                                                   cAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAg                                                                   <  1:3462964/150‑1 (MQ=255)

                                                                                   cAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAg                                                                   <  2:150133/150‑1 (MQ=255)

                                                                                    aGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAgg                                                                  >  2:4461513/1‑150 (MQ=255)

                                                                                    aGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAgg                                                                  <  2:55317/150‑1 (MQ=255)

                                                                                        gcgcCATATCGGACAGACTGCGACCCTGCTCGTCGCAGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCtt                                                              >  2:4967284/1‑150 (MQ=255)

                                                                                        gcgcCATATCGGACAGACTGCGACCCTGCTCGTCGCAGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCtt                                                              >  2:2397913/1‑150 (MQ=255)

                                                                                         cgcCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTc                                                             <  1:2234241/150‑1 (MQ=255)

                                                                                           ccATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTcgc                                                           <  2:1956851/150‑1 (MQ=255)

                                                                                           ccATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTcgc                                                           <  2:1979352/150‑1 (MQ=255)

                                                                                           ccATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTcgc                                                           <  1:2354737/150‑1 (MQ=255)

                                                                                             atatCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCgcgc                                                         <  2:2628141/150‑1 (MQ=255)

                                                                                              tatCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCa                                                        >  2:2350671/1‑150 (MQ=255)

                                                                                               atCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCaa                                                       <  1:76551/150‑1 (MQ=255)

                                                                                                tCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAAc                                                      <  2:4381807/150‑1 (MQ=255)

                                                                                                tCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAAc                                                      <  2:1803131/150‑1 (MQ=255)

                                                                                                  ggACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATc                                                                                                             >  2:83763/1‑93 (MQ=255)

                                                                                                      aGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTgcag                                                <  1:4046291/150‑1 (MQ=255)

                                                                                                          tGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCaaa                                            <  1:2719812/150‑1 (MQ=255)

                                                                                                          tGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCaaa                                            >  2:1314717/1‑150 (MQ=255)

                                                                                                            cGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATg                                          <  2:4721842/150‑1 (MQ=255)

                                                                                                            cGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACGGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATg                                          <  2:4721929/150‑1 (MQ=255)

                                                                                                              aCCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTTAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATGCa                                        <  2:1226406/150‑1 (MQ=255)

                                                                                                              aCCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATGCa                                        <  1:3568642/150‑1 (MQ=255)

                                                                                                                  tGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCTGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATGCAGCCg                                    >  2:313404/1‑150 (MQ=255)

                                                                                                                  tGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATGCAGCCg                                    >  1:4574891/1‑150 (MQ=255)

                                                                                                                   gCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATGCAGCCGt                                   <  2:2557407/150‑1 (MQ=255)

                                                                                                                    cTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATGCAGCCGtt                                  <  1:4992397/150‑1 (MQ=255)

                                                                                                                    cTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATGCAGCCGtt                                  <  1:2421937/150‑1 (MQ=255)

                                                                                                                          gcgGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATGCAGCCGTTGCTCAt                            <  2:3810520/150‑1 (MQ=255)

                                                                                                                          gcgGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATGCAGCCGTTGCTCAt                            <  2:3654191/150‑1 (MQ=255)

                                                                                                                             gAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATGCAGCCGTTGCTCAtctc                         <  1:3360083/150‑1 (MQ=255)

                                                                                                                             gAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCCCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATGCAGCCGTTGCTCAtctc                         >  2:3134057/1‑150 (MQ=255)

                                                                                                                               aaGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCGCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATGCAGCCGTTGCTCATCTCGa                       >  2:2845963/1‑150 (MQ=255)

                                                                                                                                aGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATGCAGCCGTTGCTCATCTCGAc                      <  1:2245150/150‑1 (MQ=255)

                                                                                                                                   gtatctctgAGCCACTGCCAGTGGCTGAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATGCAGCCGTTGCTCATCTCGACCAg                   >  2:100784/6‑150 (MQ=255)

                                                                                                                                   gAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATGCAGCCGTTGCTCATCTCGACCAg                   >  1:3958047/1‑150 (MQ=255)

                                                                                                                                    aaGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAATCGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATGCAGCCGTTGCTCATCTCGACCAGt                  >  2:3676537/1‑150 (MQ=255)

                                                                                                                                    aaGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATGCAGCCGTTGCTCATCTCGACCAGt                  >  2:1298496/1‑150 (MQ=255)

                                                                                                                                      ggTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATGCAGCCGTTGCTCATCTCGACCAGttt                <  1:1071597/150‑1 (MQ=255)

                                                                                                                                         ctcAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATGCAGCCGTTGCTCATCTCGACCAGTTtctc             <  2:4525728/150‑1 (MQ=255)

                                                                                                                                          tcAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATGCAGCCGTTGCTCATCTCGACCAGTTTctct            >  1:810482/1‑150 (MQ=255)

                                                                                                                                            aGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATGCAGCCGTTGCTCATCTCGACCAGTTTCTCTTc          <  1:1050040/150‑1 (MQ=255)

                                                                                                                                             gCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATGCAGCCGTTGCTCATCTCGACCAGTTTCTCTTCa         <  1:3962730/150‑1 (MQ=255)

                                                                                                                                             gCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATGCAGCCGTTGCTCATCTCGACCAGTTTCTCTTCa         <  1:954679/150‑1 (MQ=255)

                                                                                                                                              ccACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATGCAGCCGTTGCTCATCTCGACCAGTTTCTCTTCAg        <  1:3757777/150‑1 (MQ=255)

                                                                                                                                               cACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATGCAGCCGTTGCTCATCTCGACCAGTTTCTCTTCAgc       <  2:4633000/150‑1 (MQ=255)

GTTTCGCCCCGGCTGGCCAGGCTGTCTGCTGCCTGATAGTCGCGCAGGAAGTTCAAGCCATCGACCACGGTGACCATGGTGTCCAGTCGCGCCATATCGGACAGACTGCGACCCTGCTCGTCGCGGAAAGTGAAGGTCTCAGCCACTGGCAGTGGCTCAGAGATACCGGTGGACTCGATCAGCAAATAATCGAAACGCCCCTCTTCGGCCAGTCGCGCAACCTCCTCGAGCAGGTCTTCGCGCAACGTGCAGCAAATGCAGCCGTTGCTCATCTCGACCAGTTTCTCTTCAGCGCGGT  >  NC_002947/6113018‑6113315

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