BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,365,444	G→A	pseudogene (260/371 nt)	yhcE →	putative uncharacterized protein YhcE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,365,444	0	G	A	100.0%	1356.7 / NA	309	pseudogene (260/371 nt)	yhcE	putative uncharacterized protein YhcE
Reads supporting (aligned to +/- strand): ref base G (0/0); new base A (171/138); total (171/138)

Only 100 of 309 total aligned reads displayed.

AATCACGCTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGGTTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGA  >  NC_000913/3365346‑3365541

aaTCACGCTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGAtt                                                                                                 >  1:5509494/1‑101 (MQ=255)

aaTCACGCTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGAtt                                                                                                 >  2:5509491/1‑101 (MQ=255)

 aTCACGCTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGATTc                                                                                                >  2:5509497/1‑101 (MQ=255)

    aCGCTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGATTCAtt                                                                                             >  1:5509500/1‑101 (MQ=255)

     cGCTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGATTCATTa                                                                                            <  1:5508913/101‑1 (MQ=255)

      gCTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGATTCATTaa                                                                                           <  1:5508854/101‑1 (MQ=255)

       cTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGATTCATTaaa                                                                                          <  1:5509346/101‑1 (MQ=255)

       cTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGATTCATTaaa                                                                                          >  2:5509504/1‑101 (MQ=255)

         ggTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGATTCATTAAATc                                                                                        >  1:5509508/1‑101 (MQ=255)

         ggTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGATTCATTAAATc                                                                                        <  1:5508990/101‑1 (MQ=255)

          gTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGATTCATTAAATCt                                                                                       <  2:5508874/101‑1 (MQ=255)

           tAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGATTCATTAAATCTg                                                                                      >  2:5509511/1‑101 (MQ=255)

           tAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGATTCATTAAATCTg                                                                                      <  1:5508938/101‑1 (MQ=255)

            aatgcaatgcGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGATTCATTAAATCTGa                                                                                     >  1:5509514/1‑101 (MQ=255)

             atgcaatgcGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGATTCATTAAATCTGAt                                                                                    <  2:5509322/101‑1 (MQ=255)

                caatgcGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGATTCATTAAATCTGATCGc                                                                                 <  1:5509060/101‑1 (MQ=255)

                caatgcGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGATTCATTAAATCTGATCGc                                                                                 >  1:5509518/1‑101 (MQ=255)

                caatgcGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGATTCATTAAATCTGATCGc                                                                                 <  1:5508969/101‑1 (MQ=255)

                  atgcGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGATTCATTAAATCTGATCGCCa                                                                               <  1:5509076/101‑1 (MQ=255)

                  atgcGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGATTCATTAAATCTGATCGCCa                                                                               >  1:5509521/1‑101 (MQ=255)

                   tgcGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGATTCATTAAATCTGATCGCCAg                                                                              >  1:5509524/1‑101 (MQ=255)

                    gcGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGATTCATTAAATCTGATCGCCAGc                                                                             >  1:5509527/1‑101 (MQ=255)

                     cGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGATTCATTAAATCTGATCGCCAg                                                                              <  2:5509016/99‑1 (MQ=255)

                     cGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGATTCATTAAATCTGATCGCCAGCg                                                                            <  2:5509230/101‑1 (MQ=255)

                      gATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGATTCATTAAATCTGATCGCCAGc                                                                             <  1:5508958/99‑1 (MQ=255)

                      gATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGATTCATTAAATCTGATCGCCAGCGa                                                                           >  2:5509528/1‑101 (MQ=255)

                        tCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGATTCATTAAATCTGATCGCCAGCGATg                                                                         >  2:5509533/1‑101 (MQ=255)

                         cAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGATTCATTAAATCTGATCGCCAGCGATGc                                                                        >  1:5509535/1‑101 (MQ=255)

                          aGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGATTCATTAAATCTGATCGCCAGCGATGCa                                                                       <  2:5509035/101‑1 (MQ=255)

                            gACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGATTCATTAAATCTGATCGCCAGCGATGCAtt                                                                     >  2:5509538/1‑101 (MQ=255)

                             aCGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGATTCATTAAATCTGATCGCCAGCGATGCAttt                                                                    <  2:5509114/101‑1 (MQ=255)

                              cGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGATTCATTAAATCTGATCGCCAGCGATGCAtttt                                                                   <  1:5509068/101‑1 (MQ=255)

                               gCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGATTCATTAAATCTGATCGCCAGCGATGCATTTTc                                                                  >  1:5509542/1‑101 (MQ=255)

                               gCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGATTCATTAAATCTGATCGCCAGCGATGCATTTTc                                                                  <  1:5509117/101‑1 (MQ=255)

                                cTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGATTCATTAAATCTGATCGCCAGCGATGCATTTTCa                                                                 <  2:5509234/101‑1 (MQ=255)

                                cTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGATTCATTAAATCTGATCGCCAGCGATGCATTTTCa                                                                 <  2:5508980/101‑1 (MQ=255)

                                cTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGATTCATTAAATCTGATCGCCAGCGATGCATTTTCa                                                                 >  2:5509546/1‑101 (MQ=255)

                                 tCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGATTCATTAAATCTGATCGCCAGCGATGCATTTTCAg                                                                >  1:5509549/1‑101 (MQ=255)

                                   atatTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGATTCATTAAATCTGATCGCCAGCGATGCATTTTCAGgt                                                              <  2:5508942/101‑1 (MQ=255)

                                     atTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGATTCATTAAATCTGATCGCCAGCGATGCATTTTCAGgtag                                                            >  2:5509554/1‑101 (MQ=255)

                                       tCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGATTCATTAAATCTGATCGCCAGCGATGCAtttt                                                                   >  1:5509557/1‑92 (MQ=255)

                                         gtgtTCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGATTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAAc                                                        >  1:5509563/1‑101 (MQ=255)

                                          tgtTCAGGTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGATTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACa                                                       >  2:5509565/1‑101 (MQ=255)

                                           gtTCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGATTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACaa                                                      <  1:5509240/101‑1 (MQ=255)

                                            ttCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGATTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAAt                                                     <  1:5509189/101‑1 (MQ=255)

                                             tCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGATTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAAtc                                                    >  1:5509569/1‑101 (MQ=255)

                                             tCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGATTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAAtc                                                    >  2:5509571/1‑101 (MQ=255)

                                              cAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGATTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAAtct                                                   <  2:5508996/101‑1 (MQ=255)

                                               aGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGATTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAAtctc                                                  >  1:5509574/1‑101 (MQ=255)

                                                gTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGATTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCg                                                 >  2:5509576/1‑101 (MQ=255)

                                                 ttGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGATTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGa                                                <  1:5508923/101‑1 (MQ=255)

                                                  tGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGATTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGAc                                               >  2:5509580/1‑101 (MQ=255)

                                                   gAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGATTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGAcc                                              >  2:5509583/1‑101 (MQ=255)

                                                     aTACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGATTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTg                                            >  1:5509587/1‑101 (MQ=255)

                                                     aTACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGATTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTg                                            >  1:5509590/1‑101 (MQ=255)

                                                      tACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGATTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTgg                                           <  1:5508694/101‑1 (MQ=255)

                                                       aCCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGATTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGc                                          <  1:5509043/101‑1 (MQ=255)

                                                       aCCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGATTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGc                                          >  1:5509593/1‑101 (MQ=255)

                                                        ccGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGATTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCg                                         >  2:5509595/1‑101 (MQ=255)

                                                         cGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGATTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGAc                                               >  1:5509601/1‑94 (MQ=255)

                                                         cGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGATTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGt                                        >  2:5509598/1‑101 (MQ=255)

                                                           cTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGATTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTa                                      <  2:5509275/101‑1 (MQ=255)

                                                            tAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGATTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAg                                     <  2:5509123/101‑1 (MQ=255)

                                                             aaTATTTCCAACAACGGATTTGGATTTGGCCCGAATGATTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGc                                    >  1:5509605/1‑101 (MQ=255)

                                                              atatTTCCAACAACGGATTTGGATTTGGCCCGAATGATTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGcc                                   >  1:5509609/1‑101 (MQ=255)

                                                                atTTCCAACAACGGATTTGGATTTGGCCCGAATGATTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAg                                 >  2:5509612/1‑101 (MQ=255)

                                                                 tttCCAACAACGGATTTGGATTTGGCCCGAATGATTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGt                                <  2:5508964/101‑1 (MQ=255)

                                                                  ttCCAACAACGGATTTGGATTTGGCCCGAATGATTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTg                               <  2:5509082/101‑1 (MQ=255)

                                                                  ttCCAACAACGGATTTGGATTTGGCCCGAATGATTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTg                               <  1:5509091/101‑1 (MQ=255)

                                                                    ccAACAACGGATTTGGATTTGGCCCGAATGATTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGc                             >  1:5509616/1‑101 (MQ=255)

                                                                      aacaacGGATTTGGATTTGGCCCGAATGATTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCaa                           >  2:5509619/1‑101 (MQ=255)

                                                                        caacGGATTTGGATTTGGCCCGAATGATTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCaaaa                         >  1:5509622/1‑101 (MQ=255)

                                                                        caacGGATTTGGATTTGGCCCGAATGATTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCaaaa                         >  2:5509623/1‑101 (MQ=255)

                                                                          acGGATTTGGATTTGGCCCGAATGATTCATTAAATCTGATCGCCAGCGATg                                                                         >  2:5509629/1‑51 (MQ=255)

                                                                          acGGATTTGGATTTGGCCCGAATGATTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAAt                       <  1:5509284/101‑1 (MQ=255)

                                                                          acGGATTTGGATTTGGCCCGAATGATTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAAt                       >  1:5509628/1‑101 (MQ=255)

                                                                          acGGATTTGGATTTGGCCCGAATGATTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAAt                       <  2:5509046/101‑1 (MQ=255)

                                                                            ggatttggatttGGCCCGAATGATTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGa                     >  1:5509635/1‑101 (MQ=255)

                                                                            ggatttggatttGGCCCGAATGATTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGa                     <  1:5509026/101‑1 (MQ=255)

                                                                            ggatttggatttGGCCCGAATGATTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGa                     >  2:5509634/1‑101 (MQ=255)

                                                                              atttggatttggCCCGAATGATTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACa                   <  2:5509372/101‑1 (MQ=255)

                                                                               tttggatttggCCCGAATGATTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACaa                  >  1:5509639/1‑101 (MQ=255)

                                                                                ttggatttggCCCGAATGATTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAAt                 <  1:5509175/101‑1 (MQ=255)

                                                                                  ggatttggCCCGAATGATTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATcc               >  2:5509644/1‑101 (MQ=255)

                                                                                  ggatttggCCCGAATGATTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATcc               >  1:5509643/1‑101 (MQ=255)

                                                                                   gatttggCCCGAATGATTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCg              >  1:5509648/1‑101 (MQ=255)

                                                                                      ttggCCCGAATGATTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGgc           <  2:5509067/101‑1 (MQ=255)

                                                                                      ttggCCCGAATGATTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGgc           <  1:5509157/101‑1 (MQ=255)

                                                                                      ttggCCCGAATGATTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGgc           <  1:5509055/101‑1 (MQ=255)

                                                                                      ttggCCCGAATGATTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGgc           >  2:5509651/1‑101 (MQ=255)

                                                                                       tggCCCGAATGATTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGgca          <  1:5509145/101‑1 (MQ=255)

                                                                                         gCCCGAATGATTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGgcagc        <  2:5509194/101‑1 (MQ=255)

                                                                                         gCCCGAATGATTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGgcagc        <  2:5509167/101‑1 (MQ=255)

                                                                                         gCCCGAATGATTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGgcagc        >  1:5509657/1‑101 (MQ=255)

                                                                                         gCCCGAATGATTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGgcagc        >  1:5509654/1‑101 (MQ=255)

                                                                                          cccGAATGATTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGgcagca       <  1:5509205/101‑1 (MQ=255)

                                                                                            cGAATGATTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGcaca     <  2:5509143/101‑1 (MQ=255)

                                                                                             gAATGATTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAg    >  1:5509664/1‑101 (MQ=255)

                                                                                             gAATGATTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAg    >  1:5509661/1‑101 (MQ=255)

                                                                                               aTGATTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGa  >  1:5509669/1‑101 (MQ=255)

AATCACGCTGGTAATGCAATGCGATCAGGACGCTCATATTCGTGTTCAGTTGAATACCGCTAATATTTCCAACAACGGATTTGGATTTGGCCCGAATGGTTCATTAAATCTGATCGCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGA  >  NC_000913/3365346‑3365541

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: `‑`