BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,435,016	C→G	5.3%	R100R (CGG→CGC)	PP_4769 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,435,016	0	C	G	5.3%	1224.4 / 40.1	359	R100R (CGG→CGC)	PP_4769	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (166/174); new base G (10/9); total (176/183)
Fisher's exact test for biased strand distribution p-value = 8.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 367 total aligned reads displayed.

GCACGTAGCGGGTGGCAGCACGGATCAGTGACAACGCCGTCAGCCCCGGCACGCCATACACCACCGCTTCGGTGCCATGGCGCTGGATTACCCGCTCCAGCAGCAAGTCGAAATCACCGTCGCCGGAGGCCAGCACCACCTGGTCGACCCGGCTGGCCGCGTCGATCACATCCAGGGTAATGCCCACATCCCAGTCACCTTTGGCCGAGCCATCACTGCGCTGGATGAACGGTTTGAGGCGCACCTCGAAACCGAGGTTGCGCAAGATCTGCTGGAACTGCTGTTGTTTGCTGTCGCC  >  NC_002947/5434868‑5435165

gCACGTAGCGGGTGGCAGCACGGATCAGTGACAACGCCGTCAGCCCCGGCACGCCATACACCACCGCTTCGGTGCCATGGCGCTGGATTACCCGCTCCAGCAGCAAGTCGAAATCACCGTCGCCGGAGGCCAGCACCACCTGGTCGAccc                                                                                                                                                      <  1:2738488/150‑1 (MQ=255)

  aCGTAGCGGGTGGCAGCACGGATCAGTGACAACGCCGTCAGCCCCGGCACGCCATACACCACCGCTTCGGTGCCATGGCGCTGGATTACCCGCTCCAGCAGCAAGTCGAAATCACCGTCGCCGGAGGCCAGCACCACCTGGTCGACCCgg                                                                                                                                                    >  2:5490576/1‑150 (MQ=255)

  aCGTAGCGGGTGGCAGCACGGATCAGTGACAACGCCGTCAGCCCCGACACGCCATACACCACCGCTTCGGTGCCATGGCGCTGGATTACCCGCTCCAGCAGCAAGTCGAAATCACCGTCGCCGGAGGCCAGCACCACCTGGTCGACCCgg                                                                                                                                                    <  2:829026/150‑1 (MQ=255)

     tAGCGGGTGGCAGCACGGATCAGTGACAACGCCGTCAGCCCCGGCACGCCATACACCACCGCTTCGGTGCCATGGCGCTGGATTACCCGCTCCAGCAGCAAGTCGAAATCACCGTCGCCGGAGGCCAGCACCACCTGGTCGACCCGGCTg                                                                                                                                                 <  2:2657266/150‑1 (MQ=255)

      aGCGGGTGGCAGCACGGATCAGTGACAACGCCGTCAGCCCCGGCACGCCATACACCACCGCTTCGGTGCCATGGCGCTGGATTACCCGCTCCAGCAGCAAGTCGAAATCACCGTCGCCGGAGGCCAGCACCACCTGGTCGACCCGGCTgg                                                                                                                                                >  1:1027065/1‑150 (MQ=255)

      aGCGGGTGGCAGCACGGATCAGTGACAACGCCGTCAGCCCCGGCACGCCATACACCACCGCTTCGGTGCCATGGCGCTGGATTACCCGCTCCAGCAGCAAGTCGAAATCACCGTCGCCGGAGGCCAGCACCACCTGGTCGACCCGGCTgg                                                                                                                                                <  1:6109066/150‑1 (MQ=255)

       gCGGGTGGCAGCACGGATCAGTGACAACGCCGTCAGCCCCGGCACGCCATACACCACCGCTTCGGTGCCATGGCGCTGGATTACCCGCTCCAGCAGCAAGTCGAAATCACCGTCGCCGGAGGCCAGCACCACCTGGTCGACCCGGCTGGc                                                                                                                                               <  1:4138136/150‑1 (MQ=255)

         gggTGGCAGCACGGATCAGTGACAACGCCGTCAGCCCCGGCACGCCATACACCACCGCTTCGGTGCCATGGCGCTGGATTACCCGCTCCAGCAGCAAGTCGAAATCACCGTCGCCGGAGGCCAGCACCACCTGGTCGACCCGGCTGGCcg                                                                                                                                             >  2:7706124/1‑150 (MQ=255)

         gggTGGCAGCACGGATCAGTGACAACGCCGTCAGCCCCGGCACGCCATACACCACCGCTTCGGTGCCATGGCGCTGGATTACCCGCTCCAGCAGCAAGTCGAAATCACCGTCGCCGGAGGCCAGCACCACCTGGTCGACCCGGCTGGCcg                                                                                                                                             >  2:183127/1‑150 (MQ=255)

          ggTGGCAGCACGGATCAGTGACAACGCCGTCAGCCCCGGCACGCCATACACCACCGCTTCGGTGCCATGGCGCTGGATTACCCGCTCCAGCAGCAAGTCGAAATCACCGTCGCCGGAGGCCAGCACCACCTGGTCGACCCGGCTGGCcgc                                                                                                                                            <  1:812365/150‑1 (MQ=255)

            tGGCAGCACGGATCAGTGACAACGCCGTCAGCCCCGGCACGCCATACACCACCGCTTCGGTGCCATGGCGCTGGATTACCCGCTCCAGCAGCAAGTCGAAATCACCGTCGCCGGAGGCCAGCACCACCTGGTCGACCCGGCTGGCCGCGt                                                                                                                                          <  1:3818773/150‑1 (MQ=255)

             ggCAGCACGGATCAGTGACAACGCCGTCAGCCCCGGCACGCCATACACCACCGCTTCGGTGCCATGGCGCTGGATTACCCGCTCCAGCAGCAAGTCGAAATCACCGTCGCCGGAGGCCAGCACCACCTGGTCGACCCGGCTGGCCGCGTc                                                                                                                                         <  1:965254/150‑1 (MQ=255)

              gcagcaCGGATCAGTGACAACGCCGTCAGCCCCGGCACGCCATACACCACCGCTTCGGTGCCATGGCGCTGGATTACCCGCTCCAGCAGCAAGTCGAAATCACCGTCGCCGGAGGCCAGCACCACCTGGTCGACCCGGCTGGCCGCGTCg                                                                                                                                        >  1:1436623/1‑150 (MQ=255)

                agcaCGGATCAGTGACAACGCCGTCAGCCCCGGCACGCCATACACCACCGCTTCGGTGCCATGGCGCTGGATTACCCGCTCCAGCAGCAAGTCGAAATCACCGTCGCCGGAGGCCAGCACCACCTGGTCGACGCGGCCAGCCGGGTCGac                                                                                                                                      <  1:3283680/150‑2 (MQ=255)

                agcaCGGATCAGTGACAACGCCGTCAGCCCCGGCACGCCATACACCACCGCTTCGGTGCCATGGCGCTGGATTACCCGCTCCAGCAGCAAGTCGAAATCACCGTCGCCGGAGGCCAGCACCACCTGGTCGACCCGGCTGGCCGCGTCGAt                                                                                                                                      <  1:4969075/150‑1 (MQ=255)

                agcaCGGATCAGTGACAACGCCGTCAGCCCCGGCACGCCATACACCACCGCTTCGGTGCCATGGCGCTGGATTACCCGCTCCAGCAGCAAGTCGAAATCACCGTCGCCGGAGGCCAGCACCACCTGATCGACGCGGCCAGCCGGGTCGac                                                                                                                                      <  2:6387077/150‑2 (MQ=255)

                 gcaCGGATCAGTGACAACGCCGTCAGCCCCGGCACGCCATACACCACCGCTTCGGTGCCATGGCGCTGGATTACCCGCTCCAGCAGCAAGTCGAAATCACCGTCGCCGGAGGCCAGCACCACCTGGTCGACCCGGCTGGCCGCGTCGATc                                                                                                                                     >  1:4887869/1‑150 (MQ=255)

                  caCGGATCAGTGACAACGCCGTCAGCCCCGGCACGCCATACACCACCGCTTCGGTGCCATGGCGCTGGATTACCCGCTCCAGCAGCAAGTCGAAATCACCGTCGCCGGAGGCCAGCACCACCTGGTCGACCCGGCTGGCCGCGTCGATca                                                                                                                                    >  2:5963165/1‑150 (MQ=255)

                  caCGGATCAGTGACAACGCCGTCAGCCCCGGCACGCCATACACCACCGCTTCGGTGCCATGGCGCTGGATTACCCGCTCCAGCAGCAAGTCGAAATCACCGTCGCCGGAGGCCAGCACCACCTGATCGACGCGGCCAGCCGGGTCGACca                                                                                                                                    <  2:2238649/150‑1 (MQ=255)

                   aCGGATCAGTGACAACGCCGTCAGCCCCGGCACGCCATACACCACCGCTTCGGTGCCATGGCGCTGGATTACCCGCTCCAGCAGCAAGTCGAAATCACCGTCGCCGGAGGCCAGCACCACCTGGTCGACCCGGCTGGCCGCGTCGATcac                                                                                                                                   <  1:7344906/150‑1 (MQ=255)

                    cGGATCAGTGACAACGCCGTCAGCCCCGGCACGCCATACACCACCGCTTCGGTGCCATGGCGCTGGATTACCCGCTCCAGCAGCAAGTCGAAATCACCGTCGCCGGAGGCCAGCACCACCTGGTCGACCCGGCTGGCCGCGTCGATcaca                                                                                                                                  <  1:5503473/150‑1 (MQ=255)

                       aTCAGTGACAACGCCGTCAGCCCCGGCACGCCATACACCACCGCTTCGGTGCCATGGCGCTGGATTACCCGCTCCAGCAGCAAGTCGAAATCACCGTCGCCGGAGGCCAGCACCACCTGGTCGACCCGGCTGGCCGCGTCGATCACATcc                                                                                                                               <  2:533590/150‑1 (MQ=255)

                        tCAGTGACAACGCCGTCAGCCCCGGCACGCCATACACCACCGCTTCGGTGCCATGGCGCTGGATTACCCGCTCCAGCAGCAAGTCGAAATCACCGTCGCCGGAGGCCAGCACCACCTGGTCGACCCGGCTGGCCGCGTCGa                                                                                                                                       <  2:5046207/141‑1 (MQ=255)

                         cAGTGACAACGCCGTCAGCCCCGGCACGCCATACACCACCGCTTCGGTGCCATGGCGCTGGATTACCCGCTCCAGCAGCAAGTCGAAATCACCGTCGCCGGAGGCCAGCACCACCTGGTCGACCCGGCTGGCCGCGTCGATCACATCCAg                                                                                                                             >  1:8299712/1‑150 (MQ=255)

                          ttgcaCAACGCCGTCAGCCCCGGCACGCCATACACCACCGCTTCGGTGCCATGGCGCTGGATTACCCGCTCCAGCAGCAAGTCGAAATCACCGTCGCCGGAGGCCAGCACCACCTGGTCGACCCGGCTGGCCGCGTCGATCACATCCAgg                                                                                                                            <  2:3302533/146‑1 (MQ=255)

                          aGTGACAACGCCGTCAGCCCCGGCACGCCATACACCACCGCTTCGGTGCCATGGCGCTGGATTACCCGCTCCAGCAGCAAGTCGAAATCACCGTCGCCGGAGGCCAGCACCACCTGGTCGACCCGGCTGGCCGCGTCGATCACATCCAgg                                                                                                                            <  1:3773364/150‑1 (MQ=255)

                           gTGACAACGCCGTCAGCCCCGGCACGCCATACACCACCGCTTCGGTGCCATGGCGCTGGATTACCCGCTCCAGCAGCAAGTCGAAATCACCGTCGCCGGAGGCCAGCACCACCTGGTCGACCCGGCTGGCCGCGTCGATCACATCCAggg                                                                                                                           >  1:531993/1‑150 (MQ=255)

                           gTGACAACGCCGTCAGCCCCGGCACGCCATACACCACCGCTTCGGTGCCATGGCGCTGGATTACCCGCTCCAGCAGCAAGTCGAAATCACCGTCGCCGGAGGCCAGCACCACCTGGTCGACCCGGCTGGCCGCGTCGATCACATCCAggg                                                                                                                           >  1:5573686/1‑150 (MQ=255)

                           gTGACAACGCCGTCAGCCCCGGCACGCCATACACCACCGCTTCGGTGCCATGGCGCTGGATTACCCGCTCCAGCAGCAAGTCGAAATCACCGTCGCCGGAGGCCAGCACCACCTGGTCGACCCGGCTGGCCGCGTCGATCACATCCAggg                                                                                                                           >  2:8412021/1‑150 (MQ=255)

                            tGACAACGCCGTCAGCCCCGGCACGCCATACACCACCGCTTCGGTGCCATGGCGCTGGATTACCCGCTCCAGCAGCAAGTCGAAATCACCGTCGCCGGAGGCCAGCACCACCTGGTCGACCCGGCTGGCCGCGTCGATCACATCCAGGGt                                                                                                                          >  1:1642878/1‑150 (MQ=255)

                               cAACGCCGTCAGCCCCGGCACGCCATACACCACCGCTTCGGTGCCATGGCGCTGGATTACCCGCTCCAGCAGCAAGTCGAAATCACCGTCGCCGGAGGCCAGCACCACCTGGTCGACCCGGCTGGCCGCGTCGATCACATCCAGGGTAAt                                                                                                                       >  2:6905244/1‑150 (MQ=255)

                                  cGCCGTCAGCCCCGGCACGCCATACACCACCGCTTCGGTGCCATGGCGCTGGATTACCCGCTCCAGCAGCAAGTCGAAATCACCGTCGCCGGAGGCCAGCACCACCTGGTCGACCCGGCTGGCCGCGTCGATCACATCCAGGGTAATGcc                                                                                                                    <  2:4668022/150‑1 (MQ=255)

                                   gCCGTCAGCCCCGGCACGCCATACACCACCGCTTCGGTGCCATGGCGCTGGATTACCCGCTCCAGCAGCAAGTCGAAATCACCGTCGCCGGAGGCCAGCACCACCTGGTCGACCCGGCTGGCCGCGTCGATCACATCCAGGGTAATGccc                                                                                                                   <  1:7173397/150‑1 (MQ=255)

                                    ccGTCAGCCCCGGCACGCCATACACCACCGCTTCGGTGCCATGGCGCTGGATTACCCGCTCCAGCAGCAAGTCGAAATCACCGTCGCCGGAGGCCAGCACCACCTGGTCGACCCGGCTGGCCGCGTCGATCACATCCAGGGTAATGCCca                                                                                                                  <  2:42167/150‑1 (MQ=255)

                                           ccccGGCACGCCATACACCACCGCTTCGGTGCCATGGCGCTGGATTACCCGCTCCAGCAGCAAGTCGAAATCACCGTCGCCGGAGGCCAGCACCACCTGGTCGACCCGGCTGGCCGCGTCGATCACATCCAGGGTAATGCCCACATCCCa                                                                                                           <  1:2581636/150‑1 (MQ=255)

                                           ccccGGCACGCCATACACCACCGCTTCGGTGCCATGGCGCTGGATTACCCGCTCCAGCAGCAAGTCGAAATCACCGTCGCCGGAGGCCAGCACCACCTGGTCGACCCGGCTGGCCGCGTCGATCACATCCAGGGTAATGCCCACATCCCa                                                                                                           <  1:7626849/150‑1 (MQ=255)

                                           ccccGGCACGCCATACACCACCGCTTCGGTGCCATGGCGCTGGATTACCCGCTCCAGCAGCAAGTCGAAATCACCGTCGCCGGAGGCCAGCACCACCTGGTCGACCCGGCTGGCCGCGTCGATCACATCCAGGGTAATGCCCACATCCCa                                                                                                           <  1:8745774/150‑1 (MQ=255)

                                           ccccGGCACGCCATACACCACCGCTTCGGTGCCATGGCGCTGGATTACCCGCTCCAGCAGCAAGTCGAAATCACCGTCGCCGGAGGCCAGCACCACCTGGTCGACCCGGCTGGCCGCGTCGATCACATCCAGGGTAATGCCCACATCCCa                                                                                                           <  2:5680557/150‑1 (MQ=255)

                                           ccccGGCACGCCATACACCACCGCTTCGGTGCCATGGCGCTGGATTACCCGCTCCAGCAGCAAGTCGAAATCACCGTCGCCGGAGGCCAGCACCACCTGGTCGACCCGGCTGGCCGCGTCGATCACATCCAGGGTAATGCCCACATCCCa                                                                                                           <  1:4834403/150‑1 (MQ=255)

                                             ccGGCACGCCATACACCACCGCTTCGGTGCCATGGCGCTGGATTACCCGCTCCAGCAGCAAGTCGAAATCACCGTCGCCGGAGGCCAGCACCACCTGGTCGACCCGGCTGGCCGCGTCGATCACATCCAGGGTAATGCCCACATCCCAGt                                                                                                         >  2:7811979/1‑150 (MQ=255)

                                             ccGGCACGCCATACACCACCGCTTCGGTGCCATGGCGCTGGATTACCCGCTCCAGCAGCAAGTCGAAATCACCGTCGCCGGAGGCCAGCACCACCTGGTCGACCCGGCTGGCCGCGTCGATCACATCCAGGGTAATGCCCACATCCCAGt                                                                                                         >  2:3371313/1‑150 (MQ=255)

                                                gCACGCCATACACCACCGCTTCGGTGCCATGGCGCTGGATTACCCGCTCCAGCAGCAAGTCGAAATCACCGTCGCCGGAGGCCAGCACCACCTGGTCGACCCGGCTGGCCGCGTCGATCACATCCAGGGTAATGCCCACATCCCAGTCAc                                                                                                      >  2:6713515/1‑150 (MQ=255)

                                                  aCGCCATACACCACCGCTTCGGTGCCATGGCGCTGGATTACCCGCTCCAGCAGCAAGTCGAAATCACCGTCGCCGGGGGCCAGCACCACCTGGTCGACCCGGCTGGCCGCGTCGATCACATCCAGGGTAATGCCCACATCCCAGTCACCt                                                                                                    <  1:6500261/150‑1 (MQ=255)

                                                  aCGCCATACACCACCGCTTCGGTGCCATGGCGCTGGATTACCCGCTCCAGCAGCAAGTCGAAATCACCGTCGCCGGAGGCCAGCACCACCTGGTCGACCCGGCTGGCCGCGTCGATCACATCCAGGGTAATGCCCACATCCCAGTCACCt                                                                                                    <  1:1942370/150‑1 (MQ=255)

                                                  aCGCCATACACCACCGCTTCGGTGCCATGGCGCTGGATTACCCGCTCCAGCAGCAAGTCGAAATCACCGTCGCCGGAGGCCAGCACCACCTGGTCGACCCGGCTGGCCGCGTCGATCACATCCAGGGTAATGCCCACATCCCAGTCACCt                                                                                                    >  2:7284044/1‑150 (MQ=255)

                                                     ccATACACCACCGCTTCGGTGCCATGGCGCTGGATTACCCGCTCCAGCAGCAAGTCGAAATCACCGTCGCCGGAGGCCAGCACCACCTGGTCGACCCGGCTGGCCGCGTCGATCACATCCAGGGTAATGCCCACATCCCAGTCACCTTTg                                                                                                 <  1:4182201/150‑1 (MQ=255)

                                                     ccATACACCACCGCTTCGGTGCCATGGCGCTGGATTACCCGCTCCAGCAGCAAGTCGAAATCACCGTCGCCGGAGGCCAGCACCACCTGGTCGACCCGGCTGGCCGCGTCGATCACATCCAGGGTAATGCCCACATCCCAGTCACCTTTg                                                                                                 <  1:1506826/150‑1 (MQ=255)

                                                     ccATACACCACCGCTTCGGTGCCATGGCGCTGGATTACCCGCTCCAGCAGCAAGTCGAAATCACCGTCGCCGGAGGCCAGCACCACCTGGTCGACCCGGCTGGCCGCGTCGATCACATCCAGGGTAATGCCCACATCCCAGTCACCTTTg                                                                                                 >  1:6983480/1‑150 (MQ=255)

                                                      cATACACCACCGCTTCGGTGCCATGGCGCTGGATTACCCGCTCCAGCAGCAAGTCGAAATCACCGTCGCCGGAGGCCAGCACCACCTGGTCGACCCGGCTGGCCGCGTCGATCACATCCAGGGTAATGCCCACATCCCAGTCACCTTTgg                                                                                                >  1:8325840/1‑150 (MQ=255)

                                                          caccacCGCTTCGGTGCCATGGCGCTGGATTACCCGCTCCAGCAGCAAGTCGAAATCACCGTCGCCGGAGGCCAGCACCACCTGGTCGACCCGGCTGGCCGCGTCGATCACATCCAGGGTAATGCCCACATCCCAGTCACCTTTGGCCGa                                                                                            <  1:6139401/150‑1 (MQ=255)

                                                          caccacCGCTTCGGTGCCATGGCGCTGGATTACCCGCTCCAGCAGCAAGTCGAAATCACCGTCGCCGGAGGCCAGCACCACCTGGTCGACCCGGCTGGCCGCGTCGATCACATCCAGGGTAATGCCCACATCCCAGTCACCTTTGGCCGa                                                                                            <  2:1957704/150‑1 (MQ=255)

                                                           accaccGCTTCGGTGCCATGGCGCTGGATTACCCGCTCCAGCAGCAAGTCGAAATCACCGTCGCCGGAGGCCAGCACCACCTGGTCGACCCGGCTGGCCGCGTCGATCACATCCAGGGTAATGCCCACATCCCAGTCACCTTTGGCCGAg                                                                                           >  1:3452496/1‑150 (MQ=255)

                                                           accaccGCTTCGGTGCCATGGCGCTGGATTACCCGCTCCAGCAGCAAGTCGAAATCACCGTCGCCGGAGGCCAGCACCACCTGGTCGACCCGGCTGGCCGCGTCGATCACATCCAGGGTAATGCCCACATCCCAGTCACCTTTGGCCGAg                                                                                           >  2:3165779/1‑150 (MQ=255)

                                                             caccGCTTCGGTGCCATGGCGCTGGATTACCCGCTCCAGCAGCAAGTCGAAATCACCGTCGCCGGAGGCCAGCACCACCTGGTCGACCCGGCTGGCCGCGTCGATCACATCCAGGGTAATGCCCACATCCCa                                                                                                           >  2:5562785/1‑132 (MQ=255)

                                                             caccGCTTCGGTGCCATGGCGCTGGATTACCCGCTCCAGCAGCAAGTCGAAATCACCGTCGCCGGAGGCCAGCACCACCTGGTCGACCCGGCTGGCCGCGTCGATCACATCCAGGGTAATGCCCACATCCCAGTCACCTTTGGCCGAGcc                                                                                         >  2:1180643/1‑150 (MQ=255)

                                                               ccGCTTCGGTGCCATGGCGCTGGATTACCCGCTCCAGCAGCAAGTCGAAATCACCGTCGCCGGAGGCCAGCACCACCTGGTCGACCCGGCTGGCCGCGTCGATCACATCCAGGGTAATGCCCACATCCCAGTCACCTTTGGCCGAGCCAt                                                                                       <  1:3036926/150‑1 (MQ=255)

                                                                     cGGTGCCATGGCGCTGGATTACCCGCTCCAGCAGCAAGTCGAAATCACCGTCGCCGGAGGCCAGCACCACCTGGTCGACCCGGCTGGCCGCGTCGATCACATCCAGGGTAATGCCCACATCCCAGTCACCTTTGGCCGAGCCATCACTgc                                                                                 <  1:314639/150‑1 (MQ=255)

                                                                       gTGCCATGGCGCTGGATTACCCGCTCCAGCAGCAAGTCGAAATCACCGTCGCCGGAGGCCAGCACCACCTGGTCGACCCGGCTGGCCGCGTCGATCACATCCAGGGTAATGCCCACATCCCAGTCACCTTTGGCCGAGCCATCACTgcgc                                                                               >  2:4267513/1‑150 (MQ=255)

                                                                              ggCGCTGGATTACCCGCTCCAGCAGCAAGTCGAAATCACCGTCGCCGGAGGCCAGCACCACCTGGTCGACCCGGCTGGCCGCGTCGATCACATCCAGGGTAATGCCCACATCCCAGTCACCTTTGGCCGAGCCATCACTGCGCTGGATGa                                                                        >  2:8013154/1‑150 (MQ=255)

                                                                                  cTGGATTACCCGCTCCAGCAGCAAGTCGAAATCACCGTCGCCGGAGGCCAGCACCACCTGGTCGACCCGGCTGGCCGCGTCGATCACATCCAGGGTAATGCCCACATCCCAGTCACCTTTGGCCGAGCCATCACTGCGCTGGATGAACgg                                                                    >  1:2028807/1‑150 (MQ=255)

                                                                                    ggATTACCCGCTCCAGCAGCAAGTCGAAATCACCGTCGCCGGAGGCCAGCACCACCTGGTCGACCCGGCTGGCCGCGTCGATCACATCCAGGGTAATGCCCACATCCCAGTCACCTTTGGCCGAGCCATCACTGCGCTGGTTGAACGGtt                                                                  <  2:3884325/150‑1 (MQ=255)

                                                                                    ggATTACCCGCTCCAGCAGCAAGTCGAAATCACCGTCGCCGGAGGCCAGCACCACCTGGTCGACCCGGCTGGCCGCGTCGATCACATCCAGGGTAATGCCCACATCCCAGTCACCTTTGGCCGAGCCATCACTGCGCTGGATGAACGGtt                                                                  <  1:4251077/150‑1 (MQ=255)

                                                                                    ggATTACCCGCTCCAGCAGCAAGTCGAAATCACCGTCGCCGGAGGCCAGCACCACCTGGTCGACCCGGCTGGCCGCGTCGATCACATCCAGGGTAATGCCCACATCCCAGTCACCTTTGGCCGAGCCATCACTGCGCTGGATGAACGGtt                                                                  <  2:8572430/150‑1 (MQ=255)

                                                                                       ttACCCGCTCCAGCAGCAAGTCGAAATCACCGTCGCCGGAGGCCAGCACCACCTGGTCGACCCGGCTGGCCGCGTCGATCACATCCAGGGTAATGCCCACATCCCAGTCACCTTTGGCCGAGCCATCACTGCGCTGGATGAACGGTTTGa                                                               <  1:4064531/150‑1 (MQ=255)

                                                                                        tACCCGCTCCAGCAGCAAGTCGAAATCACCGTCGCCGGAGGCCAGCACCACCTGGTCGACCCGGCTGGCCGCGTCGATCACATCCAGGGTAATGCCCACATCCCAGTCACCTTTGGCCGAGCCATCACTGCGCTGGATGAACGGTTTGAg                                                              <  1:1180643/150‑1 (MQ=255)

                                                                                        tACCCGCTCCAGCAGCAAGTCGAAATCACCGTCGCCGGAGGCCAGCACCACCTGGTCGACCCGGCTGGCCGCGTCGATCACATCCAGGGTAATGCCCACATCCCAGTCACCTTTGGCCGAGCCATCACTGCGCTGGATGAACGGTTTGAg                                                              <  1:637751/150‑1 (MQ=255)

                                                                                        tACCCGCTCCAGCAGCAAGTCGAAATCACCGTCGCCGGAGGCCAGCACCACCTGGTCGACCCGGCTGGCCGCGTCGATCACATCCAGGGTAATGCCCACATCCCAGTCACCTTTGGCCGAGCCATCACTGCGCTGGATGAACGGTTTGAg                                                              <  1:352714/150‑1 (MQ=255)

                                                                                          cccGCTCCAGCAGCAAGTCGAAATCACCGTCGCCGGAGGCCAGCACCACCTGGTCGACCCGGCTGGCCGCGTCGATCACATCCAGGGTAATGCCCACATCCCAGTCACCTTTGGCCGAGCCATCACTGCGCTGGATGAACGGTTTGAGgc                                                            >  2:8173137/1‑150 (MQ=255)

                                                                                           ccGCTCCAGCAGCAAGTCGAAATCACCGTCGCCGGAGGCCAGCACCACCTGGTCGACCCGGCTGGCCGCGTCGATCACATCCAGGGTAATGCCCACATCCCAGTCACCTTTGGCCGAGCCATCACTGCGCTGGATGAACGGTTTGAGgcg                                                           <  1:6039569/150‑1 (MQ=255)

                                                                                            cGCTCCAGCAGCAAGTCGAAATCACCGTCGCCGGAGGCCAGCACCACCTGGTCGACCCGGCTGGCCGCGTCGATCACATCCAGGGTAATGCCCACATCCCAGTCACCTTTGGCCGAGCCATCACTGCGCTGGATGAACGGTTTGAGgcgc                                                          >  1:8591713/1‑150 (MQ=255)

                                                                                                 cagcagCAAGTCGAAATCACCGTCGCCGGAGGCCAGCACCACCTGGTCGACCCGGCTGGCCGCGTCGATCACATCCAGGGTAATGCCCACATCCCAGTCACCTTTGGCCGAGCCATCACTGCGCTGGATGAACGGTTTGAGGCGCACCTc                                                     <  2:4887869/150‑1 (MQ=255)

                                                                                                    cagcaAGTCGAAATCACCGTCGCCGGAGGCCAGCACCACCTGGTCGACCCGGCTGGCCGCGTCGATCACATCCAGGGTAATGCCCACATCCCAGTCACCTTTGGCCGAGCCATCACTGCGCTGGATGAACGGTTTGAGGCGCACCTCGaa                                                  <  2:1712111/150‑1 (MQ=255)

                                                                                                        aaGTCGAAATCACCGTCGCCGGAGGCCAGCACCACCTGGTCGACCCGGCTGGCCGCGTCGATCACATCCAGGGTAATGCCCACATCCCAGTCACCTTTGGCCGAGCCATCACTGCGCTGGATGAACGGTTTGAGGCGCACCTCGAAACCg                                              >  1:8409401/1‑150 (MQ=255)

                                                                                                        aaGTCGAAATCACCGTCGCCGGAGGCCAGCACCACCTGGTCGACCCGGCTGGCCGCGTCGATCACATCCAGGGTAATGCCCACATCCCAGTCACCTTTGGCCGAGCCATCACTGCGCTGGATGAACGGTTTGAGGCGCACCTCGAAACCg                                              >  1:5095649/1‑150 (MQ=255)

                                                                                                        aaGTCGAAATCACCGTCGCCGGAGGCCAGCACCACCTGGTCGACCCGGCTGGCCGCGTCGATCACATCCAGGGTAATGCCCACATCCCAGTCACCTTTGGCCGAGCCATCACTGCGCTGGATGAACGGTTTGAGGCGCACCTCGAAACCg                                              >  2:6133780/1‑150 (MQ=255)

                                                                                                         aGTCGAAATCACCGTCGCCGGAGGCCAGCACCACCTGGTCGACCCGGCTGGCCGCGTCGATCACATCCAGGGTAATGCCCACATCCCAGTCACCTTTGGCCGAGCCATCACTGCGCTGGATGAACGGTTTGAGGCGCACCTCGAAACCGa                                             <  1:776367/150‑1 (MQ=255)

                                                                                                         aGTCGAAATCACCGTCGCCGGAGGCCAGCACCACCTGGTCGACCCGGCTGGCCGCGTCGATCACATCCAGGGTAATGCCCACATCCCAGTCACCTTTGGCCGAGCCATCACTGCGCTGGATGAACGGTTTGAGGCGCACCTCGAAACCGa                                             <  2:7611891/150‑1 (MQ=255)

                                                                                                          gTCGAAATCACCGTCGCCGGAGGCCAGCACCACCTGGTCGACCCGGCTGGCCGCGTCGATCACATCCAGGGTAATGCCCACATCCCAGTCACCTTTGGCCGAGCCATCACTGCGCTGGATGAACGGTTTGAGGCGCACCTCGAAACCGAg                                            <  2:6990854/150‑1 (MQ=255)

                                                                                                           tctactacACCGTCGCCGGAGGCCAGCACCACCTGGTCGACCCGGCTGGCCGCGTCGATCACATCCAGGGTAATGCCCACATCCCAGTCACCTTTGGCCGAGCCATCACTGCGCTGGATGAACGGTTTGAGGCGCACCTCGAAACCGAgg                                           >  2:457405/8‑150 (MQ=255)

                                                                                                              aaaTCACCGTCGCCGGAGGCCAGCACCACCTGGTCGACCCGGCTGGCCGCGTCGATCACATCCAGGGTAATGCCCACATCCCAGTCACCTTTGGCCGAGCCATCACTGCGCTGGATGAACGGTTTGAGGCGCACCTCGAAACCGAGGTTg                                        >  1:1124283/1‑150 (MQ=255)

                                                                                                               aaTCACCGTCGCCGGAGGCCAGCACCACCTGGTCGACCCGGCTGGCCGCGTCGATCACATCCAGGGTAATGCCCACATCCCAGTCACCTTTGGCCGAGCCATCACTGCGCTGGATGAACGGTTTGAGGCGCACCTCGAAACCGAGGTTgc                                       >  1:7991403/1‑150 (MQ=255)

                                                                                                                aTCACCGTCGCCGGAGGCCAGCACCACCTGGTCGACCCGGCTGGCCGCGTCGATCACATCCAGGGTAATGCCCACATCCCAGTCACCTTTGGCCGAGCCATCACTGCGCTGGATGAACGGTTTGAGGCGCACCTCGAAACCGAGGTTgcg                                      >  2:2969273/1‑150 (MQ=255)

                                                                                                                    ccGTCGCCGGAGGCCAGCACCACCTGGTCGACCCGGCTGGCCGCGTCGATCACATCCAGGGTAATGCCCACATCCCAGTCACCTTTGGCCGAGCCATCACTGCGCTGGATGAACGGTTTGAGGCGCACCTCGAAACCGAGGTTGCGCAAg                                  >  2:6801073/1‑150 (MQ=255)

                                                                                                                         gCCGGAGGCCAGCACCACCTGGTCGACCCGGCTGGCCGCGTCGATCACATCCAGGGTAATGCCCACATCCCAGTCACCTTTGGCCGAGCCATCACTGCGCTGGATGAACGGTTTGAGGCGCACCTCGAAACCGAGGTTGCGCAAGATctg                             >  2:149791/1‑150 (MQ=255)

                                                                                                                          ccGGAGGCCAGCACCACCTGGTCGACCCGGCTGGCCGCGTCGATCACATCCAGGGTAATGCCCACATCCCAGTCACCTTTGGCCGAGCCATCACTGCGCTGGATGAACGGTTTGAGGCGCACCTCGAAACCGAGGTTGCGCAAGATctgc                            >  2:2842898/1‑150 (MQ=255)

                                                                                                                             gAGGCCAGCACCACCTGGTCGACCCGGCTGGCCGCGTCGATCACATCCAGGGTAATGCCCACATCCCAGTCACCTTTGGCCGAGCCATCACTGCGCTGGATGAACGGTTTGAGGCGCACCTCGAAACCGAGGTTGCGCAAGATCTGCTgg                         >  1:8804691/1‑150 (MQ=255)

                                                                                                                              aGGCCAGCACCACCTGGTCGACCCGGCTGGCCGCGTCGATCACATCCAGGGTAATGCCCACATCCCAGTCACCTTTGGCCGAGCCATCACTGCGCTGGATGAACGGTTTGAGGCGCACCTCGAAACCGAGGTTGCGCAAGATCTGCTGGa                        <  2:7336634/150‑1 (MQ=255)

                                                                                                                                gCCAGCACCACCTGGTCGACCCGGCTGGCCGCGTCGATCACATCCAGGGTAATGCCCACATCCCAGTCACCTTTGGCCGAGCCATCACTGCGCTGGATGAACGGTTTGAGGCGCACCTCGAAACCGAGGTTGCGCAAGATCTGCTGGAAc                      <  2:2451938/150‑1 (MQ=255)

                                                                                                                                  cAGCACCACCTGGTCGACCCGGCTGGCCGCGTCGATCACATCCAGGGTAATGCCCACATCCCAGTCACCTTTGGCCGAGCCATCACTGCGCTGGATGAACGGTTTGAGGCGCACCTCGAAACCGAGGTTGCGCAAGATCTGCTGGAActg                    >  1:2326936/1‑150 (MQ=255)

                                                                                                                                  cAGCACCACCTGGTCGACCCGGCTGGCCGCGTCGATCACATCCAGGGTAATGCCCACATCCCAGTCACCTTTGGCCGAGCCATCACTGCGCTGGATGAACGGTTTGAGGCGCACCTCGAAACCGAGGTTGCGCAAGATCTGCTGGAActg                    <  1:7706124/150‑1 (MQ=255)

                                                                                                                                    gCACCACCTGGTCGACCCGGCTGGCCGCGTCGATCACATCCAGGGTAATGCCCACATCCCAGTCACCTTTGGCCGAGCCATCACTGCGCTGGATGAACGGTTTGAGGCGCACCTCGAAACCGAGGTTGCGCAAGATCTGCTGGAActgct                  <  1:2887111/150‑1 (MQ=255)

                                                                                                                                          ccTGGTCGACCCGGCTGGCCGCGTCGATCACATCCAGGGTAATGCCCACATCCCAGTCACCTTTGGCCGAGCCATCACTGCGCTGGATGAACGGTTTGAGGCGCACCTCGAAACCGAGGTTGCGCAAGATCTGCTGGAACTGCTgttgtt            <  2:8299770/150‑1 (MQ=255)

                                                                                                                                              gTCGACCCGGCTGGCCGCGTCGATCACATCCAGGGTCATGCCCACATCCCAGTCACCTTTGGCCGAGCCATCACTGCGCTGGATGAACGGTTTGAGGCGCACCTCGAAACCGAGGTTGCGCAAGATCTGCTGGAACTGCTGTTGTTTGCt        <  2:4745352/150‑1 (MQ=255)

                                                                                                                                              gTCGACCCGGCTGGCCGCGTCGATCACATCCAGGGTAATGCCCACATCCCAGTCACCTTTGGCCGAGCCATCACTGCGCTGGATGAACGGTTTGAGGCGCACCTCGAAACCGAGGTTGCGCAAGATCTGCTGGAACTGCTGTTGTTTGCt        <  1:6807355/150‑1 (MQ=255)

                                                                                                                                              gTCGACCCGGCTGGCCGCGTCGATCACATCCAGGGTAATGCCCACATCCCAGTCACCTTTGGCCGAGCCATCACTGCGCTGGATGAACGGTTTGAGGCGCACCTCGAAACCGAGGTTGCGCAAGATCTGCTGGAACTGCTGTTGTTTGCt        <  2:8072673/150‑1 (MQ=255)

                                                                                                                                              gTCGACCCGGCTGGCCGCGTCGATCACATCCAGGGTAATGCCCACATCCCAGTCACCTTTGGCCGAGCCATCACTGCGCTGGATGAACGGTTTGAGGCGCACCTCGAAACCGAGGTTGCGCAAGATCTGCTGGAACTGCTGTTGTTTGCt        >  1:4548104/1‑150 (MQ=255)

                                                                                                                                              atCGACGCGGCCAGCCGGGTCGACCACATCCAGGGTAATGCCCACATCCCAGTCACCTTTGGCCGAGCCATCACTGCGCTGGATGAACGGTTTGAGGCGCACCTCGAAACCGAGGTTGCGCAAGATCTGCTGGAACTGCTGTTGTTTGCt        >  1:8039471/2‑150 (MQ=255)

                                                                                                                                               tCGACGCGGCCAGCCGGGTCGATCACATCCAGGGTAATGCCCACATCCCAGTCACCTTTGGCCGAGCCATCACTGCGCTGGATGAACGGTTTGAGGCGCACCTCGAAACCGAGGTTGCGCAAGATCTGCTGGAACTGCTGTTGTTTGCTg       >  2:4337274/1‑150 (MQ=255)

                                                                                                                                                   cccGGCTGGCCGCGTCGATCACATCCAGGGTAATGCCCACATCCCAGTCACCTTTGGCCGAGCCATCACTGCGCTGGATGAACGGTTTGAGGCGCACCTCGAAACCGAGGTTGCGCAAGATCTGCTGGAACTGCTGTTGTTTGCTGTCGc   >  1:8491462/1‑150 (MQ=255)

                                                                                                                                                    ccGGCTGGCCGCGTCGATCACATCCAGGGTAATGCCCACATCCCAGTCACCTTTGGCCGAGCCATCACTGCGCTGGATGAACGGTTTGAGGCGCACCTCGAAACCGAGGTTGCGCAAGATCTGCTGGAACTGCTGTTGTTTGCTGTCGcc  >  1:1805157/1‑150 (MQ=255)

GCACGTAGCGGGTGGCAGCACGGATCAGTGACAACGCCGTCAGCCCCGGCACGCCATACACCACCGCTTCGGTGCCATGGCGCTGGATTACCCGCTCCAGCAGCAAGTCGAAATCACCGTCGCCGGAGGCCAGCACCACCTGGTCGACCCGGCTGGCCGCGTCGATCACATCCAGGGTAATGCCCACATCCCAGTCACCTTTGGCCGAGCCATCACTGCGCTGGATGAACGGTTTGAGGCGCACCTCGAAACCGAGGTTGCGCAAGATCTGCTGGAACTGCTGTTGTTTGCTGTCGCC  >  NC_002947/5434868‑5435165

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