BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	991,560	T→C	6.7%	I60T (ATC→ACC)	der →	GTPase Der

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	991,560	0	T	C	6.7%	597.6 / 20.6	193	I60T (ATC→ACC)	der	GTPase Der
Reads supporting (aligned to +/- strand): ref base T (100/80); new base C (7/6); total (107/86)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

Only 100 of 195 total aligned reads displayed.

CCCGAACGTCGGCAAATCCACCATGTTCAACCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCGGCAGACCAGATGATCGCCGAGCACCTGCGCAAGCG  >  NC_002947/991411‑991707

cccGAACGTCGGCAAATCCACCATGTTCAACCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTAt                                                                                                                                                     <  1:784547/150‑1 (MQ=255)

  cGAACGTCGGCAAATCCACCATGTTCAACCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCa                                                                                                                                                   <  1:1244738/150‑1 (MQ=255)

     aCGTCGGCAAATCCACCATGTTCAACCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCg                                                                                                                                                <  1:5147684/150‑1 (MQ=255)

        tCGGCAAATCCACCATGTTCAACCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTg                                                                                                                                             >  1:686670/1‑150 (MQ=255)

         cGGCAAATCCACCATGTTCAACCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGa                                                                                                                                            <  1:4303466/150‑1 (MQ=255)

         cGGCAAATCCACCATGTTCAACCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGATACCACCGGTgt                                                                                                                                            <  2:5039236/150‑2 (MQ=255)

           gCAAATCCACCATGTTCAACCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACg                                                                                                                                          >  2:2082834/1‑150 (MQ=255)

             aaaTCCACCATGTTCAACCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGTCACCGGTGATACCACCGGTGTCGat                                                                                                                                        <  2:3375586/150‑2 (MQ=255)

             aaaTCCACCATGTTCAACCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGaa                                                                                                                                        >  1:1890117/1‑150 (MQ=255)

               aTCCACCATGTTCAACCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGTCACCGGTGATACCACCGGTGTCGatca                                                                                                                                      <  1:1411518/150‑4 (MQ=255)

                tCCACCATGTTCAACCGCCTGACCCAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGGGACCGCCAGGATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGGATCACCGGTGACGAAGTg                                                                                                                                     >  1:3851780/1‑150 (MQ=255)

                tCCACCATGTTCAACCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCCGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTg                                                                                                                                     >  1:1514588/1‑150 (MQ=255)

                tCCACCATGTTCAACCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTg                                                                                                                                     >  2:2945502/1‑150 (MQ=255)

                  caccaTGTTCAACCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTggg                                                                                                                                   >  2:4313335/1‑150 (MQ=255)

                  caccaTGTTCAACCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTggg                                                                                                                                   >  1:12989/1‑150 (MQ=255)

                   accaTGTTCAACCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGc                                                                                                                                  >  2:4987576/1‑150 (MQ=255)

                   accaTGTTCAACCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGc                                                                                                                                  >  1:3296819/1‑150 (MQ=255)

                   accaTGTTCAACCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGc                                                                                                                                  >  1:2848286/1‑150 (MQ=255)

                     caTGTTCAACCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCAt                                                                                                                                >  1:4634739/1‑150 (MQ=255)

                       tGTTCAACCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATgg                                                                                                                              >  1:3926113/1‑150 (MQ=255)

                          tCAACCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACg                                                                                                                           >  1:3390701/1‑150 (MQ=255)

                            aaCCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGaa                                                                                                                         >  2:5202239/1‑150 (MQ=255)

                            aaCCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGaa                                                                                                                         >  1:1004858/1‑150 (MQ=255)

                            aaCCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGaa                                                                                                                         >  1:4530654/1‑150 (MQ=255)

                             aCCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGGGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGaaa                                                                                                                        <  2:1738525/150‑1 (MQ=255)

                                gCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCCGGATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTCTCACCGGTGACGAAGTGGGCAGGGACGaaaaaa                                                                                                                     >  2:3393834/1‑150 (MQ=255)

                                 ccTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAaataat                                                                                                                    >  2:2825555/1‑150 (MQ=255)

                                 ccTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAAt                                                                                                                    <  2:3283751/150‑1 (MQ=255)

                                 ccTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAAt                                                                                                                    <  1:1023838/150‑1 (MQ=255)

                                    gACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTCTGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGc                                                                                                                 >  2:5135482/1‑150 (MQ=255)

                                    gACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGc                                                                                                                 >  2:2137950/1‑150 (MQ=255)

                                       cAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGa                                                                                                              <  1:4827572/150‑1 (MQ=255)

                                        aaGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAg                                                                                                             <  1:4990338/150‑1 (MQ=255)

                                         aGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGc                                                                                                            >  2:753998/1‑150 (MQ=255)

                                          gACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCa                                                                                                           <  1:4148751/150‑1 (MQ=255)

                                            cccGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCCTGGACGAAAAAATGGCCGAGCAGt                                                                                                         >  1:30223/1‑150 (MQ=255)

                                             ccGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGGATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTc                                                                                                        >  1:937163/1‑150 (MQ=255)

                                               gcgATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGATACCACCGGt                                                                                                                                              >  2:1540978/1‑110 (MQ=255)

                                                  aTGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCa                                                                                                   >  2:4936430/1‑150 (MQ=255)

                                                  aTGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCa                                                                                                   <  1:220333/150‑1 (MQ=255)

                                                      cATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGc                                                                                               >  1:1713935/1‑150 (MQ=255)

                                                       aTCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATgg                                                                                                                              >  2:401490/1‑118 (MQ=255)

                                                       aTCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATgg                                                                                                                              <  1:401490/118‑1 (MQ=255)

                                                       aTCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGcc                                                                                              >  1:1529606/1‑150 (MQ=255)

                                                        tcgtcgGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCAAGGCCa                                                                                             >  1:1826147/1‑150 (MQ=255)

                                                          gtcgGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCAtt                                                                                           <  2:1004858/150‑1 (MQ=255)

                                                           tcgGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTg                                                                                          >  2:705781/1‑150 (MQ=255)

                                                             ggTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTgaa                                                                                        <  1:4154111/150‑1 (MQ=255)

                                                             ggTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTgaa                                                                                        <  1:1918091/150‑1 (MQ=255)

                                                              gTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACG‑AAAAATGGCCGAGCAGTCGCTCATGGCCATTgaaga                                                                                      <  1:3214021/150‑1 (MQ=255)

                                                                gACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTgaagaa                                                                                     <  2:2681261/150‑1 (MQ=255)

                                                                      tCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGAc                                                                               <  2:1657918/150‑1 (MQ=255)

                                                                         ggCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCCTTGAAGAAGCCGACTAt                                                                            >  2:3195093/1‑150 (MQ=255)

                                                                           ccTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTc                                                                                                    <  2:378252/124‑1 (MQ=255)

                                                                           ccTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTAtgt                                                                          <  2:1903812/150‑1 (MQ=255)

                                                                              gACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCt                                                                       <  1:3041780/150‑1 (MQ=255)

                                                                                cccGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGTACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGt                                                                     <  1:817948/150‑1 (MQ=255)

                                                                                 ccGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGtt                                                                    >  1:1804860/1‑150 (MQ=255)

                                                                                 ccGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGtt                                                                    >  1:4958061/1‑150 (MQ=255)

                                                                                    tGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCt                                                                 >  1:4365243/1‑150 (MQ=255)

                                                                                     gACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTg                                                                <  1:3946707/150‑1 (MQ=255)

                                                                                      aCCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTgg                                                               <  1:2930528/150‑1 (MQ=255)

                                                                                        cGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGtt                                                             <  2:3926177/150‑1 (MQ=255)

                                                                                        cGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGtt                                                             <  1:631389/150‑1 (MQ=255)

                                                                                        cGCCAGTATGGTGACGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGTTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGtt                                                             <  1:861170/150‑1 (MQ=255)

                                                                                         gCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTg                                                            >  1:3455583/1‑150 (MQ=255)

                                                                                           cAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGCCTATGTGCTGTTCCTGGTTGAt                                                          <  1:3584388/150‑1 (MQ=255)

                                                                                           cAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGAt                                                          <  2:5364070/150‑1 (MQ=255)

                                                                                             gTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGc                                                        <  2:2329474/150‑1 (MQ=255)

                                                                                                 ggTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGt                                                    >  1:2234621/1‑150 (MQ=255)

                                                                                                 ggTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGt                                                    >  1:2234825/1‑150 (MQ=255)

                                                                                                   tGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGc                                                  >  2:1577493/1‑150 (MQ=255)

                                                                                                    gATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCg                                                 <  1:3336226/150‑1 (MQ=255)

                                                                                                     aTGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCgg                                                <  2:356072/150‑1 (MQ=255)

                                                                                                         cAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAg                                                                                                          <  1:2816946/88‑1 (MQ=255)

                                                                                                         cAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTAt                                            >  2:1251975/1‑150 (MQ=255)

                                                                                                         cAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGAAGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGAATATGTGCTGTTCCTGGTTGATGCCCGTGAGGGTAt                                            >  2:2529151/1‑150 (MQ=255)

                                                                                                           gCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGa                                          <  1:166731/150‑1 (MQ=255)

                                                                                                           gCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGa                                          >  2:1175119/1‑150 (MQ=255)

                                                                                                                       gTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAACCCCACTATGTGCTGTTCCTGGTTGATGCCCGAGCGGGTATGACTGCGGCAGAcc                              >  1:148030/1‑150 (MQ=255)

                                                                                                                        ttCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCGGCAGACCa                             >  1:1974723/1‑150 (MQ=255)

                                                                                                                          ccTTCATCCTGATCGACCCCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAACAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTCTGTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCGGCAGACCAga                           <  2:3461891/150‑1 (MQ=255)

                                                                                                                           cTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCGGCAGACCAgat                          <  2:686670/150‑1 (MQ=255)

                                                                                                                            ttCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCGGCAGACCAgatg                         >  1:1281651/1‑150 (MQ=255)

                                                                                                                            ttCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCGGCAGACCAgatg                         >  2:3315499/1‑150 (MQ=255)

                                                                                                                             tCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCGGCAGACCAgatga                        <  1:3315499/150‑1 (MQ=255)

                                                                                                                             tCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCGGCAGACAAgatga                        >  1:5417208/1‑150 (MQ=255)

                                                                                                                              cATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCGGCAGACCAggatcg                       <  2:4530654/150‑6 (MQ=255)

                                                                                                                                tCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCGGCAGACCAGATGATCg                     <  2:3040249/150‑1 (MQ=255)

                                                                                                                                    gATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGc                                                                                                      >  2:2031295/1‑65 (MQ=255)

                                                                                                                                    gATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGGTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCGGGTGGATGCCCGGGCGGGTCTGACGGCGGCAGACCAGATGATCGCCgc                 >  2:3844126/1‑149 (MQ=255)

                                                                                                                                       cGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCGGCAGACCAGATGATCGCCGAGCa              >  1:3109594/1‑150 (MQ=255)

                                                                                                                                         acacCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCGGCAGAcc                              >  1:2028450/1‑132 (MQ=255)

                                                                                                                                         acacCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCGGCAGACCAGATGATCGCCGAGCAcc            >  2:1962036/1‑150 (MQ=255)

                                                                                                                                         acacCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCGGGTTGATGCCCGTGCGGGTATGACTGCGGCAGACCAGATGATCGCCGAGCAcc            >  1:1855397/1‑150 (MQ=255)

                                                                                                                                          cacCGGTGGTATCACCGGTGACGAAGGGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGCAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCGGCAGACCAGATGATCGCCGAGCACCt           <  1:2121715/150‑1 (MQ=255)

                                                                                                                                             cGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTg                                                                <  2:1935691/94‑1 (MQ=255)

                                                                                                                                              ggtggtATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCGGCAGACCAGATGATCGCCGAGCACCTgcgc       >  1:2490235/1‑150 (MQ=255)

                                                                                                                                               gtgatACCACCGGTGTCGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCCTGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCGGCAGACCAGATGATCGCCGAGCACCTGCGCa      >  2:4210838/1‑150 (MQ=255)

                                                                                                                                                  gtATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCGGCAGACCAGATGATCGCCGAGCACCTGCGCAAGc   <  2:3109594/150‑1 (MQ=255)

CCCGAACGTCGGCAAATCCACCATGTTCAACCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGATCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCGGCAGACCAGATGATCGCCGAGCACCTGCGCAAGCG  >  NC_002947/991411‑991707

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 11 ≤` `ATCG/ATCG` `< 25 ≤` `ATCG/ATCG` `< 37 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`