BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	570,683	T→C	intergenic (+255/‑210)	ybcK → / → ybcL	DLP12 prophage; putative recombinase/DLP12 prophage; periplasmic protein YbcL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	570,683	0	T	C	100.0%	1370.0 / NA	324	intergenic (+255/‑210)	ybcK/ybcL	DLP12 prophage; putative recombinase/DLP12 prophage; periplasmic protein YbcL
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (165/159); total (165/159)

Only 100 of 324 total aligned reads displayed.

CATTTGTTAGGTGAGAACTCTTGGTCGCAGATTCAAATACTGAAAATACGTGACAAATTATTATGAGCAAAATGGTGTATGTCACGTATTTTGAATGGTAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAA  >  NC_000913/570585‑570780

 aTTTGTTAGGTGAGAACTCTTGGTCGCAGATTCAAATACTGAAAATACGTGACAAATTATTATGAGCAAAATGGTGTATGTCACGTATTTTGAATGGCAgg                                                                                                >  1:941026/1‑101 (MQ=255)

   ttGTTAGGTGAGAACTCTTGGTCGCAGATTCAAATACTGAAAATACGTGACAAATTATTATGAGCAAAATGGTGTATGTCACGTATTTTGAATGGCAGGtt                                                                                              <  2:940604/101‑1 (MQ=255)

   ttGTTAGGTGAGAACTCTTGGTCGCAGATTCAAATACTGAAAATACGTGACAAATTATTATGAGCAAAATGGTGTATGTCACGTATTTTGAATGGCAGGtt                                                                                              >  2:941028/1‑101 (MQ=255)

     gTTAGGTGAGAACTCTTGGTCGCAGATTCAAATACTGAAAATACGTGACAAATTATTATGAGCAAAATGGTGTATGTCACGTATTTTGAATGGCAGGTTaa                                                                                            >  2:941032/1‑101 (MQ=255)

        aGGTGAGAACTCTTGGTCGCAGATTCAAATACTGAAAATACGTGACAAATTATTATGAGCAAAATGGTGTATGTCACGTATTTTGAATGGCAGGTTaaaaa                                                                                         >  2:941036/1‑101 (MQ=255)

        aGGTGAGAACTCTTGGTCGCAGATTCAAATACTGAAAATACGTGACAAATTATTATGAGCAAAATGGTGTATGTCACGTATTTTGAATGGCAGGTTaaaaa                                                                                         >  1:941035/1‑101 (MQ=255)

         ggTGAGAACTCTTGGTCGCAGATTCAAATACTGAAAATACGTGACAAATTATTATGAGCAAAATGGTGTATGTCACGTATTTTGAATGGCAGGTTaaaaaa                                                                                        >  2:941041/1‑101 (MQ=255)

          gTGAGAACTCTTGGTCGCAGATTCAAATACTGAAAATACGTGACAAATTATTATGAGCAAAATGGTGTATGTCACGTATTTTGAATGGCAGGTTAAAAAAt                                                                                       <  1:940346/101‑1 (MQ=255)

            gagaACTCTTGGTCGCAGATTCAAATACTGAAAATACGTGACAAATTATTATGAGCAAAATGGTGTATGTCACGTATTTTGAATGGCAGGTTAAAAAATaa                                                                                     >  2:941044/1‑101 (MQ=255)

             agaACTCTTGGTCGCAGATTCAAATACTGAAAATACGTGACAAATTATTATGAGCAAAATGGTGTATGTCACGTATTTTGAATGGCAGGTTAAAAAATAac                                                                                    >  1:941045/1‑101 (MQ=255)

             agaACTCTTGGTCGCAGATTCAAATACTGAAAATACGTGACAAATTATTATGAGCAAAATGGTGTATGTCACGTATTTTGAATGGCAGGTTAAAAAATAac                                                                                    <  2:940405/101‑1 (MQ=255)

               aaCTCTTGGTCGCAGATTCAAATACTGAAAATACGTGACAAATTATTATGAGCAAAATGGTGTATGTCACGTATTTTGAATGGCAGGTTAAAAAATAacac                                                                                  >  1:941049/1‑101 (MQ=255)

                aCTCTTGGTCGCAGATTCAAATACTGAAAATACGTGACAAATTATTATGAGCAAAATGGTGTATGTCACGTATTTTGAATGGCAGGTTAAAAAATAACAcc                                                                                 <  1:940563/101‑1 (MQ=255)

                aCTCTTGGTCGCAGATTCAAATACTGAAAATACGTGACAAATTATTATGAGCAAAATGGTGTATGTCACGTATTTTGAATGGCAGGTTAAAAAATAACAcc                                                                                 <  2:940680/101‑1 (MQ=255)

                aCTCTTGGTCGCAGATTCAAATACTGAAAATACGTGACAAATTATTATGAGCAAAATGGTGTATGTCACGTATTTTGAATGGCAGGTTAAAAAATAACAcc                                                                                 <  1:940477/101‑1 (MQ=255)

                 ctctTGGTCGCAGATTCAAATACTGAAAATACGTGACAAATTATTATGAGCAAAATGGTGTATGTCACGTATTTTGAATGGCAGGTTAAAAAATAACACCg                                                                                <  2:940801/101‑1 (MQ=255)

                  tctTGGTCGCAGATTCAAATACTGAAAATACGTGACAAATTATTATGAGCAAAATGGTGTATGTCACGTATTTTGAATGGCAGGTTAAAAAATAACACCGa                                                                               >  2:941055/1‑101 (MQ=255)

                  tctTGGTCGCAGATTCAAATACTGAAAATACGTGACAAATTATTATGAGCAAAATGGTGTATGTCACGTATTTTGAATGGCAGGTTAAAAAATAACACCGa                                                                               <  1:940704/101‑1 (MQ=255)

                   ctTGGTCGCAGATTCAAATACTGAAAATACGTGACAAATTATTATGAGCAAAATGGTGTATGTCACGTATTTTGAATGGCAGGTTAAAAAATAACACCGAc                                                                              >  1:941057/1‑101 (MQ=255)

                    ttGGTCGCAGATTCAAATACTGAAAATACGTGACAAATTATTATGAGCAAAATGGTGTATGTCACGTATTTTGAATGGCAGGTTAAAAAATAACACCGACt                                                                             <  1:940726/101‑1 (MQ=255)

                    ttGGTCGCAGATTCAAATACTGAAAATACGTGACAAATTATTATGAGCAAAATGCTGTATGTCACGTATTTTGAATGGCAGGTTAAAAAATAACACCGACt                                                                             >  2:941059/1‑101 (MQ=255)

                      ggTCGCAGATTCAAATACTGAAAATACGTGACAAATTATTATGAGCAAAATGGTGTATGTCACGTATTTTGAATGGCAGGTTAAAAAATAACCCCGACttt                                                                           <  2:940490/101‑1 (MQ=255)

                      ggTCGCAGATTCAAATACTGAAAATACGTGACAAATTATTATGAGCAAAATGGTGTATGTCACGTATTTTGAATGGCAGGTTAAAAAATAACACCGACttt                                                                           <  2:940569/101‑1 (MQ=255)

                      ggTCGCAGATTCAAATACTGAAAATACGTGACAAATTATTATGAGCAAAATGGTGTATGTCACGTATTTTGAATGGCAGGTTAAAAAATAACACCGACttt                                                                           >  2:941063/1‑101 (MQ=255)

                       gTCGCAGATTCAAATACTGAAAATACGTGACAAATTATTATGAGCAAAATGGTGTATGTCACGTATTTTGAATGGCAGGTTAAAAAATAACACCGACTTTc                                                                          <  1:940396/101‑1 (MQ=255)

                         cGCAGATTCAAATACTGAAAATACGTGACAAATTATTATGAGCAAAATGGTGTATGTCACGTATTTTGAATGGCAGGTTAAAAAATAACACCGACTTTCGt                                                                        >  1:941067/1‑101 (MQ=255)

                          gCAGATTCAAATACTGAAAATACGTGACAAATTATTATGAGCAAAATGGTGTATGTCACGTATTTTGAATGGCAGGTTAAAAAATAACACCGACTTTCGTa                                                                       >  2:941068/1‑101 (MQ=255)

                          gCAGATTCAAATACTGAAAATACGTGACAAATTATTATGAGCAAAATGGTGTATGTCACGTATTTTGAATGGCAGGTTAAAAAATAACACCGACTTTCGTa                                                                       <  1:940634/101‑1 (MQ=255)

                          gCAGATTCAAATACTGAAAATACGTGACAAATTATTATGAGCAAAATGGTGTATGTCACGTATTTTGAATGGCAGGTTAAAAAATAACACCGACTTTCGTa                                                                       <  2:940754/101‑1 (MQ=255)

                          gCAGATTCAAATACTGAAAATACGTGACAAATTATTATGAGCAAAATGGTGTATGTCACGTATTTTGAATGGCAGGTTAAAAAATAACACCGACTTTCGTa                                                                       <  1:940573/101‑1 (MQ=255)

                           cAGATTCAAATACTGAAAATACGTGACAAATTATTATGAGCAAAATGGTGTATGTCACGTATTTTGAATGGCAGGTTAAAAAATAACACCGACTTTCGTAg                                                                      <  1:940305/101‑1 (MQ=255)

                           cAGATTCAAATACTGAAAATACGTGACAAATTATTATGAGCAAAATGGTGTATGTCACGTATTTTGAATGGCAGGTTAAAAAATAACACCGACTTTCGTAg                                                                      <  1:940522/101‑1 (MQ=255)

                            aGATTCAAATACTGAAAATACGTGACAAATTATTATGAGCAAAATGGTGTATGTCACGTATTTTGAATGGCAGGTTAAAAAATAACACCGACTTTCGTAgg                                                                     <  1:940577/101‑1 (MQ=255)

                            aGATTCAAATACTGAAAATACGTGACAAATTATTATGAGCAAAATGGTGTATGTCACGTATTTTGAATGGCAGGTTAAAAAATAACACCGACTTTCGTAgg                                                                     >  1:941073/1‑101 (MQ=255)

                            aGATTCAAATACTGAAAATACGTGACAAATTATTATGAGCAAAATGGTGTATGTCACGTATTTTGAATGGCAGGTTAAAAAATAACACCGACTTTCGTAgg                                                                     >  2:941075/1‑101 (MQ=255)

                                tCAAATACTGAAAATACGTGACAAATTATTATGAGCAAAATGGTGTATGTCACGTATTTTGAATGGCAGGTTAAAAAATAACACCGACTTTCGTAGGTGtt                                                                 >  1:941079/1‑101 (MQ=255)

                                 cAAATACTGAAAATACGTGACAAATTATTATGAGCAAAATGGTGTATGTCACGTATTTTGAATGGCAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTa                                                                <  2:940540/101‑1 (MQ=255)

                                  aaaTACTGAAAATACGTGACAAATTATTATGAGCAAAATGGTGTATGTCACGTATTTTGAATGGCAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTAc                                                               <  1:940605/101‑1 (MQ=255)

                                   aaTACTGAAAATACGTGACAAATTATTATGAGCAAAATGGTGTATGTCACGTATTTTGAATGGCAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACt                                                              >  1:941083/1‑101 (MQ=255)

                                     tACTGAAAATACGTGACAAATTATTATGAGCAAAATGGTGTATGTCACGTATTTTGAATGGCAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACtaa                                                            >  1:941088/1‑101 (MQ=255)

                                     tACTGAAAATACGTGACAAATTATTATGAGCAAAATGGTGTATGTCACGTATTTTGAATGGCAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACtaa                                                            >  2:941087/1‑101 (MQ=255)

                                        tGAAAATACGTGACAAATTATTATGAGCAAAATGGTGTATGTCACGTATTTTGAATGGCAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACtaataa                                                         <  2:940708/101‑1 (MQ=255)

                                        tGAAAATACGTGACAAATTATTATGAGCAAAATGGTGTATGTCACGTATTTTGAATGGCAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACtaataa                                                         <  1:940046/101‑1 (MQ=255)

                                         gAAAATACGTGACAAATTATTATGAGCAAAATGGTGTATGTCACGTATTTTGAATGGCAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACtaataat                                                        <  1:940642/101‑1 (MQ=255)

                                         gAAAATACGTGACAAATTATTATGAGCAAAATGGTGTATGTCACGTATTTTGAATGGCAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACtaataat                                                        <  2:940649/101‑1 (MQ=255)

                                          aaaaTACGTGACAAATTATTATGAGCAAAATGGTGTATGTCACGTATTTTGAATGGCAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACtaataata                                                       <  1:940829/101‑1 (MQ=255)

                                            aaTACGTTACAAATTATTATGAGCAAAATGGTGTATGTCACGTATTTTGAATGGCAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATaaa                                                     >  1:941092/1‑101 (MQ=255)

                                             aTACGTGACAAATTATTATGAGCAAAATGGTGTATGTCACGTATTTTGAATGGCAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAg                                                    <  1:940782/101‑1 (MQ=255)

                                              tACGTGACAAATTATTATGAGCAAAATGGTGTATGTCACGTATTTTGAATGGCAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGc                                                   >  2:941096/1‑101 (MQ=255)

                                               aCGTGACAAATTATTATGAGCAAAATGGTGTATGTCACGTATTTTGAATGGCAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCa                                                  >  1:941100/1‑101 (MQ=255)

                                                 gTGACAAATTATTATGAGCAAAATGGTGTATGTCACGTATTTTGAATGGCAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCaga                                                >  2:941101/1‑101 (MQ=255)

                                                     cAAATTATTATGAGCAAAATGGTGTATGTCACGTATTTTGAATGGCAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGttt                                            <  1:940759/101‑1 (MQ=255)

                                                      aaaTTATTATGAGCAAAATGGTGTATGTCACGTATTTTGAATGGCAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGtttt                                           <  1:940670/101‑1 (MQ=255)

                                                      aaaTTATTATGAGCAAAATGGTGTATGTCACGTATTTTGAATGGCAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGtttt                                           >  2:941110/1‑101 (MQ=255)

                                                      aaaTTATTATGAGCAAAATGGTGTATGTCACGTATTTTGAATGGCAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGtttt                                           >  2:941107/1‑101 (MQ=255)

                                                       aaTTATTATGAGCAAAATGGTGTATGTCACGTATTTTGAATGGCAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGttttt                                          >  1:941112/1‑101 (MQ=255)

                                                         ttattaTGAGCAAAATGGTGTATGTCACGTATTTTGAATGGCAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAg                                        >  2:941115/1‑101 (MQ=255)

                                                             tatGAGCAAAATGGTGTATGTCACGTATTTTGAATGGCAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAg                                    >  2:941119/1‑101 (MQ=255)

                                                              atGAGCAAAATGGTGTATGTCACGTATTTTGAATGGCAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGt                                   >  1:941120/1‑101 (MQ=255)

                                                               tGAGCAAAATGGTGTATGTCACGTATTTTGAATGGCAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTa                                  <  2:940620/101‑1 (MQ=255)

                                                               tGAGCAAAATGGTGTATGTCACGTATTTTGAATGGCAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTa                                  >  1:941123/1‑101 (MQ=255)

                                                                 aGCAAAATGGTGTATGTCACGTATTTTGAATGGCAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATc                                >  1:941128/1‑101 (MQ=255)

                                                                  gCAAAATGGTGTATGTCACGTATTTTGAATGGCAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCa                               <  1:940719/101‑1 (MQ=255)

                                                                   cAAAATGGTGTATGTCACGTATTTTGAATGGCAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCaa                              >  1:941133/1‑101 (MQ=255)

                                                                    aaaaTGGTGTATGTCACGTATTTTGAATGGCAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAAt                             >  2:941135/1‑101 (MQ=255)

                                                                     aaaTGGTGTATGTCACGTATTTTGAATGGCAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAAtg                            >  1:941139/1‑101 (MQ=255)

                                                                     aaaTGGTGTATGTCACGTATTTTGAATGGCAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAAtg                            <  1:940501/101‑1 (MQ=255)

                                                                      aaTGGTGTATGTCACGTATTTTGAATGGCAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAAtgt                           <  1:940746/101‑1 (MQ=255)

                                                                      aaTGGTGTATGTCACGTATTTTGAATGGCAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAAtgt                           >  2:941141/1‑101 (MQ=255)

                                                                       aTGGTGTATGTCACGTATTTTGAATGGCAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAAtgtg                          <  2:940465/101‑1 (MQ=255)

                                                                       aTGGTGTATGTCACGTATTTTGAATGGCAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAAtgtg                          <  1:940933/101‑1 (MQ=255)

                                                                        tGGTGTATGTCACGTATTTTGAATGGCAGGTTAAAAAATAACACCCACttgc                                                                          <  1:940535/52‑3 (MQ=255)

                                                                         ggTGTATGTCACGTATTTTGAATGGCAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCt                        >  1:941145/1‑101 (MQ=255)

                                                                         ggTGTATGTCACGTATTTTGAATGGCAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCt                        <  1:940805/101‑1 (MQ=255)

                                                                         ggTGTATGTCACGTATTTTGAATGGCAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCt                        >  2:941146/1‑101 (MQ=255)

                                                                           tgtATGTCACGTATTTTGAATGGCAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCttt                      <  2:940663/101‑1 (MQ=255)

                                                                           tgtATGTCACGTATTTTGAATGGCAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCttt                      <  1:940770/101‑1 (MQ=255)

                                                                              aTGTCACGTATTTTGAATGGCAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTtgtg                   >  1:941150/1‑101 (MQ=255)

                                                                               tGTCACGTATTTTGAATGGCAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTgtgt                  <  2:940843/101‑1 (MQ=255)

                                                                               tGTCACGTATTTTGAATGGCAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTgtgt                  <  1:940796/101‑1 (MQ=255)

                                                                                gTCACGTATTTTGAATGGCAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGta                 <  1:940441/101‑1 (MQ=255)

                                                                                gTCACGTATTTTGAATGGCAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGta                 <  1:940625/101‑1 (MQ=255)

                                                                                 tCACGTATTTTGAATGGCAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGtat                >  1:941154/1‑101 (MQ=255)

                                                                                  cACGTATTTTGAATGGCAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGtata               >  2:941162/1‑101 (MQ=255)

                                                                                  cACGTATTTTGAATGGCAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGtata               >  2:941157/1‑101 (MQ=255)

                                                                                     gTATTTTGAATGGCAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATtg            >  2:941165/1‑101 (MQ=255)

                                                                                      tATTTTGAATGGCAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATtgt           >  1:941167/1‑101 (MQ=255)

                                                                                      tATTTTGAATGGCAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATtgt           <  1:940661/101‑1 (MQ=255)

                                                                                       aTTTTGAATGGCAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATtgtg          >  1:941171/1‑101 (MQ=255)

                                                                                       aTTTTGAATGGCAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATtgtg          <  1:940684/101‑1 (MQ=255)

                                                                                        ttttGAATGGCAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTgg         >  1:941174/1‑101 (MQ=255)

                                                                                         tttGAATGGCAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGc        >  2:941175/1‑101 (MQ=255)

                                                                                         tttGAATGGCAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGc        <  1:940593/101‑1 (MQ=255)

                                                                                          ttGAATGGCAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCa       >  2:941179/1‑101 (MQ=255)

                                                                                             aaTGGCAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAaat    <  2:940363/101‑1 (MQ=255)

                                                                                              aTGGCAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAaata   <  2:940892/101‑1 (MQ=255)

                                                                                              aTGGCAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAaata   >  1:941185/1‑101 (MQ=255)

                                                                                              aTGGCAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAaata   <  2:940696/101‑1 (MQ=255)

                                                                                              aTGGCAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAaata   <  2:940724/101‑1 (MQ=255)

                                                                                               tGGCAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAaataa  >  1:941188/1‑101 (MQ=255)

CATTTGTTAGGTGAGAACTCTTGGTCGCAGATTCAAATACTGAAAATACGTGACAAATTATTATGAGCAAAATGGTGTATGTCACGTATTTTGAATGGTAGGTTAAAAAATAACACCGACTTTCGTAGGTGTTACTAATAATAAAGCAGAGTTTTTAGATAGTATCAATGTGCTTTGTGTATATTGTGGCAAATAA  >  NC_000913/570585‑570780

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