BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,765,413	T→C	T74A (ACC→GCC)	menI ←	1,4‑dihydroxy‑2‑naphthoyl‑CoA hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,765,413	0	T	C	100.0%	1397.5 / NA	333	T74A (ACC→GCC)	menI	1,4‑dihydroxy‑2‑naphthoyl‑CoA hydrolase
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (151/182); total (151/182)

Only 100 of 336 total aligned reads displayed.

ATGCAACGGTTTGCATACGCCGCGCACCCGCCCTTCTCGTGCCGAGCGGACGTGGTTAGCATTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGTACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTACCACGGATGCTCCTCCATGCAGCAACCCGAAAGGCTGCTTTGTCCGCGAGTCTACT  >  NC_000913/1765315‑1765510

   cAACGGTTTGCATACGCCGCGCACCCGCCCTTCTCGTGCCGAGCGGACGTGGTTAGCATTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGCACata                                                                                              <  2:3385182/101‑1 (MQ=255)

   cAACGGTTTGCATACGCCGCGCACCCGCCCTTCTCGTGCCGAGCGGACGTGGTTAGCATTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGCACata                                                                                              <  2:3385390/101‑1 (MQ=255)

   cAACGGTTTGCATACGCCGCGCACCCGCCCTTCTCGTGCCGAGCGGACGTGGTTAGCATTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGCACata                                                                                              <  1:3385195/101‑1 (MQ=255)

    aaCGGTTTGCATACGCCGCGCACCCGCCCTTCTCGTGCCGAGCGGACGTGGTTAGCATTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGCACataa                                                                                             <  1:3384783/101‑1 (MQ=255)

    aaCGGTTTGCATACGCCGCGCACCCGCCCTTCTCGTGCCGAGCGGACGTGGTTAGCATTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGCACataa                                                                                             >  2:3385716/1‑101 (MQ=255)

    aaCGGTTTGCATACGCCGCGCACCCGCCCTTCTCGTGCCGAGCGGACGTGGTTAGCATTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGCACataa                                                                                             >  1:3385717/1‑101 (MQ=255)

     aCGGTTTGCATACGCCGCGCACCCGCCCTTCTCGTGCCGAGCGGACGTGGTTAGCATTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGCACataaa                                                                                            <  1:3385147/101‑1 (MQ=255)

        gTTTGCATACGCCGCGCACCCGCCCTTCTCGTGCCGAGCGGACGTGGTTAGCATTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGCACataaataa                                                                                         <  2:3385300/101‑1 (MQ=255)

        gTTTGCATACGCCGCGCACCCGCCCTTCTCGTGCCGAGCGGACGTGGTTAGCATTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGCACataaataa                                                                                         >  1:3385723/1‑101 (MQ=255)

         tttGCATACGCCGCGCACCCGCCCTTCTCGTGCCGAGCGGACGTGGTTAGCATTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGCACataaataa                                                                                         >  2:3385725/1‑100 (MQ=255)

         tttGCATACGCCGCGCACCCGCCCTTCTCGTGCCGAGCGGACGTGGTTAGCATTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGCACATAAATAAc                                                                                        <  2:3385355/101‑1 (MQ=255)

          ttGCATACGCCGCGCACCCGCCCTTCTCGTGCCGAGCGGACGTGGTTAGCATTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGCACATAAATAAcc                                                                                       <  1:3384939/101‑1 (MQ=255)

          ttGCATACGCCGCGCACCCGCCCTTCTCGTGCCGAGCGGACGTGGTTAGCATTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGCACATAAATAAcc                                                                                       <  2:3385248/101‑1 (MQ=255)

          ttGCATACGCCGCGCACCCGCCCTTCTCGTGCCGAGCGGACGTGGTTAGCATTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGCACATAAATAAcc                                                                                       <  1:3385232/101‑1 (MQ=255)

            gCATACGCCGCGCACCCGCCCTTCTCGTGCCGAGCGGACGTGGTTAGCATTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGCACATAAATAACCgg                                                                                     >  1:3385727/1‑101 (MQ=255)

             cATACGCCGCGCACCCGCCCTTCTCGTGCCGAGCGGACGTGGTTAGCATTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGCACATAAATAACCGGc                                                                                    <  2:3385009/101‑1 (MQ=255)

               tACGCCGCGCACCCGCCCTTCTCGTGCCGAGCGGACGTGGTTAGCATTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGCACATAAATAACCGGCaa                                                                                  <  1:3385069/101‑1 (MQ=255)

               tACGCCGCGCACCCGCCCTTCTCGTGCCGAGCGGACGTGGTTAGCATTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGCACATAAATAACCGGCaa                                                                                  >  2:3385730/1‑101 (MQ=255)

                aCGCCGCGCACCCGCCCTTCTCGTGCCGAGCGGACGTGGTTAGCATTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGCACATAAATAACCGGCAAc                                                                                 <  2:3385203/101‑1 (MQ=255)

                aCGCCGCGCACCCGCCCTTCTCGTGCCGAGCGGACGTGGTTAGCATTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGCACATAAATAACCGGCAAc                                                                                 <  1:3385424/101‑1 (MQ=255)

                aCGCCGCGCACCCGCCCTTCTCGTGCCGAGCGGACGTGGTTAGCATTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGCACATAAATAACCGGCAAc                                                                                 <  2:3385074/101‑1 (MQ=255)

                aCGCCGCGCACCCGCCCTTCTCGTGCCGAGCGGACGTGGTTAGCATTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGCACATAAATAACCGGCAAc                                                                                 <  2:3385222/101‑1 (MQ=255)

                   ccgcGCACCCGCCCTTCTCGTGCCGAGCGGACGTGGTTAGCATTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGCACATAAATAACCGGCAACGGa                                                                              <  2:3385097/101‑1 (MQ=255)

                    cgcgCACCCGCCCTTCTCGTGCCGAGCGGACGTGGTTAGCATTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGCACATAAATAACCGGCAACGGaa                                                                             >  1:3385735/1‑101 (MQ=255)

                     gcgcACCCGCCCTTCTCGTGCCGAGCGGACGTGGTTAGCATTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGCACATAAATAACCGGCAACGGAAc                                                                            >  2:3385737/1‑101 (MQ=255)

                       gcACCCGCCCTTCTCGTGCCGAGCGGACGTGGTTAGCATTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGCACATAAATAACCGGCAACGGAACCg                                                                          <  2:3385038/101‑1 (MQ=255)

                       gcACCCGCCCTTCTCGTGCCGAGCGGACGTGGTTAGCATTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGCACATAAATAACCGGCAACGGAACCg                                                                          <  1:3385025/101‑1 (MQ=255)

                        cACCCGCCCTTCTCGTGCCGAGCGGACGTGGTTAGCATTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGCACATAAATAACCGGCAACGGAACCGa                                                                         >  2:3385742/1‑101 (MQ=255)

                        cACCCGCCCTTCTCGTGCCGAGCGGACGTGGTTAGCATTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGCACATAAATAACCGGCAACGGAACCGa                                                                         <  1:3385247/101‑1 (MQ=255)

                         aCCCGCCCTTCTCGTGCCGAGCGGACGTGGTTAGCATTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGCACATAAATAACCGGCAACGGAACCGAt                                                                        >  1:3385743/1‑101 (MQ=255)

                          cccGCCCTTCTCGTGCCGAGCGGACGTGGTTAGCATTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGCACATAAATAACCGGCAACGGAAc                                                                            <  2:3384986/96‑1 (MQ=255)

                          cccGCCCTTCTCGTGCCGAGCGGACGTGGTTAGCATTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGCACATAAATAACCGGCAACGGAACCGATa                                                                       <  1:3385053/101‑1 (MQ=255)

                           ccGCCCTTCTCGTGCCGAGCGGACGTGGTTAGCATTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGCACATAAATAACCGGCAACGGAACCGATAc                                                                      >  1:3385748/1‑101 (MQ=255)

                           ccGCCATTCTCGTGCCGAGCGGACGTGGTTAGCATTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGCACataaataa                                                                                         <  1:3384981/82‑1 (MQ=255)

                            cGCCCTTCTCGTGCCGAGCGGACGTGGTTAGCATTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGCACATAAATAACCGGCAACGGAACCGATACt                                                                     >  2:3385749/1‑101 (MQ=255)

                              cccTTCTCGTGCCGAGCGGACGTGGTTAGCATTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGCACATAAATAACCGGCAACGGAACCGATACttt                                                                   >  1:3385755/1‑101 (MQ=255)

                                  tctcGTGCCGAGCGGACGTGGTTAGCATTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGCACATAAATAACCGGCAACGGAACCGATACTTTCGGc                                                               >  1:3385759/1‑101 (MQ=255)

                                   ctcGTGCCGAGCGGACGTGGTTAGCATTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGCACATAAATAACCGGCAACGGAACCGATACTTTCGGcc                                                              >  2:3385761/1‑101 (MQ=255)

                                   ctcGTGCCGAGCGGACGTGGTTAGCATTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGCACATAAATAACCGGCAACGGAACCGATACTTTCGGcc                                                              <  1:3385259/101‑1 (MQ=255)

                                    tcGTGCCGAGCGGACGTGGTTAGCATTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGCACATAAATAACCGGCAACGGAACCGATACTTTCGGCCa                                                             >  1:3385763/1‑101 (MQ=255)

                                       tGCCGAGCGGACGTGGTTAGCATTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGCACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTa                                                          >  1:3385767/1‑101 (MQ=255)

                                       tGCCGAGCGGACGTGGTTAGCATTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGCACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTa                                                          <  1:3385185/101‑1 (MQ=255)

                                        gCCGAGCGGACGTGGTTAGCATTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGCACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTAc                                                         >  2:3385771/1‑101 (MQ=255)

                                         ccGAGCGGACGTGGTTAGCATTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGCACATAAATAACCg                                                                                      >  2:3385775/1‑71 (MQ=255)

                                          cGAGCGGACGTGGTTAGCATTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGCACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTACCa                                                       >  1:3385776/1‑101 (MQ=255)

                                          cGAGCGGACGTGGTTAGCATTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGCACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTACCa                                                       <  1:3385269/101‑1 (MQ=255)

                                           gAGCGGACGTGGTTAGCATTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGCACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTACCAc                                                      >  1:3385779/1‑101 (MQ=255)

                                           gAGCGGACGTGGTTAGCATTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGCACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTACCAc                                                      <  2:3385286/101‑1 (MQ=255)

                                            aGCGGACGTGGTTAGCATTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGCACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTACCACg                                                     <  2:3385334/101‑1 (MQ=255)

                                             gCGGACGTGGTTAGCATTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGCACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTACCACgg                                                    <  1:3385202/101‑1 (MQ=255)

                                              cGGACGTGGTTAGCATTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGCACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTACCACGGa                                                   >  2:3385784/1‑101 (MQ=255)

                                              cGGACGTGGTTAGCATTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGCACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTACCACGGa                                                   <  1:3385318/101‑1 (MQ=255)

                                               ggACGTGGTTAGCATTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGCACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTACCACGGAt                                                  >  1:3385788/1‑101 (MQ=255)

                                               ggACGTGATTAGCATTGATTTCCAGCCCAACCACTTTTTGCTCACCTTCGGCACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTACCACGGAt                                                  <  2:3385278/101‑1 (MQ=255)

                                                 aCGTGGTTAGCATTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGCACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTACCACGGATGc                                                >  1:3385792/1‑101 (MQ=255)

                                                  cGTGGTTAGCATTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGCACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTACCACGGATGct                                               <  1:3385516/101‑1 (MQ=255)

                                                  cGTGGTTAGCATTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGCACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTACCACGGATGct                                               >  2:3385794/1‑101 (MQ=255)

                                                   gTGGTTAGCATTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGCACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTACCACGGATGctc                                              <  1:3385351/101‑1 (MQ=255)

                                                    tGGTTAGCATTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGCACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTACCACGGATGctcc                                             >  1:3385796/1‑101 (MQ=255)

                                                     ggTTAGCATTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGCACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTACCACGGATGctcct                                            <  1:3385104/101‑1 (MQ=255)

                                                         aGCATTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGCACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTACCACGGATGCTCCTCCAt                                        <  1:3385596/101‑1 (MQ=255)

                                                         aGCATTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGCACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTACCACGGATGCTCCTCCAt                                        <  2:3385322/101‑1 (MQ=255)

                                                          gCATTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGCACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTACCACGGATGCTCCTCCATg                                       >  1:3385803/1‑101 (MQ=255)

                                                          gCATTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGCACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTACCACGGATGCTCCTCCATg                                       <  2:3385193/101‑1 (MQ=255)

                                                           cATTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGCACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTACCACGGATGCTCCTCCATgc                                      >  2:3385804/1‑101 (MQ=255)

                                                            aTTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGCACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTACCACGGATGCTCCTCCATgca                                     <  1:3385371/101‑1 (MQ=255)

                                                            aTTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGCACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTACCACGGATGCTCCTCCATgca                                     <  2:3385132/101‑1 (MQ=255)

                                                             ttGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGCACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTACCACGGATGct                                               >  2:3385809/1‑90 (MQ=255)

                                                              tGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGCACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTACCACGGATGCTCCTCCATgcagc                                   <  1:3385279/101‑1 (MQ=255)

                                                               gATTTCCAGACCAACCACTTTTTGCTCACCTTCGGCACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTACCACGGATGCTCCTCCATgcagca                                  >  2:3385812/1‑101 (MQ=255)

                                                               gATTTCCAGACCAACCACTTTTTGCTCACCTTCGGCACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTACCACGGATGCTCCTCCATgcagca                                  <  2:3385508/101‑1 (MQ=255)

                                                                aTTTCCAGACCAACCACTTTTTGCTCACCTTCGGCACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTACCACCGATGCTCCTCCATGCAGCaa                                 <  1:3385228/101‑1 (MQ=255)

                                                                  ttCCAGACCAACCACTTTTTGCTCACCTTCGGCACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTACCACGGATGCTCCTCCATGCAGCAAcc                               >  1:3385814/1‑101 (MQ=255)

                                                                   tCCAGACCAACCACTTTTTGCTCACCTTCGGCACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTACCACGGATGCTCCTCCATGCAGCAAccc                              <  2:3385260/101‑1 (MQ=255)

                                                                    ccAGACCAACCACTTTTTGCTCACCTTCGGCACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTACCACGGATGCTCCTCCATGCAGCAACCCg                             <  2:3385444/101‑1 (MQ=255)

                                                                    ccAGACCAACCACTTTTTGCTCACCTTCGGCACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTACCACGGATGCTCCTCCATGCAGCAACCCg                             <  1:3385139/101‑1 (MQ=255)

                                                                     cAGACCAACCACTTTTTGCTCACCTTCGGCACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTACCACGGATGCTCCTCCATGCAGCAACCCGa                            >  2:3385818/1‑101 (MQ=255)

                                                                       gACCAACCACTTTTTGCTCACCTTCGGCACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTACCACGGATGCTCCTCCATGCAGCAACCCGaaa                          >  1:3385822/1‑101 (MQ=255)

                                                                        accaaccaCTTTTTGCTCACCTTCGGCACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTACCACGGATGCTCCTCCATGCAGCAACCCGAAAg                         >  2:3385824/1‑101 (MQ=255)

                                                                        accaaccaCTTTTTGCTCACCTTCGGCACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTACCACGGATGCTCCTCCATGCAGCAACCCGAAAg                         <  1:3385396/101‑1 (MQ=255)

                                                                          caaccaCTTTTTGCTCACCTTCGGCACATAAATAACCGGCAACgg                                                                               >  1:3385828/1‑45 (MQ=255)

                                                                             ccaCTTTTTGCTCACCTTCGGCACATAAATAACCGGCAACgg                                                                               <  1:3385092/42‑1 (MQ=255)

                                                                             ccaCTTTTTGCTCACCTTCGGCACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTACCACGGATGCTCCTCCATGCAGCAACCCGAAAGgctgc                    >  2:3385831/1‑101 (MQ=255)

                                                                               aCTTTTTGCTCACCTTCGGCACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTACCACGGATGCTCCTCCATGCAGCAACCCGAAAGGCTGCtt                  <  1:3385457/101‑1 (MQ=255)

                                                                                cTTTTTGCTCACCTTCGGCACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTACCACGGATGCTCCTCCATGCAGCAACCCGa                            <  2:3385211/90‑1 (MQ=255)

                                                                                cTTTTTGCTCACCTTCGGCACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTACCACGGATGCTCCTCCATGCAGCAACCCGAAAGGCTGCttt                 <  2:3385618/101‑1 (MQ=255)

                                                                                cTTTTTGCTCACCTTCGGCACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTACCACGGATGCTCCTCCATGCAGCAACCCGAAAGGCTGCttt                 >  1:3385836/1‑101 (MQ=255)

                                                                                 tttttGCTCACCTTCGGCACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTACCACGGATGCTCCTCCATGCAGCAACCCGAAAGGCTGCTTTg                <  1:3385157/101‑1 (MQ=255)

                                                                                 tttttGCTCACCTTCGGCACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTACCACGGATGCTCCTCCATGCAGCAACCCGAAAGGCTGCTTTg                <  1:3385306/101‑1 (MQ=255)

                                                                                  ttttGCTCACCTTCGGCACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTACCACGGATGCTCCTCCATGCAGCAACCCGAAAGGCTGCTTTGt               >  2:3385839/1‑101 (MQ=255)

                                                                                   tttGCTCACCTTCGGCACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTACCACGGATGCTCCTCCATGCAGCAACCCGAAAGGCTGCTTTGTc              >  1:3385841/1‑101 (MQ=255)

                                                                                    ttGCTCACCTTCGGCACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTACCACGGATGCTCCTCCATGCAGCAACCCGAAAGGCTGCTTTGTcc             >  2:3385845/1‑101 (MQ=255)

                                                                                      gCTCACCTTCGGCACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTACCACGGATGCTCCTCCATGCAGCAACCCGAAAGGCTGCTTTGTCcgc           >  1:3385850/1‑101 (MQ=255)

                                                                                      gCTCACCTTCGGCACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTACCACGGATGCTCCTCCATGCAGCAACCCGAAAGGCTGCTTTGTCcgc           >  2:3385849/1‑101 (MQ=255)

                                                                                       cTCACCTTCGGCACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTACCACGGATGCTCCTCCATGCAGCAACCCGAAAGGCTGCTTTGTCcgcg          <  2:3385153/101‑1 (MQ=255)

                                                                                         cACCTTCGGCACATAAATAACCGGCAACGGAACCGATACTTTCGGCCa                                                             <  2:3385160/48‑1 (MQ=255)

                                                                                          aCCTTCGGCACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTACCACGGATGCTCCTCCATGCAGCAACCCGAAAGGCTGCTTTGTCCGCGAGt       >  1:3385854/1‑101 (MQ=255)

                                                                                            cTTCGGCACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTACCACGGATGCTCCTCCATGCAGCAACCCGAAAGGCTGCTTTGTCCGCGAGTCt     >  2:3385858/1‑101 (MQ=255)

                                                                                             ttCGGCACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTACCACGGATGCTCCTCCATGCAGCAACCCGAAAGGCTGCTTTGTCCGCGAGTCTa    >  2:3385861/1‑101 (MQ=255)

                                                                                               cGGCACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTACCACGGATGCTCCTCCATGCAGCAACCCGAAAGGCTGCTTTGTCCGCGAGTCTACt  <  1:3385116/101‑1 (MQ=255)

ATGCAACGGTTTGCATACGCCGCGCACCCGCCCTTCTCGTGCCGAGCGGACGTGGTTAGCATTGATTTCCAGACCAACCACTTTTTGCTCACCTTCGGTACATAAATAACCGGCAACGGAACCGATACTTTCGGCCAGTACCACGGATGCTCCTCCATGCAGCAACCCGAAAGGCTGCTTTGTCCGCGAGTCTACT  >  NC_000913/1765315‑1765510

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