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%	667.7 / NA	212	R44R (CGT→CGC)	yjaB	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (106/106); total (106/106)

Only 100 of 253 total aligned reads displayed.

ATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCAACGAGTTC  >  NC_000913/4211481‑4211660

aTCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGgagcg                                                                                         >  2:2144584/1‑93 (MQ=255)

aTCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGgagcg                                                                                         >  2:2834593/1‑93 (MQ=255)

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

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

  ccATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGgagcgaa                                                                                       <  2:693858/93‑1 (MQ=255)

  ccATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGgagcgaa                                                                                       <  2:3693012/93‑1 (MQ=255)

   cATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGgagcgaac                                                                                      <  2:1322600/93‑1 (MQ=255)

      atGCTTCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAg                                                                                   <  2:760181/93‑1 (MQ=255)

       tgctgcCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAgg                                                                                  >  1:4100257/1‑93 (MQ=255)

        gctgcCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGt                                                                                 >  2:2757269/1‑93 (MQ=255)

          tgcCCACTTAGCAACATAAATCCACCCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGtcc                                                                               <  2:4540309/93‑1 (MQ=255)

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

          tgcCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGtcc                                                                               <  2:1251229/93‑1 (MQ=255)

            cccACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGtcctc                                                                             >  2:2700846/1‑93 (MQ=255)

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

              cACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCa                                                                           >  1:4232339/1‑93 (MQ=255)

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

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

              cACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGCCCTCCa                                                                           >  1:480357/1‑93 (MQ=255)

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

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

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

                    gCAACATAAATCAAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCgg                                                                     <  1:797995/93‑1 (MQ=255)

                     cAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCggg                                                                    >  2:3057819/1‑92 (MQ=255)

                     cAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGt                                                                    <  2:4309315/93‑1 (MQ=255)

                       aCATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTcc                                                                  <  1:3328693/93‑1 (MQ=255)

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

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

                        cATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCg                                                                 >  2:3587370/1‑93 (MQ=255)

                        cATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCg                                                                 >  2:1785101/1‑93 (MQ=255)

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

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

                             aTCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATAct                                                            >  1:3586934/1‑93 (MQ=255)

                              tCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATActc                                                           <  1:1973109/93‑1 (MQ=255)

                              tCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATActc                                                           <  2:2577200/93‑1 (MQ=255)

                               ccAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATActct                                                          <  2:3855570/93‑1 (MQ=255)

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

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

                                 aaCCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTctgc                                                        <  2:4426787/93‑1 (MQ=255)

                                   ccGGCTGATCCCGCTCATTAACCGCGACCCACAGCGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTctgctg                                                      <  1:1083213/93‑1 (MQ=255)

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

                                        tGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTTATAGa                                                 >  1:948739/1‑93 (MQ=255)

                                        tGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGa                                                 <  2:3124333/93‑1 (MQ=255)

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

                                        tGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGa                                                 <  1:1199918/93‑1 (MQ=255)

                                          aTCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTGGGTCCGATACTCTGGTGATAGaga                                               >  2:1161582/1‑91 (MQ=255)

                                           tCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGaaaa                                              <  2:4240871/93‑1 (MQ=255)

                                           acccGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGaaaa                                              <  1:1364490/92‑1 (MQ=255)

                                            cccGCTCATTACCCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAAt                                             <  2:1030612/93‑1 (MQ=255)

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

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

                                              cGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCg                                           <  2:3337586/93‑1 (MQ=255)

                                                  cATTACCCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAg                                       <  2:1053916/93‑1 (MQ=255)

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

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

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

                                                       aCCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCa                                  >  1:2957294/1‑93 (MQ=255)

                                                        ccGCGACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCAt                                 >  2:1444316/1‑93 (MQ=255)

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

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

                                                         cgcgACCCACAACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATc                                >  1:4177225/1‑93 (MQ=255)

                                                              cccACAACGGCGCTTCCGGCAGCAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAg                           >  2:3169808/1‑93 (MQ=255)

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

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

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

                                                                 ccaACGGCGCTTCGGGCAGGAAGGAGCGAACTAGGTCCTCGAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGTCATCGACAGAAc                        <  2:1730351/92‑1 (MQ=255)

                                                                 ccaACGGCGCTTCCGACAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAAc                        <  1:1738365/92‑1 (MQ=255)

                                                                 acaACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAAc                        >  1:2138766/1‑93 (MQ=255)

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

                                                                  caACGGCGCTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACg                       >  2:4069597/1‑93 (MQ=255)

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

                                                                       gcgcTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACAcc                  >  2:2820411/1‑93 (MQ=255)

                                                                       gcgcTTCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACAcc                  >  1:2013120/1‑93 (MQ=255)

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

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

                                                                           ttCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGCACGACACCAAAt              >  1:4133699/1‑93 (MQ=255)

                                                                           ttCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAAt              >  2:420111/1‑93 (MQ=255)

                                                                           ttCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAAt              >  1:2385043/1‑93 (MQ=255)

                                                                           ttCCGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAAt              <  1:866796/93‑1 (MQ=255)

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

                                                                              cGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGc           <  1:1865787/93‑1 (MQ=255)

                                                                              cGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGc           <  1:614556/93‑1 (MQ=255)

                                                                              cGGCAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGc           >  1:4582043/1‑93 (MQ=255)

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

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

                                                                                 cAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCAAc        <  1:2790881/93‑1 (MQ=255)

                                                                                 cAGGAAGGAGCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCAAc        <  2:2426818/93‑1 (MQ=255)

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

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

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

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

                                                                                    gaaggagCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCAACGAg     <  1:3885817/93‑1 (MQ=255)

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

                                                                                    gaaggagCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCAACGAg     <  1:3187684/93‑1 (MQ=255)

                                                                                     aaggagCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCAACGAGt    <  2:2101631/93‑1 (MQ=255)

                                                                                     aaggagCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCAACGAGt    <  2:1877602/93‑1 (MQ=255)

                                                                                     aaggagCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCAACGAGt    <  1:3562418/93‑1 (MQ=255)

                                                                                      aggagCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCAACGAGtt   >  1:4480423/1‑93 (MQ=255)

                                                                                      aggagCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCAACGAGtt   >  2:1102671/1‑93 (MQ=255)

                                                                                       ggagCGAACTAGGTCCTCCAGCTCGGTCCGATACTCTGCTGATAGAAAATCGTGAGTGGCATCGACAGAACGACACCAAATCGCAACGAGTtc  <  1:3960628/93‑1 (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: `‑`