BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,729,098	Δ1 bp	8.9%	coding (346/645 nt)	PP_3292 →	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,729,098	0	A	.	8.9%	468.1 / 40.0	146	coding (346/645 nt)	PP_3292	membrane protein
Reads supporting (aligned to +/- strand): ref base A (64/68); new base . (6/7); total (71/75)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.71e-01

Only 100 of 153 total aligned reads displayed.

CGGTGCAGCGGGGCGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCGGCAAACCGGTGTGCACGTGCCTGGCGCAGACCTT  >  NC_002947/3728949‑3729246

cGGTGCAGCGGGGCGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGca                                                                                                                                                      <  2:107740/150‑1 (MQ=255)

  gTGCAGCGGGGCGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGTCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGcaca                                                                                                                                                    <  2:1009558/150‑1 (MQ=255)

  gTGCAGCGGGGCGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGcaca                                                                                                                                                    <  1:2131948/150‑1 (MQ=255)

    gCAGCGGGGCGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAgg                                                                                                                                                  <  1:5376329/150‑1 (MQ=255)

     cAGCGGGGCGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGc                                                                                                                                                 >  1:5201546/1‑150 (MQ=255)

      aGCGGGGCGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCa                                                                                                                                                >  2:3960725/1‑150 (MQ=255)

      aGCGGGGCGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCa                                                                                                                                                >  1:2428453/1‑150 (MQ=255)

       gCGGGGCGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGTGCCCAGCAGGGC‑CAGGCAc                                                                                                                                               <  2:69487/150‑1 (MQ=255)

       gCGGGGCGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGTGCCCAGCAGGGC‑CAGGCAc                                                                                                                                               <  2:69322/150‑1 (MQ=255)

           ggCGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGTGCCCAGCAGGGC‑CAGGCAcaaca                                                                                                                                           <  2:4869124/150‑5 (MQ=255)

           ggCGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGTGCCCAGCAGGGC‑CAGGCAcaaca                                                                                                                                           <  2:4164075/150‑5 (MQ=255)

              gTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTtgg                                                                                                                                        <  1:324179/150‑1 (MQ=255)

              gTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTtgg                                                                                                                                        <  2:1137040/150‑1 (MQ=255)

                  aGCGCACCGTTGGCAAAGCCATTTTCAACCGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTtggctgg                                                                                                                                    <  2:3462829/150‑1 (MQ=255)

                   gcgcACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAgg                                                                                                                                                  >  1:1261860/1‑135 (MQ=255)

                    cgcACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCa                                                                                                                                  >  1:5161272/1‑150 (MQ=255)

                    cgcACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCa                                                                                                                                  <  1:2230867/150‑1 (MQ=255)

                       aCCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATTGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAggg                                                                                                                               <  2:1455506/150‑1 (MQ=255)

                       aCCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAggg                                                                                                                               >  2:3795671/1‑150 (MQ=255)

                         cGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGgcg                                                                                                                             >  1:2881103/1‑150 (MQ=255)

                              gCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAg                                                                                                                        >  2:3082398/1‑150 (MQ=255)

                                aaaGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTg                                                                                                                      >  1:2634973/1‑150 (MQ=255)

                                  aGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGAc                                                                                                                    <  2:4572362/150‑1 (MQ=255)

                                   gCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACg                                                                                                                   <  1:2353177/150‑1 (MQ=255)

                                       ttttCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTtg                                                                                                                                         >  2:3818170/1‑124 (MQ=255)

                                       ttttCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTtg                                                                                                                                         <  1:3818170/124‑1 (MQ=255)

                                           cAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGt                                                                                                           <  1:2366623/150‑1 (MQ=255)

                                             aCGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTgg                                                                                                         <  1:956997/150‑1 (MQ=255)

                                             aCGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTgg                                                                                                         <  1:5109632/150‑1 (MQ=255)

                                              cGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTATTGGCTGGCAGGGCGCTTAGTGACGCAGATCGTGGc                                                                                                        <  2:3716722/150‑1 (MQ=255)

                                                  ggaacgccAGTCTGGCCTGCCAGCCTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAAc                                                                                                    >  1:4097973/8‑150 (MQ=255)

                                                  cAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAAc                                                                                                    <  2:3924011/150‑1 (MQ=255)

                                                      agaCAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCg                                                                                                <  1:3691667/150‑1 (MQ=255)

                                                        aCAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAg                                                                                              <  1:1811937/150‑1 (MQ=255)

                                                            gCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCa                                                                                          <  2:5201546/150‑1 (MQ=255)

                                                                    gCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAg                                                                                  >  1:3193619/1‑150 (MQ=255)

                                                                          cTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCaa                                                                            >  2:3105241/1‑150 (MQ=255)

                                                                          cTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCaa                                                                            >  2:789626/1‑150 (MQ=255)

                                                                           tGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAg                                                                           <  2:2130249/150‑1 (MQ=255)

                                                                            ggCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGt                                                                          <  1:2547633/150‑1 (MQ=255)

                                                                             gCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTTGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTg                                                                         >  1:5198647/1‑150 (MQ=255)

                                                                             gCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGa                                                                                                      <  2:1642444/121‑1 (MQ=255)

                                                                             gCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTg                                                                         <  2:338118/150‑1 (MQ=255)

                                                                             gCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTg                                                                         <  1:1333330/150‑1 (MQ=255)

                                                                             gCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTg                                                                         >  1:5198611/1‑150 (MQ=255)

                                                                             gCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCATTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGa                                                                                                      >  1:1642444/1‑121 (MQ=255)

                                                                              cGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTgg                                                                        <  1:1119440/150‑1 (MQ=255)

                                                                               ggggagtCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGa                                                                       <  2:371999/144‑1 (MQ=255)

                                                                               gTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGa                                                                       <  1:3245278/150‑1 (MQ=255)

                                                                                ttCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACAGCGAGGCCAACCTGCAGGCTCAAGTGGat                                                                      >  1:1410356/1‑150 (MQ=255)

                                                                                 tCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGata                                                                     <  1:2036810/150‑1 (MQ=255)

                                                                                   cTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGAtata                                                                   >  1:4379166/1‑150 (MQ=255)

                                                                                   cTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCGTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACAGCGAGGCCAACCTGCAGGCTCAAGTGGAtata                                                                   >  1:1923585/1‑150 (MQ=255)

                                                                                     tCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATAcc                                                                 <  2:12374/150‑1 (MQ=255)

                                                                                      cGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATAccc                                                                <  1:794382/150‑1 (MQ=255)

                                                                                         tCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGaa                                                             <  1:342808/150‑1 (MQ=255)

                                                                                              ggCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGa                                                        <  1:4642561/150‑1 (MQ=255)

                                                                                                cTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGaaa                                                      <  2:564189/150‑1 (MQ=255)

                                                                                                cTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGaaa                                                      <  2:2103019/150‑1 (MQ=255)

                                                                                                cTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCATTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGaaa                                                      <  1:1046462/150‑1 (MQ=255)

                                                                                                 tCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGa                                                                                               <  2:1937808/108‑1 (MQ=255)

                                                                                                 tCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCCACCTGCAGGCTCAAGTGGATATACCCGAAGATGaaaa                                                     >  1:2382540/1‑150 (MQ=255)

                                                                                                  cGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGaaaaa                                                    <  1:4049804/150‑1 (MQ=255)

                                                                                                        tCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACCGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATAGTGGCGAACTGCGAGGCCACCCTGCAGGCTCAAGTGGAGATACCCGTAGATGAAAAACTCAcc                                              >  1:4811309/1‑150 (MQ=255)

                                                                                                        tCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAACAACTCAcc                                              >  1:4811235/1‑150 (MQ=255)

                                                                                                         cGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACcg                                             <  1:3651338/150‑1 (MQ=255)

                                                                                                          gAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACcgc                                            >  2:2706197/1‑150 (MQ=255)

                                                                                                          gAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACcgc                                            >  2:3355628/1‑150 (MQ=255)

                                                                                                           aaCATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCCACCTGCCGGCTCAAGTGGATATACCCGAAGATGAAAAAATCACcgcg                                           >  1:4250771/1‑150 (MQ=255)

                                                                                                           aaCATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATAGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATAACCGAAGATGAAAAACTCACcgcg                                           >  1:4250565/1‑150 (MQ=255)

                                                                                                            aCATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCgcgc                                          <  1:86639/150‑1 (MQ=255)

                                                                                                                 ctgctgGTGGTGCTGTTGTGCCTG‑GCCCGGCTGGGCACAGGCACCGTTTGGCTGGGAGGGCGCTCAGTGACGCAGATCGTGGCGCACTGCGAGGCCAACCAGCAGGCTCAAGTGGATATACCCGCAGATGAAAAACTCACCGCGCGtctc                                     >  1:2172/1‑150 (MQ=255)

                                                                                                                   gctgGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCCACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCgg                                   >  1:4159457/1‑150 (MQ=255)

                                                                                                                   gctgGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCgg                                   <  2:130702/150‑1 (MQ=255)

                                                                                                                    ctgGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGata                                                                     <  2:591797/115‑1 (MQ=255)

                                                                                                                    ctgGTGGTGCTGTTGTGCCTG‑GCCCAGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGata                                                                     >  1:591797/1‑115 (MQ=255)

                                                                                                                          gtgCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCGGCGAACCg                            <  2:4159457/150‑1 (MQ=255)

                                                                                                                           ccaaaaggtgCCTGTGCCCAGCAGGGC‑CAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCGGCAAACCgg                           >  2:4058839/8‑150 (MQ=255)

                                                                                                                            gCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCGGCAAACCGgt                          <  2:1403450/150‑1 (MQ=255)

                                                                                                                              tgttgtGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCGGCAAACCGgtgt                        >  1:1534787/1‑150 (MQ=255)

                                                                                                                              tgttgtGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCGGCAAACCGgtgt                        <  2:2145051/150‑1 (MQ=255)

                                                                                                                                ttgtGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCGGCAAACCGGTGTGc                      >  2:2815843/1‑150 (MQ=255)

                                                                                                                                  gtgCCTGTGCCCAGCAGGGC‑CAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCGGCAAACCGGTGTGCAc                    <  2:501630/150‑1 (MQ=255)

                                                                                                                                   tgCCTGTGCCCAGCAGGGC‑CAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCGGCAAACCGGTGTGCACg                   >  1:2410738/1‑150 (MQ=255)

                                                                                                                                   tgCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCGGCAAACCGGTGTGCACg                   >  1:1621604/1‑150 (MQ=255)

                                                                                                                                    gcctgtgCCCAGCAGGGC‑CAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCGGCAAACCGGTGTGCACGt                  >  1:3518759/7‑150 (MQ=255)

                                                                                                                                     cctgtgCCCAGCAGGGC‑CAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCGGCAAACCGGTGTGCAc                    <  1:684936/142‑1 (MQ=255)

                                                                                                                                     cctgtgCCCAGCAGGGC‑CAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCGGCAAACCGGTGTGCACGTg                 >  2:580974/6‑150 (MQ=255)

                                                                                                                                     cctgtgCCCAGCAGGGC‑CAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCGGCAAACCGGTGTGCACGTg                 >  2:2827190/6‑150 (MQ=255)

                                                                                                                                     ccTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCGGCAAACCGGTGTGCACGTg                 >  1:4217684/1‑150 (MQ=255)

                                                                                                                                          gCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCGGCAAACCGGTGTGCACGTGCCTg             >  2:651591/1‑150 (MQ=255)

                                                                                                                                          gCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCGGCAAACCGGTGTGCACGTGCCTg             >  2:651929/1‑150 (MQ=255)

                                                                                                                                           cccTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCGGCAAACCGGTGTGCACGTGCCTgg            <  1:21197/150‑1 (MQ=255)

                                                                                                                                            ccTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCGGCAAACCGGTGTGCACGTGCCTGgc           >  2:5333247/1‑150 (MQ=255)

                                                                                                                                            ccTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCGGCAAACCGGTGTGCACGTGCCTGgc           >  2:2762046/1‑150 (MQ=255)

                                                                                                                                                 tgGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCGGCAAACCGGTGTGCACGTGCCTGGCGCAGa      >  1:1742777/1‑150 (MQ=255)

                                                                                                                                                 tgGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCGGCAAACCGGTGTGCACGTGCCTGGCGCAGa      >  1:1742720/1‑150 (MQ=255)

                                                                                                                                                 tgGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCGGCAAACCGGTGTGCACGTGCCTGGCGCAGa      >  2:3543374/1‑150 (MQ=255)

                                                                                                                                                     cacaGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCGGCAAACCGGTGTGCACGTGCCTGGCGCAGACCtt  >  1:2089284/1‑150 (MQ=255)

                                                                                                                                                     cacaGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCGGCAAACCGGTGTGCACGTGCCTGGCGCAGACCtt  >  2:1746587/1‑150 (MQ=255)

CGGTGCAGCGGGGCGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTG‑GCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCGGCAAACCGGTGTGCACGTGCCTGGCGCAGACCTT  >  NC_002947/3728949‑3729246

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