BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,392,471	A→C	5.2%	intergenic (+36/‑1)	glyV → / → glyX	tRNA‑Gly/tRNA‑Gly

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,392,471	0	A	C	5.2%	1264.6 / 15.1	439	intergenic (+36/‑1)	glyV/glyX	tRNA‑Gly/tRNA‑Gly
Reads supporting (aligned to +/- strand): ref base A (140/266); new base C (2/31); total (142/297)
Fisher's exact test for biased strand distribution p-value = 3.53e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 442 total aligned reads displayed.

AGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGC‑ACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCA  >  NC_000913/4392322‑4392619

aGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGC‑ACAGACCACCCaa                                                                                                                                                      <  1:1548655/150‑1 (MQ=255)

    aTGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGC‑ACAGACCACCCAAGCgg                                                                                                                                                  <  1:2376373/150‑1 (MQ=255)

     tGCGAAATTTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGC‑ACAGACCACCCAAGCggg                                                                                                                                                 <  1:7000082/150‑1 (MQ=255)

     tGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGC‑ACAGACCACCCAAGCggg                                                                                                                                                 <  1:10002338/150‑1 (MQ=255)

     tGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGC‑ACAGACCACCCAAGCggg                                                                                                                                                 <  2:3382883/150‑1 (MQ=255)

     tGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGC‑ACAGACCACCCAAGCggg                                                                                                                                                 <  2:9597304/150‑1 (MQ=255)

      gCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATGTGAAAAGTGATGCAAAGC‑ACAGACCACCCAAGCGGGa                                                                                                                                                >  1:11035548/1‑150 (MQ=255)

       cGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGC‑ACAGACCACCCAAGCGGGaa                                                                                                                                               <  2:7447083/150‑1 (MQ=255)

       cGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTGGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGC‑ACAGACCACCCAAGCGGGaa                                                                                                                                               <  2:9931064/150‑1 (MQ=255)

        gAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTG‑AAAGTGCTGTAAGGC‑ACCTACCCCCCAAGCGGGAATa                                                                                                                                             >  2:5358042/1‑150 (MQ=14)

        gAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATATGAAAAGTGATGCAAAGC‑ACAGACCACCCAAGCGGGAAt                                                                                                                                              >  2:8113893/1‑150 (MQ=255)

         aaaaTTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGC‑ACAGACCACCCAAGCGGGAATa                                                                                                                                             <  1:1168111/150‑1 (MQ=255)

         aaaaTTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGC‑ACAGACCACCCAAGCGGGAATa                                                                                                                                             <  1:5855099/150‑1 (MQ=255)

         aaaaTTACGAAAGCAAAATTAAGAAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAACGTCGGGGTCGCGAGTTCGAGTCTCGATTCCCGCTCAAAAATTTGAAAAGTGCTGCAAAGC‑ACAGACCACCCAAGCGGGAATa                                                                                                                                             <  1:10292130/150‑1 (MQ=255)

          aaaTTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGC‑ACAGACCACCCAAGCGGGAATAg                                                                                                                                            <  1:3506606/150‑1 (MQ=255)

            aTTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGC‑ACAGACCACCCAAGCGGGAATAGCt                                                                                                                                          <  2:8355204/150‑1 (MQ=255)

            aTTAAGAAAGCAAAATTAAGTAGTACGCGTGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGC‑ACAGACCACGCAAGCGGGAATAGCt                                                                                                                                          <  2:6610974/150‑1 (MQ=255)

             ttACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGC‑ACAGACCACCCAAGCGGGAATAGCTc                                                                                                                                         <  1:9511573/150‑1 (MQ=255)

              tACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGC‑ACAGACCACCCAAGCGGGAATAGCTCa                                                                                                                                        <  1:8488244/150‑1 (MQ=255)

               aCGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTTAAAAGTGATGAAAAGC‑ACATACCACCCAAGCGGGAATAGCTCAg                                                                                                                                       >  2:4749353/1‑150 (MQ=255)

                    aGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGC‑ACAGACCACCCAAGCGGGAATAGCTCAGTTGGt                                                                                                                                  <  2:3005919/150‑1 (MQ=255)

                    aGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGC‑ACAGACCACCCAAGCGGGAATAGCTCAGTTGGt                                                                                                                                  <  1:1382577/150‑1 (MQ=255)

                     gCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGC‑ACAGACCACCCAAGCGGGAATAGCTCAGTTGGTa                                                                                                                                 <  2:6253774/150‑1 (MQ=255)

                     gCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGC‑ACAGACCACCCAAGCGGGAATAGCTCAGTTGGTa                                                                                                                                 <  1:5805570/150‑1 (MQ=255)

                     gCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGC‑ACAGACCACCCAAGCGGGAATAGCTCAGTTGGTa                                                                                                                                 <  2:3744231/150‑1 (MQ=255)

                       aaaaTTAAGTAGTACGCGGGAATAGCCCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGC‑ACAGACCACCCAAGCGGGAATAGCTCAGTTGGTaga                                                                                                                               <  1:2177565/150‑1 (MQ=255)

                        aatttAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAATGTGCTGCAAAGC‑ACAGACCACCCAAGCGGGACTAGCTCAGTGGGTagtg                                                                                                                              >  1:6324693/4‑148 (MQ=255)

                        aaaTTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGC‑ACAGACCACCCAAGCGGGAATAGCTCAGTTGGTagag                                                                                                                              <  1:8009523/150‑1 (MQ=255)

                         aaTTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGC‑ACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGc                                                                                                                             <  2:4063238/150‑1 (MQ=255)

                                 agtaCGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGC‑ACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCtt                                                                                                                     <  1:7774881/150‑1 (MQ=255)

                                  gtaCGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGC‑ACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTg                                                                                                                    <  1:8761756/150‑1 (MQ=255)

                                   taCGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATGTGAAAACAGCAGCAAAGCAACAG‑CCACCCAAGCGGGAATAGCTCGGTTGGTATAGCACGAccaacc                                                                                                                   >  2:10714453/1‑146 (MQ=14)

                                     cgcgGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGC‑ACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCa                                                                                                                 <  1:9676978/150‑1 (MQ=40)

                                      gcgGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGC‑ACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCaa                                                                                                                <  1:10888738/150‑1 (MQ=39)

                                        gggAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGC‑ACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAgg                                                                                                              >  2:425154/1‑150 (MQ=39)

                                         ggAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGC‑ACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGt                                                                                                             <  1:5209035/150‑1 (MQ=39)

                                            aTAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGC‑ACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCgg                                                                                                          <  1:9316040/150‑1 (MQ=39)

                                             tAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGC‑ACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCggg                                                                                                         <  1:6194553/150‑1 (MQ=39)

                                              aGCTCAGTTGGTAGAGCACGACCTTGCCATGGTCGGGGTCGCGTGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGC‑ACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCgggg                                                                                                        <  2:10006266/150‑1 (MQ=255)

                                                    gTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGC‑ACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGa                                                                                                  <  2:1931997/150‑1 (MQ=39)

                                                    gTTGGGAGAGCACGCCCGGGCCCAGGGCGGGGTCGCGAGTGCGAGTCTCGGTTCCCGCTCCAAACCTGGAAAAGTGCTGCCACGC‑ACAGACCCCCCAAGCGGGAAGAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGa                                                                                                  <  1:10576985/150‑1 (MQ=11)

                                                     ttGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGC‑ACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAg                                                                                                 <  1:1866207/150‑1 (MQ=39)

                                                       ggTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGC‑ACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGtt                                                                                               <  2:6160535/150‑1 (MQ=39)

                                                          tgtggtcGAGCTTGCCTATTTCGGGGTCGCTTGTTCGAGTCTCGTTTCCCTCTCCTAAATTTGAAAAGTGCTGCAAAGC‑ACTGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGa                                                                                            <  2:2645689/144‑1 (MQ=14)

                                                           gagCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGC‑ACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTAGCCAAGGTCGGGGTCGCGAGATCGAg                                                                                           >  2:10291024/1‑150 (MQ=255)

                                                            agCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTCCTGCAAAGC‑ACAGACCACCCAAGCGGGAATAGCTCCGTTGGTAGAGCACGACCTTGCCCAGGTCGGGGTCGCGAGTTCGAGt                                                                                          >  1:8557534/1‑150 (MQ=255)

                                                               aCGACCTTGCCTAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTACAAAATTTGAAAAGTGCTGCAAAGC‑ACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGtctc                                                                                       <  2:10887251/150‑1 (MQ=255)

                                                               aCGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGC‑ACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGtctc                                                                                       <  1:5640384/150‑1 (MQ=39)

                                                                     ttGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGC‑ACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTcc                                                                                 <  1:4050591/150‑1 (MQ=39)

                                                                        ccAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGC‑ACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGc                                                                              >  1:5370551/1‑150 (MQ=39)

                                                                          aaGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGC‑ACAGACCACCCACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTc                                                                            <  1:10194354/150‑1 (MQ=255)

                                                                          aaGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGC‑ACAGACCACCCAAGCGGGAATAGCTCAGTTGGGAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTc                                                                            <  1:3141354/150‑1 (MQ=255)

                                                                             gTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGC‑ACAGACCTCCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCaa                                                                         <  2:2559008/150‑1 (MQ=255)

                                                                             gTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGATGAAAAGC‑ACAGACCACCCAATCGGGAATAGCTCGGTTGTTATAGCACTACCTTGCCAATGTCGGGGTAGCGAGTTCGAGACTCGTTTCCCGCTCCaa                                                                         >  2:1117311/1‑150 (MQ=11)

                                                                                 gggTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGC‑ACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAAtt                                                                     <  1:9831060/150‑1 (MQ=255)

                                                                                 gggTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGC‑ACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAAtt                                                                     >  2:2010693/1‑150 (MQ=255)

                                                                                  ggTCGCGATTTCGAGTATCGTTTCCCTCTCCAAAATTTGAAAAGTGCTGCAAAGA‑ACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAAttt                                                                    <  2:4613739/150‑1 (MQ=255)

                                                                                    tCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGC‑ACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGa                                                                  >  1:4643409/1‑150 (MQ=255)

                                                                                     cgcgAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGC‑ACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGaa                                                                 <  1:2649402/150‑1 (MQ=255)

                                                                                         aGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGAAAAGC‑ACAGACCACCCAAGCGGGAATAGCACAGTTGGTAGAGCAAGACCTAGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAgag                                                             >  1:481852/1‑148 (MQ=255)

                                                                                           ttCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGAAAAGC‑ACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGaaataact                                                           >  2:5740630/1‑145 (MQ=255)

                                                                                              gAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGC‑ACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTa                                                        >  2:1850719/1‑150 (MQ=255)

                                                                                               cctgacGTTTCCGTCTCATAAATTCGAAAAGTGCTGCAACTC‑ACTGACGACCCACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTaa                                                       <  2:3870183/145‑1 (MQ=17)

                                                                                                gTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGC‑ACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAg                                                      <  2:789341/150‑1 (MQ=255)

                                                                                                  ctcGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGC‑ACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGc                                                    <  2:5690731/150‑1 (MQ=255)

                                                                                                   tcGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGC‑ACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATATTAAATAGCTTTaatgac                                                   >  2:6442073/1‑146 (MQ=14)

                                                                                                   tcGTTTACGGCTCCCAAATTAGAACAGTTCTTCATATC‑ACTGACGACGCAAGCGGGAATAGCTCAGTTGGTAGAGCAAGACATTGCCAAGGTCGGGGTCTCGTGTTCGAGTCTAGTTTCCCGCTCCAAACTTTGAAAGTGCTGTAAGGca                                                   <  2:2228232/150‑1 (MQ=17)

                                                                                                        tCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGC‑ACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGac                                              <  2:10603920/150‑1 (MQ=255)

                                                                                                          ccGCTCCAAAATTTGAAAAGTGCTGCAAAGC‑ACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGacca                                            <  1:1287362/150‑1 (MQ=255)

                                                                                                          ccGCTCCAAAATTTGAAAAGTGCTGCAAAGC‑ACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGacca                                            <  1:7492134/150‑1 (MQ=255)

                                                                                                           cGCTCCTAAATTTGAAAAGTGCTGCAAAGC‑ACAGACCACCCAAGCGGGAATAGCTCAGTTGGCAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGaccac                                           <  2:7205562/150‑1 (MQ=255)

                                                                                                             cTCCAAAATTTGAAAAGTGCTGCAAAGC‑ACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCAccc                                         <  1:4252792/150‑1 (MQ=255)

                                                                                                                cAAAATTTGAAAAGTGCTGCAAAGC‑ACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGt                                                                                     <  2:708761/103‑1 (MQ=255)

                                                                                                                cAAAATTTGAAAAGTGCTGCAAAGC‑ACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATAAGAAAGTGCTTTAAGACACAGACCACCCaat                                      >  1:1769984/1‑149 (MQ=255)

                                                                                                                 aaaaTTTGAAAAGTGCTGCAAAGC‑ACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGc                                     >  2:8010807/1‑150 (MQ=255)

                                                                                                                 aaaaTTTGAAAAGTGCTGCAAAGC‑ACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGc                                     >  1:2980210/1‑150 (MQ=255)

                                                                                                                  aaaTTTGAAAAGTGCTGCAAAGC‑ACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCg                                    <  2:1182452/150‑1 (MQ=255)

                                                                                                                   aaTTTGAAAAGTGCTGCAAAGC‑ACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCgg                                   <  1:9026168/150‑1 (MQ=255)

                                                                                                                     tttGAAAAGTGCTGCAAAGC‑ACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTcc                                                                           >  2:2363549/1‑108 (MQ=255)

                                                                                                                      ttGAAAAGTGCTGCAAAGC‑ACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGaa                                <  1:648552/150‑1 (MQ=255)

                                                                                                                      ttGAAAAGTGCTGCAAAGC‑ACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATATGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGaa                                >  2:1472226/1‑150 (MQ=255)

                                                                                                                         aaaaGTGCTGCAAAGC‑ACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGATAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAg                             >  1:5422182/1‑150 (MQ=255)

                                                                                                                         aaaaGTGCTGCAAAGC‑ACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAg                             <  2:5236389/150‑1 (MQ=255)

                                                                                                                           aaGTGCTGCAAAGC‑ACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCATCCAAGCGGGAATAGCt                           <  2:926383/150‑1 (MQ=255)

                                                                                                                            aGTGCTGCAAAGC‑ACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAGTCTGAACGAGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTc                          >  2:9000178/1‑150 (MQ=255)

                                                                                                                               gctgcAAAGC‑ACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGt                       <  1:1082997/150‑1 (MQ=255)

                                                                                                                                 tgcAAAGC‑ACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTg                     >  2:6798671/1‑150 (MQ=255)

                                                                                                                                   cAAAGC‑ACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAATCCACAGACCACCCAAGCGGGAATAGGGCCGTTGGt                   >  2:8544756/1‑150 (MQ=17)

                                                                                                                                   cAAAGC‑ACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCGGTTGGt                   >  1:4431095/1‑150 (MQ=255)

                                                                                                                                   cAAAGC‑ACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGt                   >  1:6889983/1‑150 (MQ=255)

                                                                                                                                   cAAAGC‑ACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGATGTAAGGCACAGACCGCCCAATCGGGTATAGGTCAGTTGGt                   >  1:2814274/1‑150 (MQ=17)

                                                                                                                                       gC‑ACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATGTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTagag               >  2:3148100/1‑150 (MQ=255)

                                                                                                                                        c‑acaGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATGTGACAGTGCTGTAAGGCACAGACCACCCAAGCGGGACTAGCTCCGTGGGTCGAGc              >  2:8849903/1‑150 (MQ=17)

                                                                                                                                             gACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGa          <  2:9414465/150‑1 (MQ=255)

                                                                                                                                              accaccCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGAc         <  1:10418643/150‑1 (MQ=255)

                                                                                                                                               ccaccCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGAcc        <  1:1988152/150‑1 (MQ=255)

                                                                                                                                                  cccAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCCAGCGGGAATAGCTCAGTTGGTAGAGCACGACCGTg     >  1:3698571/1‑150 (MQ=255)

                                                                                                                                                   ccAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGc    >  1:8122847/1‑150 (MQ=255)

                                                                                                                                                   ccAAGCGGGAATAGCGCAGTTGGTAGAGCACGACCTTGCCAAGGGCGGGGTCGCGAGGTCGAGTCTCGTTTCCCGCGCCAAAATCGGAAAGGGCTGTAAGGCCCAGACCGCCCCCGCGCGAAGAGCGCCGTGGGTCGAGCACGAccgtcc    >  2:5094853/1‑146 (MQ=14)

                                                                                                                                                    cAAGCGGGAATAGCACAGTTGGTAGAGCACGACCTTGCCATGGTCGGGGTCGCTAGTTCGAGTCTCGTTTCCCGCCCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGcc   <  2:1725809/150‑1 (MQ=255)

AGAAATGCGAAAATTACGAAAGCAAAATTAAGTAGTACGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAAGTGCTGCAAAGC‑ACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCAAGGTCGGGGTCGCGAGTTCGAGTCTCGTTTCCCGCTCCAAAATTTGAAAGTGCTGTAAGGCACAGACCACCCAAGCGGGAATAGCTCAGTTGGTAGAGCACGACCTTGCCA  >  NC_000913/4392322‑4392619

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