BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	144	A→G	intergenic (–/‑46)	– / → thrL	–/thr operon leader peptide

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	144	0	A	G	100.0%	574.0 / NA	162	intergenic (–/‑46)	–/thrL	–/thr operon leader peptide
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (96/66); total (96/66)

Only 100 of 164 total aligned reads displayed.

GAGTGTCTGATAGCAGCTTCTGAACTGGTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTAC  >  NC_000913/54‑233

  gtgtCTGATAGCAGCTTCTGAACTGGTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCa                                                                                       >  1:5846666/1‑93 (MQ=255)

  gtgtCTGATAGCAGCTTCTGAACTGGTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCa                                                                                       >  1:4985531/1‑93 (MQ=255)

  gtgtCTGATAGCAGCTTCTGAACTGGTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCa                                                                                       >  1:4235497/1‑93 (MQ=255)

   tgtCTGATAGCAGCTTCTGAACTGGTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCAt                                                                                      >  2:5203767/1‑93 (MQ=255)

   tgtCTGATAGCAGCTTCTGAACTGGTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCAt                                                                                      <  1:157616/93‑1 (MQ=255)

   tgtCTGATAGCAGCTTCTGAACTGGTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCAt                                                                                      <  1:3339515/93‑1 (MQ=255)

     tCTGATAGCAGCTTCTGAACTGGTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAg                                                                                    <  1:881901/93‑1 (MQ=255)

        gATAGCAGCTTCTGAACTGGTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAgcgc                                                                                 >  2:5967284/1‑93 (MQ=255)

        gATAGCAGCTTCTGAACTGGTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAgcgc                                                                                 <  1:3418542/93‑1 (MQ=255)

           agcagcTTCTGAACTGGTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGcaca                                                                              <  1:1252640/93‑1 (MQ=255)

           agcagcTTCTGAACTGGTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGcaca                                                                              >  2:5122341/1‑93 (MQ=255)

           agcagcTTCTGAACTGGTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGcaca                                                                              >  1:5915863/1‑93 (MQ=255)

           agcagcTTCTGAACTGGTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGcaca                                                                              <  1:1205829/93‑1 (MQ=255)

            gcagcTTCTGAACTGGTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCacag                                                                             >  2:4468709/1‑93 (MQ=255)

             cagcTTCTGAACTGGTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCacaga                                                                            <  2:6201099/93‑1 (MQ=255)

              agcTTCTGAACTGGTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCacagac                                                                           >  1:6038888/1‑93 (MQ=255)

              agcTTCTGAACTGGTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCacagac                                                                           <  2:1843651/93‑1 (MQ=255)

               gcTTCTGAACTGGTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCacagaca                                                                          <  1:4806017/93‑1 (MQ=255)

                cTTCTGAACTGTTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCacagacag                                                                         <  1:366590/93‑1 (MQ=255)

                cTTCTGAACTGGTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCacagacag                                                                         <  1:4456743/93‑1 (MQ=255)

                 ttCTGAACTGGTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCAcagacaga                                                                        <  1:2851576/93‑1 (MQ=255)

                 ttCTGAACTGGTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCAcagacaga                                                                        >  1:5798361/1‑93 (MQ=255)

                 ttCTGAACTGGTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCAcagacaga                                                                        <  2:329499/93‑1 (MQ=255)

                  tCTGAACTGGTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGAt                                                                       >  2:2581926/1‑93 (MQ=255)

                  tCTAAACTGGTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGAt                                                                       <  2:6196281/93‑1 (MQ=255)

                   cTGAACTGGTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATa                                                                      >  1:4755492/1‑93 (MQ=255)

                     gAACTGGTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATaaa                                                                    >  2:707473/1‑93 (MQ=255)

                     gAACTGGTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATaaa                                                                    >  1:5230972/1‑93 (MQ=255)

                       aCTGGTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATaaaaa                                                                  >  2:275472/1‑93 (MQ=255)

                       aCTGGTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATaaaaa                                                                  >  2:430346/1‑93 (MQ=255)

                       aCTGGTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATaaaaa                                                                  <  1:195818/93‑1 (MQ=255)

                         tGGTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAAtt                                                                >  2:4744351/1‑93 (MQ=255)

                          ggTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTa                                                               >  1:5968539/1‑93 (MQ=255)

                           gTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGTGCACAGACAGATAAAAATTAc                                                              >  2:2334180/1‑93 (MQ=255)

                           gTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTAc                                                              >  2:1316825/1‑93 (MQ=255)

                           gTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTAc                                                              <  1:5862504/93‑1 (MQ=255)

                             tACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACag                                                            <  1:4186419/93‑1 (MQ=255)

                             tACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACag                                                            <  2:3261250/93‑1 (MQ=255)

                              aCCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACaga                                                           <  2:5702378/93‑1 (MQ=255)

                              aCCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACaga                                                           >  1:938480/1‑93 (MQ=255)

                               ccTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACagag                                                          >  2:4901089/1‑93 (MQ=255)

                                cTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGt                                                         >  1:4930803/1‑93 (MQ=255)

                                 tGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTa                                                        >  1:3752874/1‑93 (MQ=255)

                                   ccGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACTGACAGATAAAAATTACAGAGTaca                                                      <  1:6226590/93‑1 (MQ=255)

                                   ccGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTaca                                                      <  1:5452782/93‑1 (MQ=255)

                                    cGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTacac                                                     <  1:620530/93‑1 (MQ=255)

                                     gTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTAcaca                                                    >  2:903555/1‑93 (MQ=255)

                                      tGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACacaa                                                   >  1:3314969/1‑93 (MQ=255)

                                       gAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACacaac                                                  >  2:1427563/1‑93 (MQ=255)

                                       gAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACacaac                                                  >  1:4022913/1‑93 (MQ=255)

                                       gAGTAAATGAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACacaac                                                  >  1:6318825/1‑93 (MQ=255)

                                        aGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACacaaca                                                 <  2:3248223/93‑1 (MQ=255)

                                         gTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACAt                                                >  2:4395687/1‑93 (MQ=255)

                                          tAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATc                                               >  2:2287176/1‑93 (MQ=255)

                                            aaTTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCa                                             >  1:2893201/1‑93 (MQ=255)

                                             aTTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCAt                                            >  1:790755/1‑93 (MQ=255)

                                             aTTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCAt                                            >  1:1949310/1‑93 (MQ=255)

                                              ttAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCAta                                           >  2:5719682/1‑92 (MQ=255)

                                              ttAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATg                                           <  1:1370977/93‑1 (MQ=255)

                                               tAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGa                                          <  2:3231746/93‑1 (MQ=255)

                                                 aaaTTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGaaa                                        >  1:4972374/1‑93 (MQ=255)

                                                 aaaTTTTATTGACTGAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGaaa                                        <  2:2881221/93‑1 (MQ=255)

                                                  aaTTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCCTGAAAc                                       >  1:2608431/1‑93 (MQ=255)

                                                  aaTTTTATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAAc                                       >  2:3074177/1‑93 (MQ=255)

                                                       tATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCAtt                                  <  2:1942136/93‑1 (MQ=255)

                                                       tATTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCAtt                                  <  1:134480/93‑1 (MQ=255)

                                                        aTTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTa                                 >  1:2304429/1‑93 (MQ=255)

                                                        aTTGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTa                                 <  1:3424005/93‑1 (MQ=255)

                                                         ttGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAg                                <  2:3909697/93‑1 (MQ=255)

                                                         ttGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAg                                >  2:3635831/1‑93 (MQ=255)

                                                          tGACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGc                               <  2:5426776/93‑1 (MQ=255)

                                                           gACTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCA‑TAGcac                             >  1:5084450/1‑93 (MQ=255)

                                                            aCTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGcac                             >  1:3883316/1‑93 (MQ=255)

                                                             cTTAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCGTGAAACGCATTAGcacc                            <  2:3816415/93‑1 (MQ=255)

                                                               tAGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGcaccac                          <  2:6338231/93‑1 (MQ=255)

                                                                aGGTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAACATTACAGAGTACACAACATCCATGAAACGCATTAGCaccacc                         <  2:3432341/93‑1 (MQ=255)

                                                                 ggTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCAccacca                        <  1:5475000/93‑1 (MQ=255)

                                                                  gTCACTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCAt                       <  2:1276431/93‑1 (MQ=255)

                                                                     aCTAAATACTTTAACCAATATGGGCATAGCGCACAAACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCAATAGCACCACCATTac                    >  1:1425148/1‑93 (MQ=255)

                                                                     aCTAAATACTTTAACCAATATGGGCATAGAGCACAGACAGATAAAAATTACAGCGTACACAACATCCATGAAACGCATTAGCACCACCATTac                    >  1:3227518/1‑93 (MQ=255)

                                                                      cTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTacc                   >  2:1168257/1‑93 (MQ=255)

                                                                      cTAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTacc                   >  2:2334219/1‑93 (MQ=255)

                                                                       tAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTacca                  >  2:1137168/1‑93 (MQ=255)

                                                                       tAAATACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTacca                  >  2:4095834/1‑93 (MQ=255)

                                                                        aaaTACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTaccac                 >  1:2619919/1‑93 (MQ=255)

                                                                        aaaTACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTaccac                 >  1:3490375/1‑93 (MQ=255)

                                                                         aaTACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTaccacc                >  2:4109543/1‑93 (MQ=255)

                                                                         aaTACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTaccacc                >  2:5224742/1‑93 (MQ=255)

                                                                         aaTACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTaccacc                >  2:519870/1‑93 (MQ=255)

                                                                          aTACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTaccacca               <  1:4231783/93‑1 (MQ=255)

                                                                           tACTTTAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTaccaccac              >  1:1836208/1‑93 (MQ=255)

                                                                              tttAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACcaccat           >  2:2537181/1‑93 (MQ=255)

                                                                               ttAACCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACcaccatc          >  2:860893/1‑93 (MQ=255)

                                                                                 aaCCAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACcaccatcac        >  1:5230085/1‑93 (MQ=255)

                                                                                    cAATATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACcaccatcaccat     <  1:1559649/93‑1 (MQ=255)

                                                                                     aaTATGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCAtt    >  2:1141755/1‑93 (MQ=255)

                                                                                      atatGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTa   >  2:5682944/1‑93 (MQ=255)

                                                                                      atatGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTa   >  1:3145696/1‑93 (MQ=255)

                                                                                      atatGGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTa   >  2:4873684/1‑93 (MQ=255)

                                                                                       tatgGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTAc  >  2:3944763/1‑93 (MQ=255)

GAGTGTCTGATAGCAGCTTCTGAACTGGTTACCTGCCGTGAGTAAATTAAAATTTTATTGACTTAGGTCACTAAATACTTTAACCAATATAGGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTAC  >  NC_000913/54‑233

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