BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,818,046	A→G	5.0%	intergenic (‑186/+13)	argQ ← / ← argZ	tRNA‑Arg/tRNA‑Arg

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,818,046	0	A	G	5.0%	596.5 / 12.1	188	intergenic (‑186/+13)	argQ/argZ	tRNA‑Arg/tRNA‑Arg
Reads supporting (aligned to +/- strand): ref base A (104/69); new base G (3/6); total (108/80)
Fisher's exact test for biased strand distribution p-value = 1.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 191 total aligned reads displayed.

GCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTA  >  NC_000913/2817897‑2818192

    ctaAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTAGCCAAGAGGCTCGCCCTTCGGGCCGTTGCTAAGGCAACGATATCATCCCTAGTGCCTGCGGATAATTCCGCTACCTCTTCGATACCTTGATTGCTGAAGATTTc                                                                                                                                                >  1:3466237/1‑150 (MQ=17)

      aaaGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAATCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCgg                                                                                                                                              <  2:1342762/150‑1 (MQ=255)

      aaaGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGAGGCTCGCCCTTCGGGCCGTTGCTAACGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCgg                                                                                                                                              <  1:2318802/150‑1 (MQ=255)

      aaaGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTCGCGGCTCGCGCTTCGGGGCGTTTCTAAAGCACCGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCgg                                                                                                                                              <  2:308712/150‑1 (MQ=18)

       aaGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGga                                                                                                                                             >  1:2766268/1‑150 (MQ=255)

          ccATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCACTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGaa                                                                                                                                          <  2:2851557/150‑1 (MQ=255)

           cATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAat                                                                                                                                         <  1:688884/150‑1 (MQ=255)

              ttgatttAAGAGTATGTGCCTCAGGGTGTTTTAGTCGCCCTGAGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATg                                                                                                                                      <  1:3949596/145‑1 (MQ=18)

               tCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATgg                                                                                                                                     >  2:3889945/1‑150 (MQ=255)

               tCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCGAAAGCCCCGGTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTGCGATACCTTCATTGCTGAAGATGTCGGAGAATATgg                                                                                                                                     >  2:2811674/1‑150 (MQ=18)

               tCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCAAAAGAAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATgg                                                                                                                                     >  2:757932/1‑150 (MQ=255)

                 aaTTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTg                                                                                                                                   >  2:3511376/1‑150 (MQ=255)

                      agagATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATcc                                                                                                                              >  1:3607769/1‑150 (MQ=255)

                        agaTGGTGCATCCGGGAGGATTACTCGCGTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGTAGATTTCGGAGAATATGGTGCATCCgg                                                                                                                            >  1:2416005/1‑150 (MQ=255)

                        agaTGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCgg                                                                                                                            >  2:1824137/1‑150 (MQ=255)

                          aTGGTGCCTCCGGGAGTCTTAGTCTCCTTGCGGCTCGCCCTTCGGGCCGTTGCTTTAGCAACGTTATCCTCCCTCGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGgga                                                                                                                          <  1:1296307/150‑1 (MQ=17)

                          aTGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGgga                                                                                                                          >  1:3955326/1‑150 (MQ=255)

                          aTGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGgga                                                                                                                          <  1:515833/150‑1 (MQ=255)

                           tGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGggag                                                                                                                         >  1:1001049/1‑150 (MQ=255)

                            ggTGCATCCGGGAGGATTACTCGCCTTGCGCCTCCCCCTGCGCGCCGTGGCTAAACAAACGTTCTCATCCCTAGTCCATGCGGTTAATTCCGCTCCCGCTTAGATACCTTCAGGGCTGTAGAGTGCGGTGAAGATGGTGCATCCGggagg                                                                                                                        >  2:1414463/1‑150 (MQ=12)

                             gggccgccGGGCGGCTGACTCGCCTTGCGGCTCGCCCGTCGGGCCGTTGCTCACGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCGTCGCTACCGTCATTGCGGAAGATTTCGGAGAATATGGTGCATCCGggagga                                                                                                                       <  1:2096010/144‑1 (MQ=18)

                              tGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCAGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGAt                                                                                                                      >  2:1978683/1‑150 (MQ=255)

                                caaccGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTCAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGATACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTc                                                                                                                    <  2:3693599/147‑1 (MQ=255)

                                 ctCCGGGAGGATTCCTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCg                                                                                                                   <  2:3380906/149‑1 (MQ=255)

                                 aTCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCg                                                                                                                   >  1:623337/1‑150 (MQ=255)

                                  tCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGa                                                                                                                  <  2:1767185/150‑1 (MQ=255)

                                    cGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAAc                                                                                                                <  1:558460/150‑1 (MQ=255)

                                    cGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAAc                                                                                                                <  1:2380489/150‑1 (MQ=255)

                                      ggaggaTTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCt                                                                                                              >  1:2072103/1‑150 (MQ=255)

                                       gaggaTTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCACTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACATTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTAGAACCTc                                                                                                             >  2:1108566/1‑150 (MQ=255)

                                       gaggaTTACTCGCCTTGCGGCTCGCCCATCGGGCCGTTGCTAAAGCAACGGTATCATCCCGAGTGCCTGCGGTTAATTCCGCTACCACTTCGATACCGTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGCGGATTCGAACCTc                                                                                                             >  1:335041/1‑150 (MQ=18)

                                          gaTTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGa                                                                                                          >  2:3337534/1‑150 (MQ=255)

                                               cTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCt                                                                                                     <  2:782674/150‑1 (MQ=255)

                                                tCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAGCCTCCGACCGCTc                                                                                                    >  1:3357412/1‑150 (MQ=255)

                                                tCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTc                                                                                                    >  1:240953/1‑150 (MQ=255)

                                                tCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTc                                                                                                    >  2:2034035/1‑150 (MQ=255)

                                                tCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTc                                                                                                    >  2:3799446/1‑150 (MQ=255)

                                                tCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTc                                                                                                    >  2:178960/1‑150 (MQ=255)

                                                tCGCCTTGCGGCTCGCCCGTCGGGCCGTTGCTAAAGCAACCGTATCAACCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGAAACCAGCAATGCTGAAGATTTAGGAGCAAATGGTGCAACCCGGAGGATTCGAACCTCCGACCGcgc                                                                                                    >  1:3166228/1‑148 (MQ=18)

                                                  gCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCAACCGGGAGGATTCGAACCTCCGACCGCTCgg                                                                                                  >  1:3049697/1‑150 (MQ=255)

                                                  gCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGAGACCTTCATTGCTGAAGATTTCGGAGAGGCGGGTGCATCCGGGGGGATTCGAACCTCCGACCGCTCgg                                                                                                  >  2:1042412/1‑150 (MQ=18)

                                                         ggCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAg                                                                                           >  2:951497/1‑150 (MQ=255)

                                                            tCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCg                                                                                        >  2:590172/1‑150 (MQ=255)

                                                                ccTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTa                                                                                    <  1:590172/150‑1 (MQ=255)

                                                                ccTTCGGGCCGTTGCTAAAGCAACGTTATCATACCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCATAATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCACCGACCGCTCGGTTCGTAGCCGAGTa                                                                                    >  1:1705597/1‑150 (MQ=255)

                                                                       gCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCAGCGGTTAATTACGCTACCACTTCGATACCATCAATGCTGAAGATTTCGGAGAATATGGTGCAACCGGGAGGATTCGAACCTCCGACCGCACGGTTCGTAGCCGAGTACTCTATc                                                                             >  1:1594567/1‑150 (MQ=17)

                                                                         cGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCa                                                                           >  1:3378912/1‑150 (MQ=255)

                                                                               tAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAg                                                                     >  1:1599091/1‑150 (MQ=255)

                                                                               tAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAg                                                                     >  1:705797/1‑150 (MQ=255)

                                                                               tAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAg                                                                     >  1:3083123/1‑150 (MQ=255)

                                                                                 aaGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCt                                                                   >  1:2523521/1‑150 (MQ=255)

                                                                                   gCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTAc                                                                 >  1:671733/1‑150 (MQ=255)

                                                                                      aCGTTATCCGCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGTTTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGa                                                              <  1:2104733/150‑1 (MQ=255)

                                                                                      aCGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGa                                                              >  1:1921966/1‑150 (MQ=255)

                                                                                      aCGTTATCATCCCTAGTGCCAGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGa                                                              <  1:1978683/150‑1 (MQ=255)

                                                                                        gTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGTTTTGCGTGAATTTGGAGCATCCGGGAGGATTCGAACCTCGGACCTATCTGTTCGTGGCCGACTACTCTTTCCAGCTGAGCTAGGGATg                                                            >  2:2189063/1‑150 (MQ=14)

                                                                                        gTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATg                                                            >  1:1081118/1‑150 (MQ=255)

                                                                                        gTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATg                                                            <  1:3630266/150‑1 (MQ=255)

                                                                                        gTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATg                                                            <  1:2170960/150‑1 (MQ=255)

                                                                                        gTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATg                                                            >  1:2232277/1‑150 (MQ=255)

                                                                                              atcCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTAGAACCTCAGACCGCTCGGTTCGTAGCCTAGTACTCTATCCAGCAGAGCTACGGATGCATCgg                                                      >  2:3485574/1‑150 (MQ=255)

                                                                                                  cTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGaaa                                                  <  2:2880403/150‑1 (MQ=255)

                                                                                                  cTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGaaa                                                  >  2:228090/1‑150 (MQ=255)

                                                                                                  cTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTCCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGaaa                                                  >  2:3626781/1‑150 (MQ=255)

                                                                                                        ccTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAacttact                                            >  2:305704/1‑148 (MQ=255)

                                                                                                            cGGTTAAGTCCGCTACCTCTTCGATACCGTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTAc                                        <  1:1042668/150‑1 (MQ=255)

                                                                                                                tAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTActgct                                    <  2:671733/150‑1 (MQ=255)

                                                                                                                caagtCCGCTACCACTTCGATACCGTCACTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCTGTTCTGTTCCGAGGTCTGCACCCCGCTGAGCTACGGATGCATCGGGAAATTTATTTTActgct                                    <  1:3357702/146‑1 (MQ=14)

                                                                                                                   ttCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAACTTACTTTACTGCAGat                                 >  1:775181/1‑150 (MQ=255)

                                                                                                                     ccGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGatat                               >  1:2470501/1‑150 (MQ=255)

                                                                                                                      cGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGCGTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGTGGCATCGGGAAAATTAGTTTACTGCGTATAtt                              >  2:2650760/1‑150 (MQ=14)

                                                                                                                        cTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGa                            <  1:3002250/150‑1 (MQ=255)

                                                                                                                             tctTCGATACCTTCATTGCTGAAGATTTCGGAGAatat                                                                                                                                       <  2:2919424/38‑1 (MQ=255)

                                                                                                                             tctTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTCCGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCg                       >  2:742816/1‑150 (MQ=255)

                                                                                                                             tctTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCg                       >  1:2798706/1‑150 (MQ=255)

                                                                                                                              ctTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGc                      >  2:446052/1‑150 (MQ=255)

                                                                                                                               ttCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGa                                                                                                                  <  2:327071/57‑1 (MQ=255)

                                                                                                                                tCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGTGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGcta                    <  2:1086191/150‑1 (MQ=255)

                                                                                                                                 cGATACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCTACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGctac                   <  2:1595755/150‑1 (MQ=16)

                                                                                                                                  gATACCTTCATTGCTGAAGTTTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCg                                                       <  1:2973324/113‑1 (MQ=17)

                                                                                                                                  gATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAACTTATTTTACTGCTGATATTGATACCGctaca                  >  2:850450/1‑149 (MQ=255)

                                                                                                                                    tACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTaa                >  2:1158428/1‑150 (MQ=255)

                                                                                                                                    tACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTaa                >  2:1786300/1‑150 (MQ=255)

                                                                                                                                       cTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGc             <  1:3659817/150‑1 (MQ=18)

                                                                                                                                       cTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGc             <  2:335152/150‑1 (MQ=255)

                                                                                                                                       cTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGc             >  1:1722845/1‑150 (MQ=255)

                                                                                                                                       cTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGc             >  1:929892/1‑150 (MQ=255)

                                                                                                                                       cTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGc             <  2:1594567/150‑1 (MQ=255)

                                                                                                                                          cATTGCTGAAGATTTCGGAGAATTTGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGccgt          <  2:1693179/150‑3 (MQ=255)

                                                                                                                                          cATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTa                                                                                            <  1:814026/68‑1 (MQ=40)

                                                                                                                                           aTTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTAct                                       >  1:2199713/1‑120 (MQ=19)

                                                                                                                                             tGCTGAAGTTTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATc       <  2:2041203/150‑1 (MQ=255)

                                                                                                                                              gCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCa      <  2:2773526/150‑1 (MQ=255)

                                                                                                                                               cTGTAGATTTCGGAGATTAGTGTGCATTCGCGAGGTTTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCaa     <  1:3316492/150‑1 (MQ=14)

                                                                                                                                               cTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCaa     <  2:3818718/150‑1 (MQ=255)

                                                                                                                                                tGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATAACGCTACTAAAGACATATCAAg    >  1:3472259/1‑150 (MQ=255)

                                                                                                                                                 gAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTAAAAGCGGATATTGATACCGCTACTAAAGCCATATAAAGt   >  2:2982112/1‑150 (MQ=11)

                                                                                                                                                 gAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTa                                                                                            >  2:216174/1‑61 (MQ=40)

                                                                                                                                                 gAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTa                                                                                            <  1:216159/61‑1 (MQ=40)

                                                                                                                                                  aaGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAGTTAGTTTACTGCTGATATTGATACCGCTACTCAAGCCATATCAAGTa  >  1:2366134/1‑150 (MQ=255)

GCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTA  >  NC_000913/2817897‑2818192

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