BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,286,944	T→G	5.3%	intergenic (+418/+122)	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,286,944	0	T	G	5.3%	428.6 / 6.7	155	intergenic (+418/+122)	narI/rttR	nitrate reductase 1, gamma (cytochrome b(NR)) subunit/rtT sRNA, processed from tyrT transcript
Reads supporting (aligned to +/- strand): ref base T (41/96); new base G (9/2); total (57/98)
Fisher's exact test for biased strand distribution p-value = 1.01e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 155 total aligned reads displayed.

TCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCA  >  NC_000913/1286796‑1287089

      cTTGGTCTTAGCTTCTCATCCTTCCCCGCTTGGGCAGATTATTTGATTGCGGATTAGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCaa                                                                                                                                            <  2:3453689/150‑1 (MQ=255)

          gTCTAAGCTTATCATCCTTCCCCGCTTGGGCATATTATTTGATTGCGTATTCGTTTTAGATTTCGGGTTTTTTTGAAAGTTATGGTGGTGGGGTATGTTTTACTCAGCGCTGCGCGCTTCGCCCTTCGGTTCGTTGCCGCAGGCAACgac                                                                                                                                        <  2:706259/150‑3 (MQ=2)

          gTCGAAGCTTCTCATCCTGCCCCGCTTGGGCAGAATTTTTGTTTGCGGATTCGGTTGAGAGTTCCGGGTCTTTTGAATGTGATGGTGGTGGGGGATGGATTACTCAGCGCGGCGCGCTTCGCCCTTCGGGTCTTTGCCGCAGGCAACgac                                                                                                                                        <  2:1386309/150‑3 (MQ=11)

           tCGAAGCTTCTCATACTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACAAAGCGCTGCGCGCTACGCCCTACGCGTCCATGCCTGCGGCACCGctc                                                                                                                                       >  1:2379852/1‑150 (MQ=14)

            cGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCgg                                                                                                                                               >  1:3929375/1‑141 (MQ=255)

            cGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGctct                                                                                                                                      >  1:2980229/1‑150 (MQ=255)

             gAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGGCTGGGGGAACGctctg                                                                                                                                     >  2:1935693/1‑149 (MQ=255)

             gAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTAATCAGAGCTGCGCGCTTCGCCCTTCGGGTCTTTGCCGCCGGGAACGcgaac                                                                                                                                     >  2:603528/1‑146 (MQ=14)

             gAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTAATCAGAGCAGCGAGCTTCGCCCTTCGTGTCGATGCCTCCGGCAACgaaaac                                                                                                                                     >  2:333194/1‑145 (MQ=14)

             gAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATAGTAGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGAGGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCGTCGCCCTTCGGGTCGTTGCCTGCGGCCCCGcgctc                                                                                                                                     >  2:822093/1‑150 (MQ=14)

               aGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAAGTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCtctctc                                                                                                                                   <  1:3628215/150‑1 (MQ=255)

               aGCTTCTCATCCTTCACCGCTTGGGCAGAATCTTTCATTGCGGATTCTTTTGAGAAGTGCTGGGCTTTTGAAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCtctctc                                                                                                                                   <  1:2422304/150‑1 (MQ=14)

                gCTTCTCATCCGTCCCCGCTTGGGCAGAATTTTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCg                                                                                                                                  <  2:3377111/150‑1 (MQ=255)

                gCTTCTCAGCCGTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTAGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCg                                                                                                                                  <  2:621667/150‑1 (MQ=255)

                ctgagtaaTCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCg                                                                                                                                  >  2:4016346/8‑150 (MQ=255)

                 cTTCTCATCCGTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGCAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGc                                                                                                                                 <  1:3833537/150‑1 (MQ=255)

                    ctccgccTTCCGCTCTTGGGCAGAATACTTGATTCCGGATTAGTTTGAGAATTCCGGGTCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGa                                                                                                                              <  2:1102504/145‑1 (MQ=14)

                    ctcATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTAGCAGa                                                                                                                              >  1:537029/1‑150 (MQ=255)

                      cctCCGTCCCCGCTTGGGCAGAATTTTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACg                                                                                                                            <  2:3153642/148‑1 (MQ=255)

                         cgttCCCCGCTTGGGCATAATGTTTCATTTCGTATGAGTTCGAGTATTCCTGTGCTTTTGTAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTc                                                                                                                         <  2:167485/148‑1 (MQ=14)

                           ttCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGa                                                                                                                       >  1:311251/1‑150 (MQ=255)

                           ttCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACGCTCACTCCCCGAACCTCGa                                                                                                                       >  1:495906/1‑150 (MQ=14)

                           gtCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGa                                                                                                                       <  1:3735948/149‑1 (MQ=255)

                             ccgctcTTTGGCAGATTTTTTGATTGGGTATTCGTTTGAGTATTCGGGTGCTTTTGTAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGc                                                                                                                     <  2:3602235/145‑1 (MQ=14)

                                cGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGa                                                                                                                  <  2:64251/150‑1 (MQ=255)

                                cGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGa                                                                                                                  <  2:64490/150‑1 (MQ=255)

                                 gCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCACTTCGCCCTTCGGGTCGTTGCCGGCGGCAACCCACGCGCGCAGACGCACGAGCCGac                                                                                                                 >  1:2016777/1‑149 (MQ=14)

                                 gCTTGGGCAGAATATTTGATTGCGGATTCGTTGGAGAATTCCGGGGCGTTTGAAAGGGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCGCGCTGACGCTCCAGCCGaa                                                                                                                 >  1:281445/1‑150 (MQ=14)

                                  cTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAAc                                                                                                                <  1:3815642/150‑1 (MQ=255)

                                     tggCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGCCTTTAGAAAGTGATGGAGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCtt                                                                                                             <  2:1047845/149‑1 (MQ=255)

                                        cAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGt                                                                                                          <  2:321336/150‑1 (MQ=255)

                                        cAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCACCGCTGCGCGCTTCGCCCATCGGGTCGTTGCCAGCGGCCACCCTCACACCCAGACCCGCCAGCCCAAcccaaga                                                                                                          >  1:3667924/1‑145 (MQ=14)

                                        cAGAATATTTGATTGCGGATTCGGTTGAGAATTCCGGGGCTTTAGACAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTGCGGGTCGTTGCCTGCGGCCACGCACTCTCCCTCACGCGCGCGCCGCACCGTAGt                                                                                                          >  2:2718712/1‑150 (MQ=14)

                                         aGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGGCTGCGGCAACGCTCTCACGCAGACGCACGAGCCCCACCATAGTc                                                                                                         >  1:378401/1‑150 (MQ=14)

                                               tttGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCt                                                                                                   >  2:4095663/1‑150 (MQ=255)

                                               tttGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCATCGTCGAAGCt                                                                                                   >  2:2082937/1‑150 (MQ=255)

                                                ttGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGAGCGCTTCGCCCTTCGGGTCGTTGCCAGCGGCAACCCTCTATAGCTGACGCACGAGCCGAACCTTAGTCGAAgggt                                                                                                  >  1:2569263/1‑147 (MQ=14)

                                                    ttGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCa                                                                                              <  1:2993444/150‑1 (MQ=255)

                                                    ttGCGGATTCGTTAGAGTATTCGGGGTCTTTTGCAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCa                                                                                              <  1:227717/150‑1 (MQ=14)

                                                     tggggATTCGTTTGAGTATTCCTGGGCTTTTGTAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCAt                                                                                             <  2:3971464/147‑1 (MQ=255)

                                                     tGCGGATTCGTTTGAGAATTCCGTGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCAt                                                                                             <  2:3647983/150‑1 (MQ=255)

                                                     tGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCAt                                                                                             <  2:2629589/150‑1 (MQ=255)

                                                     tGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCAt                                                                                             <  2:2869412/150‑1 (MQ=255)

                                                      gCGGATTCGTTTGAGAATTCCGGGTCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATc                                                                                            <  1:2656956/150‑1 (MQ=255)

                                                       gggATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATcc                                                                                           <  1:721219/149‑1 (MQ=255)

                                                       cGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATcc                                                                                           <  2:2485117/150‑1 (MQ=255)

                                                        ggATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCACCGCTGCGAGCTTCTCCCTTCGGGTCGTTGCCTGCCCCAAAGCTCTCTCGCTGACCCTCGAGGCGAACATGGGTCGAAGCTACTCATCCt                                                                                          >  2:239307/1‑150 (MQ=14)

                                                          ctgcGTTTGAGTATTGCTGTGCTTTTGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTc                                                                                        <  2:2542935/147‑1 (MQ=14)

                                                                   agtaTTCCGGGGCTTTAGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATgg                                                                               <  2:3607757/147‑1 (MQ=255)

                                                                    gaATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATggg                                                                              <  1:1935693/150‑1 (MQ=255)

                                                                     aaTTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGc                                                                             <  2:2228922/150‑1 (MQ=255)

                                                                          cGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCACCTCTGCGACCTTCGCCCATCTGTTCGTTGACAGGGCCACCGCACTCTCGCGAACGCTAGAGCCGAACCGTATACGAATCTTATCAACCATCCCCGATAGGGCAGAat                                                                        >  2:4101242/1‑150 (MQ=9)

                                                                             ggCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATAtt                                                                     <  2:493324/150‑1 (MQ=38)

                                                                             ggCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCACTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATAtt                                                                     <  1:2320100/150‑1 (MQ=38)

                                                                              gCTTTTGACAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATAttt                                                                    <  2:4220083/150‑1 (MQ=38)

                                                                              gCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATAttt                                                                    <  2:1870855/150‑1 (MQ=38)

                                                                              gCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCACTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTGCTCATCCTTCCCCGCATGGGCCGAAGAttt                                                                    >  2:3538743/1‑150 (MQ=255)

                                                                                ttttGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGa                                                                  <  1:2446283/150‑1 (MQ=38)

                                                                                 tttGTAAGGGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGAGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCTCTGACGCTCGAGTAGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGAt                                                                 <  1:2770446/150‑1 (MQ=11)

                                                                                 tttGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAAAGCGCTCTCGCTGACCCACGAGCCCAACCTTTGTCGAAGCTTCTCATCCTACCCCCCAAGGGCAGAatggttaat                                                                 >  2:2177758/1‑143 (MQ=11)

                                                                                  ttGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACCCACAAACGCAGACGCACCAGCCGAACCATAGTAGAAGCATCAAAAACGTACCCCCCTGGGCAGAATATTTTAtt                                                                >  1:126612/1‑150 (MQ=2)

                                                                                    gAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTc                                                                                                                         <  1:2764426/91‑1 (MQ=33)

                                                                                    gAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGc                                                              <  1:3647888/150‑1 (MQ=39)

                                                                                     aaaGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCg                                                             <  1:4016346/150‑1 (MQ=39)

                                                                                       aGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCAGCGGCAACCCACACTCGCTGAAGCGCGAGCCCAACCATATTCGAAGCTTCTCATCCATCCCCCCAAGGGCCGAAAAATTGTTTGCGGa                                                           >  1:701660/1‑150 (MQ=14)

                                                                                         cgATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGAtt                                                         <  2:293826/149‑1 (MQ=39)

                                                                                          gATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTTCCGGACGCAACGCACTCTCGCTGACGCGCCAGCCGAACCTTAGTCGAAGCTTCTAATACTTCACCGCATGGGCAGAATATTTAATTGCGGATTc                                                        >  2:758737/1‑150 (MQ=11)

                                                                                          gATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTc                                                        <  1:2136074/150‑1 (MQ=39)

                                                                                          gATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTc                                                        <  1:2136024/150‑1 (MQ=39)

                                                                                              gtggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCtt                                                    <  2:240291/150‑1 (MQ=255)

                                                                                                ggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGCCGCTCCAGCCGAACCTTCGTCGAAGCTTCTCATCCTTCCCCCCCTGGGCCGAATATTTGATTGCGGCTTCGCTTga                                                  >  2:1020160/1‑150 (MQ=11)

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

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

                                                                                                ggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCACACTCGCAGACGCTCGCGCCGAACCCTAGTCGAAGCTTCACAACCTACACCCCATGGGCAGAATAAGTGATTGCGGAGTCgttggg                                                  >  1:2460224/1‑145 (MQ=11)

                                                                                                 gggggggAAGGATTACTCGGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTgag                                                 <  1:3703409/148‑1 (MQ=255)

                                                                                                     gggAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGtt                                             <  1:1020160/150‑1 (MQ=255)

                                                                                                     gggAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGtt                                             <  1:4095663/150‑1 (MQ=255)

                                                                                                     gggAAGGATGACGCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGtt                                             <  1:485584/150‑1 (MQ=255)

                                                                                                       gAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGACGCTTCTCATCCTTCCCCGCATGGGCAGACTATTTGATTGCGGATTCGGTGGagaggtgc                                           >  2:2433230/1‑146 (MQ=14)

                                                                                                       gAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCgttg                                                    >  2:2433151/1‑138 (MQ=37)

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

                                                                                                          ggATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAggg                                        <  1:821571/150‑1 (MQ=255)

                                                                                                           gATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGa                                                                                                                       >  2:3675925/1‑70 (MQ=35)

                                                                                                                 tCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACg                                                                                                                            >  2:3661739/1‑59 (MQ=35)

                                                                                                                   aGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCTAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGaa                               <  2:2104264/150‑1 (MQ=255)

                                                                                                                   aGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCGTGAGAGTTCAGGGACTTTTGaa                               >  2:2657284/1‑150 (MQ=255)

                                                                                                                       cTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTg                           <  2:1292372/150‑1 (MQ=255)

                                                                                                                         gcgcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGGTTGGGAGTTCAGGGACTTTTGAAAGTGAt                         >  2:2409534/1‑150 (MQ=255)

                                                                                                                            cgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGCCGAACCTTCGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGAtggt                      <  2:3792163/150‑1 (MQ=255)

                                                                                                                            cgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGAtggt                      <  2:545547/150‑1 (MQ=255)

                                                                                                                             gcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGAtggtg                     <  1:1899545/150‑1 (MQ=255)

                                                                                                                               ttCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGAtggtggt                   <  1:2657284/150‑1 (MQ=255)

                                                                                                                                   cccTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTgggg               <  1:2612001/150‑1 (MQ=255)

                                                                                                                                    ccTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTggggg              <  1:2278312/150‑1 (MQ=255)

                                                                                                                                     cTTCGGGTCGTTGCCTGCGGCAACGCTCTAGCGCTGACGCTCGGGCCGAACCTTAGTCGAAGCATATCAAACTTACCCCCAAGGGCAGAATATTTGAATTCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGa             >  2:3314638/1‑150 (MQ=17)

                                                                                                                                      ttCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCATAGTCGAAGCTTCTCATCCTTCCCCCCAGGGGCAGAATATTTGATTGCGGATTCGCTTGTGAGTTCAGGGACTTTTGAAACTGAAGGTGGTGGGGGac            >  2:1262537/1‑149 (MQ=14)

                                                                                                                                        gagagcGTTGCCGCAGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAgg          >  1:81497/6‑150 (MQ=255)

                                                                                                                                         gggTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCAGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGa         <  1:333194/150‑1 (MQ=255)

                                                                                                                                                tGCCTGCGGCAACGCTCTCTCGCTGCCGCTCGAGCCGACCCTTCGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCa  <  2:2970813/150‑1 (MQ=255)

                                                                                                                                                tGCCTGCGGCAACGCTCTCTCGCTGCCGCTCCCGCCGAACCGGGGTCGTAGCGTCTGATCGTTCCCCGGATGGGCAGATTATTTGATTTCGTATTCGCTTGAGCGTTCCGGGACTTTTGAAAGTGATAGTGGTGGGGGAAGGATTACTCa  <  2:1147713/150‑1 (MQ=11)

TCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCA  >  NC_000913/1286796‑1287089

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