BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,327,860	T→C	L24P (CTT→CCT)	yhbY →	ribosome assembly factor YhbY

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,327,860	0	T	C	100.0%	1579.7 / NA	375	L24P (CTT→CCT)	yhbY	ribosome assembly factor YhbY
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (184/191); total (184/191)

Only 100 of 375 total aligned reads displayed.

TTCAATCCCAACGTAAGCAAAAAATACGATGAATCTGAGTACTAAACAAAAACAGCACCTGAAAGGTCTGGCACATCCGCTCAAGCCAGTTGTTCTGCTTGGCAGTAATGGTTTGACCGAAGGGGTGCTGGCCGAGATTGAACAAGCGTTAGAGCACCATGAACTCATCAAGGTGAAAATCGCCACCGAAGATCGC  >  NC_000913/3327762‑3327957

  cAATCCCAACGTAAGCAAAAAATACGATGAATCTGAGTACTAAACAAAAACAGCACCTGAAAGGTCTGGCACATCCGCTCAAGCCAGTTGTTCTGCCTGGc                                                                                               <  1:5625090/101‑1 (MQ=255)

  cAATCCCAACGTAAGCAAAAAATACGATGAATCTGAGTACTAAACAAAAACAGCACCTGAAAGGTCTGGCACATCCGCTCAAGCCAGTTGTTCTGCCTGGc                                                                                               <  1:5624941/101‑1 (MQ=255)

   aaTCCCAACGTAAGCAAAAAATACGATGAATCTGAGTACTAAACAAAAACAGCACCTGAAAGGTCTGGCACATCCGCTCAAGCCAGTTGTTCTGCCTGGCa                                                                                              <  2:5625069/101‑1 (MQ=255)

   aaTCCCAACGTAAGCAAAAAATACGATGAATCTGAGTACTAAACAAAAACAGCACCTGAAAGGTCTGGCACATCCGCTCAAGCCAGTTGTTCTGCCTGGCa                                                                                              >  2:5625618/1‑101 (MQ=255)

     tCCCAACGTAAGCAAAAAATACGATGAATCTGAGTACTAAACAAAAACAGCACCTGAAAGGTCTGGCACATCCGCTCAAGCCAGTTGTTCTGCCTGGCAGt                                                                                            <  1:5624997/101‑1 (MQ=255)

      cccAACGTAAGCAAAAAATACGATGAATCTGAGTACTAAACAAAAACAGCACCTGAAAGGTCTGGCACATCCGCTCAAGCCAGTTGTTCTGCCTGGCAGTa                                                                                           <  2:5625080/101‑1 (MQ=255)

      cccAACGTAAGCAAAAAATACGATGAATCTGAGTACTAAACAAAAACAGCACCTGAAAGGTCTGGCACATCCGCTCAAGCCAGTTGTTCTGCCTGGCAGTa                                                                                           >  1:5625622/1‑101 (MQ=255)

        cAACGTAAGCAAAAAATACGATGAATCTGAGTACTAAACAAAAACAGCACCTGAAAGGTCTGGCACATCCGCTCAAGCCAGTTGTTCTGCCTGGCAGTAAt                                                                                         <  1:5625143/101‑1 (MQ=255)

        cAACGTAAGCAAAAAATACGATGAATCTGAGTACTAAACAAAAACAGCACCTGAAAGGTCTGGCACATCCGCTCAAGCCAGTTGTTCTGCCTCGCAGTAAt                                                                                         <  1:5625192/101‑1 (MQ=255)

         aaCGTAAGCAAAAAATACGATGAATCTGAGTACTAAACAAAAACAGCACCTGAAAGGTCTGGCACATCCGCTCAAGCCAGTTGTTCTGCCTGGCAGTAATg                                                                                        >  2:5625626/1‑101 (MQ=255)

          aCGTAAGCAAAAAATACGATGAATCTGAGTACTAAACAAAAACAGCACCTGAAAGGTCTGGCACATCCGCTCAAGCCAGTTGTTCTGCCTGGCAGTAATgg                                                                                       <  1:5625059/101‑1 (MQ=255)

          aCGTAAGCAAAAAATACGATGAATCTGAGTACTAAACAAAAACAGCACCTGAAAGGTCTGGCACATCCGCTCAAGCCAGTTGTTCTGCCTGGCAGTAATgg                                                                                       >  2:5625629/1‑101 (MQ=255)

          aCGTAAGCAAAAAATACGATGAATCTGAGTACTAAACAAAAACAGCACCTGAAAGGTCTGGCACATCCGCTCAAGCCAGTTGTTCTGCCTGGCAGTAATgg                                                                                       <  1:5625274/101‑1 (MQ=255)

            gTAAGCAAAAAATACGATGAATCTGAGTACTAAACAAAAACAGCACCTGAAAGGTCTGGCACATCCGCTCAAGCCAGTTGTTCTGCCTGGCAGTAATGGtt                                                                                     >  1:5625634/1‑101 (MQ=255)

               aGCAAAAAATACGATGAATCTGAGTACTAAACAAAAACAGCACCTGAAAGGTCTGGCACATCCGCTCAAGCCAGTTGTTCTGCCTGGCAGTAATGGTTTGa                                                                                  <  2:5625015/101‑1 (MQ=255)

               aGCAAAAAATACGATGAATCTGAGTACTAAACAAAAACAGCACCTGAAAGGTCTGGCACATCCGCTCAAGCCAGTTGTTCTGCCTGGCAGTAATGGTTTGa                                                                                  >  1:5625638/1‑101 (MQ=255)

                  aaaaaaTACGATGAATCTGAGTACTAAACAAAAACAGCACCTGAAAGGTCTGGCACATCCGCTCAAGCCAGTTGTTCTGCCTGGCAGTAATGGTTTGACCg                                                                               <  1:5624867/101‑1 (MQ=255)

                  aaaaaaTACGATGAATCTGAGTACTAAACAAAAACAGCACCTGAAAGGTCTGGCACATCCGCTCAAGCCAGTTGTTCTGCCTGGCAGTAATGGTTTGACCg                                                                               >  2:5625643/1‑101 (MQ=255)

                  aaaaaaTACGATGAATCTGAGTACTAAACAAAAACAGCACCTGAAAGGTCTGGCACATCCGCTCAAGCCAGTTGTTCTGCCTGGCAGTAATGGTTTGACCg                                                                               <  1:5625020/101‑1 (MQ=255)

                    aaaaTACGATGAATCTGAGTACTAAACAAAAACAGCACCTGAAAGGTCTGGCACATCCGCTCAAGCCAGTTGTTCTGCCTGGCAGTAATGGTTTGACCGaa                                                                             >  2:5625649/1‑101 (MQ=255)

                    aaaaTACGATGAATCTGAGTACTAAACAAAAACAGCACCTGAAAGGTCTGGCACATCCGCTCAAGCCAGTTGTTCTGCCTGGCAGTAATGGTTTGACCGaa                                                                             <  1:5625112/101‑1 (MQ=255)

                    aaaaTACGATGAATCTGAGTACTAAACAAAAACAGCACCTGAAAGGTCTGGCACATCCGCTCAAGCCAGTTGTTCTGCCTGGCAGTAATGGTTTGACCGaa                                                                             >  1:5625647/1‑101 (MQ=255)

                     aaaTACGATGAATCTGAGTACTAAACAAAAACAGCACCTGAAAGGTCTGGCACATCCGCTCAAGCCAGTTGTTCTGCCTGGCAGTAATGGTTTGACCGAAg                                                                            <  2:5625298/101‑1 (MQ=255)

                       aTACGATGAATCTGAGTACTAAACAAAAACAGCACCTGAAAGGTCTGGCACATCCGCTCAAGCCAGTTGTTCTGCCTGGCAGTAATGGTTTGACCg                                                                               >  1:5625655/1‑96 (MQ=255)

                        tACGATGAATCTGAGTACTAAACAAAAACAGCACCTGAAAGGTCTGGCACATCCGCTCAAGCCAGTTGTTCTGCCTGGCAGTAATGGTTTGACCGAAgggg                                                                         >  2:5625656/1‑101 (MQ=255)

                          cGATGAATCTGAGTACTAAACAAAAACAGCACCTGAAAGGTCTGGCACATCCGCTCAAGCCAGTTGTTCTGCCTGGCAGTAATGGTTTGACCGAAgggg                                                                         >  1:5625662/1‑99 (MQ=255)

                          cGATGAATCTGAGTACTAAACAAAAACAGCACCTGAAAGGTCTGGCACATCCGCTCAAGCCAGTTGTTCTGCCTGGCAGTAATGGTTTGACCGAAGGGGTg                                                                       <  2:5625044/101‑1 (MQ=255)

                          cGATGAATCTGAGTACTAAACAAAAACAGCACCTGAAAGGTCTGGCACATCCGCTCAAGCCAGTTGTTCTGCCTGGCAGTAATGGTTTGACCGAAGGGGTg                                                                       <  1:5625400/101‑1 (MQ=255)

                            aTGAATCTGAGTACTAAACAAAAACAGCACCTGAAAGGTCTGGCACATCCGCTCAAGCCAGTTGTTCTGCCTGGCAGTAATGGTTTGACCGAAGGGGt                                                                        <  2:5625178/98‑1 (MQ=255)

                            aTGAATCTGAGTACTAAACAAAAACAGCACCTGAAAGGTCTGGCACATCCGCTCAAGCCAGTTGTTCTGCCTGGCAGTAATGGTTTGACCGAAGGGGTGCt                                                                     <  1:5625162/101‑1 (MQ=255)

                            aTGAATCTGAGTACTAAACAAAAACAGCACCTGAAAGGTCTGGCACATCCGCTCAAGCCAGTTGTTCTGCCTGGCAGTAATGGTTTGACCGAAGGGGTGCt                                                                     >  2:5625666/1‑101 (MQ=255)

                              gAATCTGAGTACTAAACAAAAACAGCACCTGAAAGGTCTGGCACATCCGCTCAAGCCAGTTGTTCTGCCTGGCAGTAATGGTTTGACCGAAGGGGTGCTgg                                                                   >  1:5625669/1‑101 (MQ=255)

                               aaTCTGAGTACTAAACAAAAACAGCACCTGAAAGGTCTGGCACATCCGCTCAAGCCAGTTGTTCTGCCTGGCAGTAATGGTTTGACCGAAGGGGTGCTGGc                                                                  >  1:5625673/1‑101 (MQ=255)

                                aTCTGAGTACTAAACAAAAACAGCACCTGAAAGGTCTGGCACATCCGCTCAAGCCAGTTGTTCTGCCTGGCAGTAATGGTTTGACCGAAGGGGTGCTGGcc                                                                 <  2:5625170/101‑1 (MQ=255)

                                aTCTGAGTACTAAACAAAAACAGCACCTGAAAGGTCTGGCACATCCGCTCAAGCCAGTTGTTCTGCCTGGCAGTAATGGTTTGACCGAAGGGGTGCTGGcc                                                                 >  2:5625677/1‑101 (MQ=255)

                                  cTGAGTACTAAACAAAAACAGCACCTGAAAGGTCTGGCACATCCGCTCAAGCCAGTTGTTCTGCCTGGCAGTAATGGTTTGACCGAAGGGGTGCTGGCCga                                                               >  1:5625681/1‑101 (MQ=255)

                                    gAGTACTAAACAAAAACAGCACCTGAAAGGTCTGGCACATCCGCTCAAGCCAGTTGTTCTGCCTGGCAGTAATGGTTTGACCGAAGGGGTGCTGGCCgaga                                                             >  2:5625685/1‑101 (MQ=255)

                                     aGTACTAAACAAAAACAGCACCTGAAAGGTCTGGCACATCCGCTCAAGCCAGTTGTTCTGCCTGGCAGTAATGGTTTGACCGAAGGGGTGCTGGCCGAGAt                                                            >  1:5625687/1‑101 (MQ=255)

                                       tACTAAACAAAAACAGCACCTGAAAGGTCTGGCACATCCGCTCAAGCCAGTTGTTCTGCCTGGCAGTAATGGTTTGACCGAAGGGGTGCTGGCCGAGATTg                                                          <  2:5625136/101‑1 (MQ=255)

                                        aCTAAACAAAAACAGCACCTGAAAGGTCTGGCACATCCGCTCAAGCCAGTTGTTCTGCCTGGCAGTAATGGTTTGACCGAAGGGGTGCTGGCCGAGATTGa                                                         <  2:5625097/101‑1 (MQ=255)

                                        aCTAAACAAAAACAGCACCTGAAAGGTCTGGCACATCCGCTCAAGCCAGTTGTTCTGCCTGGCAGTAATGGTTTGACCGAAGGGGTGCTGGCCGAGATTGa                                                         <  2:5625309/101‑1 (MQ=255)

                                        aCTAAACAAAAACAGCACCTGAAAGGTCTGGCACATCCGCTCAAGCCAGTTGTTCTGCCTGGCAGTAATGGTTTGACCGAAGGGGTGCTGGCCGAGATTGa                                                         >  2:5625692/1‑101 (MQ=255)

                                          tAAACAAAAACAGCACCTGAAAGGTCTGGCACATCCGCTCAAGCCAGTTGTTCTGCCTGGCAGTAATGGTTTGACCGAAGGGGTGCTGGCCGAGATTGAAc                                                       <  2:5625050/101‑1 (MQ=255)

                                          tAAACAAAAACAGCACCTGAAAGGTCTGGCACATCCGCTCAAGCCAGTTGTTCTGCCTGGCAGTAATGGTTTGACCGAAGGGGTGCTGGCCGAGATTGAAc                                                       >  1:5625695/1‑101 (MQ=255)

                                           aaaCAAAAACAGCACCTGAAAGGTCTGGCACATCCGCTCAAGCCAGTTGTTCTGCCTGGCAGTAATGGTTTGA‑CGAAgggg                                                                         >  1:5625699/1‑81 (MQ=255)

                                            aaCAAAAACAGCACCTGAAAGGTCTGGCACATCCGCTCAAGCCAGTTGTTCTGCCTGGCAGTAATGGTTTGACCGAAGGGGTGCTGGCCGAGATTGAACaa                                                     <  1:5624980/101‑1 (MQ=255)

                                             aCAAAAACAGCACCTGAAAGGTCTGGCACATCCGCTCAAGCCAGTTGTTCTGCCTGGCAGTAATGGTTTGACCGAAGGGGTGCTGGCCGAGATTGAACAAg                                                    >  2:5625704/1‑101 (MQ=255)

                                             aCAAAAACAGCACCTGAAAGGTCTGGCACATCCGCTCAAGCCAGTTGTTCTGCCTGGCAGTAATGGTTTGACCGAAGGGGTGCTGGCCGAGATTGAACAAg                                                    <  1:5625032/101‑1 (MQ=255)

                                              cAAAAACAGCACCTGAAAGGTCTGGCACATCCGCTCAAGCCAGTTGTTCTGCCTGGCAGTAATGGTTTGACCGAAGGGGTGCTGGCCGAGATTGAACAAGc                                                   >  1:5625706/1‑101 (MQ=255)

                                               aaaaaCAGCACCTGAAAGGTCTGGCACATCCGCTCAAGCCAGTTGTTCTGCCTGGCAGTAATGGTTTGACCGAAGGGGTGCTGGCCGAGATTGAACAAGCg                                                  >  1:5625711/1‑101 (MQ=255)

                                               aaaaaCAGCACCTGAAAGGTCTGGCACATCCGCTCAAGCCAGTTGTTCTGCCTGGCAGTAATGGTTTGACCGAAGGGGTGCTGGCCGAGATTGAACAAGCg                                                  <  1:5625205/101‑1 (MQ=255)

                                                aaaaCAGCACCTGAAAGGTCTGGCACATCCGCTCAAGCCAGTTGTTCTGCCTGGCAGTAATGGTTTGACCGAAGGGGTGCTGGCCGAGATTGAACAAGCGt                                                 <  2:5625190/101‑1 (MQ=255)

                                                aaaaCAGCACCTGAAAGGTCTGGCACATCCGCTCAAGCCAGTTGTTCTGCCTGGCAGTAATGGTTTGACCGAAGGGGTGCTGGCCGAGATTGAACAAGCGt                                                 >  2:5625715/1‑101 (MQ=255)

                                                 aaaCAGCACCTGAAAGGTCTGGCACATCCGCTCAAGCCAGTTGTTCTGCCTGGCAGTAATGGTTTGACCGAAGGGGTGCTGGCCGAGATTGAACAAGCGtt                                                >  1:5625716/1‑101 (MQ=255)

                                                  aaCAGCACCTGAAAGGTCTGGCACATCCGCTCAAGCCAGTTGTTCTGCCTGGCAGTAATGGTTTGACCGAAGGGGTGCTGGCCGAGATTGAACAAGCGtt                                                >  1:5625721/1‑100 (MQ=255)

                                                  aaCAGCACCTGAAAGGTCTGGCACATCCGCTCAAGCCAGTTGTTCTGCCTGGCAGTAATGGTTTGACCGAAGGGGTGCTGGCCGAGATTGAACAAGCGTTa                                               <  2:5624870/101‑1 (MQ=255)

                                                     aGCACCTGAAAGGTCTGGCACATCCGCTCAAGCCAGTTGTTCTGCCTGGCAGTAATGGTTTGACCGAAGGGGTGCTGGCCGAGATTGAACAAGCGTTagag                                            <  2:5625347/101‑1 (MQ=255)

                                                      gCACCTGAAAGGTCTGGCACATCCGCTCAAGCCAGTTGTTCTGCCTGGCAGTAATGGTTTGACCGAAGGGGTGCTGGCCGAGATTGAACAAGCGTTAGAGc                                           <  2:5625491/101‑1 (MQ=255)

                                                       cACCTGAAAGGTCTGGCACATCCGCTCAAGCCAGTTGTTCTGCCTGGCAGTAATGGTTTGACCGAAGGGGTGCTGGCCGAGATTGAACAAGCGTTAGAGCa                                          >  2:5625726/1‑101 (MQ=255)

                                                       cACCTGAAAGGTCTGGCACATCCGCTCAAGCCAGTTGTTCTGCCTGGCAGTAATGGTTTGACCGAAGGGGTGCTGGCCGAGATTGAACAAGCGTTAGAGCa                                          >  1:5625728/1‑101 (MQ=255)

                                                       cACCTGAAAGGTCTGGCACATCCGCTCAAGCCAGTTGTTCTGCCTGGCAGTAATGGTTTGACCGAAGGGGTGCTGGCCGAGATTGAACAAGCGTTAGAGCa                                          <  1:5625242/101‑1 (MQ=255)

                                                          cTGAAAGGTCTGGCACATCCGCTCAAGCCAGTTGTTCTGCCTGGCAGTAATGGTTTGACCGAAGGGGTGCTGGCCGAGATTGAACAAGCGTTAGAGCACCa                                       <  2:5625386/101‑1 (MQ=255)

                                                           tGAAAGGTCTGGCACATCCGCTCAAGCCAGTTGTTCTGCCTGGCAGTAAt                                                                                         >  2:5625734/1‑50 (MQ=255)

                                                            gAAAGGTCTGGCACATCCGCTCAAGCCAGTTGTTCTGCCTGGCAGTAATGGTTTGACCGAAGGGGTGCTGGCCGAGATTGAACAAGCGTTAGAGCACCATg                                     >  1:5625736/1‑101 (MQ=255)

                                                            gAAAGGTCTGGCACATCCGCTCAAGCCAGTTGTTCTGCCTGGCAGTAATGGTTTGACCGAAGGGGTGCTGGCCGAGATTGAACAAGCGTTAGAGCACCATg                                     <  2:5625315/101‑1 (MQ=255)

                                                             aaaGGTCTGGCACATCCGCTCAAGCCAGTTGTTCTGCCTGGCAGTAATGGTTTGACCGAAGGGGTGCTGGCCGAGATTGAACAAGCGTTAGAGCACCATGa                                    >  1:5625740/1‑101 (MQ=255)

                                                             aaaGGTCTGGCACATCCGCTCAAGCCAGTTGTTCTGCCTGGCAGTAATGGTTTGACCGAAGGGGTGCTGGCCGAGATTGAACAAGCGTTAGAGCACCATGa                                    <  2:5625215/101‑1 (MQ=255)

                                                               aGGTCTGGCACATCCGCTCAAGCCAGTTGTTCTGCCTGGCAGTAATGGTTTGACCGAAGGGGTGCTGGCCGAGATTGAACAAGCGTTAGAGCACCATGAAc                                  <  2:5625230/101‑1 (MQ=255)

                                                               aGGTCTGGCACATCCGCTCAAGCCAGTTGTTCTGCCTGGCAGTAATGGTTTGACCGAAGGGGTGCTGGCCGAGATTGAACAAGCGTTAGAGCACCATGAAc                                  >  1:5625745/1‑101 (MQ=255)

                                                                ggTCTGGCACATCCGCTCAAGCCAGTTGTTCTGCCTGGCAGTAATGGTTTGACCGAAGGGGTGCTGGCCGAGATTGAACAAGCGTTAGAGCACCATGAACt                                 >  2:5625748/1‑101 (MQ=255)

                                                                  tCTGGCACATCCGCTCAAGCCAGTTGTTCTGCCTGGCAGTAATGGTTTGACCGAAGGGGTGCTGGCCGAGATTGAACAAGCGTTAGAGCACCATGAACtca                               >  1:5625752/1‑101 (MQ=255)

                                                                   cTGGCACATCCGCTCAAGCCAGTTGTTCTGCCTGGCAGTAATGGTTTGACCGAAGGGGTGCTGGCCGAGATTGAACAAGCGTTAGAGCACCATGAACtcat                              <  1:5625502/101‑1 (MQ=255)

                                                                    tGGCACATCCGCTCAAGCCAGTTGTTCTGCCTGGCAGTAATGGTTTGACCGAAGGGGTGCTGGCCGAGATTGAACAAGCGTTAGAGCACCATGAACtcatc                             <  2:5625203/101‑1 (MQ=255)

                                                                    tGGCACATCCGCTCAAGCCAGTTGTTCTGCCTGGCAGTAATGGTTTGACCGAAGGGGTGCTGGCCGAGATTGAACAAGCGTTAGAGCACCATGAACtcatc                             >  2:5625758/1‑101 (MQ=255)

                                                                    tGGCACATCCGCTCAAGCCAGTTGTTCTGCCTGGCAGTAATGGTTTGACCGAAGGGGTGCTGGCCGAGATTGAACAAGCGTTAGAGCACCATGAACtcatc                             >  1:5625759/1‑101 (MQ=255)

                                                                          aTCCGCTCAAGCCAGTTGTTCTGCCTGGCAGTAATGGTTTGACCGAAGGGGTGCTGGCCGAGATTGAACAAGCGTTAGAGCACCATGAACTCATCAAGGTg                       <  1:5625335/101‑1 (MQ=255)

                                                                          aTCCGCTCAAGCCAGTTGTTCTGCCTGGCAGTAATGGTTTGACCGAAGGGGTGCTGGCCGAGATTGAACAAGCGTTAGAGCACCATGAACTCATCAAGGTg                       >  2:5625765/1‑101 (MQ=255)

                                                                          aTCCGCTCAAGCCAGTTGTTCTGCCTGGCAGTAATGGTTTGACCGAAGGGGTGCTGGCCGAGATTGAACAAGCGTTAGAGCACCATGAACTCATCAAGGTg                       >  1:5625766/1‑101 (MQ=255)

                                                                           tCCGCTCAAGCCAGTTGTTCTGCCTGGCAGTAATGGTTTGACCGAAGGGGTGCTGGCCGAGATTGAACAAGCGTTAGAGCACCATGAACTCATCAAGGTGa                      <  2:5624950/101‑1 (MQ=255)

                                                                             cGCTCAAGCCAGTTGTTCTGCCTGGCAGTAATGGTTTGACCGAAGGGGTGCTGGCCGAGATTGAACAAGCGTTAGAGCACCATGAACTCATCAAGGTGaaa                    >  1:5625771/1‑101 (MQ=255)

                                                                              gCTCAAGCCAGTTGTTCTGCCTGGCAGTAATGGTTTGACCGAAGGGGTGCTGGCCGAGATTGAACAAGCGTTAGAGCACCATGAACTCATCAAGGTGaaaa                   >  2:5625775/1‑101 (MQ=255)

                                                                              gCTCAAGCCAGTTGTTCTGCCTGGCAGTAATGGTTTGACCGAAGGGGTGCTGGCCGAGATTGAACAAGCGTTAGAGCACCATGAACTCATCAAGGTGaaaa                   <  2:5624912/101‑1 (MQ=255)

                                                                               cTCAAGCCAGTTGTTCTGCCTGGCAGTAATGGTTTGACCGAAGGGGTGc                                                                      <  2:5624714/49‑1 (MQ=255)

                                                                                tCAAGCCAGTTGTTCTGCCTGGCAGTAATGGTTTGACCGAAGGGGTGCTGGCCGAGATTGAACAAGCGTTAGAGCACCATGAACTCATCAAGGTGAAAATc                 <  2:5625412/101‑1 (MQ=255)

                                                                                tCAAGCCAGTTGTTCTGCCTGGCAGTAATGGTTTGACCGAAGGGGTGCTGGCCGAGATTGAACAAGCGTTAGAGCACCATGAACTCATCAAGGTGAAAATc                 <  1:5625219/101‑1 (MQ=255)

                                                                                tCAAGCCAGTTGTTCTGCCTGGCAGTAATGGTTTGACCGAAGGGGTGCTGGCCGAGATTGAACAAGCGTTAGAGCACCATGAACTCATCAAGGTGAAAATc                 <  2:5625109/101‑1 (MQ=255)

                                                                                 cAAGCCAGTTGTTCTGCCTGGCAGTAATGGTTTGACCGAAGGGGTGCTGGCCGAGATTGAACAAGCGTTAGAGCACCATGAACTCATCAAGGTGAAAATCg                <  2:5625636/101‑1 (MQ=255)

                                                                                  aaGCCAGTTGTTCTGCCTGGCAGTAATGGTTTGACCGAAGGGGTGCTGGCCGAGATTGAACAAGCGTTAGAGCACCATGAACTCATCAAGGTGAAAATCGc               >  1:5625779/1‑101 (MQ=255)

                                                                                   aGCCAGTTGTTCTGCCTGGCAGTAATGGTTTGACCGAAGGGGTGCTGGCCGAGATTGAACAAGCGTTAGAGCACCATGAACTCATCAAGGTGAAAATCGcc              >  2:5625782/1‑101 (MQ=255)

                                                                                   aGCCAGTTGTTCTGCCTGGCAGTAATGGTTTGACCGAAGGGGTGCTGGCCGAGATTGAACAAGCGTTAGAGCACCATGAACTCATCAAGGTGAAAATCGcc              <  1:5625128/101‑1 (MQ=255)

                                                                                     ccAGTTGTTCTGCCTGGCAGTAATGGTTTGACCGAAGGGGTGCTGGCCGAGATTGAACAAGCGTTAGAGCACCATGAACTCATCAAGGTGAAAATCGCCAc            <  2:5625265/101‑1 (MQ=255)

                                                                                      cAGTTGTTCTGCCTGGCAGTAATGGTTTGACCGAAGGGGTGCTGGCCGAGATTGAACAAGCGTTAGAGCACCATGAACTCATCAAGGTGAAAATCGCCAcc           <  1:5624716/101‑1 (MQ=255)

                                                                                      cAGTTGTTCTGCCTGGCAGTAATGGTTTGACCGAAGGGGTGCTGGCCGAGATTGAACAAGCGTTAGAGCACCATGAACTCATCAAGGTGAAAATCGCCAcc           >  2:5625789/1‑101 (MQ=255)

                                                                                      cAGTTGTTCTGCCTGGCAGTAATGGTTTGACCGAAGGGGTGCTGGCCGAGATTGAACAAGCGTTAGAGCACCATGAACTCATCAAGGTGAAAATCGCCAcc           >  1:5625788/1‑101 (MQ=255)

                                                                                       aGTTGTTCTGCCTGGCAGTAATGGTTTGACCGAAgg                                                                           <  2:5625152/36‑1 (MQ=37)

                                                                                         ttgttCTGCCTGGCAGTAATGGTTTGACCGAAGGGGTGCTGGCCGAGATTGAACAAGCGTTAGAGCACCATGAACTCATCAAGGTGAAAATCGCCACCGaa        <  1:5625261/101‑1 (MQ=255)

                                                                                          tgttCTGCCTGGCAGTAATGGTTTGACCGAAGGGGTGCTGGCCGAGATTGAACAAGCGTTAGAGCACCATGAACTCATCAAGGTGAAAATCGCCACCGAAg       >  1:5625792/1‑101 (MQ=255)

                                                                                          tgttCTGCCTGGCAGTAATGGTTTGACCGAAGGGGTGCTGGCCGAGATTGAACAAGCGTTAGAGCACCATGAACTCATCAAGGTGAAAATCGCCACCGAAg       <  1:5625067/101‑1 (MQ=255)

                                                                                           gttCTGCCTGGCAGTAATGGTTTGACCGAAGGGGTGCTGGCCGAGATTGAACAAGCGTTAGAGCACCATGAACTCATCAAGGTGAAAATCGCCACCGAAGa      >  1:5625797/1‑101 (MQ=255)

                                                                                             tCTGCCTGGCAGTAATGGTTTGACCGAAGGGGTGCTGGCCGAGATTGAACAAGCGTTAGAGCACCATGAACTCATCAAGGTGAAAATCGCCACCGAAGATc    <  2:5625431/101‑1 (MQ=255)

TTCAATCCCAACGTAAGCAAAAAATACGATGAATCTGAGTACTAAACAAAAACAGCACCTGAAAGGTCTGGCACATCCGCTCAAGCCAGTTGTTCTGCTTGGCAGTAATGGTTTGACCGAAGGGGTGCTGGCCGAGATTGAACAAGCGTTAGAGCACCATGAACTCATCAAGGTGAAAATCGCCACCGAAGATCGC  >  NC_000913/3327762‑3327957

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