BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,287,041	A→G	6.3%	intergenic (+515/+25)	narI → / ← rttR	nitrate reductase 1, gamma (cytochrome b(NR)) subunit/rtT sRNA, processed from tyrT transcript

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,287,041	0	A	G	6.3%	434.6 / 11.0	156	intergenic (+515/+25)	narI/rttR	nitrate reductase 1, gamma (cytochrome b(NR)) subunit/rtT sRNA, processed from tyrT transcript
Reads supporting (aligned to +/- strand): ref base A (69/69); new base G (5/4); total (81/75)
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 162 total aligned reads displayed.

GGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCA  >  NC_000913/1286892‑1287187

ggtggGGGATGTATTACTCAGCGCTGCGCGCTTCGCCGTTCGTGTGGTTGCCTGCGGCAACGCTCTCTCGCTGCCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTga                                                                                                                                                    <  2:4156826/150‑1 (MQ=14)

ggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTga                                                                                                                                                    <  1:3155721/150‑1 (MQ=255)

ggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTga                                                                                                                                                    <  1:2363660/150‑1 (MQ=255)

  tggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTgaga                                                                                                                                                  <  1:341086/150‑1 (MQ=255)

       gAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCa                                                                                                                                             <  1:663963/150‑1 (MQ=255)

       gAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCa                                                                                                                                             <  1:3850456/150‑1 (MQ=255)

        taGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAg                                                                                                                                            <  1:3437156/149‑1 (MQ=255)

        aaGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAg                                                                                                                                            <  1:3437130/150‑1 (MQ=255)

                cTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACGCTCTCTCGCTGACGCGCGAGCCGAACCTTAGTCGAAGCGTCTCATCCTTCCCCGCATGGGCAGAATATTTTATTGCGGATTCGCTTGAGAGTTCAGGGACtttt                                                                                                                                    >  2:742577/1‑150 (MQ=255)

                 tCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGGCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTGCCCGCATGGGCAGTATATTTTATTGGGGATTGGGTTGTGAGTTCAGGGTCTTTTg                                                                                                                                   >  2:3635612/1‑150 (MQ=14)

                 tCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCCGAATATTTGATTGCGGAGTCGGTGGGGAGTTCAGGGACTTTTg                                                                                                                                   >  2:2462474/1‑150 (MQ=11)

                 tCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGGTGGGGAGTTCAGGGACTTTTg                                                                                                                                   >  2:2959826/1‑150 (MQ=255)

                 tCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTg                                                                                                                                   <  1:722839/150‑1 (MQ=255)

                 tCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCGTCTCATCCTTCCCCGCATGGGCCGAATATTTGATTGCGGATTCGCTGGAGAGTTCCGGGACGTTTg                                                                                                                                   >  2:3013316/1‑150 (MQ=14)

                  cAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGa                                                                                                                                  <  1:2011556/150‑1 (MQ=255)

                   aGCGCTGCGCGTTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGaa                                                                                                                                 <  1:2099884/150‑1 (MQ=255)

                    gcgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGaaa                                                                                                                                <  1:3124109/150‑1 (MQ=255)

                    gcgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGaaa                                                                                                                                <  1:3600731/150‑1 (MQ=255)

                    gcgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGaaa                                                                                                                                <  1:3664495/150‑1 (MQ=255)

                    gcgcTGCGCGCTTCGCCCTTCGGGTCGGTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGGCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCCTGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGaaa                                                                                                                                <  1:1583218/150‑1 (MQ=255)

                      gcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCTGCATCGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGt                                                                                                                              <  1:2606700/150‑1 (MQ=255)

                        tGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGA                                                                                                                            <  2:2819197‑M1/150‑1 (MQ=255)

                            cgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGT                                                                                                                        <  1:2744684‑M1/150‑1 (MQ=255)

                               ttCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCGTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGT                                                                                                                     <  1:579786‑M1/150‑1 (MQ=255)

                               ttCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATTTTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGT                                                                                                                     <  1:4079417‑M1/150‑1 (MQ=255)

                               ttCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGT                                                                                                                     <  1:2447010‑M1/150‑1 (MQ=255)

                                  gCCCTTCGTGTCGTTGCCTGCGGCATCGCTCTCTCGCTGTCGCTCGAGCGTAACCTTTGTCGATGCTTCTCATCCTTCCCCTCTTGGGCTGAATATTTGATTGGGGATTCGATTGAGTTTTCTGGTACTTTTGAAAGTGATGGTGGTggg                                                                                                                  <  2:1872430/150‑1 (MQ=11)

                                  gCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGG                                                                                                                  <  2:1719397‑M1/150‑1 (MQ=255)

                                  gCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGG                                                                                                                  <  1:2691248‑M1/150‑1 (MQ=255)

                                       tcgggtcgggGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAG                                                                                                             <  1:303576‑M1/150‑1 (MQ=255)

                                       tCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCGTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAG                                                                                                             <  1:1483612‑M1/150‑1 (MQ=255)

                                       tCGGGTCGGTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAG                                                                                                             <  1:303587‑M1/150‑1 (MQ=255)

                                        gagagcgttgcCGCAGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTTAAAGTGATGGTGGTGGGGGAAGG                                                                                                            >  1:3543814‑M1/6‑150 (MQ=255)

                                          ggTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCACATCCTTCCCCGCCTGGGCAGGATAGTTGATTGCGGATTCGCTGGGGAGGTCAGGGACGTTGGGAAGGGGTGGTGGTGGGGGAAGGAT                                                                                                          >  2:2475567‑M1/1‑150 (MQ=255)

                                            tCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTa                                                                                                        <  2:3064663/150‑1 (MQ=255)

                                            tCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTa                                                                                                        <  1:2304185/150‑1 (MQ=255)

                                               ttGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTc                                                                                                     >  2:788759/1‑150 (MQ=255)

                                                 gCCTGCGGCAACGCTCTCTCGCTGCCGCTCGCGCCGAACGTGCGGGGTAGCGTCTCCTCCGTCCCCGCCTGGGCAGATTATTTGATTGCGGATTCGCTTGAGCGTTCCGTGACTTTTGCAAGTGATGGTGGTGGGGGAAGGATTACTCAg                                                                                                   <  2:3963305/150‑1 (MQ=14)

                                                  ccTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCATAGTCGAAGCTTCTCAAACTTCCCCGCCTGGGCAGAATATTTGATTGCGGATTCGCTTGTGATGTAAGGGACTGTTGAAATTGATGGTGGTGGGGGAAGGATTAATCAgc                                                                                                  >  1:2266291/1‑150 (MQ=14)

                                                  ccTGCGGCAACGCTCTCACGCTGACGCTCGAGCCCAACCATAGTCGAAGCATATCAACCTTCCCCGCATGGGCAGAATATTTTATTGCGGATTCGCTTGTGAGTTAAGGGACTTTTGAAATTGTTTGTGGTGGGGGAAGGATTAATAAgc                                                                                                  >  1:390794/1‑150 (MQ=14)

                                                        gCAACGATCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTgcg                                                                                            <  2:3211537/150‑1 (MQ=255)

                                                          aaCGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTcgaaccttcgaagtc                                                                                          <  1:2157057‑M1/150‑16 (MQ=255)

                                                           aCGCTCTCTCGCTGACGCTCGAGCCGAACCTCCGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTgcgcgc                                                                                         <  2:1835526/150‑1 (MQ=255)

                                                             gCTCTCTCGCTGACGCTCGAGCCGCACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGTCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCtt                                                                                       <  1:367628/150‑1 (MQ=255)

                                                             gCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGTGGGAAGGATTACTCACAGCTGCGCCCtt                                                                                       >  2:1116549/1‑150 (MQ=255)

                                                             gCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCACCTCTGCGCCCtt                                                                                       >  2:1784163/1‑150 (MQ=255)

                                                               tctctcGCTGCCGCTCCAGCCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGCTTGAGACTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCg                                                                                     <  2:1166525/150‑1 (MQ=11)

                                                               tctctcGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCg                                                                                     >  2:1315412/1‑150 (MQ=255)

                                                                  ctcGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTcgaaccttcgaagtcgatgacgg                                                                                  <  2:2029600‑M1/150‑24 (MQ=255)

                                                                     gccgccGCTCGAGCCGAACCTTCGACGAAGCTTCTCATCCCTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGATTTCAGGGTGTTTTGATTGTTTTGGTGTTGTGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTc                                                                               <  1:3561698/145‑1 (MQ=14)

                                                                     gCTGTCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTc                                                                               <  1:1470740/150‑1 (MQ=255)

                                                                     gCTGCCGCTCGAGCCGACCCTTCGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTc                                                                               <  1:1470688/150‑1 (MQ=255)

                                                                     gCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTc                                                                               >  2:2232018/1‑150 (MQ=255)

                                                                     gCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTc                                                                               >  2:113840/1‑150 (MQ=255)

                                                                     gCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTc                                                                               >  2:1445884/1‑150 (MQ=255)

                                                                           gCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCg                                                                         >  1:119820/1‑150 (MQ=255)

                                                                           gCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCg                                                                         >  1:3725170/1‑150 (MQ=255)

                                                                           gCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCg                                                                         >  1:2455309/1‑150 (MQ=255)

                                                                           gCTCGAGCCGAACCCTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCATCGGGTCg                                                                         >  1:119780/1‑150 (MQ=255)

                                                                             tCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGtt                                                                       <  2:1744753/150‑1 (MQ=255)

                                                                             tCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTACGCGCTGCGCGCATCTCCCTTCGGGACGtt                                                                       >  2:1602086/1‑150 (MQ=17)

                                                                              ctaGCCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTg                                                                      <  2:2871147/148‑1 (MQ=255)

                                                                               gAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGGTGGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGAGCGGCGAGCTTCGCCCGTCGTGTAGATGc                                                                     >  2:308361/1‑150 (MQ=14)

                                                                                    gAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTTCGCGCTTCGCCCTTCGGGTCGTTGcaggcg                                                                >  1:2659202/1‑145 (MQ=255)

                                                                                    gAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCGTCGGGTCGTTGCCTGCg                                                                >  1:955547/1‑150 (MQ=255)

                                                                                    gAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGGTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTGGCCTGCg                                                                >  1:3842930/1‑150 (MQ=255)

                                                                                        cTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCaa                                                            >  1:67883/1‑150 (MQ=255)

                                                                                           aGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGc                                                         <  1:788759/150‑1 (MQ=255)

                                                                                              cgcaGCTGCTCATCCTTCCCCGCATGGGCAGAATATTGGAGTGCGGACTCGCTTGAGTGTTCCGGGACTTTTGCAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGctct                                                      <  2:3329077/147‑1 (MQ=11)

                                                                                               gAAGCTTCTCATCTTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGctctc                                                     >  2:4162295/1‑150 (MQ=255)

                                                                                               gAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGGAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGctctc                                                     >  2:3298251/1‑150 (MQ=255)

                                                                                               gAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGctctc                                                     >  2:530220/1‑150 (MQ=255)

                                                                                               ccaGCTTCCCACCCCCCCCCGCATGGGCAGAATATTCGATTGCGGATTCGCTTGAGAGTTCAGGGTCTTTTGAAAGTGATGGTTGTTGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGctctc                                                     <  1:1788174/148‑1 (MQ=14)

                                                                                                 aGCGTCTCACCCGTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCtctctc                                                   <  2:676224/150‑1 (MQ=255)

                                                                                                    ttCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGGTTGTGAGTTCAGGGAGTTTTGTAAGTGATGGTGGTGGGGGAAGTTTTTATCAGGGCTGCGAGCTTGGGGGTTGGGGTCGTTGCGTGCGGCAAGGctctctgggt                                                >  2:2735852/1‑146 (MQ=14)

                                                                                                    ttCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGGGAGTTCAGGGACTATGGAAAGTGAGGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTGCGCCCTTCGGGGCGTTGCCTGCGGCAACGCtcgctcgcg                                                >  2:4007458/1‑149 (MQ=14)

                                                                                                     tctcATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTg                                               >  2:1754531/1‑150 (MQ=255)

                                                                                                     tctcATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTg                                               >  2:1754212/1‑150 (MQ=255)

                                                                                                     tctcATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTg                                               >  1:1294738/1‑150 (MQ=255)

                                                                                                     tctcATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCGTCGGGTCGTGGCCTGGGGCAACGCGCTCTCGcgg                                               >  2:2234502/1‑148 (MQ=255)

                                                                                                         aTCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTTGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGGGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGgcgc                                           >  2:1332551/1‑150 (MQ=255)

                                                                                                            cTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCg                                        <  2:3592460/150‑1 (MQ=255)

                                                                                                              tCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAg                                      >  1:556918/1‑150 (MQ=255)

                                                                                                                 ccGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCTGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCg                                   >  2:3108789/1‑150 (MQ=255)

                                                                                                                     aTGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCAGCCGCAACCCACTCACGCACGCGCACCAGTCGAAcc                               >  2:3564687/1‑150 (MQ=14)

                                                                                                                         tcAGAATATTTGATTGCGGATTCGATTGAGATTTCAGGGTCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTgg                           <  1:1157647/149‑1 (MQ=255)

                                                                                                                         gCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTgg                           <  1:1157840/150‑1 (MQ=255)

                                                                                                                         gCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGAGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCGTgg                           >  2:1469605/1‑150 (MQ=255)

                                                                                                                           aGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCCACGCTCTCTCGCTGGCGCTCGAGTCGAACCGTGGTc                         >  1:2991642/1‑150 (MQ=255)

                                                                                                                           aGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACCCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTc                         >  1:2652432/1‑150 (MQ=255)

                                                                                                                             aaTATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCTCTGCGCGCTTCGCCCTACGGGTCGTTGCCTGAGGCAACGCACTCTCGCTGTCGCTCGAGTCGAACCTTGGTCGa                       >  2:322231/1‑150 (MQ=14)

                                                                                                                                atTTGATTGCGGATTCGCTTGAGAGTTCAGGGTCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGc                    <  2:1222281/150‑1 (MQ=255)

                                                                                                                                       tGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTcgaaccttcgaagtcgatgacggtagatttacagtctgctccctttggccgctcgggaaccccaccacggggtaatgctttttactggcctg             <  2:380306‑M1/150‑93 (MQ=255)

                                                                                                                                       tGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCAGCGCGCTTCGCCCATCGGGTCGTTGCCAGCGGCCACCCACACACCCAGGCGCACCACTCCAACCAAGGGCGAAGcaacacaa             >  1:900717/1‑143 (MQ=9)

                                                                                                                                           gATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCtt         <  1:77715/150‑1 (MQ=14)

                                                                                                                                           gATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCGTGGTCGAAGCTTCTCATCCtt         >  2:2879722/1‑150 (MQ=255)

                                                                                                                                                gCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGGTCTCTCGCTTGCGCTCGATTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCg    >  1:469582/1‑150 (MQ=255)

                                                                                                                                                gCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGACGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAATCGAACCATGGTCGAAGCTTCTCAACCGTccccc    >  1:630394/1‑149 (MQ=14)

                                                                                                                                                  ttGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCATGGTCGAACCTTCTCATCCATacccccc  >  1:2923236/1‑147 (MQ=14)

                                                                                                                                                  ttGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCACTCTCGCTGGCGCTCGAGTCGCACCATGGTCGAAGCTTCTCAACCTTACCCCCa  >  1:2923590/1‑150 (MQ=14)

GGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCA  >  NC_000913/1286892‑1287187

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