BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,545,558	C→G	6.7%	S9S (TCC→TCG)	rnz →	ribonuclease Z

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,545,558	0	C	G	6.7%	490.0 / 19.4	148	S9S (TCC→TCG)	rnz	ribonuclease Z
Reads supporting (aligned to +/- strand): ref base C (69/69); new base G (6/4); total (75/73)
Fisher's exact test for biased strand distribution p-value = 7.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 152 total aligned reads displayed.

GACACCAGCAAGCTCATCTGGCGACAGGATAACTAGCGGGCGGTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCTTTCGATCCGCGACCTGCGCGCGATCTT  >  NC_002947/4545409‑4545707

gACACCAGCAAGCTCATCTGGCGACAGGATAACTAGCGGGCGGTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTcc                                                                                                                                                       <  1:3989839/150‑1 (MQ=255)

     cAGCAAGCTCATCTGGCGACAGGATAACTAGCGGGCGGTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCgg                                                                                                                                                  <  2:2719540/150‑1 (MQ=255)

     cAGCAAGCTCATCTGGCGACAGGATAACTAGCGGGCGGTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCgg                                                                                                                                                  >  2:2097833/1‑150 (MQ=255)

      aGCAAGCTCATCTGGCGACAGGATAACTAGCGGGCGGTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCggg                                                                                                                                                 <  2:2188301/150‑1 (MQ=255)

         aaGCTCATCTGGCGACAGGATAACTAGCGGGCGGTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTg                                                                                                                                              >  2:1868048/1‑150 (MQ=255)

             tCATCTGGCGACAGGATAACTAGCGGGCGGTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCCCGGCGGAGGTGccaa                                                                                                                                          <  2:4680256/150‑3 (MQ=255)

             tCATCTGGCGACAGGATAACTAGCGGGCGGTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCCCGGCGGAGGTGccaa                                                                                                                                          <  2:431463/150‑3 (MQ=255)

                  tGGCGACAGGATAACTAGCGGGCGGTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAg                                                                                                                                     >  1:4955006/1‑150 (MQ=255)

                     cGACAGGATAACTAGCGGGCGGTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCa                                                                                                                                  <  2:1403742/150‑1 (MQ=255)

                      gACAGGATAACTAGCGGGCGGTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCAc                                                                                                                                 >  2:4043751/1‑150 (MQ=255)

                      gACAGGATAACTAGCGGGCGGTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCAc                                                                                                                                 <  1:1956856/150‑1 (MQ=255)

                      gACAGGATAACTAGCGGGCGGTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCAc                                                                                                                                 >  1:4154233/1‑150 (MQ=255)

                              aaCTAGCGGGCGGTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTc                                                                                                                         >  1:2449143/1‑150 (MQ=255)

                              aaCTAGCGGGCGGTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTc                                                                                                                         >  1:1073073/1‑150 (MQ=255)

                               aCTAGCGGGCGGTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCa                                                                                                                        >  2:3888402/1‑150 (MQ=255)

                                 tAGCGGGCGGTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAgc                                                                                                                      >  1:3729703/1‑150 (MQ=255)

                                 tAGCGGGCGGTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAgc                                                                                                                      >  2:3229064/1‑150 (MQ=255)

                                 tAGCGGGCGGTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAgc                                                                                                                      <  1:3389329/150‑1 (MQ=255)

                                     gggcggTCAACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCca                                                                                                                  <  2:1779585/150‑1 (MQ=255)

                                            cgccggcggGTGGCAATCCGGCTACGGTGTGCGCCCCGTATCACCTTGGCCGCGTCTGGCCTGACACCGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGt                                                                                                           <  2:4409938/147‑1 (MQ=255)

                                              aCGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGa                                                                                                         <  1:1617531/150‑1 (MQ=255)

                                              aCGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGa                                                                                                         <  1:1748706/150‑1 (MQ=255)

                                               cggcggGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGAt                                                                                                        <  2:772679/150‑1 (MQ=255)

                                               cggcggGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGAt                                                                                                        <  2:1618083/150‑1 (MQ=255)

                                               cggcggGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGAt                                                                                                        <  2:2128385/150‑1 (MQ=255)

                                                ggcggGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATc                                                                                                       <  1:3993308/150‑1 (MQ=255)

                                                   gggTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGaa                                                                                                    <  2:309967/150‑1 (MQ=255)

                                                   gggTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGaa                                                                                                    <  1:3634177/150‑1 (MQ=255)

                                                    ggTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAg                                                                                                   <  2:4412041/150‑1 (MQ=255)

                                                       ggCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCa                                                                                                <  2:3337003/150‑1 (MQ=255)

                                                         cAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGc                                                                                              <  2:2043511/150‑1 (MQ=255)

                                                          aaTCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCACTGTCAGCGCCACAGCCGTGATCGAAGCCAGCg                                                                                             >  2:1235854/1‑150 (MQ=255)

                                                          aaTCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCg                                                                                             >  2:1235745/1‑150 (MQ=255)

                                                          aaTCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCg                                                                                             <  1:1873484/150‑1 (MQ=255)

                                                          aaTCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCg                                                                                             >  1:3486857/1‑150 (MQ=255)

                                                          aaTCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCg                                                                                             >  1:2883715/1‑150 (MQ=255)

                                                           aTCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCgg                                                                                            >  1:4686461/1‑150 (MQ=255)

                                                            tCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGc                                                                                           >  1:682608/1‑150 (MQ=255)

                                                             ccGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCa                                                                                          <  2:3309291/150‑1 (MQ=255)

                                                             ccGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCa                                                                                          <  1:2790437/150‑1 (MQ=255)

                                                             ccGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCa                                                                                          <  2:1539770/150‑1 (MQ=255)

                                                               ggCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGc                                                                                        <  1:5052556/150‑1 (MQ=255)

                                                                gCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGcc                                                                                       <  1:2212827/150‑1 (MQ=255)

                                                                gCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGcc                                                                                       >  1:2637586/1‑150 (MQ=255)

                                                                gCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGcc                                                                                       >  2:1885188/1‑150 (MQ=255)

                                                                  tAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCAc                                                                                     >  2:882520/1‑150 (MQ=255)

                                                                      gtgtGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGt                                                                                 >  2:1652850/1‑150 (MQ=255)

                                                                         tgCGCCACGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATc                                                                              <  1:3904335/150‑1 (MQ=255)

                                                                           cgcCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTg                                                                            <  2:1847391/150‑1 (MQ=255)

                                                                           cgcCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTg                                                                            <  1:4687177/150‑1 (MQ=255)

                                                                              cccGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTc                                                                         <  1:4870337/150‑1 (MQ=255)

                                                                                     cACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTgcg                                                                  <  2:4047429/150‑1 (MQ=255)

                                                                                         ttGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGa                                                              <  1:1664707/150‑1 (MQ=255)

                                                                                         ttGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGa                                                              <  2:1734382/150‑1 (MQ=255)

                                                                                         ttGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGa                                                              <  2:1734603/150‑1 (MQ=255)

                                                                                          tGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGaa                                                             <  2:410921/150‑1 (MQ=255)

                                                                                          tGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGaa                                                             <  2:1440102/150‑1 (MQ=255)

                                                                                          tGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGaa                                                             <  1:4375152/150‑1 (MQ=255)

                                                                                           ggCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAg                                                            <  2:4479322/150‑1 (MQ=255)

                                                                                           ggCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAg                                                            <  2:4479728/150‑1 (MQ=255)

                                                                                                cTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCAcc                                                       <  2:241896/150‑1 (MQ=255)

                                                                                                 ttCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCAccc                                                      >  2:2670525/1‑150 (MQ=255)

                                                                                                      gCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCAc                                                 >  1:3199196/1‑150 (MQ=255)

                                                                                                        cTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAATGCACCCAGCACCg                                               >  1:4835795/1‑150 (MQ=255)

                                                                                                        cTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCg                                               >  2:4329041/1‑150 (MQ=255)

                                                                                                         tGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCCCGGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGc                                              >  1:1049671/1‑150 (MQ=255)

                                                                                                         tGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGc                                              >  2:3549365/1‑150 (MQ=255)

                                                                                                         tGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGc                                              >  1:355147/1‑150 (MQ=255)

                                                                                                         tGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGc                                              <  2:4874615/150‑1 (MQ=255)

                                                                                                           aaaCAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCct                                            >  1:1658374/1‑150 (MQ=255)

                                                                                                             aCAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCctgc                                          <  2:2352892/150‑1 (MQ=255)

                                                                                                             aCAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCctgc                                          >  2:2586598/1‑150 (MQ=255)

                                                                                                               aGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCctgctg                                        >  1:2438902/1‑150 (MQ=255)

                                                                                                                ggAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCtgctgc                                       >  1:2645227/1‑150 (MQ=255)

                                                                                                                ggAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCtgctgc                                       >  2:1634018/1‑150 (MQ=255)

                                                                                                                      tGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGAtt                                                                     >  1:30331/1‑114 (MQ=255)

                                                                                                                       gAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACAcccc                                >  2:422396/1‑150 (MQ=255)

                                                                                                                       gAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACAcccc                                >  1:4381798/1‑150 (MQ=255)

                                                                                                                         gCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGc                              >  1:1896801/1‑150 (MQ=255)

                                                                                                                                    ctgTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCTTTCGATCcgc                   <  1:3901908/150‑1 (MQ=255)

                                                                                                                                    ctgTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCTTTCGATCcgc                   >  2:3535802/1‑150 (MQ=255)

                                                                                                                                     tgTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCTTTCGATCcgcg                  >  1:2072065/1‑150 (MQ=255)

                                                                                                                                      cttggtgggCACCCCGGCGGAGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCTTTCGATCCGCGa                 >  2:1219460/8‑150 (MQ=255)

                                                                                                                                      cttggtgggCACCCCGGCGGAGGTGCCAACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCTTTCGATCCGCGa                 >  2:2979606/8‑150 (MQ=255)

                                                                                                                                          gtgggCACCCCGGCGGAGGTGCCAACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCTTTCGATCCGCGACCTg             >  1:1408352/4‑150 (MQ=255)

                                                                                                                                           ttGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCTTTCGATCCGCGACCTgc            >  1:3686819/1‑150 (MQ=255)

                                                                                                                                           ttGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCTTTCGATCCGCGACCTgc            >  1:121924/1‑150 (MQ=255)

                                                                                                                                           ttGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCTTTCGATCCGCGACCTgc            <  2:1181436/150‑1 (MQ=255)

                                                                                                                                           ttGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCTTTCGATCCGCGACCTgc            >  2:4375972/1‑150 (MQ=255)

                                                                                                                                           tgggCACCCCGGCGGAGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCctgct                                         >  1:769708/3‑121 (MQ=255)

                                                                                                                                           tgggCACCCCGGCGGAGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCctgct                                         <  2:769708/119‑1 (MQ=255)

                                                                                                                                           tgggCACCCCGGCGGAGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCTTTCGATCCGCGACCTgc            >  1:3478112/3‑150 (MQ=255)

                                                                                                                                             ggCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCTTTCGATCCGCGACCTgcgc          >  2:3231079/1‑150 (MQ=255)

                                                                                                                                                aCCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCTTTCGATCCGCGACCTGcgcgcg       >  2:1157258/1‑150 (MQ=255)

                                                                                                                                                 ccTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCTTTCGATCCGCGACCTGCGCGCGa      >  1:2950430/1‑150 (MQ=255)

                                                                                                                                                  cTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCTTTCGATCCGCGACCTGCGCGCGAt     <  2:4876751/150‑1 (MQ=255)

                                                                                                                                                   tCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCTTTCGATCCGCGACCTGCGCGCGATc    <  1:4283962/150‑1 (MQ=255)

                                                                                                                                                   tCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCTGTCGATCCGCGACCTGCGCGCGATc    >  2:3790701/1‑150 (MQ=255)

                                                                                                                                                     cgccgGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCTTTCGATCCGCGACCTGCGCGCGATCtt  >  1:2702107/1‑150 (MQ=255)

                                                                                                                                                     cgccgGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCTTTCGATCCGCGACCTGCGCGCGATCtt  >  1:2273949/1‑150 (MQ=255)

GACACCAGCAAGCTCATCTGGCGACAGGATAACTAGCGGGCGGTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCTTTCGATCCGCGACCTGCGCGCGATCTT  >  NC_002947/4545409‑4545707

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