BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,176,096	T→A	5.4%	L101Q (CTG→CAG)	PP_2786 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,176,096	0	T	A	5.4%	451.6 / 7.0	138	L101Q (CTG→CAG)	PP_2786	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (67/63); new base A (2/6); total (69/69)
Fisher's exact test for biased strand distribution p-value = 2.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 139 total aligned reads displayed.

CCTGATCGACAAGCAAGGTAACAGCCGCATCCGCGAGTTCAGCTACCTGCAGAAGGATTACCCGGACAGCGACAAATTCACCATGTATTTCTC‑C‑GCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCGGGGCGTTCATGGGTAGCGACTACTCCTATGCCGACCTCGACAAGATTCGCGCCGAAGACTATCGCCAAAC  >  NC_002947/3175947‑3176243

cctgatcGACAAGCAATGTAACAGCCGCATCCGCGAGTTCAGCTACCTGCAGAAGGATTACCCGGACAGCGACAAATTCACCATGTATTTCTC‑C‑GCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCt                                                                                                                                                     >  2:3847082/1‑150 (MQ=255)

 ctgatcGACAAGCAAGGTAACAGCCGCATCCGCGAGTTCAGCTACCTGCAGAAGGATTACCCGGACAGCGACAAATTCACCATGTATTTCTC‑C‑GCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTg                                                                                                                                                    <  2:3237454/150‑1 (MQ=255)

  ttatcGACAAGCAAGGTAACAGCCTCATCCGCGAGTTCAGCTACCTGCAGAAGGATTACCCGGACAGCGACAAATTCACCATGTATTTCTC‑C‑GCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTgg                                                                                                                                                   <  1:2799041/148‑1 (MQ=255)

  tgatcGACAAGCAAGGTAACAGCCGCATCCGCGAGTTCAGCTACCTGCAGAAGGATTACCCGGACAGCGACAAATTCACCATGTATTTCTC‑C‑GCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTgg                                                                                                                                                   <  1:3058100/150‑1 (MQ=255)

   gatcGACAAGCAAGGTAACAGCCGCATCCGCGAGTTCAGCTACCTGCAGAAGGATTACCCGGACAGCGACAAATTCACCATGTATTTCTC‑C‑GCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTggg                                                                                                                                                  >  1:2309173/1‑150 (MQ=255)

       gACAATCAAGGTAACAGCCGCAGCCGCGAGTTCAGCTACCTGCAGACGGATTACCCGGACAGCGACAAATTCACCAGGTATTTCTC‑C‑GAGCCAAGTGATGTGCGTGACGTCGCCTTGCACATCGAAAACCCGCGCGAGGCCCAGGGCCTc                                                                                                                                              <  2:4751452/150‑1 (MQ=255)

        aCAAGCAAGGTAACAGCCGCATCCGCGAGTTCAGCTACCTGCAGAAGGATTACCCGGACAGCGACAAATTCACCATGTATTTCTC‑C‑GCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCg                                                                                                                                             >  2:1331098/1‑150 (MQ=255)

          aagcaagGTAACAGCCGCATCCGCGAGTTCAGCTACCTGCAGACGGATTACCCGGACAGCGACAAATTCACCCTGTATGTCTC‑C‑GCGCCAAGGGATGTGCGTGACGTCGCCTTCCACATCGACAACCCGCGCGCGGCCCTGGGCCTCGaa                                                                                                                                           <  2:5176940/150‑1 (MQ=255)

          aagcaagGTAACAGCCGCATCCGCGAGTTCAGCTACCTGCAGAAGGATTACCCGGACAGCGACAAATTCACCATGTATTTCTC‑C‑GCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGaa                                                                                                                                           <  1:3177481/150‑1 (MQ=255)

           agcaagGTAACAGCCGCATCCGCGAGTTCAGCTACCTGCAGAAGGATTACCCGGACAGCGACAAATTCACCATGTATTTCTC‑C‑GCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCAGGGCCTCgcgc                                                                                                                                          <  1:2159421/150‑4 (MQ=255)

           agcaagGTAACAGCCGCATCCGCGAGTTCAGCTACCTGCAGAAGGATTACCCGGACAGCGACAAATTCACCATGTATTTCTC‑C‑GCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCAGGGCCTCgcgc                                                                                                                                          <  1:4010386/150‑4 (MQ=255)

             caagGTAACAGCCGCATCCGCGAGTTCAGCTACCTGCAGAAGGATTACCCGGACAGCGACAAATTCACCATGTATTTCTC‑C‑GCGCCAAGTGTTGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCAGGGCCTCgcgcgg                                                                                                                                        <  2:1035548/150‑6 (MQ=255)

              aagGTAACAGCCGCATCCGCGAGTTCAGCTACCTGCAGAAGGATTACCCGGACAGCGACAAATTCACCATGTATTTCTCTCTGCGCCAAGTGATGCGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGa                                                                                                                                         <  1:2122847/150‑1 (MQ=255)

                ggTAACAGCCGCATCCGCGAGTTCAGCTACCTGCAGAAGGATTACCCGGACAGCGACAAATTCACCATGTATTTCTC‑C‑GCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGc                                                                                                                                     <  1:1847194/150‑1 (MQ=255)

                ggTAACAGCCGCATCCGCGAGTTCAGCTACCTGCAGAAGGATTACCCGGACAGCGACAAATTCACCATGTATTTCTC‑C‑GCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGc                                                                                                                                     <  2:309469/150‑1 (MQ=255)

                 gTAACAGCCGCATCCGCGAGTTCAGCTACCTGCAGAAGGATTACCCGGACAGCGACAAATTCACCATGTATTTCTC‑C‑GCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGcc                                                                                                                                    >  1:2713902/1‑150 (MQ=255)

                     cAGCCGCATCCGCGAGTTCAGCTACCTGCAGAAGGATTACCCGGACAGCGACAAATTCACCATGTATTTCTC‑C‑GCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTg                                                                                                                                <  1:1093034/150‑1 (MQ=255)

                      aGCCGCATCCGCGAGTTCAGCTACCTGCAGAAGGATTACCCGGACAGCGACAAATTCACCATGTATTTCTC‑C‑GCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATAGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTgg                                                                                                                               >  1:2273238/1‑150 (MQ=255)

                       gCCGCATCCGCGAGTTCAGCTACCTGCAGAAGGATTACCCGGACAGCGACAAATTCACCATGTATTTCTC‑C‑GCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGa                                                                                                                              <  1:144858/150‑1 (MQ=255)

                        ccGCATCCGCGAGTTCAGCTACCTGCAGAAGGATTACCCGGACAGCGACATATTCACCATGTATTTCTC‑C‑GCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGAt                                                                                                                             >  1:1787789/1‑150 (MQ=255)

                        ccGCATCCGCGAGTTCAGCTACCTGCAGAAGGATTACCCGGACAGCGACAAATTCACCATGTATTTCTC‑C‑GCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGAt                                                                                                                             <  1:5050219/150‑1 (MQ=255)

                        ccGCATCCGCGAGTTCAGCTACCTGCAGAAGGATTACCCGGACAGCGACAAATTCACCATGTATTTCTC‑C‑GCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGAt                                                                                                                             <  1:2483749/150‑1 (MQ=255)

                        ccGCATCCGCGAGTTCAGCTACCTGCAGAAGGATTACCCGGACAGCGACAAATTCACCATGTATTTCTC‑C‑GCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCACTGGAt                                                                                                                             <  2:2758945/150‑1 (MQ=255)

                              ccGCGAGTTCAGCTACCTGCAGAAGGATTACCCGGACAGCGACAAATTCACCATGTATTTCTC‑C‑GCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTAttt                                                                                                                       >  2:280411/1‑150 (MQ=255)

                                gcgAGTTCAGCTACCTGCAGAAGGATTACCCGGACAGCGACAAATTCACCATGTATTTCTC‑C‑GCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGc                                                                                                                     >  1:4974441/1‑150 (MQ=255)

                                gcgAGTTCAGCTACCTGCAGAAGGATTACCCGGACAGCGACAAATTCACCATGTATTTCTC‑C‑GCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGc                                                                                                                     <  1:4863853/150‑1 (MQ=255)

                                gcgAGTTCAGCTACCTGCAGAAGGATTACCCGGACAGCGACAAATTCACCATGTATTTCTC‑C‑GCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGc                                                                                                                     >  1:1638048/1‑150 (MQ=255)

                                 cgAGTTCAGCTACCTGCAGAAGGATTACCCGGACAGCGACAAATTCACCATGTATTTCTC‑C‑GCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGcc                                                                                                                    >  2:1454872/1‑150 (MQ=255)

                                 cgAGTTCAGCTACCTGCAGAAGGATTACCCGGACAGCGACAAATTCACCATGTATTTCTC‑C‑GCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGcc                                                                                                                    >  2:604586/1‑150 (MQ=255)

                                 cgAGTTCAGCTACCTGCAGAAGGATTACCCGGACAGCGACAAATTCACCATGTATTTCTC‑C‑GCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGcc                                                                                                                    >  2:344553/1‑150 (MQ=255)

                                 cgAGTTCAGCTACCTGCAGAAGGATTACCCGGACAGCGACAAATTCACCATGTATTTCTC‑C‑GCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGcc                                                                                                                    <  1:1681265/150‑1 (MQ=255)

                                     ttCAGCTACCTGCAGAAGGATTACCCGGACAGCGACAAATTCACCTTGTATTTCTC‑C‑GCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTc                                                                                                                <  1:3271281/150‑1 (MQ=255)

                                        aGCTACCTGCAGAAGGATTACCCGGACAGCGACAAATTCACCATGTATTTCTC‑C‑GCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGc                                                                                                             <  1:3643083/150‑1 (MQ=255)

                                         gCTACCTGCAGAAGGATTACCCGGACAGCGACAAATTCACCATGTATTTCTC‑C‑GCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGcc                                                                                                            <  1:1533055/150‑1 (MQ=255)

                                           tACCTGCAGAAGGATTACCCGGACAGCGACAAATTCACCATGTATTTCTC‑C‑GCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGa                                                                                                          <  1:2950348/150‑1 (MQ=255)

                                           tACCTGCAGAAGGATTACCCGGACAGCGACAAATTCACCATGTATTTCTC‑C‑GCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGa                                                                                                          <  1:3731297/150‑1 (MQ=255)

                                             ccTGCAGAAGGATTACCCGGACAGCGACAAATTCACCATGTATTTCTC‑C‑GCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACa                                                                                                        >  1:3587556/1‑150 (MQ=255)

                                              cTGCAGAAGGATTACCCGGACAGCGACAAATTCACCATGTATTTCTC‑C‑GCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAACCAGTGGATGTATTTGCCGGTCAGCCGACAg                                                                                                       <  2:2103682/150‑1 (MQ=255)

                                               tGCAGAAGGATTACCCGGACAGCGACAAATTCACCATGTATTTCTC‑C‑GCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGa                                                                                                      >  2:2788886/1‑150 (MQ=255)

                                                     aGGATTACCCGGACAGCGACAAATTCACCATGTATTTCTC‑C‑GCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGaa                                                                                                <  1:1480190/150‑1 (MQ=255)

                                                            cccGGACAGCGACAAATTCACCATGTATTTCTC‑C‑GCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTg                                                                                                                                                    <  2:546126/91‑1 (MQ=255)

                                                            cccGGACAGCGACAAATTCACCATGTATTTCTC‑C‑GCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTc                                                                                         <  2:344965/150‑1 (MQ=255)

                                                             ccGGACAGCGACAAATTCACCATGTATTTCTC‑C‑GCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTcc                                                                                        >  2:5062282/1‑150 (MQ=255)

                                                               ggACAGCGACAAATTCACCATGTATTTCTC‑C‑GCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTccac                                                                                      <  2:3451578/150‑1 (MQ=255)

                                                                      gACAAATTCACCATGTATTTCTC‑C‑GCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGAc                                                                               <  1:4511165/150‑1 (MQ=255)

                                                                      gACAAATTCACCATGTATTTCTC‑C‑GCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGAc                                                                               >  2:3281636/1‑150 (MQ=255)

                                                                      gACAAATTCACCATGTATTTCTC‑C‑GCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGAc                                                                               <  1:3408018/150‑1 (MQ=255)

                                                                      gACAAATTCACCATGTATTTCTC‑C‑GCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGAc                                                                               >  2:1986167/1‑150 (MQ=255)

                                                                       aCAAATTCACCATGTATTTCTC‑C‑GCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACa                                                                              >  2:3545588/1‑150 (MQ=255)

                                                                       aCAAATTCACCATGTATTTCTC‑C‑GCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACa                                                                              >  1:930220/1‑150 (MQ=255)

                                                                         aaaTTCACCATGTATTTCTC‑C‑GCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACaaa                                                                            >  1:790531/1‑150 (MQ=255)

                                                                                ccATGTATTTCTC‑C‑GCGCCAAGTGCTGTGCGTGACGTCGCCTTCCACATCGAAAACCCGTGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCgggg                                                                     >  2:3081014/1‑150 (MQ=255)

                                                                                ccATGTATTTCTC‑C‑GCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCTGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCgggg                                                                     <  2:840371/150‑1 (MQ=255)

                                                                                 cATGTATTTCTC‑C‑GCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGAc                                                                                                                                        >  1:2598344/1‑82 (MQ=255)

                                                                                 cATGTATTTCTC‑C‑GCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGAc                                                                                                                                        <  2:2598344/82‑1 (MQ=255)

                                                                                 cATGTATTTCTC‑C‑GCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCGGGGc                                                                    >  2:4113212/1‑150 (MQ=255)

                                                                                  aTGTATTTCTC‑C‑GCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCGGGGCg                                                                   >  1:3392808/1‑150 (MQ=255)

                                                                                   tGTATTTCTC‑C‑GCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCGGGGCGt                                                                  >  2:4649042/1‑150 (MQ=255)

                                                                                      aTTTCTC‑C‑GCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCGGGGCGTTCa                                                               >  1:728117/1‑150 (MQ=255)

                                                                                       tttCTC‑C‑GCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCGGGGCGTTCAt                                                              >  1:3647381/1‑150 (MQ=255)

                                                                                       tttCTC‑C‑GCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCGGGGCGTTCAt                                                              <  2:4855337/150‑1 (MQ=255)

                                                                                                       gTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCGGGGCGTTCATGGGTAGCGactact                                                <  2:4208386/150‑1 (MQ=255)

                                                                                                       gTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCGCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCGGGGCGTTCATGGGTAGCGactact                                                <  1:647333/150‑1 (MQ=255)

                                                                                                        tGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCGGGGCGTTCATGGGTAGCGACTACTc                                               <  2:2713902/150‑1 (MQ=255)

                                                                                                          aTGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCGGGGCGTTCATGGGTAGCGACTACTCCt                                             <  1:430092/150‑1 (MQ=255)

                                                                                                           tgtgCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCGGGGCGTTCATGGGTAGCGACTACTCCTa                                            >  2:1540038/1‑150 (MQ=255)

                                                                                                           tgtgCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCGGGGCGTTCATGGGTAGCGACTACTCCTa                                            >  1:4031531/1‑150 (MQ=255)

                                                                                                            gtgCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCGGGGCGTTCATGGGTAGCGACTACTCCTAt                                           <  1:455805/150‑1 (MQ=255)

                                                                                                             tgCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCGGGGCGTTCATGGGTAGCGACTACTCCTATg                                          >  1:559813/1‑150 (MQ=255)

                                                                                                              gCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCGGGGCGTTCATGGGTAGCGACTACTCCTATGc                                         >  1:3697575/1‑150 (MQ=255)

                                                                                                               cGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTGGAAGACAGCCAGTGGATGTATTTGCCTGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCGGGGCGTTCATGGGTAGCGACTACTCCTATGcc                                        <  1:3847082/150‑1 (MQ=255)

                                                                                                                     gTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCGGGGCGTTCATGGGTAGCGACTACTCCTATGCCGACCTc                                  >  1:4678657/1‑150 (MQ=255)

                                                                                                                     gTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCGGGGCGTTCATGGGTAGCGACTACTCCTATGCCGACCTc                                  >  1:4678619/1‑150 (MQ=255)

                                                                                                                          cTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCGGGGCGTTCATGGGTAGCGACTACTCCTATGCCGACCTCGACaa                             >  2:2867504/1‑150 (MQ=255)

                                                                                                                          cTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCGGGGCGTTCATGGGTAGCGACTACTCCTATGCCGACCTCGACaa                             >  2:2867571/1‑150 (MQ=255)

                                                                                                                           ttCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCGGGGCGTTCATGGGTAGCGACTACTCCTATGCCGACCTCGACAAg                            >  2:5434910/1‑150 (MQ=255)

                                                                                                                           ttCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACACGAAGGATTTCCACCACCGACAAACGCGGGGCGTTCATGGGTAGCGACAACTCCTATGCCGACCTCGACAAg                            >  2:740531/1‑150 (MQ=255)

                                                                                                                               acaTCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGaa                                                                                                <  1:3583730/78‑1 (MQ=255)

                                                                                                                               acaTCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGaa                                                                                                >  2:3583730/1‑78 (MQ=255)

                                                                                                                               acaTCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCGGGGCGTTCATGGGTAGCGACTACTCCTATGCCGACCTCGACAAGATTc                        >  2:2581940/1‑150 (MQ=255)

                                                                                                                               acaTCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCGGGGCGTTCATGGGTAGCGACTACTCCTATGCCGACCTCGACAAGATTc                        <  2:1136326/150‑1 (MQ=255)

                                                                                                                                 aTCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCGGGGCTTTCATGGGTAGCGACTACTCCTATGCCGACCTCGACAAGATTcgc                      >  2:1662081/1‑150 (MQ=255)

                                                                                                                                 aTCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCGGGGCGTTCATGGGTAGCGACTACTCCTATGCCGACCTCGACAAGATTcgc                      <  2:2309173/150‑1 (MQ=255)

                                                                                                                                   cGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCGGGGCGTTCATGGGTAGCGACTACTCCTATGCCGACCTCGACAAGATTCgcgc                    >  2:2806481/1‑150 (MQ=255)

                                                                                                                                    gAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCGGGGCGTTCATGGGTAGCGACTACTCCTATGCCGACCTCGACAAGATTCgcgc                    >  2:2816398/1‑149 (MQ=255)

                                                                                                                                    gAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCGGGGCGTTCATGGGTAGCGACTACTCCTATGCCGACCTCGACAAGATTCgcgc                    <  1:2816398/149‑1 (MQ=255)

                                                                                                                                    gAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCGGGGCGTTCATGGGTAGCGACTACTCCTATGCCGACCTCGACAAGATTCGCGcc                   <  1:5349468/150‑1 (MQ=255)

                                                                                                                                       aaCCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCGGGGCGTTCATGGGTAGCGACTACTCCTATGCCGACCTCGACAAGATTCGCGCCGaa                <  1:2723883/150‑1 (MQ=255)

                                                                                                                                        aCCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCGGGGCGTTCATGGGTAGc                                                       <  2:4949883/110‑1 (MQ=255)

                                                                                                                                        aCCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCGGGGCGTTCATGGGTAGc                                                       >  1:4949883/1‑110 (MQ=255)

                                                                                                                                          ccGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCGGGGCGTTCATGGGTAGCGACTACTCCTATGCCGACCTCGACAAGATTCGCGCCGAAGAc             >  2:1464766/1‑150 (MQ=255)

                                                                                                                                          ccGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCGGGGCGTTCATGGGTAGCGACTACTCCTATGCCGACCTCGACAAGATTCGCGCCGAAGAc             >  2:464193/1‑150 (MQ=255)

                                                                                                                                           tcttcgAGGCCCAGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCGGGGCGTTCATGGGTAGCGACTACTCCTATGCCGACCTCGACAAGATTCGCGCCGAAGACt            >  1:687777/5‑150 (MQ=255)

                                                                                                                                           cgcgcgAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCGGGGCGTTCATGGGTAGCGACTACTCCTATGCCGACCTCGACAAGATTCGCGCCGAAGACt            <  1:5434910/150‑1 (MQ=255)

                                                                                                                                            gcgcgAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCGGGGCGTTCATGGGTAGCGACTACTCCTATGCCGACCTCGACAAGATTCGCGCCGAAGACTa           <  1:4744503/150‑1 (MQ=255)

                                                                                                                                              gcgAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCGGGGCGTTCATGGGTAGCGACTACTCCTATGCCGACCTCGACAAGATTCGCGCCGAAGACTATc         <  1:2518820/150‑1 (MQ=255)

                                                                                                                                                gAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCGGGGCGTTCATGAGTAGCGACTACTCCTATGCCGACCTCGACAAGATTCGCGCCGAAGACTATCGc       >  1:2314007/1‑150 (MQ=255)

                                                                                                                                                   gCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCGGGGCGTTCATGGGTAGCGACTACTCCTATGCCGACCTCGACAAGATTCGCGCCGAAGACTATCGCCaa    >  1:3476145/1‑150 (MQ=255)

                                                                                                                                                    cccTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCGGGGCGTTCATGGGTAGCGACTACTCCTATGCCGACCTCGACAAGATTCGCGCCGAAGACTATCGCCaac   <  2:3136963/150‑2 (MQ=255)

                                                                                                                                                     ccTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCGGGGCGTTCATGGGTAGCGACTACTCCTATGCCGACCTCGACAAGATTCGCGCCGAAGACTATCGCCAAAc  >  1:4179339/1‑150 (MQ=255)

CCTGATCGACAAGCAAGGTAACAGCCGCATCCGCGAGTTCAGCTACCTGCAGAAGGATTACCCGGACAGCGACAAATTCACCATGTATTTCTC‑C‑GCGCCAAGTGATGTGCGTGACGTCGCCTTCCACATCGAAAACCCGCGCGAGGCCCTGGGCCTCGAAGACAGCCAGTGGATGTATTTGCCGGTCAGCCGACAGACCCGAAGGATTTCCACCACCGACAAACGCGGGGCGTTCATGGGTAGCGACTACTCCTATGCCGACCTCGACAAGATTCGCGCCGAAGACTATCGCCAAAC  >  NC_002947/3175947‑3176243

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