BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,760,788	C→G	5.5%	A154A (GCG→GCC)	selA ←	selenocysteine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,760,788	0	C	G	5.5%	324.8 / 6.2	104	A154A (GCG→GCC)	selA	selenocysteine synthase
Reads supporting (aligned to +/- strand): ref base C (65/33); new base G (4/2); total (69/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 112 total aligned reads displayed.

TGGTGGTCCCTACTTCGTGTAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCA  >  NC_000913/3760640‑3760936

tggtggTCCCTACTTCGTGTAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACAccgccgccg                                                                                                                                                     <  2:546526/150‑1 (MQ=255)

  gtggtCCCTACTTCGTGTAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTACGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCa                                                                                                                                                   >  2:1785272/1‑150 (MQ=255)

      tCCCTACTTCGTGTAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTGCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCAgtgt                                                                                                                                               >  1:703273/1‑150 (MQ=255)

       cccTACTTCGTGTAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGtt                                                                                                                                              >  2:450169/1‑150 (MQ=255)

       cccTACTTCGTGTAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCCCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATAGtt                                                                                                                                              >  1:1467515/1‑150 (MQ=255)

          tACTTCGTGTAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTAtt                                                                                                                                           <  2:203562/150‑1 (MQ=255)

          tACTTCGTGTAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTAtt                                                                                                                                           <  2:1146815/150‑1 (MQ=255)

           aCTTCGTGTAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTg                                                                                                                                          <  2:1596569/150‑1 (MQ=255)

             ttCGTGTAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGAc                                                                                                                                        <  1:1698089/150‑1 (MQ=255)

                gtgtTGGGTGCAGCCTGCCTGACGCATAAAATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGATCGCCGCGAGATACGACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCGGCATTGTTATTGACGAt                                                                                                                                     <  2:1387867/150‑1 (MQ=255)

                  gtAGGGTGCAGCCTGACTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATAc                                                                                                                                   <  2:390834/150‑1 (MQ=255)

                   tAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACAATGCGGCGGCGGTGTTATTGatgttggc                                                                                                                                  <  1:387636/150‑8 (MQ=38)

                         gCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATc                                                                                                                            >  1:709059/1‑150 (MQ=255)

                           aGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCtt                                                                                                                          >  1:1268222/1‑150 (MQ=255)

                             cctgcctgACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGCGATACCACCACCTCTTTTCCGCTGGCGGTGGCCGCCAACATCAATAACACCGCCGCCGCAGTGTTATTGACGATACAGGCATCTTcc                                                                                                                        >  2:668106/1‑150 (MQ=255)

                             cctgcctgACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTcc                                                                                                                        <  2:1041009/150‑1 (MQ=255)

                             cctgcctgACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTcc                                                                                                                        >  2:681580/1‑150 (MQ=255)

                             cctgcctgACGCATAACATCGGGAATACCAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTcc                                                                                                                        <  1:809737/150‑1 (MQ=255)

                               tgcctgACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGc                                                                                                                      <  1:1827043/150‑1 (MQ=255)

                               tgcctgACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGc                                                                                                                      >  2:1724054/1‑150 (MQ=255)

                               tgcctgACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGc                                                                                                                      <  2:43458/150‑1 (MQ=255)

                                    gACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCg                                                                                                                 <  1:1025626/150‑1 (MQ=255)

                                        cATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCCCCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCAGTGTTTTTGACGATACAGGCATCTTCCGCACCCGTAAt                                                                                                             >  2:187442/1‑150 (MQ=255)

                                        cATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAAt                                                                                                             >  1:1462972/1‑150 (MQ=255)

                                            aCATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACgg                                                                                                         <  1:1302636/150‑1 (MQ=255)

                                               tCGGGAATACGAAACGCGCCGCCAATCTCCACCAGATCGCCGCGAGATACCACCACCTCTTATCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGcac                                                                                                      <  2:1392892/150‑1 (MQ=255)

                                                cGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGcaca                                                                                                     >  1:938953/1‑150 (MQ=255)

                                                 gggAATATGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTATCCCCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCAcag                                                                                                    >  2:831961/1‑150 (MQ=255)

                                                  ggAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGGGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCAcagc                                                                                                   <  1:1177084/150‑1 (MQ=255)

                                                     aTACGAAACGCGCCGCCAATCTCCACCAGATCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACagcagc                                                                                                >  2:1004995/1‑150 (MQ=255)

                                                       aCGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGGCGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCCGCTg                                                                                              >  2:1336519/1‑150 (MQ=255)

                                                         gAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACCGGCATCTTCCGCCCCCGGGATACGGCACAGCAGATgcg                                                                                            >  2:319770/1‑150 (MQ=255)

                                                          aaaCGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTgcgc                                                                                           >  2:780344/1‑150 (MQ=255)

                                                               cgccgaCAATCTCCACCACTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGTCACAGCAGCTGCGCAAgcg                                                                                      <  1:283945/150‑1 (MQ=255)

                                                                gccgccAATCTCCACAAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAgcgc                                                                                     >  1:931036/1‑150 (MQ=255)

                                                                 ccgccAATCTCCCCCAGTTCGCCGCGAGCTACCACCCCCTCTTTTCCGCTGGCAGCGGCCGCCAACATCAATAACGCCGCCGCCGCATTGTAATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCCCCTGCGCCAGCGcc                                                                                    >  2:91438/1‑150 (MQ=255)

                                                                 ccgccAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCAGCCAACATCAATAACACCGCCGCCGCTTTGTTATTGCCGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGcc                                                                                    >  2:1111686/1‑150 (MQ=255)

                                                                  cgccAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCGTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACATTCATCTTCCGCCCCCGTAATACGGCACATCAGCTGCGCCAGCGccc                                                                                   >  2:854096/1‑150 (MQ=255)

                                                                     cAATCTCCACCAGTTCGCCGCCAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCACCATCAATAACACCGCCGCCGCATTGTTATTGCCGATACAGGCATCTTCCGCCCCCGTAATACGGCACCGCAGCTGCGCCAGCGCCCgga                                                                                >  2:1519243/1‑148 (MQ=255)

                                                                       aTCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATcg                                                                              >  1:544636/1‑150 (MQ=255)

                                                                        tctcCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATAAGGCACAGCAGCTGCGCCAGCGCCCGATcgc                                                                             >  2:942453/1‑150 (MQ=255)

                                                                         ctcCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATcgcg                                                                            <  2:168354/150‑1 (MQ=255)

                                                                         ctcCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATcgcg                                                                            >  1:983927/1‑150 (MQ=255)

                                                                           ccaccaGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGAt                                                                          >  2:945901/1‑150 (MQ=255)

                                                                            caccaGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATg                                                                         >  2:597928/1‑150 (MQ=255)

                                                                             accaGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGt                                                                        >  2:452060/1‑150 (MQ=255)

                                                                                 gTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACCGCAGCTGCGCCCGCGCCCGATCGCGATGTCcgc                                                                    >  2:668090/1‑150 (MQ=255)

                                                                                 gTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCcgc                                                                    >  1:1581633/1‑150 (MQ=255)

                                                                                  gtCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACAGCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCcgcg                                                                   <  1:269819/149‑1 (MQ=255)

                                                                                   tCGCCGCGAGAAACCACCACCTCTATTCCGCTCGCAGTGGACGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATAATCCACAGCAGCTGCACCAGCGCCCGATCGCGATCTCcgcct                                                                  >  2:1738718/1‑148 (MQ=255)

                                                                                     gccgcGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTcc                                                                >  2:559955/1‑150 (MQ=255)

                                                                                      ccgcGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCg                                                               >  2:1604868/1‑150 (MQ=255)

                                                                                      ccgcGAGATACCACCACCTCTTATCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCAGAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCg                                                               <  1:216045/150‑1 (MQ=255)

                                                                                        gcgAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGc                                                             <  1:1727152/150‑1 (MQ=255)

                                                                                             aTACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGcgtcgg                                                        >  1:1520347/1‑149 (MQ=255)

                                                                                               accaccaccTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTcc                                                      >  2:1652916/1‑150 (MQ=255)

                                                                                               accaccaccTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTcc                                                      >  1:72479/1‑150 (MQ=255)

                                                                                                ccaccaccTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCa                                                     >  2:229953/1‑150 (MQ=255)

                                                                                                 caccaccTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCAGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGTCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAg                                                    >  1:863262/1‑150 (MQ=255)

                                                                                                  accaccTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGa                                                   >  1:1609962/1‑150 (MQ=255)

                                                                                                   ccaccTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACCCCGCCGCCGCAGTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCCGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGAt                                                  <  2:866488/150‑1 (MQ=255)

                                                                                                   ccaccTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGAt                                                  <  1:681580/150‑1 (MQ=255)

                                                                                                   ccaccTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGAt                                                  >  2:689918/1‑150 (MQ=255)

                                                                                                         ctTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCCGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATaag                                            >  2:1545416/1‑148 (MQ=255)

                                                                                                          ttttCCGCTGGCAGTGGCCGCCAACATCAATAACCCCGCCGCCGCATTGTTATTGACGATACAGGCAGCTTCCGCTCCCGTAATACGGGACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATAccc                                           <  1:237108/150‑3 (MQ=255)

                                                                                                           tttCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCg                                          >  1:1677629/1‑150 (MQ=255)

                                                                                                              ccGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAgg                                       >  1:292332/1‑150 (MQ=255)

                                                                                                               cGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAggg                                      >  1:1250627/1‑150 (MQ=255)

                                                                                                               cGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAggg                                      >  1:1817774/1‑150 (MQ=255)

                                                                                                                  tGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTcc                                   >  1:1242939/1‑149 (MQ=255)

                                                                                                                     cAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTg                                <  1:540414/150‑1 (MQ=255)

                                                                                                                       gTGGCCGCCAACATCAATAACCCCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGc                              <  2:1268222/150‑1 (MQ=255)

                                                                                                                         ggCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGa                            >  1:425278/1‑150 (MQ=255)

                                                                                                                          gccgccAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGaa                           >  2:220641/1‑150 (MQ=255)

                                                                                                                           ccgccAACTTCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCTCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAAc                          <  2:1138839/150‑1 (MQ=255)

                                                                                                                           ccgccAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACCGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAAc                          >  1:1197485/1‑150 (MQ=255)

                                                                                                                           ccgccAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAAc                          <  1:491736/150‑1 (MQ=255)

                                                                                                                           ccgccAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAAc                          >  2:1365241/1‑150 (MQ=255)

                                                                                                                            cgccAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACg                         <  1:1264188/150‑1 (MQ=255)

                                                                                                                             gccAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGc                        <  1:1724054/150‑1 (MQ=255)

                                                                                                                                aaCATCAATAACACCGCCGCCGCATTGTTATTGCCGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATa                     <  2:1467515/150‑1 (MQ=255)

                                                                                                                                 tcgtCAATAACAATGCGGCGGCGGTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACCCCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTACAGATCATAATCGAGGGTCACTGGCGAACGCATAg                    >  2:1570169/4‑150 (MQ=255)

                                                                                                                                 tcgtCAATAACAATGCGGCGGCGGTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACGCGAGGGTCACTGGCGAACGCATAg                    >  2:1357691/4‑150 (MQ=255)

                                                                                                                                   gtCAATAACAATGCGGCGGCGGTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGcc                  >  1:425721/2‑150 (MQ=255)

                                                                                                                                   aTCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGCTCATACTCGAGGGTCACTGGCGAACGCATAGcc                  >  1:148727/1‑150 (MQ=255)

                                                                                                                                     cAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGTACGCATAGCCTg                >  2:748234/1‑150 (MQ=255)

                                                                                                                                     cAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGcccg                <  2:544053/150‑3 (MQ=255)

                                                                                                                                     cAATAACAATGCGGCGGCGGTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTg                >  1:659957/1‑150 (MQ=255)

                                                                                                                                      aaTAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGCACGCATAGCCTgc               >  2:415956/1‑150 (MQ=255)

                                                                                                                                        tAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTgcgc             >  1:905638/1‑150 (MQ=255)

                                                                                                                                          acacCGCCGCCGCATTGTTATTGGCGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCCCTGGCGAACGCATAGCCTGCGCaa           <  1:1336519/150‑1 (MQ=255)

                                                                                                                                          acacAGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCACGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCaa           >  2:1797846/1‑150 (MQ=255)

                                                                                                                                            acCGCCGCCTCATTGTTATTGACGATACAGTCATATTCCGCCCCCGTAATACGGCACACCAGCTGCGCCCGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATAATACTAGAGGGTCACTCGCGAACGCATAGCCTGCGCAAcc         >  1:462841/1‑149 (MQ=255)

                                                                                                                                             ccgccgccgCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGGCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACgg        <  1:450169/150‑1 (MQ=255)

                                                                                                                                               gccgccgcATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCt      <  1:945901/150‑1 (MQ=255)

                                                                                                                                                ccgccgcATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGACCGCATAGCCTGCGCAACGGCtt     >  2:1502300/1‑150 (MQ=255)

                                                                                                                                                ccgccgcATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCtt     >  1:1779548/1‑150 (MQ=255)

                                                                                                                                                ccgccgcATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCtt     <  2:1686025/150‑1 (MQ=255)

                                                                                                                                                  gccgcATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCCTAGCCTGCGCACCGGCTTcc   >  1:1821145/1‑150 (MQ=255)

                                                                                                                                                   ccgcATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTccc  >  2:1747316/1‑149 (MQ=255)

TGGTGGTCCCTACTTCGTGTAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCA  >  NC_000913/3760640‑3760936

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 14 ≤` `ATCG/ATCG` `< 15 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 36 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`