BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,504,950	A→T	5.4%	S94C (AGC→TGC)	speB →	agmatinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,504,950	0	A	T	5.4%	662.2 / 17.1	198	S94C (AGC→TGC)	speB	agmatinase
Reads supporting (aligned to +/- strand): ref base A (86/101); new base T (7/4); total (93/105)
Fisher's exact test for biased strand distribution p-value = 3.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 199 total aligned reads displayed.

GGCCACCAGTAACCGCCCCGGCGCGCGCTTCGGCCCGCGGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCACGCCCGCCGCCATGGCCCGCTGGCG  >  NC_002947/2504802‑2505099

ggCCACCAGTAACCGCCCCGGCGCGCGCTTCGGCCCGCGGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAg                                                                                                                                                      >  2:3631558/1‑150 (MQ=255)

    accaGTAACCGCCCCGGCGCGCGCTTCGGCCCGCGGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGc                                                                                                                                                  >  1:4816007/1‑150 (MQ=255)

    accaGTAACCGCCCCGGCGCGCGCTTCGGCCCGCGGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGc                                                                                                                                                  >  1:438829/1‑150 (MQ=255)

     ccaGTAACCGCCCCGGCGCGCGCTTCGGCCCGCGGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCa                                                                                                                                                 <  2:3385157/150‑1 (MQ=255)

         tAACCGCCCCGGCGCGCGCTTCGGCCCGCGGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCAcgcc                                                                                                                                             <  1:5143620/150‑1 (MQ=255)

          aaCCGCCCCGGCGCGCGCTTCGGCCCGCGGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCAcgccg                                                                                                                                            <  2:5143640/150‑1 (MQ=255)

            ccGCCCCGGCGCGCGCTTCGGCCCGCGGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCa                                                                                                                                          <  1:1693792/150‑1 (MQ=255)

            ccGCCCCGGCGCGCGCTTCGGCCCGCGGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCa                                                                                                                                          >  2:3063664/1‑150 (MQ=255)

             cGCCCCGGCGCGCGCTTCGGCCCGCGGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAg                                                                                                                                         >  1:3008816/1‑150 (MQ=255)

              gCCCCGGCGCGCGCTTCGGCCCGCGGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACTGCGGCGTGCCGCTGt                                                                                                                                        <  2:4796950/150‑1 (MQ=255)

              gCCCCGGCGCGCGCTTCGGCCCGCGGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCg                                                                                                                                                    >  1:687013/1‑138 (MQ=255)

               ccccGGCGCGCGCTTCGGCCCGCGGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTc                                                                                                                                       <  2:3477704/150‑1 (MQ=255)

                cccGGCGCGCGCTTCGGCCCGCGGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACTGCGGCGTGCCGCTGTCg                                                                                                                                      <  2:779817/150‑1 (MQ=255)

                cccGGCGCGCGCTTCGGCCCGCGGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCg                                                                                                                                      >  2:4312744/1‑150 (MQ=255)

                    gcgcgcgcTTCGGCCCGCGGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTAc                                                                                                                                  <  1:4010111/150‑1 (MQ=255)

                     cgcgcgcTTCGGCCCGCGGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGG‑AACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCg                                                                                                                                <  2:2610875/150‑1 (MQ=255)

                      gcgcgcTTCGGCCCGCGGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCg                                                                                                                                <  2:4143812/150‑1 (MQ=255)

                        gcgcTTCGGCCCGCGGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCg                                                                                                                                                    <  1:647087/128‑1 (MQ=255)

                         cgcTTCGGCCCGCGGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGAc                                                                                                                             <  1:2796928/150‑1 (MQ=255)

                                gCCCGCGGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTTGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCg                                                                                                                      >  2:4024340/1‑150 (MQ=255)

                                gCCCGCGGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCg                                                                                                                      <  2:2976350/150‑1 (MQ=255)

                                    gcgGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGc                                                                                                                  <  2:446903/150‑1 (MQ=255)

                                      gggCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATc                                                                                                                       >  1:1108933/1‑143 (MQ=255)

                                      gggCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATc                                                                                                                       <  2:1108933/143‑1 (MQ=255)

                                      gggCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATc                                                                                                                       <  2:3844160/143‑1 (MQ=255)

                                        gCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCAcgccg                                                                                                                                            >  1:1806934/1‑120 (MQ=255)

                                           aTCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGcc                                                                                                           <  2:4756886/150‑1 (MQ=255)

                                            tCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCg                                                                                                          >  2:1446530/1‑150 (MQ=255)

                                             ccGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGa                                                                                                         <  1:2384465/150‑1 (MQ=255)

                                               gTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGc                                                                                                       >  1:2485260/1‑150 (MQ=255)

                                                 gctgcCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGcac                                                                                                     >  1:526525/1‑150 (MQ=255)

                                                 gctgcCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGcac                                                                                                     >  2:3203589/1‑150 (MQ=255)

                                                 gctgcCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGcac                                                                                                     >  2:3952444/1‑150 (MQ=255)

                                                     ccTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGc                                                                                                       <  1:409598/144‑1 (MQ=255)

                                                     ccTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGc                                                                                                       >  2:409598/1‑144 (MQ=255)

                                                     ccTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTc                                                                                                 <  2:1240087/150‑1 (MQ=255)

                                                     ccTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTc                                                                                                 <  2:894962/150‑1 (MQ=255)

                                                      cTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCt                                                                                                <  2:697065/150‑1 (MQ=255)

                                                       tCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTg                                                                                               <  2:731740/150‑1 (MQ=255)

                                                       tCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGGTCGACAACCTGTCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTg                                                                                               >  1:2210255/1‑150 (MQ=255)

                                                           tCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAgg                                                                                           >  1:369995/1‑150 (MQ=255)

                                                            ccAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGc                                                                                          <  1:3663844/150‑1 (MQ=255)

                                                             cagcagGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCg                                                                                         <  2:5005762/150‑1 (MQ=255)

                                                             cagcagGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCg                                                                                         <  2:4001946/150‑1 (MQ=255)

                                                              agcagGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCgg                                                                                        <  1:1256357/150‑1 (MQ=255)

                                                              agcagGCCTGGGCCCGCCACGGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCgg                                                                                        <  2:665378/150‑1 (MQ=255)

                                                                   gCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTgc                                                                                   <  1:1318594/150‑1 (MQ=255)

                                                                    ccTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTgcg                                                                                  >  1:5234389/1‑150 (MQ=255)

                                                                      tGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGcc                                                                                >  2:3138778/1‑150 (MQ=255)

                                                                          cccGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGc                                                                            >  1:784012/1‑150 (MQ=255)

                                                                             gCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGa                                                                         <  1:4659407/150‑1 (MQ=255)

                                                                             gCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGa                                                                         >  2:1993259/1‑150 (MQ=255)

                                                                                  tctccGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCTTCGTCTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCtt                                                                    <  2:2534944/147‑1 (MQ=255)

                                                                                  tGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCtt                                                                    >  1:595201/1‑150 (MQ=255)

                                                                                  tGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCtt                                                                    >  1:595302/1‑150 (MQ=255)

                                                                                   ggCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTtg                                                                   <  2:1862999/150‑1 (MQ=255)

                                                                                   ggCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTtg                                                                   <  1:5162411/150‑1 (MQ=255)

                                                                                      cGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTtggtg                                                                >  1:1385358/1‑150 (MQ=255)

                                                                                      cGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCAACGACTATGGCGACTGTCCCAACGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCAAGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTtgaga                                                                >  1:2534944/1‑147 (MQ=255)

                                                                                        tGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGc                                                              >  2:1941145/1‑150 (MQ=255)

                                                                                        tGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGc                                                              >  2:1097733/1‑150 (MQ=255)

                                                                                        tGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGc                                                              <  2:2038955/150‑1 (MQ=255)

                                                                                        tGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGc                                                              <  2:2056191/150‑1 (MQ=255)

                                                                                        tGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGc                                                              <  1:1719424/150‑1 (MQ=255)

                                                                                          ggCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGa                                                            >  2:2737972/1‑150 (MQ=255)

                                                                                             gttcgacccgttcgaccACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCa                                                         >  2:854872/1‑150 (MQ=255)

                                                                                             gttcgacccgttcgaccACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTGGGTGGCGACCa                                                         <  2:266275/150‑1 (MQ=255)

                                                                                             gttcgacccgttcgaccACCTGGCAGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCa                                                         <  1:4225615/150‑1 (MQ=255)

                                                                                               tcgacccgttcgaccACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACt                                                       <  2:1486810/150‑1 (MQ=255)

                                                                                                   cccgttcgaccACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTtcat                                                   <  2:4949950/150‑1 (MQ=255)

                                                                                                   cccgttcgaccACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTtcat                                                   <  1:1259569/150‑1 (MQ=255)

                                                                                                      gttcgaccACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAg                                                >  2:4398648/1‑150 (MQ=255)

                                                                                                      gttcgaccACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTGCATCAg                                                >  2:3271594/1‑150 (MQ=255)

                                                                                                       ttcgaccACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGc                                               >  2:170893/1‑150 (MQ=255)

                                                                                                          gaccACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTAc                                            <  1:3033364/150‑1 (MQ=255)

                                                                                                                cTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCgctg                                      <  2:2483993/150‑1 (MQ=255)

                                                                                                                     ggTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCaa                                 <  2:2822553/150‑1 (MQ=255)

                                                                                                                      gTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAg                                <  1:1977200/150‑1 (MQ=255)

                                                                                                                         aTCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCt                             >  1:26981/1‑150 (MQ=255)

                                                                                                                         aTCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCt                             >  1:4548289/1‑150 (MQ=255)

                                                                                                                          tCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTc                            >  1:2678855/1‑150 (MQ=255)

                                                                                                                           cGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCa                           <  2:3698625/150‑1 (MQ=255)

                                                                                                                            gACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCAc                          <  1:5398268/150‑1 (MQ=255)

                                                                                                                               tATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCt                                              >  2:2230104/1‑127 (MQ=255)

                                                                                                                                aTGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCAcgccc                      >  2:3547109/1‑150 (MQ=255)

                                                                                                                                  ggCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCt                                                                                                <  1:3831724/74‑1 (MQ=255)

                                                                                                                                    cGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCACGCccgccg                  >  1:325301/1‑150 (MQ=255)

                                                                                                                                     gACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCACGCCcgccgc                 >  2:4206641/1‑150 (MQ=255)

                                                                                                                                     gACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCACGCCcgccgc                 >  2:531217/1‑150 (MQ=255)

                                                                                                                                      aCTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCACGCCCgccgcc                <  1:5334672/150‑1 (MQ=255)

                                                                                                                                          tCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCa                                                         >  1:4846648/1‑105 (MQ=255)

                                                                                                                                           ggtaccGACTGCGGCGTGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCACGCCCGCCGCCATGGc           >  1:4467177/6‑150 (MQ=255)

                                                                                                                                            ccTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCACGCCCGCCGCCATGGcc          <  1:3570830/150‑1 (MQ=255)

                                                                                                                                            ccTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCACGCCCGCCGCCATGGcc          >  1:2887172/1‑150 (MQ=255)

                                                                                                                                            ccTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCACGCCCGCCGCCATGGcc          >  1:5179695/1‑150 (MQ=255)

                                                                                                                                             cTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCATGCCCGCCGCCATGGccc         >  1:3076659/1‑150 (MQ=255)

                                                                                                                                               ccGACTGCGGCGTGCCGCTGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCACGCCCGCCGCCATGGCCCGc       >  2:1043225/2‑150 (MQ=255)

                                                                                                                                                cGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCACGCCCGCCGCCATGGCCCGCt      <  1:1652798/150‑1 (MQ=255)

                                                                                                                                                 gACTGCGGCGTGCCGCTGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGa                                                                         >  1:2164976/1‑82 (MQ=255)

                                                                                                                                                  aCAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCGTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCACGCCCGCCGCCATGGCCCGCTgg    >  2:559274/1‑150 (MQ=255)

GGCCACCAGTAACCGCCCCGGCGCGCGCTTCGGCCCGCGGGCGATCCGTGCTGCCTCGGTCCAGCAGGCCTGGGCCCGCCACTGGCCGTGGGCGTTCGACCCGTTCGACCACCTGGCGGTCATCGACTATGGCGACTGTCCCTTCGACAGCGGCACGCCGCAGTCGGTACCGGACAGCATCGAGGCCCATGCCGAGCACATTCTGCAGGCGGGTTGCGCCATGCTGACCCTTGGTGGCGACCACTTCATCAGCTACCCGCTGCTCAAGGCTCACGCCCGCCGCCATGGCCCGCTGGCG  >  NC_002947/2504802‑2505099

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