BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,666,629	T→C	T64A (ACG→GCG)	yhjA ←	putative cytochrome C peroxidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,666,629	0	T	C	100.0%	1451.3 / NA	459	T64A (ACG→GCG)	yhjA	putative cytochrome C peroxidase
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (226/233); total (226/233)

Only 100 of 485 total aligned reads displayed.

TTAAAAGATTTATATCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGTGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATCGGCCTGTTTAC  >  NC_000913/3666538‑3666718

 tAAAAGATTTATATCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCgcgt                                                                                         >  1:6993280/1‑93 (MQ=255)

  aaaaGATTTATATCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGcgtg                                                                                        <  1:5410624/93‑1 (MQ=255)

   aaaGATTTATATCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTgg                                                                                       >  2:1615182/1‑93 (MQ=255)

      gATTTATATCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCaa                                                                                    >  1:1070598/1‑93 (MQ=255)

       aTTTATATCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAAt                                                                                   >  1:4661805/1‑93 (MQ=255)

       aTTTATATCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAAt                                                                                   >  2:3977147/1‑93 (MQ=255)

       aTTTATATCCAAGATTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAAt                                                                                   <  2:6932878/93‑1 (MQ=255)

        tttATATCCAAGCTTACTGGCGTAATCCATCACCTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATa                                                                                  <  2:2695859/93‑1 (MQ=255)

         ttATATCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAg                                                                                 <  2:1901994/93‑1 (MQ=255)

          tataTCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAga                                                                                >  1:4423425/1‑92 (MQ=255)

           atatCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTc                                                                               >  1:1841957/1‑93 (MQ=255)

           atatCCAAGCTTAATGTCGTAATCCATCAACTGTTGCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTc                                                                               <  1:4837590/93‑1 (MQ=255)

            tatCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCg                                                                              <  2:6249348/93‑1 (MQ=255)

             atCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGc                                                                             >  2:3355049/1‑93 (MQ=255)

              tCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCa                                                                            >  1:7115756/1‑93 (MQ=255)

              tCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCa                                                                            >  2:223880/1‑93 (MQ=255)

              tCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCa                                                                            <  1:1691885/93‑1 (MQ=255)

                cAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATc                                                                          >  2:6673769/1‑93 (MQ=255)

                   gCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCtt                                                                       >  1:4549840/1‑93 (MQ=255)

                     ttAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGTGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCtttt                                                                     <  2:6066501/93‑1 (MQ=255)

                     ttAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCtttt                                                                     >  1:5233386/1‑93 (MQ=255)

                           tCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCgg                                                               >  2:4915433/1‑93 (MQ=255)

                           tCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCgg                                                               <  2:4740669/93‑1 (MQ=255)

                             gTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGgag                                                             >  1:2092522/1‑93 (MQ=255)

                              tAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGgaga                                                            <  2:2520602/93‑1 (MQ=255)

                              tAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGgaga                                                            <  1:990185/93‑1 (MQ=255)

                              tAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGgaga                                                            <  1:5082598/93‑1 (MQ=255)

                               aaTCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGaa                                                           <  2:731901/93‑1 (MQ=255)

                               aaTCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGaa                                                           <  1:4155955/93‑1 (MQ=255)

                                aTCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGaaa                                                          >  2:362758/1‑93 (MQ=255)

                                aTCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGaaa                                                          <  1:7393712/93‑1 (MQ=255)

                                aTCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGaaa                                                          <  2:4614819/93‑1 (MQ=255)

                                  ccATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGc                                                        <  1:5686445/93‑1 (MQ=255)

                                  ccATCAACTGTTTCGCGCAAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGc                                                        <  2:5902135/93‑1 (MQ=255)

                                   cATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGcc                                                       <  1:1760632/93‑1 (MQ=255)

                                     tCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTa                                                     >  1:5922406/1‑93 (MQ=255)

                                     tCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTa                                                     >  2:4869718/1‑93 (MQ=255)

                                     acAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTa                                                     <  2:5697709/92‑1 (MQ=255)

                                      cAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTaa                                                    <  2:4105209/93‑1 (MQ=255)

                                      cAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTaa                                                    >  2:211924/1‑93 (MQ=255)

                                      cAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTaa                                                    <  1:2983346/93‑1 (MQ=255)

                                       aaCTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTaaa                                                   >  1:1932825/1‑93 (MQ=255)

                                       aaCTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTaaa                                                   <  2:6811018/93‑1 (MQ=255)

                                       aaCTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTaaa                                                   >  2:7000866/1‑93 (MQ=255)

                                           gTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCt                                               >  1:5548474/1‑93 (MQ=255)

                                           gTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCt                                               <  2:1161754/93‑1 (MQ=255)

                                            tttCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCtt                                              >  1:6204338/1‑93 (MQ=255)

                                                gcgcCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCttatta                                          <  2:888630/93‑1 (MQ=255)

                                                 cgcCAGGAATATAATAATAGGCGTGTAAGTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTtattat                                         <  2:7561952/93‑1 (MQ=255)

                                                 cgcCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTtattat                                         <  1:5734937/93‑1 (MQ=255)

                                                 cgcCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTtattat                                         >  2:1932653/1‑93 (MQ=255)

                                                 cgcCAGAAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTtattat                                         <  1:83391/93‑1 (MQ=255)

                                                   ccAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTAttt                                       <  1:3861889/93‑1 (MQ=255)

                                                    cAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTAtttt                                      >  2:6614383/1‑93 (MQ=255)

                                                     aGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTc                                     >  2:2357846/1‑93 (MQ=255)

                                                      ggAATATAATAATGGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCa                                    <  1:2289138/93‑1 (MQ=255)

                                                      ggAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGCAAGCCTAAAACCTTATTATTTTCa                                    >  2:1795489/1‑93 (MQ=255)

                                                      ggAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCa                                    >  1:4073325/1‑93 (MQ=255)

                                                      ggAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCa                                    >  1:7553148/1‑93 (MQ=255)

                                                       gAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCAc                                   <  2:4261666/93‑1 (MQ=255)

                                                       gAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCAc                                   >  1:3698587/1‑93 (MQ=255)

                                                       gAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCAc                                   >  1:7282678/1‑93 (MQ=255)

                                                       gAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCAc                                   <  2:1124458/93‑1 (MQ=255)

                                                        aaTATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACt                                  >  2:541129/1‑93 (MQ=255)

                                                        aaTATAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACt                                  <  1:1513221/93‑1 (MQ=255)

                                                         atatAATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTg                                 <  2:6447130/93‑1 (MQ=255)

                                                          tataATAATAGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGa                                >  2:6346869/1‑93 (MQ=255)

                                                            taataataGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGAca                              <  2:2827833/93‑1 (MQ=255)

                                                               taataGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGAcagca                           >  2:2928494/1‑93 (MQ=255)

                                                                aataGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGAcagcag                          >  1:6663776/1‑93 (MQ=255)

                                                                 ataGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGa                         >  1:1185091/1‑93 (MQ=255)

                                                                 ataGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGa                         <  1:3235847/93‑1 (MQ=255)

                                                                 ataGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGa                         <  2:5103943/93‑1 (MQ=255)

                                                                 ataGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGa                         <  1:2369382/93‑1 (MQ=255)

                                                                 ataGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGa                         <  1:5288983/93‑1 (MQ=255)

                                                                 ataGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGa                         <  1:7832969/93‑1 (MQ=255)

                                                                  taGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGAt                        <  1:251239/93‑1 (MQ=255)

                                                                   aGGCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATg                       >  1:13083/1‑93 (MQ=255)

                                                                     gCGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGcc                     <  1:2821833/93‑1 (MQ=255)

                                                                      cGGGTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCt                    <  2:3127849/93‑1 (MQ=255)

                                                                       gggTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTg                   >  1:7688746/1‑93 (MQ=255)

                                                                        ggTAATTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGa                  >  1:6240740/1‑93 (MQ=255)

                                                                           aaTTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACa               >  2:1356213/1‑93 (MQ=255)

                                                                           aaTTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACa               >  1:641550/1‑93 (MQ=255)

                                                                            aTTCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACAt              >  1:6830982/1‑93 (MQ=255)

                                                                             ttCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATc             <  2:5354884/93‑1 (MQ=255)

                                                                             ttCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATc             >  1:3489674/1‑93 (MQ=255)

                                                                             ttCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATc             >  2:1232391/1‑93 (MQ=255)

                                                                             ttCTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATc             <  2:1716711/93‑1 (MQ=255)

                                                                               cTGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATCgg           <  1:3097911/93‑1 (MQ=255)

                                                                                tGCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATCGGc          <  2:559760/93‑1 (MQ=255)

                                                                                 gCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATCGGcc         <  2:7062750/93‑1 (MQ=255)

                                                                                 gCCGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATCGGcc         >  2:7176868/1‑93 (MQ=255)

                                                                                  ccGAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATCGGCCt        <  2:1461935/93‑1 (MQ=255)

                                                                                    gAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATCGGCCTGt      >  1:7877007/1‑93 (MQ=255)

                                                                                    gAAGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATCGGCCTGt      <  2:7389448/93‑1 (MQ=255)

                                                                                     aaGGCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATCGGCCTGtt     <  1:6887018/93‑1 (MQ=255)

                                                                                       ggCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATCGGCCTGTTTa   <  2:3829536/93‑1 (MQ=255)

                                                                                       ggCGCGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATAATTTTCACTGACAGCAGATGCCTGAACATCGGCCTGtttc   >  1:5629697/1‑92 (MQ=255)

                                                                                        gcgcGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATCGGCCTGTTTAc  >  1:269440/1‑93 (MQ=255)

TTAAAAGATTTATATCCAAGCTTAATGTCGTAATCCATCAACTGTTTCGCGCCAGGAATATAATAATAGGCGGGTAATTCTGCCGAAGGCGTGTGGCAATAGTCGCATCCTTTTTCGCGGAGAAAGCCTAAAACCTTATTATTTTCACTGACAGCAGATGCCTGAACATCGGCCTGTTTAC  >  NC_000913/3666538‑3666718

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