BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,110,313	G→C	6.7%	W23C (TGG→TGC)	PP_5358 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,110,313	0	G	C	6.7%	609.7 / 26.2	180	W23C (TGG→TGC)	PP_5358	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (80/88); new base C (4/8); total (84/96)
Fisher's exact test for biased strand distribution p-value = 3.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

Only 100 of 182 total aligned reads displayed.

CGCAGCGCAGTATGAACTTGATTGCTACCTGTACGACAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATGCCAATGGCCGCAAGCGCCTGGCACGCGTCTA  >  NC_002947/6110165‑6110462

cgcagcgcagTATGAACTTGAGTGCTACCTGTACGACAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGc                                                                                                                                                      <  2:801538/150‑1 (MQ=255)

    gcgcagTATGAACTTGATTGCTACCTGTACGACAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACaa                                                                                                                                                  <  2:2133628/150‑1 (MQ=255)

        agTATGAACTTGATTGCTACCTGTACGACAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGCCACAACCag                                                                                                                                              <  1:765512/150‑2 (MQ=255)

        agTATGAACTTGATTGCTACCTGTACGACAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGCCACAACCag                                                                                                                                              <  2:1872034/150‑2 (MQ=255)

        agTATGAACTTGATTGCTACCTGTACGACAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGCCACAACCag                                                                                                                                              <  1:2883076/150‑2 (MQ=255)

         gTATGAACTTGATTGCTACCTGTACGACAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATg                                                                                                                                             <  1:4419235/150‑1 (MQ=255)

         gTATGAACTTGATTGCTACCTGTACGACAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATg                                                                                                                                             <  1:3747325/150‑1 (MQ=255)

         gTATGAACTTGATTGCTACCTGTACGACAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATg                                                                                                                                             <  2:2327647/150‑1 (MQ=255)

          tATGAACTTGATTGCTACCTGTACGACAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCATGGTTGTGCCACAACCaggc                                                                                                                                            <  2:1922381/150‑4 (MQ=255)

              aaCTTGATTGCTACCTGTACGACAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTg                                                                                                                                        >  1:1057399/1‑150 (MQ=255)

              aaCTTGATTGCTACCTGTACGACAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTg                                                                                                                                        >  1:1816718/1‑150 (MQ=255)

              aaCTTGATTGCTACCTGTACGACAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTg                                                                                                                                        >  2:2902625/1‑150 (MQ=255)

                 ttGATTGCTACCTGTACGACAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGt                                                                                                                                     >  1:416141/1‑150 (MQ=255)

                 ttGATTGCTACCTGTACGACAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATCGGTTGCGt                                                                                                                                     >  1:4248395/1‑150 (MQ=255)

                     ttGCTACCTGTACGACAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGaaa                                                                                                                                 <  2:1061139/150‑1 (MQ=255)

                        cTACCTGTACGACAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAgg                                                                                                                              >  1:4183544/1‑150 (MQ=255)

                         tACCTGTACGACAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATCCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGc                                                                                                                             >  2:3798646/1‑150 (MQ=255)

                          aCCTGTACGACAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGcc                                                                                                                            <  1:4333610/150‑1 (MQ=255)

                           ccTGTACGACAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCGGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGACAAGGccc                                                                                                                           >  1:1570690/1‑150 (MQ=255)

                            cTGTACGACAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCt                                                                                                                          <  1:2848396/150‑1 (MQ=255)

                             tGTACGACAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTg                                                                                                                         <  2:4036931/150‑1 (MQ=255)

                              gTACGACAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTgg                                                                                                                        >  2:280592/1‑150 (MQ=255)

                                aCGACAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGac                                                                                                                      >  2:2506750/1‑149 (MQ=255)

                                  gACAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACg                                                                                                                    <  1:3981405/150‑1 (MQ=255)

                                  gACAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGATTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACg                                                                                                                    <  2:2658760/150‑1 (MQ=255)

                                   aCAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACgg                                                                                                                   <  1:4456751/150‑1 (MQ=255)

                                           ggAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACa                                                                                                           <  1:471457/150‑1 (MQ=255)

                                                   gTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTgcg                                                                                                   >  1:1577592/1‑150 (MQ=255)

                                                                  ccTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTgg                                                                                    >  2:1396040/1‑150 (MQ=255)

                                                                   cTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGa                                                                                   >  1:432430/1‑150 (MQ=255)

                                                                   cTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGa                                                                                   >  1:3187970/1‑150 (MQ=255)

                                                                   cTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGa                                                                                   >  1:584625/1‑150 (MQ=255)

                                                                    tGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAg                                                                                  >  2:4192841/1‑150 (MQ=255)

                                                                     gTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGt                                                                                 <  1:2264604/150‑1 (MQ=255)

                                                                      tttGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGtt                                                                                <  1:3867864/150‑1 (MQ=255)

                                                                        tGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGc                                                                              >  2:1126422/1‑150 (MQ=255)

                                                                        tGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGc                                                                              <  2:470964/150‑1 (MQ=255)

                                                                         ggAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCt                                                                             <  1:2985772/150‑1 (MQ=255)

                                                                          gAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCtt                                                                            <  1:3381971/150‑1 (MQ=255)

                                                                          gAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCtt                                                                            <  1:3094520/150‑1 (MQ=255)

                                                                          gAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCtt                                                                            >  1:2686120/1‑150 (MQ=255)

                                                                             cATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGAc                                                                         <  2:3444934/150‑1 (MQ=255)

                                                                               tatGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGa                                                                       <  2:2554025/150‑1 (MQ=255)

                                                                               tatGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGa                                                                       >  1:2260636/1‑150 (MQ=255)

                                                                                atGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGAt                                                                      <  2:753199/150‑1 (MQ=255)

                                                                                 tGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATg                                                                     >  1:2787110/1‑150 (MQ=255)

                                                                                 tGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATg                                                                     >  1:3054537/1‑150 (MQ=255)

                                                                                 tGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATg                                                                     <  1:2234427/150‑1 (MQ=255)

                                                                                  gCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGc                                                                    <  2:1694956/150‑1 (MQ=255)

                                                                                  gCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGc                                                                    >  1:1075303/1‑150 (MQ=255)

                                                                                    aaaTGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCg                                                                  >  2:2108477/1‑150 (MQ=255)

                                                                                    aaaTGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCg                                                                  >  1:3103804/1‑150 (MQ=255)

                                                                                        gTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAgc                                                              >  1:1064735/1‑150 (MQ=255)

                                                                                          tGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAgcgc                                                            >  2:3514833/1‑150 (MQ=255)

                                                                                              ccTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCaa                                                        <  2:3945762/150‑1 (MQ=255)

                                                                                               cTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAg                                                       >  1:2512574/1‑150 (MQ=255)

                                                                                                tGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGc                                                      <  1:2427541/150‑1 (MQ=255)

                                                                                                tGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGc                                                      <  2:3327687/150‑1 (MQ=255)

                                                                                                tGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGc                                                      <  1:4183080/150‑1 (MQ=255)

                                                                                                 ggAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCg                                                     >  1:1905246/1‑150 (MQ=255)

                                                                                                   agaACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTa                                                   <  2:3668284/150‑1 (MQ=255)

                                                                                                       ccTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGc                                               >  1:2125900/1‑150 (MQ=255)

                                                                                                           ttCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCttt                                           <  1:2392216/150‑1 (MQ=255)

                                                                                                           ttCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCttt                                           >  2:4122526/1‑150 (MQ=255)

                                                                                                              gTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTc                                        >  2:1848945/1‑150 (MQ=255)

                                                                                                                ccTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCg                                      <  2:3311939/150‑1 (MQ=255)

                                                                                                                ccTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCg                                      >  2:1026806/1‑150 (MQ=255)

                                                                                                                 cTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCgg                                     <  2:2429971/150‑1 (MQ=255)

                                                                                                                 cTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCATTTGTCCgg                                     >  1:4502665/1‑150 (MQ=255)

                                                                                                                  tGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCggg                                    >  2:951660/1‑150 (MQ=255)

                                                                                                                   gATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGa                                   >  2:3300643/1‑150 (MQ=255)

                                                                                                                   gATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGa                                   >  1:1236656/1‑150 (MQ=255)

                                                                                                                       cTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATGcc                               >  2:3469380/1‑150 (MQ=255)

                                                                                                                        tGGTAGCCACGGCAGGCGCCT‑GTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATGCCaa                             >  1:81740/1‑150 (MQ=255)

                                                                                                                         ggTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATGCCaa                             >  1:2908178/1‑150 (MQ=255)

                                                                                                                          gTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATGCCAAt                            <  1:4101100/150‑1 (MQ=255)

                                                                                                                              ccACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATGCCAATGGcc                        <  2:3453493/150‑1 (MQ=255)

                                                                                                                                aCGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATGCCAATGGCCGc                      >  1:122424/1‑150 (MQ=255)

                                                                                                                                aCGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATGCCAATGGCCGc                      >  2:3098356/1‑150 (MQ=255)

                                                                                                                                aCGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATGCCAATGGCCGc                      >  1:3633711/1‑150 (MQ=255)

                                                                                                                                aCGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATGCCAATGGCCGc                      >  2:1342199/1‑150 (MQ=255)

                                                                                                                                  ggCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATGCCAATGGCCGCaa                    <  2:60514/150‑1 (MQ=255)

                                                                                                                                  ggCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATGCCAATGGCCGCaa                    <  1:261080/150‑1 (MQ=255)

                                                                                                                                  ggCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATGCCAATGGCCGCaa                    <  2:1598260/150‑1 (MQ=255)

                                                                                                                                    cAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATGCCAATGGCCGCAAgc                  >  2:4045088/1‑150 (MQ=255)

                                                                                                                                     aGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATGCCAATGGCCGCAAgcg                 >  1:3439690/1‑150 (MQ=255)

                                                                                                                                     aGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATGCCAATGGCCGCAAgcg                 <  1:498343/150‑1 (MQ=255)

                                                                                                                                       gcgcCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATGCCAATGGCCGCAAGCGcc               <  2:2599296/150‑1 (MQ=255)

                                                                                                                                        cccatGGTTGTGCCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATGCCAATGGCCGCAAGCGCCt              >  2:282775/5‑150 (MQ=255)

                                                                                                                                        cccatGGTTGTGCCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATGCCAATGGCCGCAAGCGCCt              >  2:282766/5‑150 (MQ=255)

                                                                                                                                         gcCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATGCCAATGGCCGCAAGCGCCTg             <  2:461650/150‑1 (MQ=255)

                                                                                                                                         gcCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATGCCAATGGCCGCAAGCGCCTg             <  1:1754690/150‑1 (MQ=255)

                                                                                                                                         gcCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATGCCAATGGCCGCAAGCGCCTg             <  1:3661134/150‑1 (MQ=255)

                                                                                                                                         gcCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGAAGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATGCCAATGGCCGCAAGCGCCTg             >  2:3878586/1‑150 (MQ=255)

                                                                                                                                          catGGTTGTGCCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATGCCAATGGCCGCAAGCGCCTgg            >  2:472721/3‑150 (MQ=255)

                                                                                                                                          catGGTTGTGCCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATGCCAATGGCCGCAAGCGCCAgg            >  1:3345833/3‑150 (MQ=255)

                                                                                                                                              gTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATc                                                 >  2:3599780/1‑109 (MQ=255)

                                                                                                                                               ttGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATGCCAATGGCCGCAAGCGCCTGGCAcgc       <  1:2055669/150‑1 (MQ=255)

                                                                                                                                               ttGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATGCCAATGGCCGCAAGCGCCTGGCAcgc       <  2:1799632/150‑1 (MQ=255)

                                                                                                                                                    gCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATGCCAATGGCCGCAAGCGCCTGGCACGCGTCTa  >  2:439535/1‑150 (MQ=255)

CGCAGCGCAGTATGAACTTGATTGCTACCTGTACGACAGGCCAGGAACGCAGTTAAGCTGGATCTACCTGTTTGGAGCATATGCAAATGTTGACCCTGGAGAACCTCTTCGTCCTGATGCTGGTAGCCACGGCAGGCGCCTGGTTGTGGCACAACCATGGGTTGCGTGAAAAGGCCCTGGAACGGGTCAAACAGCATTGCGCCAAACTCGACCTGGAGTTGCTTGACGATGCCGTAGCGCTCAAGCGTATCGCCTTTGTCCGGGATGCCAATGGCCGCAAGCGCCTGGCACGCGTCTA  >  NC_002947/6110165‑6110462

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