BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,005,760	C→G	5.2%	T76T (ACG→ACC)	PP_2624 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,005,760	0	C	G	5.2%	1262.5 / 42.0	366	T76T (ACG→ACC)	PP_2624	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (143/204); new base G (10/9); total (153/213)
Fisher's exact test for biased strand distribution p-value = 3.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 376 total aligned reads displayed.

AGGGCGACCAGTGCTTCGCGCAGCGCCAGCAGCTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCACCTCAGC  >  NC_002947/3005611‑3005909

aGGGCGACCAGTGCTTCGCGCAGCGCCAGCAGCTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTAcc                                                                                                                                                       <  2:4953439/150‑1 (MQ=255)

        cAGTGCTTCGCGCAGCGCCAGCAGCTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTc                                                                                                                                               <  1:974940/150‑1 (MQ=255)

           tGCTTCGCGCAGCGCCAGCAGCTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGCTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAg                                                                                                                                            <  1:5535800/150‑1 (MQ=255)

            gCTTCGCGCAGCGCCAGCAGCTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAgg                                                                                                                                           <  1:1949964/150‑1 (MQ=255)

             cTTCGCGCAGCGCCAGCAGCTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGc                                                                                                                                          >  1:32643/1‑150 (MQ=255)

              ttCGCGCAGCGCCAGCAGCTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCt                                                                                                                                         >  2:334162/1‑150 (MQ=255)

              ttCGCGCAGCGCCAGCAGCTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCt                                                                                                                                         <  1:2504368/150‑1 (MQ=255)

               tCGCGCAGCGCCAGCAGCTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGTGACGGTACCGTCCAGGCTc                                                                                                                                        <  1:3827438/150‑1 (MQ=38)

               tCGCGCAGCGCCAGCAGCTTGTTCAGCCCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGTGACGGTACCGTCCAGGCTc                                                                                                                                        <  1:45261/150‑1 (MQ=255)

                 gcgcAGCGCCAGCAGCTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTcac                                                                                                                                      >  1:6991609/1‑150 (MQ=255)

                      gcgcCAGCAGCTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTgg                                                                                                                                 <  1:737354/150‑1 (MQ=255)

                       cgcCAGCAGCTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGc                                                                                                                                <  1:1486813/150‑1 (MQ=255)

                         ccAGCAGCTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGcgt                                                                                                                              <  2:3161313/150‑1 (MQ=255)

                         ccAGCAGCTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGcgt                                                                                                                              >  2:7026714/1‑150 (MQ=255)

                          cagcagCTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGcgtc                                                                                                                             <  2:1347471/150‑1 (MQ=255)

                          cagcagCTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGcgtc                                                                                                                             <  2:193200/150‑1 (MQ=255)

                          cagcagCTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGcgtc                                                                                                                             <  2:817059/150‑1 (MQ=255)

                           agcagcTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGcgtcg                                                                                                                            <  2:3529378/150‑1 (MQ=255)

                            gcagcTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGcgtcgt                                                                                                                           <  1:4448026/150‑1 (MQ=255)

                             cagcTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGtt                                                                                                                          <  1:3922439/150‑1 (MQ=255)

                             cagcTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGtt                                                                                                                          <  1:2866489/150‑1 (MQ=255)

                              agcTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTc                                                                                                                         <  1:1046322/150‑1 (MQ=255)

                                       aGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTg                                                                                                                >  2:6663077/1‑150 (MQ=255)

                                        gCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGt                                                                                                               >  1:1767605/1‑150 (MQ=255)

                                        gCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGt                                                                                                               >  2:4195359/1‑150 (MQ=255)

                                         cTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGtt                                                                                                              >  2:1144419/1‑150 (MQ=255)

                                         cTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGtt                                                                                                              >  2:1089160/1‑150 (MQ=255)

                                         cTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGtt                                                                                                              <  2:1961142/150‑1 (MQ=255)

                                         cTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGtt                                                                                                              <  1:7227320/150‑1 (MQ=255)

                                            gggCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGaa                                                                                                           <  1:4590909/150‑1 (MQ=255)

                                            gggCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGaa                                                                                                           >  1:2276760/1‑150 (MQ=255)

                                            gggCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGaa                                                                                                           <  2:7348995/150‑1 (MQ=255)

                                              gCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAatg                                                                                                         >  2:1459051/1‑150 (MQ=255)

                                              gCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAatg                                                                                                         >  2:2926582/1‑150 (MQ=255)

                                                aCCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAatgta                                                                                                       <  1:6673720/150‑1 (MQ=255)

                                                 ccTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAatgtat                                                                                                      <  1:1167875/150‑1 (MQ=255)

                                                 ccTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAatgtat                                                                                                      <  2:2543785/150‑1 (MQ=255)

                                                 ccTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAatgtat                                                                                                      <  2:2290708/150‑1 (MQ=255)

                                                  cTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAatgtctg                                                                                                     >  2:6059772/1‑150 (MQ=255)

                                                  cTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAatgtatg                                                                                                     >  2:4461851/1‑150 (MQ=255)

                                                  cTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTGGAatgtatg                                                                                                     >  1:3022283/1‑150 (MQ=255)

                                                           gATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGAt                                                                                            >  2:7215139/1‑150 (MQ=255)

                                                             tccggtccggCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATgc                                                                                          <  1:2062230/150‑1 (MQ=255)

                                                             tccggtccggCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATgc                                                                                          <  1:5875069/150‑1 (MQ=255)

                                                             tccggtccggCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATgc                                                                                          <  1:2123461/150‑1 (MQ=255)

                                                                    cggCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGAc                                                                                   <  1:4885344/150‑1 (MQ=255)

                                                                     ggCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGAcc                                                                                  <  2:1642735/150‑1 (MQ=255)

                                                                      gCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGACTGTATGTTTATCGATGCGCATGACCt                                                                                 >  2:5510495/1‑150 (MQ=255)

                                                                        ccGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGa                                                                               >  2:4796296/1‑150 (MQ=255)

                                                                        ccGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTGGAATGTCTGTTTATCGATGCGCATGACCTGa                                                                               >  2:5154183/1‑150 (MQ=255)

                                                                         cGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGAc                                                                              >  2:4756420/1‑150 (MQ=255)

                                                                          gAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACg                                                                             >  1:1655363/1‑150 (MQ=255)

                                                                             aTCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTc                                                                          <  1:2979452/150‑1 (MQ=255)

                                                                             aTCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTc                                                                          <  2:5570928/150‑1 (MQ=255)

                                                                              tCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTcc                                                                         <  1:485701/150‑1 (MQ=255)

                                                                                gTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTc                                                                       >  1:6490551/1‑150 (MQ=255)

                                                                                 ttGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCa                                                                      >  1:793806/1‑150 (MQ=255)

                                                                                  tGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCaa                                                                     >  1:7072774/1‑150 (MQ=255)

                                                                                      ggACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAgcgc                                                                 >  1:4226124/1‑150 (MQ=255)

                                                                                       gACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCa                                                                <  2:6278504/150‑1 (MQ=255)

                                                                                         ccTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACg                                                              >  1:3686796/1‑150 (MQ=255)

                                                                                         ccTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACg                                                              <  1:5330159/150‑1 (MQ=255)

                                                                                          cTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTg                                                                                                                                  >  2:5420896/1‑81 (MQ=255)

                                                                                          cTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTg                                                                                                                                  <  1:5420896/81‑1 (MQ=255)

                                                                                          cTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCt                                                                        >  2:5945959/1‑139 (MQ=255)

                                                                                          cTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCt                                                                        <  1:5945959/139‑1 (MQ=255)

                                                                                           tGAACTGCAGGTGGTCGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGTCGGTTCCGTCACGTCCAGGCTCACATTGGCGTCGTTCTACACGCTGTTGAATGTATCTTTATGGATGCGCATGACCTGACGGTCCTCTAGCGCACGGc                                                            <  1:606558/150‑1 (MQ=255)

                                                                                            gAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCt                                                           <  1:7270075/150‑1 (MQ=255)

                                                                                                tGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCa                                                       <  2:7051359/150‑1 (MQ=255)

                                                                                                tGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCa                                                       >  1:3541705/1‑150 (MQ=255)

                                                                                                   aGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATcc                                                    >  1:4074529/1‑150 (MQ=255)

                                                                                                     gTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAg                                                  >  1:2175220/1‑150 (MQ=255)

                                                                                                       ggACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGc                                                <  1:1347620/150‑1 (MQ=255)

                                                                                                         aCGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCg                                              <  1:3723190/150‑1 (MQ=255)

                                                                                                           ggCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGaa                                            <  1:5662014/150‑1 (MQ=255)

                                                                                                           ggCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGaa                                            <  1:1519536/150‑1 (MQ=255)

                                                                                                           ggCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGaa                                            >  1:4657466/1‑150 (MQ=255)

                                                                                                                   cGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCcca                                    >  2:2472977/1‑150 (MQ=255)

                                                                                                                    ggTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCccag                                   <  2:2175220/150‑1 (MQ=255)

                                                                                                                        tCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGtt                                                                                                              >  2:4340717/1‑71 (MQ=255)

                                                                                                                        tCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCccagcacc                               >  1:520354/1‑150 (MQ=255)

                                                                                                                            gggCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGa                                                                                    <  2:3880527/93‑1 (MQ=255)

                                                                                                                               ccaatgtgAGCCTGGACGTGACGGTACCGTCCAGGCTCGCATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCt                        >  2:5873841/8‑150 (MQ=21)

                                                                                                                               cGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAatgtat                                                                                                      <  1:114333/72‑1 (MQ=255)

                                                                                                                                ggTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCtt                       <  2:6969366/150‑1 (MQ=255)

                                                                                                                                   gCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTtgtg                    <  1:3371815/150‑1 (MQ=255)

                                                                                                                                    ccGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGc                   <  2:2738902/150‑1 (MQ=255)

                                                                                                                                     cGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCg                  <  2:2085340/150‑1 (MQ=255)

                                                                                                                                         ccTGGACGTGACGGTACCGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAg              >  1:4172323/1‑150 (MQ=18)

                                                                                                                                         ccTGGACGGTACCGTCACGTCCAGGCTCTCATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCGTGAAGCGCACGGCTGTCATCCTGGCCGTAGTCGCCCAGCAGCTGCATCTTGTGCGGCAg              >  2:3731670/1‑150 (MQ=21)

                                                                                                                                         ccTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAg              >  2:3809070/1‑150 (MQ=255)

                                                                                                                                           tGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCt            >  2:6179691/1‑150 (MQ=255)

                                                                                                                                            ggACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTc           <  1:2598016/150‑1 (MQ=255)

                                                                                                                                               cGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCAc        <  1:403397/150‑1 (MQ=255)

                                                                                                                                                  tACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCACCTc     >  1:6306936/1‑150 (MQ=255)

                                                                                                                                                   aCCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCACCTCa    <  2:1146598/150‑1 (MQ=255)

                                                                                                                                                   aCCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTAGGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCACCTCa    >  1:6771076/1‑150 (MQ=255)

                                                                                                                                                    ccGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGa                                                                                    <  1:5975611/69‑1 (MQ=255)

                                                                                                                                                     cGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAg                                                  <  2:3472883/102‑1 (MQ=255)

                                                                                                                                                     cGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCACCTCAGc  >  1:6422636/1‑150 (MQ=255)

AGGGCGACCAGTGCTTCGCGCAGCGCCAGCAGCTTGTTCAGCTCGGGCACCTGGCGCGCGATCCGGTCCGGCCCGAAATCGTTGATGGACCTGAACTGCAGGTGGACGGCCAGCTCGGTATCCGGGGCGGTGCCGAGCCTGGACGGTACCGTCACGTCCAGGCTCACATTGGCGTCGTTCAACACGCTGTTGAATGTATGTTTATCGATGCGCATGACCTGACGGTCCTCAAGCGCACGGCTGTCATCCAGGCCGAAGTCGCCCAGCACCAGCATCTTGTGCGGCAGCTCCACCTCAGC  >  NC_002947/3005611‑3005909

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