BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,545,561	C→G	5.4%	A10A (GCC→GCG)	rnz →	ribonuclease Z

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,545,561	0	C	G	5.4%	828.5 / 27.8	240	A10A (GCC→GCG)	rnz	ribonuclease Z
Reads supporting (aligned to +/- strand): ref base C (106/121); new base G (7/6); total (113/127)
Fisher's exact test for biased strand distribution p-value = 7.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01

Only 100 of 246 total aligned reads displayed.

ACCAGCAAGCTCATCTGGCGACAGGATAACTAGCGGGCGGTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCTTTCGATCCGCGACCTGCGCGCGATCTTC  >  NC_002947/4545412‑4545708

acCAGCAAGCTCATCTGGCGACAGGATAACTAGCGGGCGGTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTccgcc                                                                                                                                                     >  2:4148119/1‑150 (MQ=255)

     cAAGCTCATCTGGCGACAGGATAACTAGCGGGCGGTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGt                                                                                                                                                <  1:4134215/150‑1 (MQ=255)

     cAAGCTCATCTGGCGACAGGATAACTAGCGGGCGGTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGt                                                                                                                                                >  1:3753793/1‑150 (MQ=255)

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

        gCTCATCTGGCGACAGGATAACTAGCGGGCGGTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGcc                                                                                                                                             >  1:1539483/1‑150 (MQ=255)

         cTCATCTGGCGACAGGATAACTAGCGGGCGGTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGccc                                                                                                                                            <  2:681447/150‑1 (MQ=255)

         cTCATCTGGCGACAGGATAACTAGCGGGCGGTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGccc                                                                                                                                            <  2:4433747/150‑1 (MQ=255)

          tCATCTGGCGACAGGATAACTAGCGGGCGGTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCCCGGCGGAGGTGccaa                                                                                                                                           <  1:766019/150‑3 (MQ=255)

          tCATCTGGCGACAGGATAACTAGCGGGCGGTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCCCGGCGGAGGTGccaa                                                                                                                                           <  1:748183/150‑3 (MQ=255)

           cATCTGGCGACAGGATAACTAGCGGGCGGTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGccgag                                                                                                                                          <  1:2639515/150‑4 (MQ=255)

           cATCTGGCGACAGGATAACTAGCGGGCGGTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCccac                                                                                                                                          >  1:2886996/1‑150 (MQ=255)

            aTCTGGCGACAGGATAACTAGCGGGCGGTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCccacc                                                                                                                                         <  1:947194/150‑1 (MQ=255)

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

                        gATAACTAGCGGGCGGTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAAt                                                                                                                             <  2:1035768/150‑1 (MQ=255)

                          tAACTAGCGGGCGGTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGt                                                                                                                           <  2:3486024/150‑1 (MQ=255)

                              tAGCGGGCGGTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCACTGTCAgc                                                                                                                       >  1:2345495/1‑150 (MQ=255)

                               aGCGGGCGGTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAgcg                                                                                                                      >  2:3028739/1‑150 (MQ=255)

                                     cggTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAg                                                                                                                >  1:2482386/1‑150 (MQ=255)

                                     cggTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAg                                                                                                                >  2:4357064/1‑150 (MQ=255)

                                      ggTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGc                                                                                                               >  2:1078250/1‑150 (MQ=255)

                                        tCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCg                                                                                                             >  1:2772710/1‑150 (MQ=255)

                                          gACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTg                                                                                                           >  1:3223706/1‑150 (MQ=255)

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

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

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

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

                                            cggcggGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGAt                                                                                                         >  2:3739318/1‑150 (MQ=255)

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

                                            cggcggGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGAt                                                                                                         >  1:1825138/1‑150 (MQ=255)

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

                                             ggcggGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATc                                                                                                        >  1:3857829/1‑150 (MQ=255)

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

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

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

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

                                                gggTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCCTCTGGCCAGAACCAGGAGCCTGAGCATGGACCAGCTGTTCCTTGGCACCTCCGCCGCGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGaa                                                                                                     <  1:1951905/150‑1 (MQ=255)

                                                  gTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGc                                                                                                   >  2:1411356/1‑150 (MQ=255)

                                                     gCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAg                                                                                                <  2:1114623/150‑1 (MQ=255)

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

                                                        aTCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGGTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCgg                                                                                             <  2:3090164/150‑1 (MQ=255)

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

                                                           cGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAg                                                                                          >  1:2439494/1‑150 (MQ=255)

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

                                                                  ggTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTgg                                                                                   <  1:3826273/150‑1 (MQ=255)

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

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

                                                                    tgtgCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTa                                                                                 >  2:1403417/1‑150 (MQ=255)

                                                                      tgCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCCAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTg                                                                                    <  2:1628743/145‑1 (MQ=255)

                                                                          ccccGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGt                                                                           <  1:1881473/150‑1 (MQ=255)

                                                                                 tCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTgc                                                                    >  2:1581272/1‑150 (MQ=255)

                                                                                   aCCTTGGCCGCTTCTGGCCTGAAACCGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTgcgg                                                                  <  2:591785/150‑1 (MQ=255)

                                                                                   aCCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTgcgg                                                                  <  1:4069512/150‑1 (MQ=255)

                                                                                   aCCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTgcgc                                                                  <  2:1727923/150‑2 (MQ=255)

                                                                                     cTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTgcggcg                                                                >  1:3991500/1‑150 (MQ=255)

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

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

                                                                                        ggCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTCTCTGGTCGATTGCGGCGAAg                                                             >  2:932643/1‑150 (MQ=255)

                                                                                         gCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAgg                                                            <  2:4500697/150‑1 (MQ=255)

                                                                                          ccGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGc                                                           >  1:876926/1‑150 (MQ=255)

                                                                                            gCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCAc                                                         >  1:1373263/1‑150 (MQ=255)

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

                                                                                                cTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAg                                                     >  2:3348118/1‑150 (MQ=255)

                                                                                                 tGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGc                                                    <  1:3069264/150‑1 (MQ=255)

                                                                                                  ggCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCa                                                   <  2:1519246/150‑1 (MQ=255)

                                                                                                                gCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGcaca                                     <  1:2808167/150‑1 (MQ=255)

                                                                                                                   tGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACAccc                                  >  2:3047973/1‑150 (MQ=255)

                                                                                                                   tGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACAccc                                  >  2:3112304/1‑150 (MQ=255)

                                                                                                                   tGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACAccc                                  >  2:137729/1‑150 (MQ=255)

                                                                                                                    gAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCtgctgcctgcgcc                                 >  1:4366651/1‑143 (MQ=255)

                                                                                                                       cATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCt                              <  1:3552786/150‑1 (MQ=255)

                                                                                                                         tGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCttt                            >  1:3992419/1‑150 (MQ=255)

                                                                                                                         tGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCttt                            <  1:2140448/150‑1 (MQ=255)

                                                                                                                         tGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCttt                            <  1:2549301/150‑1 (MQ=255)

                                                                                                                         tGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCttt                            <  1:2588108/150‑1 (MQ=255)

                                                                                                                          ggACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCTTTc                           >  2:3284661/1‑150 (MQ=255)

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

                                                                                                                          ggACCGGCTGTTCCGTGGCACCTGCGCCGGGGTGCGCACCAAGGCACGCAATGTCTGCGCGACAGCCGTGATCGAAGCCAGCGGCAGCCTCTGGTATCTGGTCGATTGCGGCGTAGGCACCCTGCTCCGCCTGCTGCACACCCCGCTTTc                           <  2:4376334/150‑1 (MQ=255)

                                                                                                                              ctgctgTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCTTTCGATc                       <  1:1898119/150‑1 (MQ=255)

                                                                                                                                gctgTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCTTTCGATCcg                     >  1:1672051/1‑150 (MQ=255)

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

                                                                                                                                  tgTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCTGTCGATCcgcg                   >  2:1617579/1‑150 (MQ=255)

                                                                                                                                    ttggtgggCACCCCGGCGGAGGTG‑CCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCTTTCGATCCGCGAcc                >  2:871189/7‑150 (MQ=255)

                                                                                                                                    ttCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCTTTCGATCCGCGAc                 >  2:3988445/1‑150 (MQ=255)

                                                                                                                                       cTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCTTTCGATCCGCGACCTg              >  1:3132636/1‑150 (MQ=255)

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

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

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

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

                                                                                                                                        tgggCACCCCGGCGGAGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCCGCACCGCCTGCTGCACACCCCGCTTTCGATCCGCGACCTgc             >  2:2450901/3‑150 (MQ=255)

                                                                                                                                        tgggCACCCCGGCGGAGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCAATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCTTTCGATCCGCGACCTgc             >  2:2200503/3‑150 (MQ=255)

                                                                                                                                            cACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCTTTCGATCCGCGACCTgcgcgc         <  1:3838623/150‑1 (MQ=255)

                                                                                                                                            cACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCTTTCGATCCGCGACCTgcgcgc         <  1:4406440/150‑1 (MQ=255)

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

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

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

                                                                                                                                               cTCCGACGGGGTCCCAACAAAGGCACGCAATGTCAGCGCCACAGACTTGAACGAAGCCAGCGGCAGCCACTGGTATCTGGACGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGAACACCCCGCTATCGATCCGCGACCTGCGCGCGAt      >  2:3063734/1‑150 (MQ=255)

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

                                                                                                                                                tCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCTTTCGATCCGCGACCTGCGCGCGATc     >  1:2558283/1‑150 (MQ=255)

                                                                                                                                                 ccgccgGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCTTTCGATCCGCGACCTGCGCGCGATCt    <  2:3917116/150‑1 (MQ=255)

                                                                                                                                                   gccgGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCTTTCGATCCGCGACCTGCGCGCGATCTtc  <  2:4099963/150‑1 (MQ=255)

ACCAGCAAGCTCATCTGGCGACAGGATAACTAGCGGGCGGTCGACGGCGGGTGGCAATCCGGCTAAGGTGTGCGCCCCGTATCACCTTGGCCGCTTCTGGCCTGAAACAGGAGCCTGAGCATGGACCTGCTGTTCCTTGGCACCTCCGCCGGGGTGCCCACCAAGGCACGCAATGTCAGCGCCACAGCCGTGATCGAAGCCAGCGGCAGCCACTGGTATCTGGTCGATTGCGGCGAAGGCACCCAGCACCGCCTGCTGCACACCCCGCTTTCGATCCGCGACCTGCGCGCGATCTTC  >  NC_002947/4545412‑4545708

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