BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,286,965	A→T	7.3%	intergenic (+439/+101)	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,965	0	A	T	7.3%	945.9 / 7.4	362	intergenic (+439/+101)	narI/rttR	nitrate reductase 1, gamma (cytochrome b(NR)) subunit/rtT sRNA, processed from tyrT transcript
Reads supporting (aligned to +/- strand): ref base A (136/175); new base T (10/20); total (158/204)
Fisher's exact test for biased strand distribution p-value = 3.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 367 total aligned reads displayed.

CTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGT  >  NC_000913/1286816‑1287114

ctcATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGa                                                                                                                                                       >  2:3887763/1‑150 (MQ=255)

ctcATCCGTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGa                                                                                                                                                       <  1:237163/150‑1 (MQ=255)

 tcATCCTTCCCCGCTTGGGCAGAATATTTGATGGCGGATTCGTTTGAGAATTCGGGGGCTTTTGAAAGTGTTGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCTACGCTCTCTCGCTGAc                                                                                                                                                      <  2:11371696/150‑1 (MQ=255)

       ttCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGGCTGCGGCAACGCTCTCTCGCTGACGCACgg                                                                                                                                                >  1:11622035/1‑149 (MQ=255)

       ttCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGa                                                                                                                                                >  1:9964509/1‑150 (MQ=255)

       ttCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGa                                                                                                                                                >  1:4634196/1‑150 (MQ=255)

       ttCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCGCTCGCTGACGCTCgg                                                                                                                                                >  1:5255771/1‑149 (MQ=255)

        tCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCAATAGCAGAAGCTCGAg                                                                                                                                               >  2:5573437/1‑150 (MQ=14)

        tCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACGCTCTCTCGCTGACGCTCGAg                                                                                                                                               >  1:2888583/1‑150 (MQ=255)

        tCCCCGCTGGGGCAGAATATTTGATGGCGGAGTCGTGTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAg                                                                                                                                               <  2:2352000/150‑1 (MQ=255)

         ccccGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGTCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGc                                                                                                                                              <  2:7563673/150‑1 (MQ=255)

         ccccGCGTGGGCAGAATCCTAGATTGGGGCTTCGTTTGAGAATTCCGGTGCTTTTGAAAGAGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGc                                                                                                                                              <  2:2878120/150‑1 (MQ=14)

          cccGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGcc                                                                                                                                             <  1:4971656/150‑1 (MQ=255)

           ccGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCAGCGGCAACGCACTCACGCTGACGCTCGAGCCg                                                                                                                                            >  1:5127395/1‑150 (MQ=255)

                   gCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGGCTGCGGCAACGCACTCTCGCTGACGCACGAGCCGAACCTTat                                                                                                                                    >  1:3290369/1‑149 (MQ=255)

                    cAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCCACCCGCACACCCACACGCACGAGCCCAAccaaata                                                                                                                                   >  1:1537210/1‑145 (MQ=14)

                    cAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAAACCTCTCACGCAGACCCTCGAGCCGAAccaaagt                                                                                                                                   >  1:9410435/1‑145 (MQ=14)

                     aGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGAGACCTTCGCCCTACGTGTCGATGCCGGCCCCAAAACACACTCGCAGCCGCTCGAGCCGAACCATAGTc                                                                                                                                  >  2:6010194/1‑150 (MQ=14)

                        ctatTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTGGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGaa                                                                                                                               <  1:1127973/149‑1 (MQ=255)

                        atatTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGaa                                                                                                                               >  2:12037106/1‑150 (MQ=255)

                         tatTTGATTGCGGATTCGTGTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAg                                                                                                                              <  1:12071183/150‑1 (MQ=255)

                          ttttGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGTAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGc                                                                                                                             <  1:7208691/149‑1 (MQ=255)

                             tGATTTCGGATTCGTTTGAGTATTGCGTGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTtc                                                                                                                          <  2:9668450/150‑1 (MQ=255)

                              gATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTtct                                                                                                                         >  2:7667522/1‑150 (MQ=12)

                              gATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTtct                                                                                                                         >  2:9955266/1‑150 (MQ=255)

                               atttggtaTTAGTTTGAGTATTGCTGGTGTTTAGTAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTtctc                                                                                                                        <  2:6539738/143‑1 (MQ=14)

                                ttGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCGGAGCGCTTCGCCCAGCGGGTCGTTGCAGGCGGCAAACCCCGCTCGCTGACTCTCGAGCCGAACCTTAGTCGAAGCTTCTCa                                                                                                                       >  2:10183315/1‑150 (MQ=14)

                                 tGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGATCGCATCGCACTACGGGTCGTTGCCTGCGGCAACGCTATATCGCTGACGCTCCAGCCGAAACATAGTCGAAGCTTCTCAt                                                                                                                      >  2:3784762/1‑150 (MQ=14)

                                    ggATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGGCGAACCTTAGTCGAAGCTTCTCATCCt                                                                                                                   >  2:338461/1‑150 (MQ=255)

                                           tttGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGACGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGc                                                                                                            >  2:2515834/1‑150 (MQ=255)

                                           tttGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCGCTCTCGCGGACGCACGAGCCGAACCTTAGTCGCAGCGTCTCATCCATCCCCGc                                                                                                            >  1:9661602/1‑150 (MQ=14)

                                           tttGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACCCACTATAGCTGACCCTCCAGCCGAACCTTAGTCGAAGCATCTCAAACTTcccccc                                                                                                            >  1:10822026/1‑148 (MQ=14)

                                           tttGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTAACAGCTGCGACCTTCTCCCTTCGGGTCCTTGCCGGCGGCAACGATCTCTCGCTGGCGCTCGAGCCGAACCGTAGTCGCAGCTTCTCATCCATCCCCGc                                                                                                            >  2:11060602/1‑150 (MQ=14)

                                             tGAGAATTCCGGGTCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCAt                                                                                                          <  1:5564228/150‑1 (MQ=255)

                                             tGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCAt                                                                                                          <  1:2223015/150‑1 (MQ=255)

                                                gaATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATggg                                                                                                       <  2:6126154/150‑1 (MQ=255)

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

                                                   ttCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGa                                                                                                                                                >  2:6669961/1‑106 (MQ=40)

                                                   ttCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGctg                                                                                                    >  2:9587201/1‑148 (MQ=255)

                                                   ttCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAg                                                                                                    <  2:10891024/150‑1 (MQ=255)

                                                       ggggCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTTCCTGCGGCAACGCTCTCTAGCGGAAGCACGACCCGAAACATAGTCGAAGCTTAACAAACATACCAGCAAGGGCAGaaaa                                                                                                >  1:7848346/1‑148 (MQ=14)

                                                       ggggCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCCACCCTCTCACGCAGACGCACGACACGAAACAGAGTCGAAGCATATCAAACATCCCCGCATGGGCCGaaga                                                                                                >  1:11196764/1‑148 (MQ=11)

                                                         ttcTTTAGCAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCACTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATAtt                                                                                              <  2:2913522/148‑1 (MQ=255)

                                                         ggCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTGGTCGAAg                                                                                                                              >  2:9820466/1‑118 (MQ=21)

                                                          gCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACCCTCTCTAGCTGACGCTCGAGCCGAACCATAGTCGAAGCTTCTCATCCGTACCCGCATGGGCAGAATAttt                                                                                             >  1:6975088/1‑150 (MQ=11)

                                                            ttaaggaaGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGa                                                                                           <  1:176537/144‑1 (MQ=255)

                                                            ttaaggaaGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGa                                                                                           <  2:1058380/144‑1 (MQ=255)

                                                            ttaaggaaGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCGTCCCCGCATGGGCAGAATATTTGa                                                                                           <  1:372513/144‑1 (MQ=255)

                                                                       aTGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCg                                                                                <  1:10237906/150‑1 (MQ=39)

                                                                         ggtggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCt                                                                              <  1:8825970/150‑1 (MQ=40)

                                                                         ggtggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCt                                                                              <  2:118368/150‑1 (MQ=40)

                                                                              tggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGTGTCGTTGCCTTCGGGAAAGCTCTCTCGCGTTCGCGCGGGGCGAACCATAGTCGAAGCTTATCATGCGTCCCCGCCGGGGCCGAATATTTGATTGCGggttgggttgggg                                                                         >  2:4298792/1‑147 (MQ=2)

                                                                               gggggAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCAGGGGCAGACTATATGGTTGGGGATTGGCGTGagag                                                                        >  2:1251898/1‑150 (MQ=11)

                                                                                ggggAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCACACTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCGTCTCAACCATCCCCGCCAGGGCAGAAAATGTGATTGCGGAGTCGCtggccagt                                                                       >  2:6751256/1‑143 (MQ=14)

                                                                                  ggAAGGCTTACTCCGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGAGTCGCTTGAGAGTTc                                                                     <  1:9850073/150‑1 (MQ=255)

                                                                                   gAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCCGAATATTTGATTGCGGATTCGCGTGGGAGTTCa                                                                    >  2:8819239/1‑150 (MQ=255)

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

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

                                                                                   gAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCGCGAGCCGAACCTGAGTCGAAGCTTCTCATCCTTCCCCGCCTGGGCAGAATAGTTGATTGCGGAGTCGCGGGTGAGTTcc                                                                    >  2:8396390/1‑149 (MQ=14)

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

                                                                                      ggATTACTCAGCGCTGCGCGCTTCTCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGACGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAggg                                                                 <  1:7699495/150‑1 (MQ=255)

                                                                                       gATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGGTGGGGAGTTCAGGGa                                                                >  2:8958899/1‑150 (MQ=255)

                                                                                         ttACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGACTATTTGATTGCGGCGTCGGTTGGGAGTTCAgggccg                                                              >  2:10056609/1‑147 (MQ=14)

                                                                                          tACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACtt                                                             <  1:10376491/150‑1 (MQ=255)

                                                                                            cTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACtttt                                                           <  1:10078220/150‑1 (MQ=255)

                                                                                               aGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGaa                                                        >  2:7904124/1‑150 (MQ=255)

                                                                                               aGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGaa                                                        >  2:10786065/1‑150 (MQ=255)

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

                                                                                                    tGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCCAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGa                                                   <  2:1431843/150‑1 (MQ=255)

                                                                                                     gcgcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGAt                                                  <  1:10970984/150‑1 (MQ=255)

                                                                                                      cgcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCATGGGCCGAATATTTGATTGCGGATTCGGTTGCGAGTTCAGGGACTTTTGAAAGTGAtg                                                 >  2:5932064/1‑150 (MQ=255)

                                                                                                         gcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCGTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGGTTGGGAGTTCAGGGACTTTTGAAAGTGAtggtg                                              >  2:8507263/1‑150 (MQ=255)

                                                                                                            tCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGAtggtggtg                                           <  1:6002209/150‑1 (MQ=255)

                                                                                                              gCCCTTCGGGTCGTTGCCTGCGGCAACGCTCCCTCGCTTACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCAACCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTggg                                         >  2:1724254/1‑150 (MQ=255)

                                                                                                               cccTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCTAGTCGAACCTTTGTCGAAGATTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTgggg                                        <  2:6931291/150‑1 (MQ=11)

                                                                                                                ccTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATGTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTggggg                                       >  2:4123798/1‑150 (MQ=255)

                                                                                                                  ttCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGCCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGaa                                     <  1:2993313/150‑1 (MQ=255)

                                                                                                                  ttCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGaa                                     <  2:1492183/150‑1 (MQ=255)

                                                                                                                  ttCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGaa                                     >  2:800236/1‑150 (MQ=255)

                                                                                                                  gtCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGaa                                     <  1:4333152/149‑1 (MQ=255)

                                                                                                                   tCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAg                                    <  1:3824660/150‑1 (MQ=255)

                                                                                                                    gagagcGTTGCCGCAGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAAGATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAgg                                   >  2:5848855/6‑150 (MQ=17)

                                                                                                                    cGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGCCGCTCCAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAgg                                   <  2:396641/150‑1 (MQ=255)

                                                                                                                       gtagTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGAtt                                <  2:2156516/147‑1 (MQ=255)

                                                                                                                       agcGTTGCCGCAGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGACGCTTCTCATCCTTCCCCGCATGGGCAGAATATTGGATTGCGGATTCGCTGGAGAGTTCAGGGACTGGTGAAAGTGCTGGTGGTGGGGGCAGGagt                                >  2:10328354/3‑148 (MQ=14)

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

                                                                                                                               cTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCATAGTCGAAGCTTCTCATCCATCCCCGCAAGGGCAGAATATTTGATTGCGGCTTCGCTTGAGAGTTCAGGGACTTTTGAAAGTTATGGTGGTGGGGGAAGGATTACTCAgcg                        >  1:6427729/1‑150 (MQ=17)

                                                                                                                                tGCGGCAACGCTCTCACGCTGACGCTCGAGCCGAACCATAGTCGAAGCTTCAAAACCATACCCGCATGGGCAGAAAATTTTATTGCGGTTTCGCTTGTGAGTTCAGGGACTTTTGTAAGTTATGGTGGTGGGGGCAGGATTACTCAgctc                       >  1:480053/1‑148 (MQ=255)

                                                                                                                                     cAACGCTCTCTCGCTGTCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTgcgc                  <  2:2833346/150‑1 (MQ=255)

                                                                                                                                       aCGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATGGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTgcgcgc                <  2:7097558/150‑1 (MQ=255)

                                                                                                                                         gCTCTCTCGCTGTCGCTCTAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCtt              <  1:8444570/150‑1 (MQ=255)

                                                                                                                                         gCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCtt              >  2:5416022/1‑150 (MQ=255)

                                                                                                                                         gCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTACCCGCATGGGCAGAATAATTAATTGCGGATTAGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTAATAAGCTCAgcgaccat              >  1:11967528/1‑145 (MQ=14)

                                                                                                                                             tctcGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCAGCCTTCCCCGCATGGGCAGATTCTTTTATTTCGGATTCGCTTGAGAGTTCAGGTACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGcc          <  2:7391249/150‑1 (MQ=17)

                                                                                                                                                 gCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCTCCCTTc      >  1:3905359/1‑150 (MQ=255)

                                                                                                                                                 gCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGAGGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTc      >  2:5272957/1‑150 (MQ=255)

                                                                                                                                                   tGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCgg    >  2:4502930/1‑150 (MQ=255)

                                                                                                                                                    gACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCggg   <  2:3567718/150‑1 (MQ=255)

                                                                                                                                                     aCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCACTTCGGGt  >  1:7973847/1‑150 (MQ=255)

                                                                                                                                                     aCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTGTTGAGAGAGAACGTGTTGGGGGAAGGATAACTCAGCGCTGCGCGCTTCGCCCTTCGGGt  >  1:6796056/1‑150 (MQ=17)

CTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGT  >  NC_000913/1286816‑1287114

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