BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	720,640	C→A	15.3%	H20N (CAC→AAC)	ybfK →	uncharacterized protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	720,640	0	C	A	15.3%	375.2 / 46.4	118	H20N (CAC→AAC)	ybfK	uncharacterized protein
Reads supporting (aligned to +/- strand): ref base C (47/53); new base A (5/13); total (52/66)
Fisher's exact test for biased strand distribution p-value = 1.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01

Only 100 of 120 total aligned reads displayed.

TCCTCAATTTCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTAGACTAACGGGCATTGAAGAAGTGAAAGTCGAAGCTATTGCGATGAGCAAACTTCATAATATGA  >  NC_000913/720492‑720778

tCCTCAATTTCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGc                                                                                                                                            >  3:488627/1‑149 (MQ=255)

tCCTCAATTTCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGc                                                                                                                                            >  1:396282/1‑149 (MQ=255)

   tCAATTTCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACt                                                                                                                                         >  1:78210/1‑149 (MQ=255)

   tCAATTTCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACt                                                                                                                                         <  1:559385/149‑1 (MQ=255)

   tCAATTTCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACt                                                                                                                                         >  2:559385/1‑149 (MQ=255)

   tCAATTTCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACt                                                                                                                                         >  4:571040/1‑149 (MQ=255)

   tCAATTTCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACt                                                                                                                                         >  4:883789/1‑149 (MQ=255)

     aaTTTCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGc                                                                                                                                       >  4:860216/1‑149 (MQ=255)

       tttCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCgg                                                                                                                                     >  1:177944/1‑149 (MQ=255)

         tctcGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGaa                                                                                                                                   <  1:316508/149‑1 (MQ=255)

         tctcGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGaa                                                                                                                                   <  2:177944/149‑1 (MQ=255)

         tctcGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGaa                                                                                                                                   <  4:488627/149‑1 (MQ=255)

          ctcGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGAACTGCGGAAt                                                                                                                                  >  2:308381/1‑149 (MQ=255)

           tcGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATa                                                                                                                                 >  3:391258/1‑149 (MQ=255)

           tcGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATa                                                                                                                                 >  4:555672/1‑149 (MQ=255)

             gAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGAACTGCGGAATACa                                                                                                                               >  1:206288/1‑149 (MQ=255)

              aaTGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAg                                                                                                                              <  1:480165/149‑1 (MQ=255)

              aaTGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAg                                                                                                                              <  1:82822/149‑1 (MQ=255)

                  aggaaTCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCat                                                                                                                          >  1:512481/4‑149 (MQ=255)

                  agaAATCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCat                                                                                                                          >  2:493606/1‑149 (MQ=255)

                   gaAATCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCata                                                                                                                         >  1:456653/1‑149 (MQ=255)

                    aaaTCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCatat                                                                                                                        >  4:826036/1‑149 (MQ=255)

                       tcgatcgTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGcc                                                                                                                     >  1:19837/1‑149 (MQ=255)

                        cgatcgTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCt                                                                                                                    <  2:433501/149‑1 (MQ=255)

                         gatcgTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTg                                                                                                                   >  2:869054/1‑149 (MQ=255)

                         gatcgTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACGGCGGAATACAGGCATATGCCTg                                                                                                                   <  4:692848/149‑1 (MQ=255)

                          atcgTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGt                                                                                                                  <  1:503584/149‑1 (MQ=255)

                             gTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTatt                                                                                                               <  4:770009/149‑3 (MQ=255)

                              tGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGAACTGCGGAATACAGGCATATGCCTGTAATa                                                                                                              <  1:308381/149‑1 (MQ=255)

                                ccGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACa                                                                                                            <  2:19837/149‑1 (MQ=255)

                                ccGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGAACTGCGGAATACAGGCATATGCCTGTAATACa                                                                                                            <  1:763995/149‑1 (MQ=255)

                                ccGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGAACTGCGGAATACAGGCATATGCCTGTAATACa                                                                                                            >  4:166407/1‑149 (MQ=255)

                                  ggAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACaga                                                                                                          <  4:595810/149‑1 (MQ=255)

                                        aaGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGCATGCAGAAACGCTTACCAGTCAGGTGGCTGACAATACAACGGTTAACAGAGAATGAGCATTGCGGAATACAGGCATATGCCCGTAATACAGAGAGTCa                                                                                                    <  3:575231/149‑1 (MQ=255)

                                        aaGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCa                                                                                                    <  2:512481/149‑1 (MQ=255)

                                         aGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAg                                                                                                   >  3:459170/1‑149 (MQ=255)

                                         aGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAg                                                                                                   <  2:396282/149‑1 (MQ=255)

                                          gCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAgg                                                                                                  >  1:400416/1‑149 (MQ=255)

                                            aaaGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGaa                                                                                                <  2:585351/149‑1 (MQ=255)

                                               gTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGAc                                                                                             >  1:680186/1‑149 (MQ=255)

                                               gTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGAc                                                                                             <  3:360921/149‑1 (MQ=255)

                                               gTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGAACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGAc                                                                                             <  4:591001/149‑1 (MQ=255)

                                                        cTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCt                                                                                    <  3:295309/149‑1 (MQ=255)

                                                           ataaataaATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGaa                                                                                 <  3:273101/149‑1 (MQ=255)

                                                           ataaataaATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGAACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGaa                                                                                 <  2:889275/149‑1 (MQ=255)

                                                           ataaataaATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGAACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGaa                                                                                 <  3:28413/149‑1 (MQ=255)

                                                             gaataaatCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTGAGGTGGGTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGt                                                                               >  1:814398/2‑149 (MQ=255)

                                                             aaataaatCAGGTAGATTGCTTTTTTCAGAATACAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGt                                                                               <  3:723267/149‑1 (MQ=255)

                                                                taaatCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGAACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGcc                                                                            <  2:813232/149‑1 (MQ=255)

                                                                   atCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCtt                                                                         >  3:675813/1‑149 (MQ=255)

                                                                     cAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAg                                                                       <  1:259485/149‑1 (MQ=255)

                                                                        gTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCg                                                                    >  2:571290/1‑149 (MQ=255)

                                                                           gATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGAACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTa                                                                 <  2:206288/149‑1 (MQ=255)

                                                                           gATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGAACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTa                                                                 <  1:696074/149‑1 (MQ=255)

                                                                              tGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTAGAc                                                              <  3:555672/149‑1 (MQ=255)

                                                                                 ttttttCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAg                                                                                                                              <  3:669668/82‑1 (MQ=255)

                                                                                 ttttttCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAg                                                                                                                              >  4:669668/1‑82 (MQ=255)

                                                                                    tttCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTAGCCTAACgg                                                        <  2:145977/149‑1 (MQ=255)

                                                                                    tttCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTAGACTAACgg                                                        >  3:491442/1‑149 (MQ=255)

                                                                                    tttCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGAACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTAGACTAACgg                                                        <  1:165629/149‑1 (MQ=255)

                                                                                    tttCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGAACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTAGACTAACgg                                                        <  3:485037/149‑1 (MQ=255)

                                                                                      tCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTAGACTAACGGGc                                                      <  3:826036/149‑1 (MQ=255)

                                                                                      tCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTAGACTAACGGGc                                                      <  1:59871/149‑1 (MQ=255)

                                                                                        aGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGAACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGCGCCCTTAGTCGGTAGACTAACGGGCAt                                                    <  4:420545/149‑1 (MQ=255)

                                                                                         gAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGaa                                                                                                                                   >  4:494852/1‑69 (MQ=255)

                                                                                         gAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGaa                                                                                                                                   <  3:494852/69‑1 (MQ=255)

                                                                                         gAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTAGACTAACGGGCAtt                                                   >  2:769364/1‑149 (MQ=255)

                                                                                            tGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTAGACTAACGGGCATTgaa                                                <  3:860216/149‑1 (MQ=255)

                                                                                            tGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTAGACTAACGGGCATTgaa                                                >  4:200056/1‑149 (MQ=255)

                                                                                               aGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTAGACTAACGGGCATTgaagaa                                             <  4:542616/149‑1 (MQ=255)

                                                                                               aGAAACCCTTACCAGTCAGGTGGCTAATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTAGACTAACGGGCATTgaagaa                                             >  2:817746/1‑149 (MQ=255)

                                                                                                gAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTAGACTAACGGGCATTGaagaag                                            <  3:742293/149‑1 (MQ=255)

                                                                                                   aCCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTAGACTAACGGGCATCGAAGAAGTGa                                         >  3:675061/1‑149 (MQ=255)

                                                                                                     ccTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTAGACTAACGGGCATTGAAGAAGTg                                          >  3:236501/1‑146 (MQ=255)

                                                                                                     ccTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTAGACTAACGGGCATTGAAGAAGTg                                          <  4:236501/146‑1 (MQ=255)

                                                                                                           cAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTAGACTAACGGGCATTgaagaa                                             <  1:841363/137‑1 (MQ=255)

                                                                                                           cAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTAGACTAACGGGCATTGAAGAAGTGAAAGTCGaa                                 <  4:155484/149‑1 (MQ=255)

                                                                                                             gTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTAGACTAACGGGCATTGAAGAAGTGAAAGTCGAAGc                               >  1:869961/1‑149 (MQ=255)

                                                                                                               cAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACATGTACGCTGAAGTGCCCTTAGTCGGTAGACTAACGGGCATTGAAGAAGTGAAAGTCGAAGCTa                             <  1:450555/149‑1 (MQ=255)

                                                                                                                       tGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTAGACTAACTGGCATTGAAGAAGTGAAAGTCGAAGCTATTGCGATg                     >  4:625007/1‑149 (MQ=255)

                                                                                                                        gATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTAGACTAACGGGCATTGAAGAAGTGAAAGTCGAAGCTATTGCGATGa                    <  2:337734/149‑1 (MQ=255)

                                                                                                                        gATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTAGACTAACGGGCATTGAAGAAGTGAAAGTCGAAGCTATTGCGATGa                    <  1:742200/149‑1 (MQ=255)

                                                                                                                        gATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTAGACTAACGGGCATTGAAGAAGTGAAAGTCGAAGCTATTGCGATGa                    <  2:81725/149‑1 (MQ=255)

                                                                                                                           aataCAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTAGACTAACGGGCATTGAAGAAGTGAAAGTCGAAGCTATTGCGATGAGCa                 <  1:571290/149‑1 (MQ=255)

                                                                                                                              aCAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTAGACTAACGGGCAt                                                    >  1:292127/1‑111 (MQ=255)

                                                                                                                              aCAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTAGACTAACGGGCAt                                                    <  2:292127/111‑1 (MQ=255)

                                                                                                                               cAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTAGACTAACGGGCATTGAAGAAGTGAAAGTCGAAGCTATTGCGATGAGCAAACt             >  1:12105/1‑149 (MQ=255)

                                                                                                                               cAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTAGACTAACGGGCATTGAAGAAGTGAAAGTCGAAGCTATTGCGATGAGCAAACt             >  3:725645/1‑149 (MQ=255)

                                                                                                                                aaGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTAGAc                                                              <  3:475189/99‑1 (MQ=255)

                                                                                                                                aaGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTAGAc                                                              >  4:475189/1‑99 (MQ=255)

                                                                                                                                aaGGATTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTAGACTAACGGGCATTGAATAAGTGAAAGTCGAAGCTATTGCGATGAGCAAAttt            <  2:78210/149‑1 (MQ=255)

                                                                                                                                 aGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCGGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTAGACTAACGGGCATTGAAGAAGTGAAAGTCGAAGCTATTGCGATGAGCAAACTTc           >  4:398460/1‑149 (MQ=255)

                                                                                                                                  gggTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTAGACTAACGGGCATTGAAGAAGTGAAAGTCGAAGCTATTGCGATGAGCAAACTTCa          <  3:181142/149‑1 (MQ=255)

                                                                                                                                   ggTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTAGACTAACGGGCATTGAAGAAGTGAAAGTCGAAGCTATTGCGATGAGCAAACTTCat         >  4:842859/1‑149 (MQ=255)

                                                                                                                                        aCAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTAGACTAACGGGCATTGAAGAAGTGAAAGTCGAAGCTATTGCGATGAGCAAACTTCATAatat    >  4:31552/1‑149 (MQ=255)

                                                                                                                                         cAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTAGACTAACGGGCATTGAAGAAGTGAAAGTCGAAGCTATTGCGATGAGCAAACTTCATAATATg   <  1:493606/149‑1 (MQ=255)

                                                                                                                                          agagAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTAGACTAACGGGCATTGAAGAAGTGAAAGTCGAAGCTATTGCGATGAGCAAACTTCATAATATGa  >  1:325932/1‑149 (MQ=255)

                                                                                                                                          agagAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTAGACTAACGGGCATTGAAGAAGTGAAAGTCGAAGCTATTGCGATGAGCAAACTTCATAATATGa  <  2:12105/149‑1 (MQ=255)

                                                                                                                                          agagAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTAGACTAACGGGCATTGAAGAAGTGAAAGTCGAAGCTATTGCGATGAGCAAACTTCATAATATGa  <  4:702482/149‑1 (MQ=255)

                                                                                                                                          agagAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTAGACTAACGGGCATTGAAGAAGTGAAAGTCGAAGCTATTGCGATGAGCAAACTTCATAATATGa  >  4:873371/1‑149 (MQ=255)

TCCTCAATTTCTCGAATGAGAAATCGATCGTGCCGGAAAAAAGCAAAGTTACCCCGCTGATAAATAAATCAGGTAGATTGCTTTTTTCAGAATGCAGAAACCCTTACCAGTCAGGTGGCTGATAATACAAGGGTTAACAGAGAATGAGCACTGCGGAATACAGGCATATGCCTGTAATACAGAGAGTCAGGAAGACCTGTACGCTGAAGTGCCCTTAGTCGGTAGACTAACGGGCATTGAAGAAGTGAAAGTCGAAGCTATTGCGATGAGCAAACTTCATAATATGA  >  NC_000913/720492‑720778

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 27 ≤` `ATCG/ATCG` `< 34 ≤` `ATCG/ATCG` `< 39 ≤` `ATCG/ATCG` `< 40 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`