BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,442,590	A→T	5.6%	E123V (GAG→GTG)	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,590	0	A	T	5.6%	635.6 / 17.2	182	E123V (GAG→GTG)	PP_3060	tail sheath protein
Reads supporting (aligned to +/- strand): ref base A (79/92); new base T (3/7); total (83/99)
Fisher's exact test for biased strand distribution p-value = 5.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 187 total aligned reads displayed.

CTTCGGCCTGGCTGACGGCAACGCCTTCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGCCTGTAACCCCC  >  NC_002947/3442441‑3442737

cttcGGCCTGGCTGACGGCAACGCCTTCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGa                                                                                                                                                     <  1:1230501/150‑1 (MQ=255)

 ttcGGCCTGGCTGACGGCAACGCCTTCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAg                                                                                                                                                    <  2:4257658/150‑1 (MQ=255)

   cGGCCTGGCTGACGGCAACGCCTTCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCt                                                                                                                                                  <  2:850434/150‑1 (MQ=255)

        tGGCTGACGGCAACGCCTTCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGc                                                                                                                                             >  2:1496231/1‑150 (MQ=255)

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

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

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

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

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

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

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

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

                     cGCCTTCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGt                                                                                                                                <  1:3598479/150‑1 (MQ=255)

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

                      gCCTTCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGa                                                                                                                                                     >  2:2097195/1‑128 (MQ=255)

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

                         ttCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCAcc                                                                                                                            <  1:1435863/150‑1 (MQ=255)

                         ttCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCAcc                                                                                                                            >  1:3235025/1‑150 (MQ=255)

                         ttCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCAcc                                                                                                                            >  2:784998/1‑150 (MQ=255)

                            aaCGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACctac                                                                                                                         >  1:876661/1‑150 (MQ=255)

                              cGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACctacta                                                                                                                       <  1:178375/150‑1 (MQ=255)

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

                                 cGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACaa                                                                                                                    >  1:1410910/1‑150 (MQ=255)

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

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

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

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

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

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

                                     ttCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACAAGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGcgc                                                                                                                >  2:3513403/1‑148 (MQ=255)

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

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

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

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

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

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

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

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

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

                                                     aGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAgg                                                                                                >  2:4340520/1‑150 (MQ=255)

                                                     aGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAgg                                                                                                >  1:1389182/1‑150 (MQ=255)

                                                       aTCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTc                                                                                              >  2:2973999/1‑150 (MQ=255)

                                                       aTCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTc                                                                                              <  2:3613896/150‑1 (MQ=255)

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

                                                         ccAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCAt                                                                                            <  2:4479241/150‑1 (MQ=255)

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

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

                                                            gAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTa                                                                                         >  1:1966674/1‑150 (MQ=255)

                                                               gggCGCAACCCGTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACga                                                                                      <  2:4485984/150‑1 (MQ=255)

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

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

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

                                                                   gAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAATCTCGAAGTCGGCGGCGAGGTCATCTACGAGAtt                                                                                  <  1:2425759/150‑1 (MQ=255)

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

                                                                               gtcgtcgtcACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGc                                                                      <  1:2032983/150‑1 (MQ=255)

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

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

                                                                                      tcACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGc                                                               >  2:471603/1‑150 (MQ=255)

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

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

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

                                                                                          ccTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCgtgt                                                           <  2:1282389/150‑1 (MQ=255)

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

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

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

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

                                                                                                     cccTGGAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCa                                                >  2:1441367/1‑150 (MQ=255)

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

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

                                                                                                      ccTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCAc                                               <  2:3377015/150‑1 (MQ=255)

                                                                                                      ccTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCAc                                               >  1:638699/1‑150 (MQ=255)

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

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

                                                                                                         gAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGa                                            >  1:992407/1‑150 (MQ=255)

                                                                                                               gCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACt                                      <  2:2949151/150‑1 (MQ=255)

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

                                                                                                                   gACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGGCGTCACCTACTACAAGCTCGACGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGcc                                  <  1:1771986/150‑1 (MQ=255)

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

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

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

                                                                                                                       tGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTc                              <  1:3354254/150‑1 (MQ=255)

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

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

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

                                                                                                                        gggCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCa                             <  1:1625964/150‑1 (MQ=255)

                                                                                                                          gCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAgc                           <  1:2936607/150‑1 (MQ=255)

                                                                                                                          gCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAgc                           <  2:1946024/150‑1 (MQ=255)

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

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

                                                                                                                                     gCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTc                <  2:1716837/150‑1 (MQ=255)

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

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

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

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

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

                                                                                                                                            aCAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCCGTCAGCGCGCCGACCTCGGCCTga         >  2:4392108/1‑149 (MQ=255)

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

                                                                                                                                              aaGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGCCTGTaa       <  1:3042061/150‑1 (MQ=255)

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

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

CTTCGGCCTGGCTGACGGCAACGCCTTCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGAGCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGCCTGTAACCCCC  >  NC_002947/3442441‑3442737

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