BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,598,547	C→G	5.2%	K436N (AAG→AAC)	PP_4071 ←	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,598,547	0	C	G	5.2%	350.2 / 7.0	107	K436N (AAG→AAC)	PP_4071	membrane protein
Reads supporting (aligned to +/- strand): ref base C (52/49); new base G (2/4); total (54/53)
Fisher's exact test for biased strand distribution p-value = 4.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 111 total aligned reads displayed.

GGCGACGCAGAAGGCCAGAACGCCTACAGAAGAACCCACCCACAGCAGTTGGCGCTTGCGCCGAGCCACCCTGACGTTATCGCCGGCAAGGCCGGCCTCTTTGTAGATGACCCGCCCAAAGGCCTGTTGAATGAAGTAGGCCAACGCCTTGGCCTGCGGCTTGCCCTCAAGCAACGGGAGGCTGGTCTTGTAAGGCTGCGCCGCCTCACGGACAAAGGCGTTGCGCACGTCACCGTGCTGTACGACTGACGACCAATACACACCGCGCACCAGCGCCGGGGTGGTGAAACGGTCACT  >  NC_002947/4598400‑4598696

ggcgaCGCAGAAGGCCAGAACGCCTACAGAAGAACCCACCCACAGCAGTTGGCGCTTGCGCCGAGCCACCCTGACGTTATCGCCGGCAAGGCCGGCCTCTTTGTAGCTGACCCGCCCAAAGGCCTGTTGAATGAAGTAGGCCAACGCCtt                                                                                                                                                     <  2:1745840/150‑1 (MQ=255)

 gcgaCGCAGAAGGCCAGAACGCCTACAGAAGAACCCACCCACAGCAGTTGGCGCTTGCGCCGAGCCACCCTGACGTTATCGCCGGCAAGGCCGGCCTCTTTGTAGATGACCCGCCCAAAGGCCTGTTGAATGAAGTAGGCCAACGCCTTg                                                                                                                                                    <  1:600936/150‑1 (MQ=255)

   gaCGCAGAAGGCCAGAACGCCTACAGAAGAACCCACCCACAGCAGTTGGCGCTTGCGCCGAGCCACCCTGACGTTATCGCCGGCAAGGCCGGCCTCTTTGTAGATGACCCGCCCAAAGGCCTGTTGAATGAAGTAGGCCAACGCCTTGGc                                                                                                                                                  <  2:1968565/150‑1 (MQ=255)

   gaCGCAGAAGGCCAGAACGCCTACAGAAGAACCCACCCACAGCAGTTGGCGCTTGCGCCGAGCCACCCTGACGTTATCGCCGGCAAGGCCGGCCTCTTTGTAGATGACCCGCCCAAAGGCCTGTTGAATGAAGTAGGCCAACGCCTTGGc                                                                                                                                                  <  2:2272393/150‑1 (MQ=255)

        aGAAGGCCAGAACGCCTACAGAAGAACCCACCCACAGCAGTTGGCGCTTGCGCCGAGCCACCCTGACGTTATCGCCGGCAAGGCCGGCCTCTTTGTAGATGACCCGCCCAAAGGCCTGTTGAATGAAGTAGGCCAAGGCGTTGGCCtact                                                                                                                                             <  2:2100206/150‑4 (MQ=255)

        aGAAGGCCAGAACGCCTACAGAAGAACCCACCCACAGCAGTTGGCGCTTGCGCCGAGCCACCCTGACGTTATCGCCGGCAAGGCCGGCCTCTTTGTAGATGACCCGCCCAAAGGCCTGTTGAATGAAGTAGGCCAACGCCTTGGCCTGCg                                                                                                                                             >  1:1126608/1‑150 (MQ=255)

            ggCCAGAACGCCTACAGAAGACCCCACCCACAGCAGTTGGCGCTTGCGCTGAGCCACCCTGACGTTATCGCCGGCAAGGCCGGCCTCTTTGTAGATGACCCGCCCAAAGGCCTGTTGAATGAAGTAGGCCAACGCGTTGGCCTGCGGCtt                                                                                                                                         <  2:379859/150‑1 (MQ=255)

            ggCCAGAACGCCTACAGAAGAACCCACCCACAGCAGTTGGCGCTTGCGCCGAGCCACCCTGACGTTATCGCCGGCAAGGCCGGCCTCTTTGTAGATGACCCGCCCAAAGGCCTGTTGAATGAAGTAGGCCAACGCCTTGGCCTGCGGCtt                                                                                                                                         <  2:1126608/150‑1 (MQ=255)

            ggCCAGAACGCCTACAGAAGAACCCACCCACAGCAGTTGGCGCTTGCGCCGAGCCACCCTGACGTTATCGCCGGCAAGGCCGGCCTCTTTGTAGATGACCCGCCCAAAGGCCTGTTGAATGAAGTAGGCCAACGCCTTGGCCTGCGGCtt                                                                                                                                         <  2:1979322/150‑1 (MQ=255)

              ccAGAACGCCTACAGAAGAACCCACCCACAGCAGTTGGCGCTTGCGCCGAGCCACCCTGACGTTATCGCCGGCAAGGCCGGCCTCTTTGTAGATGACCCGCCCAAAGGCCTGTTGAATGAAGTAGGCCAACGCCTTGGCCTGCGGcgggc                                                                                                                                       >  1:571830/1‑146 (MQ=255)

                aGAACGCCTACAGAAGAACCCACCCACAGCAGTTGGCGCTTGCGCCGAGCCACCCTGACGTTATCGCCGGCAAGGCCGGCCTCTTTGTAGATGACCCGCCCAAAGGCCTGTTGAATGAAGTAGGCCAACGCCTTGGCCTGCGGCTTGccc                                                                                                                                     <  2:2147717/150‑1 (MQ=255)

                  aaCGCCTACAGAAGAACCCACCCACAGCAGTTGGCGCTTGCGCCGAGCCACCCTGACGTTATCGCCGGCAAGGCCGGCCTCTTTGTAGATGACCCGCCCAAAGGCCTGTTGAATGAAGTAGGCCAACGCCTTGGCCTGCGGCTTGCCCTc                                                                                                                                   >  1:1399786/1‑150 (MQ=255)

                    cGCCTACAGAAGAACCCACCCACAGCAGTTGTCGCTTGCGCCGAGCCACCCTGACGTTATCGCCGGCAAGGCCGGCCTCTTTGTAGATGACCCGCCCAAAGGCCTGTTGAATGAAGTAGGCCAACGCCTTGGCCTGCGGCTTGCCCTCaa                                                                                                                                 <  1:1578953/150‑1 (MQ=255)

                    cGCCTACAGAAGAACCCACCCACAGCAGTTGGCGCTTGCGCCGAGCCACCCTGACGTTATCGCCGGCAAGGCCGGCCTCTTTGTAGATGACCCGCCCAAAGGCCTGTTGAATGAAGTAGGCCAACGCCTTGGCCTGCGGCTTGCCCTCaa                                                                                                                                 <  1:1628104/150‑1 (MQ=255)

                    cGCCTACAGAAGAACCCACCCACAGCAGTTGGCGCTTGCGCCGAGCCACCCTGACGTTATCGCCGGCAAGGCCGGCCTCTTTGTAGATGACCCGCCCAAAGGCCTGTTGAATGAAGTAGGCCAACGCCTTGGCCTGCGGCTTGCCCTCaa                                                                                                                                 <  1:923823/150‑1 (MQ=255)

                    cGCCTACAGAAGAACCCACCCACAGCAGTTGGCGCTTGCGCCGAGCCACCCTGACGTTATCGCCGGCAAGGCCGGCCTCTTTGTAGATGACCCGCCCAAAGGCCTGTTGAATGAAGTAGGCCAACGCCTTGGCCTGCGGCTTGCCCTCaa                                                                                                                                 <  1:987387/150‑1 (MQ=255)

                        tACAGAAGAACCCACCCACAGCAGTTGGCGCTTGCGCCGAGCCACCCTGACGTTATCGCCGGCAAGGCCGGCCTCTTTGTAGATGACCCGCCCAAAGGCCTGTTGAATGAAGTAGGCCAACGCCTTGGCCTGCGGCTTGCCCTCAAGCaa                                                                                                                             >  2:1629324/1‑150 (MQ=255)

                           agaagaACCCACCCACAGCAGTTGGCGCTTGCGCCGAGCCACCCTGACGTTATCGCCGGCAAGGCCGGCCTCTTTGTAGATGACCCGCCCAAAGGCCTGTTGAATGAAGTAGGCCAACGCCTTGGCCTGCGGCTTGCCCTCAAGCAACgg                                                                                                                          >  2:2195534/1‑150 (MQ=255)

                            gaagaaCCCACCCACAGCAGTTGGCGCTTGCGCCGAGCCACCCTGACGTTATCGCCGGCAAGGCCGGCCTCTTTGTAGATGACCCGCCCAAAGGCCTGTTGAATGAAGTAGGCCAACGCCTTGGCCTGCGGCTTGCCCTCAAGCAACggg                                                                                                                         >  1:166516/1‑150 (MQ=255)

                             aagagCTCACTCACAGCAGTTGGCGCTTGCGCCGAGCCACCCTGACGTTATCGCCGGCAAGGCCGGCCTCTTTGTAGATGACCCGCCCAAAGGCCTGTTGAATGAAGTAGGCCAACGCCTTGGCCTGCGGCTTGCCCTCAAGCAACGGGa                                                                                                                        <  2:571830/150‑1 (MQ=255)

                             aagaaCCCACCCACAGCAGTTGGCGCTTGCGCCGAGCCACCCTGACGTTATCGCCGGCAAGGCCGGCCTCTTTGTAGATGACCCGCCCAAAGGCCTGTTGAATGAAGTAGGCCAACGCCTTGGCCTGc                                                                                                                                              <  2:272446/128‑1 (MQ=255)

                             aagaaCCCACCCACAGCAGTTGGCGCTTGCGCCGAGCCACCCTGACGTTATCGCCGGCAAGGCCGGCCTCTTTGTAGATGACCCGCCCAAAGGCCTGTTGAATGAAGTAGGCCAACGCCTTGGCCTGc                                                                                                                                              >  1:272446/1‑128 (MQ=255)

                             aagaaCCCACCCACAGCAGTTGGCGCTTGCGCCGAGCCACCCTGACGTTATCGCCGGCAAGGCCGGCCTCTTTGTAGATGACCCGCCCAAAGGCCTGTTGAATGAAGTAGGCCAACGCCTTGGCCTGc                                                                                                                                              >  1:272152/1‑128 (MQ=255)

                             aagaaCCCACCCACAGCAGTTGGCGCTTGCGCCGAGCCACCCTGACGTTATCGCCGGCAAGGCCGGCCTCTTTGTAGATGACCCGCCCAAAGGCCTGTTGAATGAAGTAGGCCAACGCCTTGGCCTGCGGCTTGCCCTCAAGCAACGGGa                                                                                                                        <  1:2195534/150‑1 (MQ=255)

                                aaCCCACCCACAGCAGTTGGCGCTTGCGCCGAGCCACCCTGACGTTATCGCCGGCAAGGCCGGCCTCTTTGTAGATGACCCGCCCAAAGGCCTGTTGAATGAAGTAGGCCAACGCCTTGGCCTGCGGCTTGCCCTCaa                                                                                                                                 >  2:386864/1‑138 (MQ=255)

                                   ccaccaacAGCAGTAAACGCTTGCGCCGAGCCACCCTGACGTTATCGCCGGCAAGGCCGGCCTCTTTGTAGATGACCCGCCCAAAGGCCTGTTGAATGAAGTAGGCCAACGCCTTGGCCTGCGGCTTGCCCTCAAGCACCGGGAGGCTgg                                                                                                                  >  2:713199/1‑150 (MQ=255)

                                     acccacAGCAGTTGGCGCTTGCGCCGAGCCACCCTGACGTTATCGCCGGCAAGGCCGGCCTCTTTGTAGATGACCCGCCCAAAGGCCTGTTGAATGAAGTAGGCCAACGCCTTGGCCTGCGGCTTGCCCTCAAGCAACGGGAgg                                                                                                                      >  2:1675482/1‑144 (MQ=255)

                                     acccacAGCAGTTGGCGCTTGCGCCGAGCCACCCTGACGTTATCGCCGGCAAGGCCGGCCTCTTTGTAGATGACCCGCCCAAAGGCCTGTTGAATGAAGTAGGCCAACGCCTTGGCCTGCGGCTTGCCCTCAAGCAACGGGAgg                                                                                                                      <  1:1675482/144‑1 (MQ=255)

                                       ccacAGCAGTTGGCGCTTGCGCCGAGCCACCCTGACGTTATCGCCGGCAAGGCCGGCCTCTTTGTAGATGACCCGCCCAAAGGCCTGTTGAATGAAGTAGGCCAACGCCTTGGCCTGCGGCTTGCCCTCAAGCAACGGGAGGc                                                                                                                     >  2:31648/1‑143 (MQ=255)

                                       ccacAGCAGTTGGCGCTTGCGCCGAGCCACCCTGACGTTATCGCCGGCAAGGCCGGCCTCTTTGTAGATGACCCGCCCAAAGGCCTGTTGAATGAAGTAGGCCAACGCCTTGGCCTGCGGCTTGCCCTCAAGCAACGGGAGGc                                                                                                                     <  1:31648/143‑1 (MQ=255)

                                             cagTTGGCGCTTGCGCCGAGCCACCCTGACGTTATCGCCGGCAAGGCCGGCCTCTTTGTAGATGACCCGCCCAAAGGCCTGTTGAATGAAGTAGGCCAACGCCTTGGCCTGCGGCTTGCCTCCAAGCAACGGGAGGCTGGTCTTGTAAgg                                                                                                        <  1:1270790/150‑1 (MQ=255)

                                             cagTTGGCCTTTGCGCCGAGCCACCCTGACGTTATCGCCGGCAAGGCCGGCCCCATTGTATGTGCAACCAGCGAAGGCCTGTTGAATGAAGTAGGCCAACGACTTGGCCTGCGGCTTCCCCTCAAGCAACGGGAGGCTGGTCTTGTAAgg                                                                                                        <  1:1376223/150‑1 (MQ=255)

                                                 tGGCGCTTGCGCCGAGCCACCCTGACGTTATCGCCGGCAAGGCCGGCCTCTTTGTAGATGACCCGCCCAAAGGCCTGTTGAATGAAGTAGGCCAACGCCTTGGCCTGCGGCTTGCCCTCAAGCAACGGGAGGCTGGTCTTGTAAGGCTgc                                                                                                    <  2:1774426/150‑1 (MQ=255)

                                                      cTTGCGCCGAGCCACCCTGACGTTATCGCCGGCAAGGCCGGCCTCTTTGTAGATGACCCGCCCAAAGGCCTGTTGAATGAAGTAGGCCAACGCCTTGGCCTGCGGCTTGCCCTCAAGCAACGGGAGGCTGGTCTTGTAAGGCTGcgccgc                                                                                               <  2:441733/150‑1 (MQ=255)

                                                       aagcgcCGAGCCACCCTGACGTTATCGCCGGCAAGGCCGGCCTCTTTGTAGATGACCCGCCCAAAGGCCTGTTGAATGAAGTAGGCCAACGCCTTGGCCTGCGGCTTGCCCTCAAGCAACGGGAGGCTGGTCTTGTAAGGCTGCgccgcc                                                                                              >  1:615096/3‑150 (MQ=255)

                                                          cgcCGAGCCACCCTGACGTTATCGCCGGCAAGGCCGGCCTCTTTGTAGATGACCCGCCCAAAGGCCTGTTGAATGAAGTAGGCCAACGCCTTGGCCTGCGGCTTGCCCTCAAGCAACGGGAGGCTGGTCTTGTAAGGCTGCGCCGCCTCa                                                                                           >  1:283867/1‑150 (MQ=255)

                                                                     ccTGACGTTATCGCCGGCAAGGCCGGCCTCTTTGTAGATGACCCGCCCAAAGGCCTGTTGAATGAAGTAGGCCAACGCCTTGGCCTGCGGCTTGCCCTCAAGCAACGGGAGGCTGGTCTTGTAAGGCTGCGCCGCCTCACGGACAAAg                                                                                  <  2:1043505/148‑1 (MQ=255)

                                                                     ccTGACGTTATCGCCGGCAAGGCCGGCCTCTTTGTAGATGACCCGCCCAAAGGCCTGTTGAATGAAGTAGGCCAACGCCTTGGCCTGCGGCTTGCCCTCAAGCAACGGGAGGCTGGTCTTGTAAGGCTGCGCCGCCTCACGGACAAAg                                                                                  >  1:1043505/1‑148 (MQ=255)

                                                                           gTTATCGCCGGCAAGGCCGGCCTCTTTGTAGATGACCCGCCCAAAGGCCTGTTGAATGAAGTAGGCCAACGCCTTGGCCTGCGGCTTGCCCTCAAGCAACGGTAGGCTGGTCTTGTAAGGCTGCGCCGCCTCACGGACAAAGGCGTTgcg                                                                          >  2:2185439/1‑150 (MQ=255)

                                                                             tATCGCCGGCAAGGCCGGCCTCTTTGTAGATGACCCGCCCAAAGGCCTGTTGAATGAAGTAGGCCAACGCCTTGGCCTGCGGCTTGCCCTCAAGCAACGGGAGGCTGGTCTTGTAAGGCTGCGCCGCCTCACGGACAAAGGCGTTGCGCa                                                                        <  2:888003/150‑1 (MQ=255)

                                                                                        aGGCCGGCCTCTTTGTAGATGACCCGCCCAAAGGCCTGTTGAATGAAGTAGGCCAACGCCTTGGCCTGCGGCTTGCCCTCAAGCAACGGGAGGCTGGTCTTGTAAGGCTGCGCCGCCTCACGGACAAAGGCGTTGCGCACGTCACCGTGc                                                             >  1:2310441/1‑150 (MQ=255)

                                                                                         ggccggccTCTTTGTAGATGACCCGCCCAAAGGCCTGTTGAATGAAGTAGGCCAACGCCTTGGCCTGCGGCTTGCCCTCAAGCAACGGGAGGCTGGTCTTGTAAGGCTGCGCCGCCTCACGGACAAAGGCGTTGCGCACGTCACCGTGCt                                                            <  2:1399786/150‑1 (MQ=255)

                                                                                          gccggccTCTTTGTAGATGACCCGCCCAAAGGCCTGTTGAATGAAGTAGGCCAACGCCTTGGCCTGCGGCTTGCCCTCAAGCAACGGGAGGCTGGTCTTGTaa                                                                                                          >  1:1710324/1‑103 (MQ=255)

                                                                                          gccggccTCTTTGTAGATGACCCGCCCAAAGGCCTGTTGAATGAAGTAGGCCAACGCCTTGGCCTGCGGCTTGCCCTCAAGCAACGGGAGGCTGGTCTTGTaa                                                                                                          <  2:1710324/103‑1 (MQ=255)

                                                                                           ccggccTCTTTGTAGATGACCCGCCCAAAGGCCTGTTGAATGAAGTAGGCCAACGCCTTGGCCTGCGGCTTGCCCTCAAGCAACGGGAGACTGGTCTTGTAAGGCTGCGCCGCCTCACGGACAAAGGCGTTGCGCACGTCACCGTGCTGt                                                          >  2:254589/1‑150 (MQ=255)

                                                                                            cggccTCTTTGTAGATGACCCGCCCAAAGGCCTGTTGAATGAAGTAGGCCAACGCCTTGGCCTGCGGCTTGCCCTCAAGCAACGGGAGGCTGGTCTTGTAAGGCTGCGCCGCCTCACGGACAAAGGCGTTGCGCACGTCACCGTGCTGTa                                                         >  1:1765738/1‑150 (MQ=255)

                                                                                               ccTCTTTGTAGATGACCCGCCCAAAGGCCTGTTGAATGAAGTAGGCCAACGCCTTGGCCTGCGGCTTGCCCTCAAGCAACGGGAGGCTGGTCTTGTAAGGCTGCGCCGCCTCACGGACAAAGGCGTTGCGCACGTCACCGTGCTGTACGa                                                      <  1:2390945/150‑1 (MQ=255)

                                                                                               ccTCTTTGTAGATGACCCGCCCAAAGGCCTGTTGAATGAAGTAGGCCAACGCCTTGGCCTGCGGCTTGCCCTCAAGCAACGGGAGGCTGGTCTTGTAAGGCTGCGCCGCCTCACGGACAAAGGCGTTGCGCACGTCACCGTGCTGTACGa                                                      <  2:2310441/150‑1 (MQ=255)

                                                                                                ctctTTGTAGATGACCCGCCCAAAGGCCTGTTGACTGAAGTAGGCCAACGCCTTGGCCTGCGGCTTGCCCTCAAGCAACGGGAGGCTGGTCTTGTAAGGCTGCGCCGCCTCACGGACAAAGGCGTTGCGCACGTCACCGTGCTTTACGAc                                                     >  2:808496/1‑150 (MQ=255)

                                                                                                ctctTTGTAGATGACCCGCCCAAAGGCCTGTTGAATGAAGTAGGCCAACGCCTTGGCCTGCGGCTTGCCCTCAAGCAACGGGAGGCTGGTCTTGTAAGGCTGCGCCGCCTCACGGACAAAGGCGTTGCGCACGTCACCGTGCTGTACGAc                                                     >  2:678746/1‑150 (MQ=255)

                                                                                                ctctTTGTAGATGACCCGCCCAAAGGCCTGTTGAATGAAGTAGGCCAACGCCTTGGCCTGCGGCTTGCCCTCAAGCAACGGGAGGCTGGTCTTGTAAGGCTGCGCCGCCTCACGGACAAAGGCGTTGCGCACGTCACCGTGCTGTACGAc                                                     >  1:144150/1‑150 (MQ=255)

                                                                                                   tttGTAGATGACCCGCCCAAAGGCCTGTTGAATGAAGTAGGCCAACGCCTTGGCCTGCGGCTTGCCCTCAAGCAACGGGAGGCTGGTCTTGTAAGGCTGCGCCGCCTCACGGACAAAGGCGTTGCGCACGTCACCGTGCTGTACGACTga                                                  >  1:2284865/1‑150 (MQ=255)

                                                                                                    ttGTAGATGACCCGCCCAAAGGCCTGTTGAATGAAGTAGGCCAACGCCTTGGCCTGCGGCTTGCCCTCAAGCAACGGGAGGCTGGTCTTGTAAGGCTGCGCCGCCTCACGGACAAAGGCGTTGCGCACGTCACCGTGCTGTACGACTgac                                                 <  2:166516/150‑1 (MQ=255)

                                                                                                      gTAGATGACCCGCCCAAAGGCCTGTTGAATGAAGTAGGCCAACGCCTTGGCCTGCGGCTTGCCCTCAAGCAACGTGAGGCTGGTCTTGTAAGGCTGCGCCGCCTCACGGACAAAGGCGTTGCGCACGTCACCGTGCTGTACGACTgacga                                               >  2:82450/1‑150 (MQ=255)

                                                                                                      gTAGATGACCCGCCCAAAGGCCTGTTGAATGAAGTAGGCCAACGCCTTGGCCTGCGGCTTGCCCTCAAGCAACGGGAGGCTGGTCTTGTAAGGCTGCGCCGCCTCACGGACAAAGGCGTTGCGCACGTCACCGTGCTGTACGACTgacga                                               <  2:1453888/150‑1 (MQ=255)

                                                                                                       tAGATGACCCGCCCAAAGGCCTGTTGAATGAAGTAGGCCAACGCCTTGGCCTGCGGCTTGCCCTCAAGCAACGGGAGGCTGGTCTTGTAAGGCTGCGCCGCCTCACGGACAAAGGCGTTGCTCACGTCACCGTGCTGTACGACTgacgac                                              <  1:1089607/150‑1 (MQ=255)

                                                                                                         gATGACCCGCCCAAAGGCCTGTTGAATGAAGTAGGCCAACGCCTTGGCCTGCGGCTTGCCCTCAAGCAACGGGAGGCTGGTCTTGTAAGGCTGCGCCGCCTCACGGACAAAGGCGTTGCGCACGTCACCGTGCTGTACGACTGACGACCa                                            >  2:471889/1‑150 (MQ=255)

                                                                                                         gATGACCCGCCCAAAGGCCTGTTGAATGAAGTAGGCCAACGCCTTGGCCTGCGGCTTGCCCTCAAGCAACGGGAGGCTGGTCTTGTAAGGCTGCGCCGCCTCACGGACAAAGGCGTTGCGCACGTCACCGTGCTGTACGACTGACGACCa                                            >  1:1042319/1‑150 (MQ=255)

                                                                                                          aTGACCCGCCCAAAGGCCTGTTGAATGAAGTAGGCCAACGCCTTGGCCTGCGGCTTGCCCTCAAGCAACGGGAGGCTGGTCTTGTAAGGCTGCGCCGCCTCACGGACAAAGCCGTTGCGCACGTCACCGTGCTGTACGACTGACGACCaa                                           <  1:1127544/150‑1 (MQ=255)

                                                                                                             aCCCGCCCAAAGGCCTGTTGAATGAAGTAGGCCAACGCCTTGGCCTGCGGCTTGCCCTCAAGCAACGGGAGGCTGGTCTTGTAAGGCTGCGCCGCCTCACgg                                                                                        <  1:2038754/102‑1 (MQ=255)

                                                                                                             aCCCGCCCAAAGGCCTGTTGAATGAAGTAGGCCAACGCCTTGGCCTGCGGCTTGCCCTCAAGCAACGGGAGGCTGGTCTTGTAAGGCTGCGCCGCCTCACgg                                                                                        >  2:2038754/1‑102 (MQ=255)

                                                                                                              cccGCCCAAAGGCCTGTTGAATGAAGTAGGCCAACGCCTTGGCCTGCGGCTTGCCCTCAAGCAACGGGAGGCTGGTCTTGTAAGGCTGCGCCGCCTCACGGACAAAGGCGTTGCGCACGTCACCGTGCTGTACGACTGACGACCaa                                           <  2:494711/146‑1 (MQ=255)

                                                                                                              cccGCCCAAAGGCCTGTTGAATGAAGTAGGCCAACGCCTTGGCCTGCGGCTTGCCCTCAAGCAACGGGAGGCTGGTCTTGTAAGGCTGCGCCGCCTCACGGACAAAGGCGTTGCGCACGTCACCGTGCTGTACGACTGACGACCaa                                           >  1:494711/1‑146 (MQ=255)

                                                                                                               ccGCCCAAAGGCCTGTTGAATGAAGTAGGCCAACGCCTTGGCCTGCGGCTTGCCCTCAAGCAACg                                                                                                                           >  2:2087503/1‑65 (MQ=255)

                                                                                                               ccGCCCAAAGGCCTGTTGAATGAAGTAGGCCAACGCCTTGGCCTGCGGCTTGCCCTCAAGCAACg                                                                                                                           <  1:2087503/65‑1 (MQ=255)

                                                                                                               ccGCCCAAAGGCCTGTTGAATGAAGTAGGCCAACGCCTTGGCCTGCGGCTTGCCCTCAAGCAACGGGAGGCTGGTCTTGTAAGGCTGCGCCGCCTCACGGACAAAGGCGTTGCGCACGTCACCGTGCTGTACGACTGACGACCAATacac                                      >  2:676621/1‑150 (MQ=255)

                                                                                                                cGCCCAAAGGCCTGTTGAATGAAGTAGGCCAACGCCTTGGCCTGCGGCTTGCCCTCAAGCAACGGGAGGCTGGTCTTGTAAGGCTGCGCCGCCTCACGGACAAAGGCGTTGCGCACGTCACCGTGCTGTa                                                         >  1:388600/1‑130 (MQ=255)

                                                                                                                cGCCCAAAGGCCTGTTGAATGAAGTAGGCCAACGCCTTGGCCTGCGGCTTGCCCTCAAGCAACGGGAGGCTGGTCTTGTAAGGCTGCGCCGCCTCACGGACAAAGGCGTTGCGCACGTCACCGTGCTGTa                                                         <  2:388600/130‑1 (MQ=255)

                                                                                                                 gCCCAAAGGCCTGTTGAATGAAGTAGGCCAACGCCTTGGCCTGCGGCTTGCCCTCAAGCAACGGGAGGCTGGTCTTGTAAGGCTGCGCCGCCTCACGGACAAAGGCGTTGCGCACGTCAcc                                                                 >  1:1543665/1‑121 (MQ=255)

                                                                                                                 gCCCAAAGGCCTGTTGAATGAAGTAGGCCAACGCCTTGGCCTGCGGCTTGCCCTCAAGCAACGGGAGGCTGGTCTTGTAAGGCTGCGCCGCCTCACGGACAAAGGCGTTGCGCACGTCACCGTGCTGTACGACTGACGACCAATacacac                                    >  2:1802246/1‑150 (MQ=255)

                                                                                                                       aGGCCTGTTGAATGAAGTAGGCCAACGCCTTGGCCTGCGGCTTGCCCTCAAACAACGGGAGGCTGGTCTTGTAAGGCTGCGCCGCCTCACGGACAAAGGCGTTGCGCACGTCACCGTGCTGTACGACTGACGACCAATACACACCGCGCa                              <  1:996943/150‑1 (MQ=255)

                                                                                                                         gCCTGTTGAATGAAGTAGGCCAACGCCTTGGCCTGCGGCTTGCCCTCAAGCAACGGGAGGCTGGTCTTGTAAGGCTGCGCCGCCTCACGGACAAAGGCGTTGCGCACGTCAc                                                                  <  1:2158896/112‑1 (MQ=255)

                                                                                                                         gCCTGTTGAATGAAGTAGGCCAACGCCTTGGCCTGCGGCTTGCCCTCAAGCAACGGGAGGCTGGTCTTGTAAGGCTGCGCCGCCTCACGGACAAAGGCGTTGCGCACGTCAc                                                                  >  2:2158896/1‑112 (MQ=255)

                                                                                                                                gaatgaaGTAGGCCAACGCCTTGGCCTGCGGCTTGCCCTCAAGCAACGGGAGGCTGGTCTTGTAAGGCTGCGCCGCCTCACGGACAAAGGCGTTGCGCACGTCACCGTGCTGTACGACTGACGACCAATACACACCGCGCACCAGCGCCg                     >  1:1972734/1‑150 (MQ=255)

                                                                                                                                gaatgaaGTAGGCCAACGCCTTGGCCTGCGGCTTGCCCTCAAGCAACGGGAGGCTGGTCTTGTAAGGCTGCGCCGCCTCACGGACAAAGGCGTTGCGCACGTCACCGTGCTGTACGACTGACGACCAATACACACCGCGCACCAGCGCCg                     <  2:663021/150‑1 (MQ=255)

                                                                                                                                   tgaaGTAGGCCAACGCCTTGGCCTGCGGCTTGCCCTCAAGCAACGGGAGGCTGGTCTTGTAAGGCTGCGCCGCCTCACGGACAAAGGCGTTGCGCACGTCACCGTGc                                                             >  1:207877/1‑107 (MQ=255)

                                                                                                                                   tgaaGTAGGCCAACGCCTTGGCCTGCGGCTTGCCCTCAAGCAACGGGAGGCTGGTCTTGTAAGGCTGCGCCGCCTCACGGACAAAGGCGTTGCGCACGTCACCGTGc                                                             <  2:207877/107‑1 (MQ=255)

                                                                                                                                   tgaaGTAGGCCAACGCCTTGGCCTGCGGCTTGCCCTCAAGCAACGGGAGGCTGGTCTTGTAAGGCTGCGCCGCCTCACGGACAAAGGCGTTGCGCACGTCACCGTGCTGTACGACTGACGACCAATACACACCGCGCACCAGCGCCgggg                  >  1:1061892/1‑150 (MQ=255)

                                                                                                                                     ccgcaGGCCAAGGCGTTGGCCTGCGGCTTGCCCTCAAGCAACGGGAGGCTGGTCTTGTAAGGCTGCGCCGCCTCACGGACAAAGGCGTTGCGCACGTCACCGTGCTGTACGACTGACGACCAATACACACCGCGCACCa                           <  2:2105557/135‑1 (MQ=255)

                                                                                                                                     ccgcaGGCCAAGGCGTTGGCCTGCGGCTTGCCCTCAAGCAACGGGAGGCTGGTCTTGTAAGGCTGCGCCGCCTCACGGACAAAGGCGTTGCGCACGTCACCGTGCTGTACGACTGACGACCAATACACACCGCGCACCa                           >  1:2105557/5‑139 (MQ=255)

                                                                                                                                     aaGTAGGCCAACGCCTTGGCCTGCGGCTTGCCCTCAAGCAACGGGAGGCTGGTCTTGTAAGGCTGCGCCGCCTCACGGACAAAGGCGTTGCGCACGTCACCGTGCTGTACGACTGACGACCAATACACACCGCGCACCAGCGCCGGggtg                >  1:799216/1‑150 (MQ=255)

                                                                                                                                     aaGTAGGCCAACGCCTTGGCCTGCGGCTTGCCCTCAAGCAACGGGAGGCTGGTCTTGTAAGGCTGCGCCGCCTCACGGACAAAGGCGTTGCGCACGTCACCGTGCTGTACGACTGACGACCAATACACACCGCGCACCAGCGCCGGggtg                <  1:471889/150‑1 (MQ=255)

                                                                                                                                       gTAGGCCAACGCCTTGGCCTGCGGCTTGCCCTCAAGCAACGGGAGGCTGGTCTTGTAAGGCTGCGCCGCCTCACGGACAAAGGCGTTGCGCACGTCACCGTGCTGTACGACTGACGACCAATACACACCGCGCACCAGCGCCGGggtggt              <  2:933422/150‑1 (MQ=255)

                                                                                                                                        tAGGCCAACGCCTTGGCCTGCGGCTTGCCCTCAAGCAACGGGAGGCTGGTCTTGTAAGGCTGCGCCGCCTCACGGACAAAGGCGTTGCGCACGTCACCGTGCTGTACGACTGACGACCAATACACACCGCGCACCAGCGCCGGGgtggtg             >  1:584473/1‑150 (MQ=255)

                                                                                                                                        caGGCCAAGGCGTTGGCCTGCGGCTTGCCCTCAAGCAACGGGAGGCTGGTCTTGTAAGGCTGCGCCGCCTCACGGACAAAGGCGTTGCGCACGTCACCGTGCTGTACGACTGACGACCAATACACACCGCGCACCAGCGCCGGggtg                <  1:2009872/146‑1 (MQ=255)

                                                                                                                                        caGGCCAAGGCGTTGGCCTGCGGCTTGCCCTCAAGCAACGGGAGGCTGGTCTTGTAAGGCTGCGCCGCCTCACGGACAAAGGCGTTGCGCACGTCACCGTGCTGTACGACTGACGACCAATACACACCGCGCACCAGCGCCGGggtg                >  2:2009872/2‑147 (MQ=255)

                                                                                                                                          ggCCAACGCCTTGGCCTGCGGCTTGCCCTCAAGCAACGGGAGGCTGGTCTTGTAAGGCTGCGCCGCCTCACGGACAAAGGCGTTGCGCACGTCACCGTGCTGTACGACTGACGACCAATACACACCGCGCACCAGCGCCGGGGTGGTGaa           >  2:695319/1‑150 (MQ=255)

                                                                                                                                          ggCCAACGCCTTGGCCTGCGGCTTGCCCTCAAGCAACGGGAGGCTGGTCTTGTAAGGCTGCGCCGCCTCACGGACAAAGGCGTTGCGCACGTCACCGTGCTGTACGACTGACGACCAATACACACCGCGCACCAGCACCGGGGTGGTGaa           >  2:2295482/1‑150 (MQ=255)

                                                                                                                                          ggCCAACGCCTTGGCCTGCGGCTTGACCTCAAGCAACGGGAGGCTGGTCTTGTAAGGCTGCGCCGCCTAACGGACAAAGACGTTGCGCACGTCACCGTGCTGTACGACTGGCGACAAAGACACGCCGCGCACCAGCGCCGGGGTGGGGaa           >  2:1395215/1‑150 (MQ=255)

                                                                                                                                           gCCAACGCCTTGGCCTGCGGCTTGCCCTCAAGCAACGGGAGGCTGGTCTTGTAAGGCTGCGCCGCCTCACGGACAAAGGCGTTGCGCACGTCACCGTGCTGTACGACTGACGACCAATACACACCGCGCACCAGCGCCGGGGTGGTGaaa          >  2:2129860/1‑150 (MQ=255)

                                                                                                                                             cAACGCCTTGGCCTGCGGCTTGCCCTCAAGCAACGGGAGGCTGGTCTTGTAAGGCTGCGCCGCCTCACGGACAAAGGCGTTGCGCACGTCACCGTGCt                                                            <  1:1676836/98‑1 (MQ=255)

                                                                                                                                             cAACGCCTTGGCCTGCGGCTTGCCCTCAAGCAACGGGAGGCTGGTCTTGTAAGGCTGCGCCGCCTCACGGACAAAGGCGTTGCGCACGTCACCGTGCt                                                            >  2:1676836/1‑98 (MQ=255)

                                                                                                                                             cAACGCCTTGGCCTGCGGCTTGCCCTCAAGCAACGGGAGGCTGGTCTTGTAAGGCTGCGCCGCCTCACGGACAAAGGCGTTGCGCACGTCACCGTGCTGTACGACTgacgac                                              <  1:2146424/112‑1 (MQ=255)

                                                                                                                                             cAACGCCTTGGCCTGCGGCTTGCCCTCAAGCAACGGGAGGCTGGTCTTGTAAGGCTGCGCCGCCTCACGGACAAAGGCGTTGCGCACGTCACCGTGCTGTACGACTgacgac                                              >  2:2146424/1‑112 (MQ=255)

                                                                                                                                             cAACGCCTTGGCCTGCGGCTTGCCCTCAAGCAACGGGAGGCTGGTCTTGTAAGGCTGCGCCGCCTCACGGACAAAGGCGTTGCGCACGTCACCGTGCTGTACGACTGACGACCAATACACACCGCGCACCAGCGCCGGGGTGGTGAAACg        >  1:943165/1‑150 (MQ=255)

                                                                                                                                              aaCGCCTTGGCCTGCGGCTTGCCCTCAAGCAACGGGAGGCTGGTCTTGTAAGGCTGCGCCGCCTCACGGACAAAGGCGTTGCGCACGTCACCGTGCTGTACGACTGACGACCAATACACACCGCGCACCAGCGCCGGGGTGGTGAAACgg       >  2:1804966/1‑150 (MQ=255)

                                                                                                                                                  ccTTGGCCTGCGGCTTGCCCTCAAGCAACGGGAGGCTGGTCTTGTAAGGCTGCGCCGCCTCACGGACAAAGGCGTTGCGCACGTCACCGTGCTGTACGACTGACGACCAATACACACCGCGCACCAGCGCCGGGGTGGTGAAACggcttc   <  1:695319/150‑5 (MQ=255)

                                                                                                                                                  ccTTGGCCTGCGGCTTGCCCTCAAGCAACGGGAGGCTGGTCTTGTAAGGCTGCGCCGCCTCACGGACAAAGGCGTTGCGCACGTCACCGTGCTGTACGACTGACGACCAATACACACCGCGCACCAGCGCCGGGGTGGTGAAACGGTCAc   >  1:311882/1‑150 (MQ=255)

                                                                                                                                                  ccTTGGCCTGCGGCTTGCCCTCAAGCAACGGGAGGCTGGTCTTGTAAGGCTGCGCCGCCTCACGGACAAAGGCGTTGCGCACGTCACCGTGCTGTACGACTGACGACCAATACACACCGCGCACCAGCACCGGGGTGGTGAAACGGTCAc   <  1:2295482/150‑1 (MQ=255)

                                                                                                                                                   cTTGGCCTGCGGCTTGCCCTCAAGCAACGGGAGGCTGGTCTTGTAAGGCTGCGCCGCCTCACGGACAAAGGCGTTGCGCACGTCACCGTGCTGTACGACTGACGACCAATACACACCGCGCACCAGCGCCGGGGTGGTGAAACGGTCACt  >  1:1889104/1‑150 (MQ=255)

GGCGACGCAGAAGGCCAGAACGCCTACAGAAGAACCCACCCACAGCAGTTGGCGCTTGCGCCGAGCCACCCTGACGTTATCGCCGGCAAGGCCGGCCTCTTTGTAGATGACCCGCCCAAAGGCCTGTTGAATGAAGTAGGCCAACGCCTTGGCCTGCGGCTTGCCCTCAAGCAACGGGAGGCTGGTCTTGTAAGGCTGCGCCGCCTCACGGACAAAGGCGTTGCGCACGTCACCGTGCTGTACGACTGACGACCAATACACACCGCGCACCAGCGCCGGGGTGGTGAAACGGTCACT  >  NC_002947/4598400‑4598696

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