BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	AP012306	826,737	+GT	coding (177/285 nt)	ihfB →	integration host factor (IHF), DNA‑binding protein, beta subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AP012306	826,735	1	.	G	96.6%	1256.8 / 30.9	323	coding (175/285 nt)	ihfB	integration host factor (IHF), DNA‑binding protein, beta subunit
Reads supporting (aligned to +/- strand): ref base . (7/4); new base G (118/193); total (125/198)
Fisher's exact test for biased strand distribution p-value = 1.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01
*	AP012306	826,735	2	.	T	96.6%	1279.4 / 29.1	324	coding (175/285 nt)	ihfB	integration host factor (IHF), DNA‑binding protein, beta subunit
Reads supporting (aligned to +/- strand): ref base . (7/4); new base T (118/195); total (125/199)
Fisher's exact test for biased strand distribution p-value = 1.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 336 total aligned reads displayed.

GACTTGCCACCCAGCAATCGCACATTCCCGCCAAGACGGTTGAAGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCAC‑‑GTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAA  >  AP012306/826586‑826874

||

gACTTGCCACCCAGCAATCGCACATTCCCGCCAAGACGGTTGAAGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGcaccac                                                                                                                                               >  2:1304584/1‑150 (MQ=255)

    tGCCACCCAGCAATCGCACATTCCCGCCAAGACGGTTGAAGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCAC‑‑Gtgt                                                                                                                                         >  2:759853/1‑148 (MQ=255)

     gCCACCCAGCAATCGCACATTCCCGCCAAGACGGTTGAAGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCAC‑‑Gtgta                                                                                                                                        >  2:4227698/1‑147 (MQ=255)

      ccccccAGCAACCGCACATTCCCGCCAAGACGGTTGAAGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCAC‑‑Gtgtac                                                                                                                                       <  2:1594183/147‑5 (MQ=255)

         cccAGCAATCGCACATTCCCGCCAAGACGGTTGAAGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCAC‑‑Gtgtaccgg                                                                                                                                    >  2:399655/1‑143 (MQ=255)

            agcacccGCACTTTCCCGCCAAGACGGTTGAAGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACg                                                                                                                                   <  2:908977/144‑1 (MQ=255)

            aGCAATCGCACATTCCCGCCAAGACGGTTGAAGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACg                                                                                                                                   >  1:581233/1‑150 (MQ=255)

              cAATCGCACATTCCCGCCAAGACGGTTGAAGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTa                                                                                                                                 <  1:2240368/150‑1 (MQ=255)

                aTCGCACATTCCCGCCAAGACGGTTGAAGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAAt                                                                                                                               >  1:1635216/1‑150 (MQ=255)

                    cacaTTCCCGCCAAGACGGTTGAAGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGa                                                                                                                           <  1:1817894/150‑1 (MQ=255)

                     acaTTCCCGCCAAGACGGTTGAAGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGaa                                                                                                                          >  1:531710/1‑150 (MQ=255)

                       aTTCCCGCCAAGACGGTTGAAGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGa                                                                                                                        >  2:1199321/1‑150 (MQ=255)

                         tCCCGCCAAGACGGTTGAAGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACt                                                                                                                      >  2:2951947/1‑150 (MQ=255)

                         tCCCGCCAAGACGGTTGAAGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACt                                                                                                                      <  2:2416279/150‑1 (MQ=255)

                          cccGCCAAGACGGTTGAAGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTg                                                                                                                     >  2:801310/1‑150 (MQ=255)

                            cGCCAAGACGGTTGAAGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGc                                                                                                                   >  1:3501239/1‑150 (MQ=255)

                            cGCCAAGACGGTTGAAGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGc                                                                                                                   <  1:2293565/150‑1 (MQ=255)

                            cGCCAAGACGGTTGAAGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGc                                                                                                                   <  1:1383612/150‑1 (MQ=255)

                                  gACGGTTGAAGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATaaa                                                                                                             >  1:2353817/1‑150 (MQ=255)

                                    cGGTTGAAGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTTGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGt                                                                                                           <  2:3508594/150‑1 (MQ=255)

                                    cGGTTGAAGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGt                                                                                                           <  2:1441018/150‑1 (MQ=255)

                                      gTTGAAGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAg                                                                                                         >  1:1737009/1‑150 (MQ=255)

                                      gTTGAAGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAg                                                                                                         <  2:1337493/150‑1 (MQ=255)

                                       ttGAAGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGa                                                                                                        <  1:3146040/150‑1 (MQ=255)

                                       ttGAAGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGa                                                                                                        <  1:1245925/150‑1 (MQ=255)

                                         gAAGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAAc                                                                                                      <  1:1142528/150‑1 (MQ=255)

                                          aaGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACt                                                                                                     <  1:774971/150‑1 (MQ=255)

                                          aaGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACt                                                                                                     <  1:3017361/150‑1 (MQ=255)

                                             aTGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTgga                                                                                                  >  2:3750839/1‑150 (MQ=255)

                                             aTGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTgga                                                                                                  <  2:1543300/150‑1 (MQ=255)

                                             aTGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTgga                                                                                                  >  2:266095/1‑150 (MQ=255)

                                                cAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTggaagg                                                                                               <  1:604776/150‑1 (MQ=255)

                                                    aaaaGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGaaaa                                                                                           >  2:4365362/1‑150 (MQ=255)

                                                    aaaaGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGaaaa                                                                                           >  2:3438510/1‑150 (MQ=255)

                                                      aaGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATa                                                                                         <  2:2919635/150‑1 (MQ=255)

                                                       agagATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATAc                                                                                        <  1:2791186/150‑1 (MQ=255)

                                                          gaTGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGtt                                                                                     <  2:2269281/150‑1 (MQ=255)

                                                           aTGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTc                                                                                    >  1:295337/1‑150 (MQ=255)

                                                            tGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTcc                                                                                   <  1:2417281/150‑1 (MQ=255)

                                                             gCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCt                                                                                  <  1:3138783/150‑1 (MQ=255)

                                                             gCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCt                                                                                  >  2:1838318/1‑150 (MQ=255)

                                                              cTGGAGCATATGGCCTCGACTCTTGCGCAGGGCAAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGGAGAACTGGAAGGAAAATACGTTCCTc                                                                                 >  2:3887448/1‑150 (MQ=255)

                                                                 gAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACt                                                                              >  1:3397245/1‑150 (MQ=255)

                                                                  aGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACtt                                                                             >  1:3206225/1‑150 (MQ=255)

                                                                  aGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACtt                                                                             >  1:1008664/1‑150 (MQ=255)

                                                                   gCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACttt                                                                            >  2:713679/1‑150 (MQ=255)

                                                                        tGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAAcc                                                                       <  1:3777374/150‑1 (MQ=255)

                                                                        tGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAAcc                                                                       >  1:2626522/1‑150 (MQ=255)

                                                                         ggCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCt                                                                      >  1:2540222/1‑150 (MQ=255)

                                                                         ggCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCt                                                                      >  2:4079617/1‑150 (MQ=255)

                                                                          gCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTg                                                                     <  1:1898496/150‑1 (MQ=255)

                                                                             tCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTa                                                                  <  1:3712854/150‑1 (MQ=255)

                                                                              cGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACg                                                                                       >  2:114775/1‑128 (MQ=255)

                                                                               gACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTaaa                                                                >  2:1682135/1‑150 (MQ=255)

                                                                               gACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTaaa                                                                <  2:2745807/150‑1 (MQ=255)

                                                                               gACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTaaa                                                                <  2:2874476/150‑1 (MQ=255)

                                                                               gACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTaaa                                                                <  2:3237721/150‑1 (MQ=255)

                                                                                aCTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAg                                                               <  1:1591073/150‑1 (MQ=255)

                                                                                 ctctTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGa                                                              >  2:1423729/1‑150 (MQ=255)

                                                                                  tctTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGaa                                                             <  1:4164773/150‑1 (MQ=255)

                                                                                  tctTGCGCAGGGCGAGCGTATGAAAACCCTCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGaa                                                             <  2:3269110/150‑1 (MQ=255)

                                                                                    ttGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACt                                                           <  2:1479573/150‑1 (MQ=255)

                                                                                    ttGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACt                                                           <  2:3189016/150‑1 (MQ=255)

                                                                                    ttGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACt                                                           <  1:1187023/150‑1 (MQ=255)

                                                                                       cgcAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTgcg                                                        >  2:2681465/1‑150 (MQ=255)

                                                                                       cgcAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTgcg                                                        <  2:2014555/150‑1 (MQ=255)

                                                                                         cAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGcgcg                                                      >  1:1302741/1‑150 (MQ=255)

                                                                                            ggCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATc                                                   <  2:3100966/150‑1 (MQ=255)

                                                                                            ggCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATc                                                   <  2:1963539/150‑1 (MQ=255)

                                                                                              cGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATcgc                                                 <  2:2802227/150‑1 (MQ=255)

                                                                                               gAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATcgcg                                                >  1:395353/1‑150 (MQ=255)

                                                                                                 gCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACt                                                                              >  2:1780934/1‑118 (MQ=255)

                                                                                                      ttGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAatat                                         <  2:1261363/150‑1 (MQ=255)

                                                                                                           aTCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACg                                    <  1:2305502/150‑1 (MQ=255)

                                                                                                           aTCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACg                                    >  1:990901/1‑150 (MQ=255)

                                                                                                              cgcgGTTTCGGCAGTTTCTCTTTGTACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGtt                                 <  2:1655068/150‑1 (MQ=255)

                                                                                                              cgcgGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGtt                                 <  2:2176141/150‑1 (MQ=255)

                                                                                                                cgGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTaa                               <  1:1490858/150‑1 (MQ=255)

                                                                                                                 ggTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAg                              <  1:1987977/150‑1 (MQ=255)

                                                                                                                  gTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGt                             >  2:520016/1‑150 (MQ=255)

                                                                                                                   tttCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGtt                            <  1:3076390/150‑1 (MQ=255)

                                                                                                                   tttCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGtt                            <  1:344736/150‑1 (MQ=255)

                                                                                                                    ttCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGttt                           <  1:434122/150‑1 (MQ=255)

                                                                                                                     tCGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGtttt                          <  1:2580003/150‑1 (MQ=255)

                                                                                                                      cGGCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGttttt                         <  1:4021195/150‑1 (MQ=255)

                                                                                                                        gCAGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTa                       >  2:409869/1‑150 (MQ=255)

                                                                                                                          aGTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACt                     <  1:945405/150‑1 (MQ=255)

                                                                                                                           gTTTCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTc                    <  2:1158936/150‑1 (MQ=255)

                                                                                                                            tttCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCa                   <  1:3899613/150‑1 (MQ=255)

                                                                                                                             ttCTCTTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCaa                  >  2:3655973/1‑150 (MQ=255)

                                                                                                                               ctctTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAAc                <  1:2841455/150‑1 (MQ=255)

                                                                                                                                tctTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACt               >  1:2201703/1‑150 (MQ=255)

                                                                                                                                 ctTTGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACtt              <  1:1084387/150‑1 (MQ=255)

                                                                                                                                  tttGCACTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTg             >  2:628704/1‑150 (MQ=255)

                                                                                                                                       aCTACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACga        <  1:424474/150‑1 (MQ=255)

                                                                                                                                         tACCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACgaga      <  2:1050377/150‑1 (MQ=255)

                                                                                                                                          aCCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACgagag     >  1:2873890/1‑150 (MQ=255)

                                                                                                                                          aCCGCGCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACgagag     <  1:2070630/150‑1 (MQ=255)

                                                                                                                                            cgcgCACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGaa   <  1:266095/150‑1 (MQ=255)

                                                                                                                                             gcgcACCACGTGTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGaaa  <  1:701755/150‑1 (MQ=255)

||

GACTTGCCACCCAGCAATCGCACATTCCCGCCAAGACGGTTGAAGATGCAGTAAAAGAGATGCTGGAGCATATGGCCTCGACTCTTGCGCAGGGCGAGCGTATTGAAATCCGCGGTTTCGGCAGTTTCTCTTTGCACTACCGCGCACCAC‑‑GTACCGGACGTAATCCGAAGACTGGCGATAAAGTAGAACTGGAAGGAAAATACGTTCCTCACTTTAAACCTGGTAAAGAACTGCGCGATCGCGCCAATATTTACGGTTAAGTTTTTTACTCAAACTTGAACGAGAGAAA  >  AP012306/826586‑826874

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 24 ≤` `ATCG/ATCG` `< 30 ≤` `ATCG/ATCG` `< 35 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`