BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,442,599	A→T	5.9%	H126L (CAC→CTC)	PP_3060 →	tail sheath protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,442,599	0	A	T	5.9%	643.2 / 19.8	185	H126L (CAC→CTC)	PP_3060	tail sheath protein
Reads supporting (aligned to +/- strand): ref base A (81/93); new base T (4/7); total (85/100)
Fisher's exact test for biased strand distribution p-value = 5.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 186 total aligned reads displayed.

GGCTGACGGCAACGCCTTCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGCCTGTAACCCCCTTTCCCCTC  >  NC_002947/3442450‑3442746

ggCTGACGGCAACGCCTTCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCa                                                                                                                                                     >  1:3270333/1‑150 (MQ=255)

 gCTGACGGCAACGCCTTCAACGGCGTATTCTGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCAc                                                                                                                                                    >  1:2407356/1‑150 (MQ=255)

  cTGACGGCAACGCCTTCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACg                                                                                                                                                   <  2:3258183/150‑1 (MQ=255)

   tGACGGCAACGCCTTCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACgg                                                                                                                                                  >  1:2479478/1‑150 (MQ=255)

   tGACGGCAACGCCTTCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACgg                                                                                                                                                  >  1:3964883/1‑150 (MQ=255)

    gACGGCAACGCCTTCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGc                                                                                                                                                 <  1:691413/150‑1 (MQ=255)

    gACGGCAACGCCTTCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGc                                                                                                                                                 <  1:343299/150‑1 (MQ=255)

       ggCAACGCCTTCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGTGCTTCAGCTCGGCCtt                                                                                                                                              <  2:2784890/150‑1 (MQ=255)

       ggCAACGCCTTCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCAtt                                                                                                                                              <  1:127191/150‑1 (MQ=255)

       ggCAACGCCTTCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCAtt                                                                                                                                              <  2:2975884/150‑1 (MQ=255)

        gCAACGCCTTCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGTGCTTCAGCTCGGCCTTg                                                                                                                                             <  2:629530/150‑1 (MQ=255)

          aaCGCCTTCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGcc                                                                                                                                           <  2:1280292/150‑1 (MQ=255)

          aaCGCCTTCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGcc                                                                                                                                           >  2:4287072/1‑150 (MQ=255)

          aaCGCCTTCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGcc                                                                                                                                           >  1:735349/1‑150 (MQ=255)

          aaCGCCTTCAACGGCGTATTCCGAGGCTCGGTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGcc                                                                                                                                           <  2:3548369/150‑1 (MQ=255)

          aaCGCCTTCAACGGCGTATTCCGAGGCTCGGTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGcc                                                                                                                                           <  2:1981181/150‑1 (MQ=255)

             gCCTTCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGGCAATGCCGTGCTTCAGCTCGGCCTTGTCgcc                                                                                                                                        <  1:1988104/150‑3 (MQ=255)

              ccTTCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCa                                                                                                                                       <  1:3545766/150‑1 (MQ=255)

               cTTCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCAc                                                                                                                                      <  1:2299077/150‑1 (MQ=255)

               cTTCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCAc                                                                                                                                      <  1:2931295/150‑1 (MQ=255)

                 tCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACct                                                                                                                                    <  2:4512208/150‑1 (MQ=255)

                   aaCGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACctac                                                                                                                                  <  1:1430645/150‑1 (MQ=255)

                      ggCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCtactac                                                                                                                               <  2:410225/150‑1 (MQ=255)

                      ggCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCtactac                                                                                                                               <  2:410338/150‑1 (MQ=255)

                           aTTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCt                                                                                                                          <  1:1008313/150‑1 (MQ=255)

                            ttCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTc                                                                                                                         >  1:304023/1‑150 (MQ=255)

                            ttCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTc                                                                                                                         <  1:801025/150‑1 (MQ=255)

                            ttCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTc                                                                                                                         <  1:3126810/150‑1 (MQ=255)

                            ttCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTc                                                                                                                         <  1:2767511/150‑1 (MQ=255)

                            ttCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTc                                                                                                                         <  1:3546598/150‑1 (MQ=255)

                            ttCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTc                                                                                                                         <  1:744612/150‑1 (MQ=255)

                              ccGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGa                                                                                                                       <  2:4181332/150‑1 (MQ=255)

                               cGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGaa                                                                                                                      <  2:238098/150‑1 (MQ=255)

                               cGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGaa                                                                                                                      <  1:528410/150‑1 (MQ=255)

                               cGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGaa                                                                                                                      >  2:4169035/1‑150 (MQ=255)

                               cGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGaa                                                                                                                      <  1:3374467/150‑1 (MQ=255)

                               cGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGaa                                                                                                                      >  1:2022338/1‑150 (MQ=255)

                               cGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGaa                                                                                                                      <  1:1205210/150‑1 (MQ=255)

                                 aGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGt                                                                                                                    >  1:2119588/1‑150 (MQ=255)

                                     tCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTcggc                                                                                                                >  1:4353405/1‑150 (MQ=255)

                                     tCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTcggc                                                                                                                >  1:2521512/1‑150 (MQ=255)

                                           aaGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAg                                                                                                          <  2:3963049/150‑1 (MQ=255)

                                               tCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCa                                                                                                      >  1:3446008/1‑150 (MQ=255)

                                                  aGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGTGCTTCAGCTCgg                                                                                                                                                  <  2:2739366/103‑1 (MQ=255)

                                                      gggCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACga                                                                                               <  1:1107977/150‑1 (MQ=255)

                                                      gggCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACga                                                                                               <  2:4485473/150‑1 (MQ=255)

                                                      gggCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACga                                                                                               <  2:3018997/150‑1 (MQ=255)

                                                          gAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGAtt                                                                                           <  1:1407320/150‑1 (MQ=255)

                                                               ccTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGAccc                                                                                      <  1:1750785/150‑1 (MQ=255)

                                                                 ttGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCgg                                                                                    >  1:3991701/1‑150 (MQ=255)

                                                                    cGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCg                                                                                 >  1:1946024/1‑150 (MQ=255)

                                                                       tcgtcgtcACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCa                                                                              <  2:2479478/150‑1 (MQ=255)

                                                                       tcgtcgtcACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCa                                                                              >  2:1920631/1‑150 (MQ=255)

                                                                            gtcACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAg                                                                         >  2:527109/1‑150 (MQ=255)

                                                                              cACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCg                                                                       >  1:449978/1‑150 (MQ=255)

                                                                              cACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCg                                                                       >  1:2704933/1‑150 (MQ=255)

                                                                               aCCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCgt                                                                      >  2:829856/1‑150 (MQ=255)

                                                                                  cTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGtc                                                                   <  2:3285816/150‑1 (MQ=255)

                                                                                       cgGCACCCAGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAATCTAGAAGTCGGAGGCGAGGTCATATACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCaa                                                              <  1:922038/150‑1 (MQ=255)

                                                                                        ggCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAAc                                                             >  1:2056475/1‑150 (MQ=255)

                                                                                           aCCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGaa                                                                                                                      >  2:3406191/1‑90 (MQ=255)

                                                                                           aCCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGc                                                          >  2:1690335/1‑150 (MQ=255)

                                                                                           aCCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGc                                                          >  2:1725374/1‑150 (MQ=255)

                                                                                             ccTGAAAGGGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCAc                                                        <  1:49376/150‑1 (MQ=255)

                                                                                             ccTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATTAAGCGTGTCATCAACGGCAc                                                        <  1:2792115/150‑1 (MQ=255)

                                                                                                gAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACTGa                                                     >  1:3701962/1‑150 (MQ=255)

                                                                                                gAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGa                                                     <  2:2123474/150‑1 (MQ=255)

                                                                                                        tGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTgg                                             >  1:3194182/1‑150 (MQ=255)

                                                                                                        tGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTgg                                             >  1:3810943/1‑150 (MQ=255)

                                                                                                         ggACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGc                                            <  2:2876627/150‑1 (MQ=255)

                                                                                                          gACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGcc                                           >  1:1493128/1‑150 (MQ=255)

                                                                                                          gACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGcc                                           <  1:1772209/150‑1 (MQ=255)

                                                                                                           aCATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCa                                          <  2:3775825/150‑1 (MQ=255)

                                                                                                           aCATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCa                                          <  1:3977054/150‑1 (MQ=255)

                                                                                                            cATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAg                                         >  1:1516681/1‑150 (MQ=255)

                                                                                                               gggCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCa                                      >  1:1533461/1‑150 (MQ=255)

                                                                                                                     cTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAgcgcgc                                <  1:3824319/150‑1 (MQ=255)

                                                                                                                      tGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGcc                               <  2:1236700/150‑1 (MQ=255)

                                                                                                                      tGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGcc                               >  2:1338442/1‑150 (MQ=255)

                                                                                                                        gAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGa                             >  2:4289671/1‑150 (MQ=255)

                                                                                                                         aaGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCCCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGAc                            <  1:1789658/150‑1 (MQ=255)

                                                                                                                             ccGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCg                        >  2:1742244/1‑150 (MQ=255)

                                                                                                                              cGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCgg                       >  2:2966000/1‑150 (MQ=255)

                                                                                                                               ggCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGc                      <  2:1127949/150‑1 (MQ=255)

                                                                                                                               ggCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGc                      <  1:642988/150‑1 (MQ=255)

                                                                                                                                  ggcAATGCCGCGCTTCAGCTCGGCCTTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGCCTg                   >  2:2613609/3‑150 (MQ=255)

                                                                                                                                  gACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGCCTg                   >  1:2183234/1‑150 (MQ=255)

                                                                                                                                   aCAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGCCTGt                  >  2:4392244/1‑150 (MQ=255)

                                                                                                                                     aaGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGCCTGTaa                >  2:2360596/1‑150 (MQ=255)

                                                                                                                                       tgCCGTGCTTCAGCTCGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGCCTGTAAcc              >  2:1953662/2‑150 (MQ=255)

                                                                                                                                       ggCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGCCTGTAAcc              <  1:3100034/150‑1 (MQ=255)

                                                                                                                                        gCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGCCTGTAAccc             >  1:1604712/1‑150 (MQ=255)

                                                                                                                                          cGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGCCTGTAAccccc           >  1:3590329/1‑150 (MQ=255)

                                                                                                                                          cGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGCCTGTAAccccc           <  2:4380065/150‑1 (MQ=255)

                                                                                                                                          cGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGCCTGTAAccccc           <  2:4380234/150‑1 (MQ=255)

                                                                                                                                             gCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGCCTGTAACCCCCttt        <  1:4169035/150‑1 (MQ=255)

                                                                                                                                              cTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGCCTGTAACCCCCTTTc       >  2:369901/1‑150 (MQ=255)

                                                                                                                                               tGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGCCTGAAACCCCCTATcc      >  1:3602163/1‑150 (MQ=255)

                                                                                                                                                gAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGCCTGTAACCCCCTTTccc     <  2:4113721/150‑1 (MQ=255)

                                                                                                                                                 aaGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGCCTGTAACCCCCTTTcccc    >  2:3099265/1‑150 (MQ=255)

GGCTGACGGCAACGCCTTCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGCCTGTAACCCCCTTTCCCCTC  >  NC_002947/3442450‑3442746

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