BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,243,659	T→A	6.4%	A58A (GCT→GCA)	PP_5701 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,243,659	0	T	A	6.4%	989.6 / 39.6	327	A58A (GCT→GCA)	PP_5701	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (142/157); new base A (6/15); total (153/174)
Fisher's exact test for biased strand distribution p-value = 1.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.00e-01

Only 100 of 332 total aligned reads displayed.

GGCACCTTGAGGCATGGCACGATGCATGAACTGGAGGCGCTGCGGCCTCTTCTTGAGTCGGCGCTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCACAGTAGTCAAGACTTGCACAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCGGACCTGCCGACAAATATCCCCGGCCATGCACCGTCGCAGCCGCCACCCCATTCGCGGGCAAGCCTGTGC  >  NC_002947/5243510‑5243808

     cTTGAGGCATGGCACGATGCATGAACTGGAGGCGCTGCGGCCTCTTCTTGAGTCGGCGCTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCC                                                                                                                                                  <  1:6666874‑M1/150‑1 (MQ=255)

         aGGCATGGCACGATGCATGAACTGGAGGCGCTGCGGCCTCTTCTTGAGTCGGCGCTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTC                                                                                                                                              <  1:156312‑M1/150‑1 (MQ=255)

         aGGCATGGCACGATGCATGAACTGGAGGCGCTGCGGCCTCTTCTTGAGTCGGCGCTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTC                                                                                                                                              <  1:6432943‑M1/150‑1 (MQ=255)

         aGGCATGGCACGATGCATGAACTGGAGGCGCTGCGGCCTCTTCTTGAGTCGGCGCTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTC                                                                                                                                              <  2:6582119‑M1/150‑1 (MQ=255)

         aGGCATGGCACGATGCATGAACTGGAGGCGCTGCGGCCTCTTCTTGAGTCGGCGCTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTC                                                                                                                                              <  1:4736966‑M1/150‑1 (MQ=255)

         aGGCATGGCACGATGCATGAACTGGAGGCGCTGCGGCCTCTTCTTGAGTCGGCGCTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTC                                                                                                                                              >  2:1907753‑M1/1‑150 (MQ=255)

            cATGGCACGATGCATGAACTGGAGGCGCTGCGGCCTCTTCTTGAGTCGGCGCTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCA                                                                                                                                           <  2:4607632‑M1/150‑1 (MQ=255)

            cATGGCACGATGCATGAACTGGAGGCGCTGCGGCCTCTTCTTGAGTCGGCGCTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCA                                                                                                                                           <  2:817066‑M1/150‑1 (MQ=255)

              tGGCACGATGCATGAACTGGAGGCGCTGCGGCCTCTTCTTGAGTCGGCGCTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACA                                                                                                                                         >  1:1428685‑M1/1‑150 (MQ=255)

              tGGCACGATGCATGAACTGGAGGCGCTGCGGCCTCTTCTTGAGTCGGCGCTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACA                                                                                                                                         >  1:4590288‑M1/1‑150 (MQ=255)

                gctgcgggagcgcctttgccaggaaccatttgcccgaaacagttgtaggcacaggcttgcccgcgaatggggtggcggctgcgacggtgcatggccggggatatttgtcggcaggtccggccctTTCGCGGGCTCGCCCGCTCCCACAGc                                                                                                                                       <  2:4421903‑M1/26‑1 (MQ=255)

                gCACGATGCATGAACTGGAGGCGCTGCGGCCTCTTCTTGAGTCGGCGCTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGc                                                                                                                                       >  2:3307354/1‑150 (MQ=255)

                 cACGATGCATGAACTGGAGGCGCTGCGGCCTCTTCTTGAGTCGGCGCTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCg                                                                                                                                      <  1:5045735/150‑1 (MQ=255)

                   cGATGCATGAACTGGAGGCGCTGCGGCCTCTTCTTGAGTCGGCGCTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCGGt                                                                                                                                    >  1:1793136/1‑150 (MQ=255)

                    gATGCATGAACTGGAGGCGCTGCGGCCTCTTCTTGAGTCGGCGCTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTc                                                                                                                                   <  1:5721948/150‑1 (MQ=255)

                     aTGCATGAACTGGAGGCGCTGCGGCCTCTTCTTGAGTCGGCGCTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTca                                                                                                                                  <  2:1153387/150‑1 (MQ=255)

                        cATGAACTGGAGGCGCTGCGGCCTCTTCTTGAGTCGGCGCTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTcaccc                                                                                                                               >  1:5432263/1‑149 (MQ=255)

                         aTGAACTGGAGGCGCTGCGGCCTCTTCTTGAGTCGGCGCTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTcaccac                                                                                                                              >  1:941596/1‑150 (MQ=255)

                          tGAACTGGAGGCGCTGCGGCCTCTTCTTGAGTCGGCGCTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACcaca                                                                                                                             <  2:4241312/150‑1 (MQ=255)

                          tGAACTGGAGGCGCTGCGGCCTCTTCTTGAGTCGGCGCTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACGTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACcaca                                                                                                                             >  1:585531/1‑150 (MQ=255)

                             aCTGGAGGCGCTGCGGCCTCTTCTTGAGTCGGCGCTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCACagta                                                                                                                          <  1:7398158/150‑1 (MQ=255)

                                 gAGGCGCTGCGGCCTCTTCTTGAGTCGGCGCTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGCGGTCACCACAGTAGTCa                                                                                                                      >  1:7033642/1‑150 (MQ=255)

                                    gcgcTGCGGCCTCTTCTTGAGTCGGCGCTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCACAGTAGTCAAGa                                                                                                                   >  1:7426596/1‑150 (MQ=255)

                                    gcgcTGCGGCCTCTTCTTGAGTCGGCGCTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCACAGTAGTCAAGa                                                                                                                   <  2:200411/150‑1 (MQ=255)

                                      gcTGCGGCCTCTTCTTGAGTCGGCGCTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGACGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCACAGTAGTCAAGACt                                                                                                                 >  2:7126518/1‑150 (MQ=255)

                                      gcTGCGGCCTCTTCTTGAGTCGGCGCTGGTTCAGGTGAAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGGTCACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCACCCACAGCGGTCACCACAGTAGTCAAGACt                                                                                                                 >  2:2524618/1‑150 (MQ=255)

                                       cTGCGGCCTCTTCTTGAGTCGGCGCTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCTGGCTCGCCCGCTCCCACAGCGGTCACCACAGTAGTCAAGACtt                                                                                                                <  1:2747107/150‑1 (MQ=255)

                                       cTGCGGCCTCTTCTTGAGTCGGCGCTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCACAGTAGTCAAGACtt                                                                                                                <  2:1555487/150‑1 (MQ=255)

                                       cTGCGGCCTCTTCTTGAGTCGGCGCTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCACAGTAGTCAAGACtt                                                                                                                <  2:267852/150‑1 (MQ=255)

                                        tGCGGCCTCTTCTTGAGTCGGCGCTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCACAGTAGTCAAGACTTg                                                                                                               <  2:5039402/150‑1 (MQ=255)

                                            gCCTCTTCTTGAGTCGGCGCTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCACAGTAGTCAAGACTTGcaca                                                                                                           <  1:6348376/150‑1 (MQ=255)

                                             ccTCTTCTTGAGTCGGCGCTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCACAGTAGTCAAGACTTGCACAg                                                                                                          >  2:3832295/1‑150 (MQ=255)

                                              ctctTCTTGAGTCGGCGCTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCACAGTAGTCAAGACTTGCACAGt                                                                                                         >  1:5330092/1‑150 (MQ=255)

                                                cttcttGAGTCGGCGCTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCACAGTAGTCAAGACTTGCACAGTAc                                                                                                       >  1:4139687/1‑150 (MQ=255)

                                                cttcttGAGTCGGCGCTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCACAGTAGTCAAGACTTGCACAGTAc                                                                                                       >  1:579419/1‑150 (MQ=255)

                                                      gAGTCGGCGCTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGTCCGGCCCTTTCGCGGGCACGCCCGCTCCCACAGgtactgtgcaagtcttgactactgtggtgaccgctgtgg                                                                                                 <  1:2249335‑M1/150‑40 (MQ=255)

                                                      gAGTCGGCGCTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCACAGTAGTCAAGACTTGCACAGTACctgtgg                                                                                                 >  2:5223166‑M2/1‑144 (MQ=255)

                                                        gTCGGCGCTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGgtactgtgcaagtcttgactactgtggtgaccgctgtggga                                                                                               <  2:4787598‑M1/150‑42 (MQ=255)

                                                        gTCGGCGCTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGgtactgtgcaagtcttgactactgtggtgaccgctgtggga                                                                                               <  2:1434842‑M1/150‑42 (MQ=255)

                                                        gTCGGCGCTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCACAGTAGTCAAGACTTGCACAGTACctgtggga                                                                                               >  1:1968775‑M2/1‑142 (MQ=255)

                                                        agttgtaggcacaggcttgcccgcgaatggggtggcggctgcgacggtgcatggccggggatatttgtcggcaggtccggcctcTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCACAGTAGTCAAGACTTGCACAGTACCTGTGGGa                                                                                               >  2:7438838‑M1/85‑150 (MQ=255)

                                                            gcgcTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCACAGTAGTCAAGACTTGCACAGTACctgtgggagcgg                                                                                           <  2:2323380‑M2/150‑13 (MQ=255)

                                                            gcgcTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCACAGTAGTCAAGACTTGCACAGTACctgtgggagcgg                                                                                           >  1:5131187‑M2/1‑138 (MQ=255)

                                                             cgcTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCACAGTAGTCAAGACTTGAACAGTACctgtgggagcggg                                                                                          >  1:1889784‑M2/1‑137 (MQ=255)

                                                               cTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCACAGTAGTCAAGACTTGCACAGTACctgtgggagcgggcg                                                                                        <  2:7181379‑M2/150‑16 (MQ=255)

                                                                tGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGgtactgtgcaagtcttgactactgtggtgaccgctgtgggagcgggcga                                                                                       <  1:1278581‑M1/150‑50 (MQ=255)

                                                                  gTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCACGCCCGCTCCCACAGgtactgtgcaagtcttgactactgtggtgaccgctgtgggagcgggcgagc                                                                                     <  1:3314968‑M1/150‑52 (MQ=255)

                                                                  acaggcttgcccgcgaatggggtggcggctgcgacggtgcatggccggggatatttgtcggcaggtccggccctTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCACAGTAGTCAAGACTTGCACAGTACCTGTGGGAGCGGGCGTGc                                                                                     >  1:4421903‑M1/75‑150 (MQ=255)

                                                                   ttCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCACAGTAGTCAAGACTTGCACAGTACCTGTGGGAGCGGGCGTGcc                                                                                    >  1:5556267/1‑150 (MQ=255)

                                                                    tCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGCAGGTCCGGCCCTTTCGCGGGCACGCCCGCTCCCACAGgtactgtgcaagtcttgactactgtggtgaccgctgtgggagcgggcgagccc                                                                                   <  2:1072639‑M1/150‑54 (MQ=255)

                                                                    tCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCACAGTAGTCAAGACTTGCACAGTACCTGTGGGAGCGGGCGTGccc                                                                                   <  2:586915/150‑1 (MQ=255)

                                                                          gCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCACAGTAGTCAAGACTTGCACAGTACCTGTGGGAGCGGGCGTGCCCGCGaaa                                                                             <  2:5556267/150‑1 (MQ=255)

                                                                          gCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCACAGTAGTCAAGACTTGCACAGTACCTGTGGGAGCGGGCGTGCCCGCGaaa                                                                             <  1:5679759/150‑1 (MQ=255)

                                                                          gCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCACAGTAGTCAAGACTTGCACAGTACCTGTGGGAGCGGGCGTGCCCGCGaaa                                                                             <  2:1783162/150‑1 (MQ=255)

                                                                           cAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCACAGTAGTCAAGACTTGCACAGTACCTGTGGGAGCGGGCATGCCCGCGAAAg                                                                            <  2:3917544/150‑1 (MQ=255)

                                                                            aGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCACAGTAGTCAAGACTTGCACAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAgg                                                                           >  1:3141030/1‑150 (MQ=255)

                                                                               tGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCACAGTAGTCAAGACTTGCACAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAGGGcc                                                                        >  1:6792774/1‑150 (MQ=255)

                                                                                gTCCGCTTAGCCCAGGAACTTCGCACTAAGGCTGGTTACCGTGCCGGACGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACCGCGGTCACCACCGTAGTAAAGACTTGCAAATTACCTGTGGGAGGGGGCGTGCCCGCGAAAGGGCCg                                                                       <  1:2491602/150‑1 (MQ=255)

                                                                                 aatggggtggcggctgcgacggtgcatggccggggatatttgtcggcaggtccggcctcTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCACAGTAGTCAAGACTTGCACAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCgg                                                                      >  1:3473295‑M1/60‑150 (MQ=255)

                                                                                  ccGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCACAGTAGTCAAGACTTGCACAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCGGa                                                                     <  2:6885095/150‑1 (MQ=255)

                                                                                    gCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCACAGTAGTCAAGACTTGCACAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCGGAcc                                                                   >  1:5740711/1‑150 (MQ=255)

                                                                                    gCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCCCAGCGGTCCCCCCAGTAGTCAAGACTTGCACAGTACCGGTGGGAGCGGGCGTGCCCGCGAAAGGGCCGGAcc                                                                   >  1:3599534/1‑150 (MQ=38)

                                                                                     gggtggcggctgcgacggtgcatggccggggatatttgtcggcagggccggcctcTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCACAGTAGTCAAGACTTGCACAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCGGACCt                                                                  >  1:1305295‑M1/56‑150 (MQ=255)

                                                                                     cTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACTTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACCGCGGTCACCACTGTAGTCAAGACTTGCACAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCGGACCt                                                                  <  1:1078079/150‑1 (MQ=255)

                                                                                       gAGCCCAGTAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCACAGTAGTCAAGACTTGCACAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCGGACCTGc                                                                <  2:3468100/150‑1 (MQ=255)

                                                                                       gAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCACAGTAGTCAAGACTTGCACCGTACCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCGGACCTGc                                                                <  2:1341844/150‑1 (MQ=255)

                                                                                        cgCCCCGGACCTTGGCACTAGGGGTGGTTAACGTGCCGGACGGGCCGGCCTCTTCGGGGGCTCGCCCGCTCCCACCGCGGTCACCACAGTAGTCAAGACTTGCACAGTACCTGTGGGTGGGGGCGTGCCCGCGAAAGGGCCGGACCTGcc                                                               <  1:3833342/149‑1 (MQ=255)

                                                                                         gCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTCAACACAGTAGTCAAGACTTGCACAGTACCTGTGGGAGCGGGCGTGCCCGCGAAACGGCCGGACCAGccc                                                              >  1:2235517/1‑149 (MQ=255)

                                                                                            cAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCACAGTAGTCAAGACTTGCACAGTACctgt                                                                                                   >  1:4203355‑M2/1‑106 (MQ=255)

                                                                                            cAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCACAGTAGTCAAGACTTGCACAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCGGACCTGCCGACa                                                           >  2:6712871/1‑150 (MQ=255)

                                                                                             aGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCACAGTAGTCAAGACTTGCACAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCGGACCTGCCGACaa                                                          >  2:5648931/1‑150 (MQ=255)

                                                                                             aGGAACTCCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCACAGTTGTCAAGACTTGCACAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCGGACCTGCCGACaa                                                          <  2:965410/150‑1 (MQ=255)

                                                                                               gAACTTCGCACTAAGGCTGGTTAACGTGTCGGCAGGTCCGGCCCTGTCGCGGGCACGCCCGCTCCCACAGgtactgtgcaagtcttgactactgtggtgaccgctgtgggagcgggcgagcccgcgaagaggccggcccttccggcacgt                                                        <  2:3408634‑M1/150‑81 (MQ=255)

                                                                                                 aCTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCACAGTAGTCAAGACTTGCACAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCGGACCTGCCGAAAAatat                                                      >  1:608407/1‑150 (MQ=255)

                                                                                                   acggtgcatggccggggatatttgtcggcaggtccggccctTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCACAGTAGTCAAGACTTGCACAGTACCTGTGGGAGCgg                                                                                           <  1:493673‑M1/70‑1 (MQ=255)

                                                                                                   acggtgcatggccggggatatttgtcggcaggtccggccctTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCACAGTAGTCAAGACTTGCACAGTACCTGTGGGAGCgg                                                                                           >  2:493673‑M1/42‑111 (MQ=255)

                                                                                                    tCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCACAGTAGTCAAGACTTGCACAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCGGACCTGCCGACAAATATccc                                                   >  1:2361386/1‑150 (MQ=255)

                                                                                                                       cGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCACAGTAGTCAAGACTTGCACAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCGGACCTGCCGACAAATATCCCCGGCCATGCACCGTCGCAg                                >  2:4132891/1‑150 (MQ=255)

                                                                                                                          gCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCACAGTAGTCAAGACTTGCACAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCGGACCTGCCGACAAATATCCCCGGCCATGCACCGTCGCAgccg                             >  2:5117749/1‑150 (MQ=255)

                                                                                                                           ccGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCACAGTAGTCAAGACTTGCACAGTACctgtgggagcgggcgtgcccgcgaaagggccggacctgccggcacgttaaccagccttagtgcgaagttcctggg                            <  2:5754901‑M2/150‑76 (MQ=255)

                                                                                                                           ccGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCACAGTAGTCAAGACTTGCACAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCGGACCTGCCGACACATATCCCCGGCCATGCACCGTCGCAgccgc                            >  2:6073414/1‑150 (MQ=255)

                                                                                                                            cggcaggtccggccctTTCGCGGGCACGCCCGCTCCCACAGCGGTCACCACAGTAGTCAAGACTTGCACAGTACCTGTgg                                                                                                 <  1:1669960‑M1/64‑1 (MQ=255)

                                                                                                                             ggAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCACAGTAGTCAAGACTTGCACAGTACCTGTGGGAGCGGGCGTGCCCCCGAAAGGGCCGGCCCTGCCGGCAAAAAAACCCCGCCAGGCAACCTCGCAccccccc                          >  1:2116075/1‑149 (MQ=21)

                                                                                                                              gatgggCCGGCCTCTTCGCTGTCTCGCCCGCTCCCACAGCGGTCACCACAGTAGTCAAGACTTGCACAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCGGACCTGCCGACAAATATCCCCGGCCATGCACCGTCGCAGCCGCCAc                         <  1:6751928/147‑1 (MQ=255)

                                                                                                                              gatgggCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCACAGTAGTCAAGACTTGCACAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCGGACCTGCCGACAAATATCCCCGGCCATGCACCGTCGCAGCCGCCAc                         <  2:3692123/147‑1 (MQ=255)

                                                                                                                              gAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCACAGTAGTCAAGACTTGCACAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCGGACCTGCCGACAAATATCCCCGGCCATGCACCGTCGCAGCCGCCAc                         <  2:6353510/150‑1 (MQ=255)

                                                                                                                                  ggccggccTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCACAGTAGTCAAGCCTTGCACAGTACctgtgggagcgggcgtgcccgcgaaagggccggacctgccggcacgttaaccagccttagtgcgaagttcctgggctcagcg                     <  2:3150291‑M2/150‑83 (MQ=255)

                                                                                                                                   gccggccTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCACAGTAGTCAAGACTTGCACAGTACCTGTGGGAGCGGGCGTGGCCCCGGAAGGGCCGGGCCGGCCGACCAAAATACCCCGCCCAGCCCCCGCGCCCCccccaccccaa                    >  1:3361283/1‑149 (MQ=18)

                                                                                                                                    ccggccTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCACAGTAGTCAAGACTTGCACAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCGGACCTGCCGACAAATATCCCCGGCCATGCACCGTCGCAGCCGCCAc                         >  2:4964684/1‑144 (MQ=255)

                                                                                                                                     cggccTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCACAGTAGTCAAGACTTGCACAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCGGACCTGCCGACAAATATCCCCGGCCATGCACCGTCGCAGCCGCCACCCCAGTc                  >  2:4313243/1‑150 (MQ=255)

                                                                                                                                         ctctTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCACAGTAGTCAAGACTTGCACAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCGGACCTGCCGACAAATATCCCCGGCCATGAACCGTCGCAGCCGCCACCCCATTCGCgg              <  2:2548036/150‑1 (MQ=255)

                                                                                                                                          tctTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCACAGTAGTCAAGACTTGCACAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCGGACCTGCCGACCAATATCCCCGGCCATGCACCGTCGCAGCCGCCACCCCATTCGCggg             >  1:7192996/1‑150 (MQ=255)

                                                                                                                                              cgggggCTCGCCCGCTCCCACAGCGGTCACCACCGTAGTCAAGCCTTGCACAGTACCTGTGGGTGGGGGCTTGCCCGCGAAAGGGCCGGACCTGCCGACAAATATCCCCGGCCATGCACCGTCGCAGCCGCCACCCCATTCGCGGGCAAg         <  1:3737378/147‑1 (MQ=38)

                                                                                                                                                 GGGCTCGCCCGCTCCCACAGCGGTCACCACAGTAGTCAAGACTTGCACAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCGGACCTGCCGACAAATATCCCCGGCCATGCACCGTCGCAGCCGCCACCCCATTCGCGGGCAAGCCt      <  1:4260548‑M1/150‑1 (MQ=255)

                                                                                                                                                  GGCTCGCCCGCTCCCACAGCGGTCACCACAGTAGTCAAGACTTGCACAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCGGACCTGCCGACAAATATCCCCGGCCATGCACCGTCGCAGCCGCCACCCCATTCGCGGGCAAGCCtg     <  2:4200764‑M1/150‑1 (MQ=255)

                                                                                                                                                   GCTCGCCCGCTCCCACAGCGGTCACCACAGTAGTCAAGACTTGCACAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCGGACCTGCCGACAAATATCCCCGGCCCTGCACCGTCGCAGCCGCCACCCCATTCGCGGGCAAGCCtgt    <  2:3727441‑M1/150‑1 (MQ=255)

                                                                                                                                                   GCTCGCCCGCTCCCACAGCGGTCACCACAGTAGTCAAGACTTGCACAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCGGACCTGCCGACAAATATCCCCGGCCATGCACCGTCGCAGCCGCCACCCCATTCGCGGGCAAGCCtgt    >  2:621458‑M1/1‑150 (MQ=255)

                                                                                                                                                    CTCGCCCGCTCCCCCAGCGGTCACCACAGTAGTCAAGACTTGCACAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCGGACCTGCCGACAAATATCCCCGGCCATGCACCGTCGCAGCCGCCACCCCATTCGCGGGCAAGCCtgtg   >  2:2814723‑M1/1‑150 (MQ=255)

                                                                                                                                                     tCGCCCGCTCCCACAGCGTTCAACACATTAGAAAAGACTTGCAAAGTACCTTTGGGAGCGGGCGTGCCCGCGAAAGGGCCGGACCTGCCGACAAATATCCCCGGCCATGCACCGTCGCAGCCGCCACCCCATTCGCGGGCAAGCCTGTGc  <  1:846102/150‑1 (MQ=255)

                                                                                                                                                     tCGCCCGCTCCCACAGCGGTCACCACAGTAGTCAAGACTTGCACAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCGGACCTGCCGACAAATATCCCCGGCCATGCACCGTCGCAGCCGCCACCCCATTCGCGGGCAAGCCTGTGc  <  1:2563075/150‑1 (MQ=255)

GGCACCTTGAGGCATGGCACGATGCATGAACTGGAGGCGCTGCGGCCTCTTCTTGAGTCGGCGCTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCACAGTAGTCAAGACTTGCACAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCGGACCTGCCGACAAATATCCCCGGCCATGCACCGTCGCAGCCGCCACCCCATTCGCGGGCAAGCCTGTGC  >  NC_002947/5243510‑5243808

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