BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,211,571	A→G	R44R (CGT→CGC)	yjaB ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,211,571	0	A	G	100.0%	701.9 / NA	219	R44R (CGT→CGC)	yjaB	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (118/101); total (118/101)

Only 100 of 250 total aligned reads displayed.

ATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCAACGAGTTC  >  NC_000913/4211481‑4211660

aTCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGgagcg                                                                                         <  2:3260851/93‑1 (MQ=255)

 tCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGgagcga                                                                                        <  2:504450/93‑1 (MQ=255)

 tCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGgagcga                                                                                        >  2:4671790/1‑93 (MQ=255)

  ccATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGgagcgaa                                                                                       >  2:5420915/1‑93 (MQ=255)

   cATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGgagcgaac                                                                                      >  2:2374468/1‑93 (MQ=255)

    atatGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACt                                                                                     >  2:783336/1‑93 (MQ=255)

     tatGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTa                                                                                    <  1:2066092/93‑1 (MQ=255)

     tatGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTa                                                                                    <  1:5539031/93‑1 (MQ=255)

        gctgcCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGt                                                                                 <  2:1607217/93‑1 (MQ=255)

        gctgcCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGt                                                                                 <  2:1357931/93‑1 (MQ=255)

         ctgcCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGtc                                                                                <  2:1802634/93‑1 (MQ=255)

          tgcCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGtcc                                                                               >  1:4569378/1‑93 (MQ=255)

           gcCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGtcct                                                                              <  2:5926686/93‑1 (MQ=255)

            cccACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGtcctc                                                                             <  1:4761758/93‑1 (MQ=255)

             ccACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGtcctcc                                                                            <  1:384360/93‑1 (MQ=255)

              cACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCa                                                                           <  2:2709971/93‑1 (MQ=255)

               aCTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAg                                                                          >  2:5855331/1‑93 (MQ=255)

               aCTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAg                                                                          >  2:5102269/1‑93 (MQ=255)

                cTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGc                                                                         >  1:5283009/1‑93 (MQ=255)

                 ttAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCt                                                                        >  1:2420493/1‑93 (MQ=255)

                 ttAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCt                                                                        <  2:2623615/93‑1 (MQ=255)

                  tAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGCCCCACCACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTc                                                                       >  2:4239811/1‑93 (MQ=255)

                    gCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCgg                                                                     >  2:2873030/1‑93 (MQ=255)

                      aaCATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTc                                                                   <  1:2219050/93‑1 (MQ=255)

                       aCATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTcc                                                                  >  2:2240254/1‑93 (MQ=255)

                         aTAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGa                                                                >  2:5575660/1‑93 (MQ=255)

                         aTAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGa                                                                <  1:1215633/93‑1 (MQ=255)

                          tAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGAt                                                               >  1:4306968/1‑93 (MQ=255)

                           aaaTCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATa                                                              <  2:1974110/93‑1 (MQ=255)

                           aaaTCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATa                                                              <  2:5718094/93‑1 (MQ=255)

                            aaTCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATAc                                                             <  1:5221295/93‑1 (MQ=255)

                            aaTCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATAc                                                             <  2:6169432/93‑1 (MQ=255)

                            aaTCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTACGATAc                                                             >  1:3436403/1‑93 (MQ=255)

                              tCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATActc                                                           >  1:5890862/1‑93 (MQ=255)

                                cAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTctg                                                         >  1:5877071/1‑93 (MQ=255)

                                   ccGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTctgctg                                                      <  1:4278043/93‑1 (MQ=255)

                                   ccGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTctgctg                                                      >  1:4095297/1‑93 (MQ=255)

                                    cGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGa                                                     >  2:6257919/1‑93 (MQ=255)

                                    cGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGa                                                     <  2:5070094/93‑1 (MQ=255)

                                    cGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGa                                                     >  1:5990485/1‑93 (MQ=255)

                                     ggCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGAt                                                    <  2:1517449/93‑1 (MQ=255)

                                     ggCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGAt                                                    <  2:1877537/93‑1 (MQ=255)

                                      gCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATa                                                   <  2:4594742/93‑1 (MQ=255)

                                        tGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGa                                                 >  2:3608056/1‑93 (MQ=255)

                                         gATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGaa                                                >  2:5684062/1‑93 (MQ=255)

                                         gATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGaa                                                >  1:4803334/1‑93 (MQ=255)

                                           tCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGaaaa                                              >  2:2676854/1‑93 (MQ=255)

                                           tCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGaaaa                                              >  2:2500357/1‑93 (MQ=255)

                                            cccGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAAt                                             >  1:2166648/1‑93 (MQ=255)

                                             ccGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATc                                            >  2:5311331/1‑93 (MQ=255)

                                               gCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGt                                          <  2:818415/93‑1 (MQ=255)

                                               gCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGt                                          >  1:4690900/1‑93 (MQ=255)

                                               gCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGt                                          >  2:2109717/1‑93 (MQ=255)

                                                cTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCNAGNTCGGTCCGATACTCTGCTGATAGAAAATCGTg                                         >  2:5172061/1‑93 (MQ=255)

                                                cTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTg                                         >  2:1438650/1‑93 (MQ=255)

                                                cTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTg                                         >  1:3053971/1‑93 (MQ=255)

                                                cTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTg                                         >  1:5803939/1‑93 (MQ=255)

                                                 tCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATAATCTGCTGATAGAAAATCGTgg                                        >  2:4966016/1‑92 (MQ=255)

                                                  cATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAg                                       <  1:2547530/93‑1 (MQ=255)

                                                  cATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAg                                       >  1:3625880/1‑93 (MQ=255)

                                                  cATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAg                                       >  1:4048720/1‑93 (MQ=255)

                                                  cATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAg                                       <  1:3894688/93‑1 (MQ=255)

                                                   aTTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGt                                      >  1:5100347/1‑93 (MQ=255)

                                                    ttAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTg                                     <  1:6081245/93‑1 (MQ=255)

                                                    ttAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTg                                     <  1:5820648/93‑1 (MQ=255)

                                                        ccGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCAt                                 <  1:1104726/93‑1 (MQ=255)

                                                        ccGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCAt                                 <  2:449478/93‑1 (MQ=255)

                                                          gcgACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCg                               >  2:887819/1‑93 (MQ=255)

                                                           cgACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGa                              <  2:4939410/93‑1 (MQ=255)

                                                           cgACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGa                              >  1:237906/1‑93 (MQ=255)

                                                            gACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGAc                             <  1:5515353/93‑1 (MQ=255)

                                                             aCCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACa                            >  2:3838825/1‑93 (MQ=255)

                                                             aCCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACa                            >  1:6225892/1‑93 (MQ=255)

                                                              cccACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAg                           >  2:1035765/1‑93 (MQ=255)

                                                               ccACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGa                          >  1:1564913/1‑93 (MQ=255)

                                                               ccACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGa                          <  1:2879805/93‑1 (MQ=255)

                                                                cacaACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGaa                         >  1:1412691/1‑93 (MQ=255)

                                                                cacaACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGaa                         <  2:4413123/93‑1 (MQ=255)

                                                                cacaACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGaa                         <  2:6271565/93‑1 (MQ=255)

                                                                cacaACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGaa                         >  2:6077064/1‑93 (MQ=255)

                                                                 acaACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAAc                        <  2:1283613/93‑1 (MQ=255)

                                                                    aCGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGac                     >  2:911252/1‑93 (MQ=255)

                                                                    aCGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGac                     >  1:161132/1‑93 (MQ=255)

                                                                    aCGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGAGGCATCGACAGAACGac                     >  2:2526596/1‑93 (MQ=255)

                                                                     cGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGaca                    <  1:1973877/93‑1 (MQ=255)

                                                                        cgcTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCa                 <  2:5513011/93‑1 (MQ=255)

                                                                          cTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCaaa               >  2:2046391/1‑93 (MQ=255)

                                                                            tCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATc             <  1:3151737/93‑1 (MQ=255)

                                                                               ggCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCa          >  1:1000414/1‑93 (MQ=255)

                                                                                gCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCaa         <  1:738934/93‑1 (MQ=255)

                                                                                gCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCaa         <  2:4494018/93‑1 (MQ=255)

                                                                                gCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCaa         <  2:717114/93‑1 (MQ=255)

                                                                                  aggaaggaGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCAACg       >  1:3521976/1‑93 (MQ=255)

                                                                                  aggaaggaGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCAACg       <  2:631498/93‑1 (MQ=255)

                                                                                   ggaaggagCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCAACGa      >  1:6254485/1‑93 (MQ=255)

                                                                                   ggaaggagCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCAACGa      >  1:818627/1‑93 (MQ=255)

                                                                                    gaaggagCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCAACGAg     >  1:3608540/1‑93 (MQ=255)

                                                                                     aaggagCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCAACGAGt    >  2:3728861/1‑93 (MQ=255)

                                                                                       ggagCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCAACGAGTtc  >  1:5065349/1‑93 (MQ=255)

                                                                                       ggagCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCAACGAGTtc  >  2:2364727/1‑93 (MQ=255)

ATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCAACGAGTTC  >  NC_000913/4211481‑4211660

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