BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,545,563	G→A	6.8%	G11E (GGG→GAG)	rnz →	ribonuclease Z

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,545,563	0	G	A	6.8%	446.8 / 14.8	146	G11E (GGG→GAG)	rnz	ribonuclease Z
Reads supporting (aligned to +/- strand): ref base G (70/66); new base A (6/4); total (76/70)
Fisher's exact test for biased strand distribution p-value = 7.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 152 total aligned reads displayed.

CAGCAAGCTCATCTGGCGACAGGATAACTAGCGGGCGGTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCTTTCGATCCGCGACCTGCGCGCGATCTTCATC  >  NC_002947/4545414‑4545711

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

cAGCAAGCTCATCTGGCGACAGGATAACTAGCGGGCGGTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCgg                                                                                                                                                      <  2:2719540/150‑1 (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)

             tGGCGACAGGATAACTAGCGGGCGGTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCGTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCGTGGCACCTCCGCCGGGGTGCCCACCAAg                                                                                                                                         >  1:895940/1‑150 (MQ=255)

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

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

                      gATAACTAGCGGGCGGTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAAt                                                                                                                                >  1:4876751/1‑150 (MQ=255)

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

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

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

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

                            tAGCGGGCGGTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAgc                                                                                                                          >  1:3729703/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:1748706/150‑1 (MQ=255)

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

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

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

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

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

                                            gcggGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCg                                                                                                          >  2:1223976/1‑150 (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                                                                                                 <  1:1873484/150‑1 (MQ=255)

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

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

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

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

                                                        ccGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCa                                                                                              <  1:2790437/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)

                                                                         cccGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTc                                                                             >  1:4742688/1‑150 (MQ=255)

                                                                            gTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGAt                                                                          <  1:4099469/150‑1 (MQ=255)

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

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

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

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

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

                                                                                     tGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGaa                                                                 <  2:1440102/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)

                                                                                                   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:1633866/1‑150 (MQ=255)

                                                                                                            gAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGca                                          <  2:1658374/150‑1 (MQ=255)

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

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

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

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

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

                                                                                                                        ggACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCTTTc                              <  1:510340/150‑1 (MQ=255)

                                                                                                                        ggACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCTTTc                              <  2:1205126/150‑1 (MQ=255)

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

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

                                                                                                                                tgTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCTTTCGATCcgcg                      >  1:2072065/1‑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:4375972/1‑150 (MQ=255)

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

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

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

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

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

                                                                                                                                           aCCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCTTTCGATCCGCGACCTGcgcgcg           <  2:1482614/150‑1 (MQ=255)

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

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

                                                                                                                                              tCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCTTTCGATCCGCGACCTGCGCGCGATc        <  1:4284115/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)

CAGCAAGCTCATCTGGCGACAGGATAACTAGCGGGCGGTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCTTTCGATCCGCGACCTGCGCGCGATCTTCATC  >  NC_002947/4545414‑4545711

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