BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,914,878	T→C	E42E (GAA→GAG)	atpD ←	membrane‑bound ATP synthase, F1 sector, beta‑subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,914,878	0	T	C	100.0%	690.0 / NA	220	E42E (GAA→GAG)	atpD	membrane‑bound ATP synthase, F1 sector, beta‑subunit
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (120/100); total (120/100)

Only 100 of 230 total aligned reads displayed.

GAGGTCTTTTACATCCAGACCGCGACGCAGACCGTCGGAGGAACCCATTGCGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACTTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCATCGTACACGCGCGGTACGGCATCCTGAGGGAATTCGACGTCAACTAC  >  NC_000913/3914788‑3914967

gAGGTCTTTTACATCCAGACCGCGACGCAGACCGTCGGAGGAACCCATTGCGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACCTc                                                                                         <  2:4379684/93‑1 (MQ=255)

 aGGTCTTTTACATCCAGACCGCGACGCAGACCGTCGGAGGAACCCATTGCGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACCTcc                                                                                        <  1:2601618/93‑1 (MQ=255)

 aGGTCTTTTACATCCAGACCGCGACGCAGACCGTCGGAGGAACCCATTGCGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACCTcc                                                                                        >  2:2340781/1‑93 (MQ=255)

  ggTCTTTTACATCCAGACCGCGACGCAGACCGTCGGAGGAACCCATTGCGATGGTACGTCCGATACCGCCGCCGAGCTGCTGCTGAACCTCCa                                                                                       <  2:2314198/93‑1 (MQ=255)

  ggTCTTTTACATCCAGACCGCGACGCAGACCGTCGGAGGAACCCATTGCGATGGTACGTACGATACCGCCGCCCAGCTGCTGCTGAACCTccg                                                                                       >  1:4263442/1‑92 (MQ=255)

    tCTTTTACATCCAGACCGCGACGCAGACCGTCGGAGGAACCCATTGCGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACCTCCAGc                                                                                     >  2:1817902/1‑93 (MQ=255)

     cTTTTACATCCAGACCGCGACGCAGACCGTCGGAGGAACCCATTGCGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACCTCCAGCa                                                                                    >  1:4066739/1‑93 (MQ=255)

     cTTTTACATCCAGACCGCGACGCAGACCGTCGGAGGAACCCATTGCGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACCTCCAGCa                                                                                    <  1:3774258/93‑1 (MQ=255)

        ttACATCCAGACCGCGACGCAGACCGTCGGAGGAACCCATTGCGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACCTCCAGCACCa                                                                                 <  1:2385512/93‑1 (MQ=255)

        ttACATCCAGACCGCGACGCAGACCGTCGGAGGAACCCATTGCGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACCTCCAGCACCa                                                                                 <  2:2725022/93‑1 (MQ=255)

        ttACATCCAGACCGCGACGCAGACCGTCGGAGGAACCCATTGCGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACCTCCAGCACCa                                                                                 >  2:2380359/1‑93 (MQ=255)

        ttACATCCAGACCGCGACGCAGACCGTCGGAGGAACCCATTGCGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACCTCCAGCACCa                                                                                 <  1:1214153/93‑1 (MQ=255)

        ttACATCCAGACCGCGACGCAGACCGTCGGAGGAACCCATTGCGATGGTACGTACGATAACGCAGCCGAGCTGCTGCTGAACCTCCAGCACCa                                                                                 >  1:2852557/1‑93 (MQ=255)

         tACATCCAGACCGCGACGCAGACCGTCGGAGGAACCCATTGCGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACCTCCAGCACCAg                                                                                >  1:2642683/1‑93 (MQ=255)

         tACATCCAGACCGCGACGCAGACCGTCGGAGGAACCCATTGCGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACCTCCAGCACCAg                                                                                >  1:1760332/1‑93 (MQ=255)

          aCATCCAGACCGCGACGCAGACCGTCGGAGGAACCCATTGCGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACCTCCAGCACCAGa                                                                               <  2:1329411/93‑1 (MQ=255)

           cATCCAGACCGCGACGCAGACCGTCGGAGGAACCCATTGCGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACCTCCAGCACCAGAc                                                                              >  1:1963457/1‑93 (MQ=255)

           cATCCAGACCGCGACGCAGACCGTCGGAGGAACCCATTGCGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACCTCCAGCACCAGAc                                                                              <  2:3449548/93‑1 (MQ=255)

            aTCCAGACCGCGACGCAGACCGTCGGAGGAACCCATTGCGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACCTCCAGCACCAGACg                                                                             >  1:2931362/1‑93 (MQ=255)

            aTCCAGACCGCGACGCAGACCGTCGGAGGAACCCATTGCGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACCTCCAGCACCAGACg                                                                             >  2:3153208/1‑93 (MQ=255)

             tCCAGACCGCGACGCAGACCGTCGGAGGAACCCATTGCGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACCTCCAGCACCAGACGc                                                                            <  1:3041617/93‑1 (MQ=255)

             tCCAGACCGCGACGCAGACCGTCGGAGGAACCCATTGCGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACCTCCAGCACCAGACGc                                                                            >  1:4338228/1‑93 (MQ=255)

               cAGACCGCGACGCAGACCGTCGGAGGAACCCATTGCGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACCTCCAGCACCAGACGCTc                                                                          >  1:2441548/1‑93 (MQ=255)

                aGACTGCGACGCAGACCGTCGGAGGAACCCATTGCGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTAAACCTCCAGCACCAGACGCTCa                                                                         >  1:502781/1‑93 (MQ=255)

                aGACCGCGACGCAGACCGTCGGAGGAACCCATTGCGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACCTCCAGCACCAGACGCTCa                                                                         <  2:1147312/93‑1 (MQ=255)

                aGACCGCGACGCAGACCGTCGGAGGAACCCATTGCGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACCTCCAGCACCAGACGCTCa                                                                         >  2:3587506/1‑93 (MQ=255)

                aGACCGCGACGCAGACCGTCGGAGGAACCCATTGCGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACCTCCAGCACCAGAAGCTCa                                                                         >  1:731680/1‑93 (MQ=255)

                 gACCGCGACGCAGACCGTCGGAGGAACCCATTGCGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACCTCCAGCACCAGACGCTCAt                                                                        >  2:1244232/1‑93 (MQ=255)

                  aCCGCGACGCAGACCGTCGGAGGAACCCATTGCGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACCTCCAGCACCAGACGCACAtt                                                                       <  2:335360/93‑1 (MQ=255)

                   ccGCGACGCAGACCGTCGGAGGAACCCATTGCGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACCTCCAGCACCAGACGCTCATTa                                                                      <  2:4118860/93‑1 (MQ=255)

                    cgcgACGCAGACCGTCGGAGGAACCCATTGCGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACCTCCAGCACCAGACGCTCATTAc                                                                     <  2:1846344/93‑1 (MQ=255)

                     gcgACGCAGACCGTCGGAGGAACCCATTGCGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACCTCCAGCACCAGACGCTCATTAcc                                                                    <  1:2803796/93‑1 (MQ=255)

                      cgACGCAGACCGTCGGAGGAACCCATTGCGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACCTCCAGCACCAGACGCTCATTACCa                                                                   <  1:229912/93‑1 (MQ=255)

                       gACGCAGACCGTCGGAGGAACCCATTGCGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACCTCCAGCACCAGACGCTCATTACCAt                                                                  <  1:918233/93‑1 (MQ=255)

                         cGCAGACCGGCGGAGGAACCCATTGCGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACCTCCAGCACCAGACGCTCATTACCAttt                                                                <  1:1350007/93‑1 (MQ=255)

                          gCAGACCGTCGGAGGAACCCATTGCGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACCTCCAGCACCAGACGCTCATTACCAtttt                                                               <  2:853889/93‑1 (MQ=255)

                          gCAGACCGTCGGAGGAACCCATTGCGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACCTCCAGCACCAGACGCTCATTACCAtttt                                                               >  1:4158365/1‑93 (MQ=255)

                          gCAGACCGTCGGAGGAACCCATTGCGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACCTCCAGCACCAGACGCTCATTACCAtttt                                                               <  1:3563887/93‑1 (MQ=255)

                           cAGACCGTCGGAGGAACCCATTGCGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACCTCCAGCACCAGACGCTCATTACCATTTTg                                                              >  2:4279616/1‑93 (MQ=255)

                           cAGACCGTCGGAGGAACCCATGGCGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAGCCTCCAGCACCAGACGCTCATTACCATTTTg                                                              <  1:942035/93‑1 (MQ=255)

                           cAGACCCTCGGAGGAACCCCTTGGGATGGTACGTACGAATCCGCCGCCGAGCTGCTTCTTAACCCCCAGCACCAGACGCTCATTACCATTTTg                                                              <  2:4178100/93‑1 (MQ=255)

                            aGACCGTCGGAGGAACCCATTGCGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACCTCCAGCACCAGACGCTCATTACCATTTTGc                                                             >  1:1549956/1‑93 (MQ=255)

                            aGACCGTCGGAGGAACCCACTGCGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACCTCCAGCACCAGACGCTCATTACCATTTTGc                                                             >  2:4562889/1‑93 (MQ=255)

                             gACCGTCGGAGGAACCCATTGCGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACCTCCAGCACCAGACGCTCATTACCATTTTGCa                                                            >  2:3866501/1‑93 (MQ=255)

                                  tCGGAGGAACCCATTGCGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACCTCCAGCACCAGACGCTCATTACCATTTTGCACCTCa                                                       >  1:1265128/1‑93 (MQ=255)

                                  tCGGAGGAACCCATTGCGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACCTCCAGCACCAGACGCCCATTACCATTTTGCACCTCa                                                       >  2:2641091/1‑93 (MQ=255)

                                   cGGAGGAACCCATTGCGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACCTCCAGCACCAGACGCTCATTACCATTTTGCACCTCaa                                                      <  2:3094189/93‑1 (MQ=255)

                                   cGGAGAAACCCATTGCGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACCTCCAGCACCAGACGCTCATTACCATTTTGCACCTCaa                                                      <  1:453773/93‑1 (MQ=255)

                                     gaggaACCCATTGCGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACCTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAaga                                                    <  1:1407797/93‑1 (MQ=255)

                                      aggaACCCATTGCGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACCTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAagag                                                   <  1:3703908/93‑1 (MQ=255)

                                      aggaACCCATTGCGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACCTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAagag                                                   <  2:2236929/93‑1 (MQ=255)

                                      aggaACCCATTGCGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACCTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAagag                                                   <  2:4251168/93‑1 (MQ=255)

                                      aggaACCCATTGCGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACCTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAagag                                                   <  1:3714223/93‑1 (MQ=255)

                                        gaACCCATTGCGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACCTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCa                                                 >  2:1189152/1‑93 (MQ=255)

                                        gaACCCATTGCGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACCTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCa                                                 <  1:1074060/93‑1 (MQ=255)

                                        gaACCCATTGCGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACCTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCa                                                 <  2:3892527/93‑1 (MQ=255)

                                          aCCCATTGCGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACCTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCATc                                               >  2:2747593/1‑93 (MQ=255)

                                          aCCCATTGCGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACCTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCATc                                               >  1:3510234/1‑93 (MQ=255)

                                          aCCCATTGCGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACCTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCATc                                               <  2:4205952/93‑1 (MQ=255)

                                            ccATTGCGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACCTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCATCGt                                             >  1:277879/1‑93 (MQ=255)

                                            ccATTGCGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACCTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCATCGt                                             >  2:3047638/1‑93 (MQ=255)

                                            ccATTGCGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACCTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCATCGt                                             >  2:2279579/1‑93 (MQ=255)

                                              aTTGCGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACCTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCATCGTac                                           >  2:1448532/1‑93 (MQ=255)

                                              aTTGCGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACCTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCATCGTac                                           <  2:2493614/93‑1 (MQ=255)

                                              aTTGCGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACCTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCATCGTac                                           >  1:1526867/1‑93 (MQ=255)

                                              aTTGCGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACCTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCATCGTac                                           <  1:2290153/93‑1 (MQ=255)

                                               ttGCGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACCTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCATCGTaca                                          <  2:365146/93‑1 (MQ=255)

                                                tGCGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACCTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCATCGTacac                                         <  1:4435267/93‑1 (MQ=255)

                                                  cGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACCTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCATCGTACAcgc                                       <  2:2574767/93‑1 (MQ=255)

                                                  cGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACCTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCATCGTACAcgc                                       >  2:1897466/1‑93 (MQ=255)

                                                   gATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACCTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCATCGTACAcgcg                                      >  1:3349130/1‑93 (MQ=255)

                                                     tGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACCTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCATCGTACAcgcgcg                                    >  1:1620625/1‑93 (MQ=255)

                                                       gtacgtacGATACCGCCGCCGAGCTGCTGCTGAACCTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCATCGTACACGCGCgga                                  >  2:2938835/1‑92 (MQ=255)

                                                            tacgATACCGCCGCCGAGCTGCTGCTGAACCTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCATCGTACACGCGCGGTACGGc                             <  2:461451/93‑1 (MQ=255)

                                                             acgATGCCGCCGCCGAGCTGCTGCGGAACCTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCATCGTACACGCGCGGTACGGCa                            <  2:3821729/93‑1 (MQ=255)

                                                             acgATACCGCCGCCGAGCTGCTGCTGAACCTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCATCGTACACGCGCGGTACGGCa                            >  1:100876/1‑93 (MQ=255)

                                                              cgATACGGCCGCCGAGCTGCTGCTGAACCTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCATCGGACACGCGCGGTACGGCAt                           >  2:2139287/1‑93 (MQ=255)

                                                              cgATACCGCCGCCGAGCTGCTGCTGAACCTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCATCGTACACGCGCGGTACGGCAt                           <  1:3805413/93‑1 (MQ=255)

                                                                aTACCGCCGCCGAGCTGCTGCTGAACCTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCATCGTACACGCGCGGTACGGCATcc                         >  1:4462863/1‑93 (MQ=255)

                                                                aTACCGCCGCCGAGCTGCTGCTGAACCTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCATCGTACACGCGCGGTACGGCATcc                         >  2:4074832/1‑93 (MQ=255)

                                                                 tACCGCCGCCGAGCTGCTGCTGAACCTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCATCGTACACGCGCGGTACGGCATCCt                        >  1:620631/1‑93 (MQ=255)

                                                                 tACCGCCGCCGAGCTGCTGCTGAACCTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCATCGTACACGCGCGGTACGGCATCCt                        >  1:2551805/1‑93 (MQ=255)

                                                                   ccgccgccgAGCTGCTGCTGAACCTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCATCGTACACGCGCGGTACGGCATCCTGa                      <  1:2066752/93‑1 (MQ=255)

                                                                    cgccgccgAGCTGCTGCTGAACCTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCATCGTACACGCGCGGTACGGCATCCTGAg                     >  2:600391/1‑93 (MQ=255)

                                                                        gccgAGCTGCTGCTGAACCTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCATCGTACACGCGCGGTACGGCATCCTGAGGGaa                 <  1:3654368/93‑1 (MQ=255)

                                                                          cgAGCTGCTCCTGAACCTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCATCGTACACGCGCGATACGGCATCCTGAGGGAAtt               <  2:3285142/93‑1 (MQ=255)

                                                                           gAGCTGCTGCTGAACCTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCATCGTACACGCGCGGGACGGCATCCTGAGGGAATTc              >  2:513399/1‑93 (MQ=255)

                                                                             gctgctgctGACCCTCAAGCATCAGACGCTCATTACCATTTTGCCCCTCAAGAGCATCGTACACGCGCTGTACGGCATCCTGAGGGAATTCGa            <  2:1641834/93‑1 (MQ=255)

                                                                             gctgctgctGAACCTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCATCGTACACGCGCGGTACGGCATCCTGAGGGAATTCGa            <  2:815845/93‑1 (MQ=255)

                                                                              ctgctgctgAACCTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCATCGTACACGCGCGGTACGGCATCCTGAGGGAATTCGAc           <  1:4622783/93‑1 (MQ=255)

                                                                               tgctgctgAACCTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCATCGTACACGCGCGGTACGGCATCCTGAGGGAATTCGACg          >  2:3086012/1‑93 (MQ=255)

                                                                                  tgctgAACCTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCATCGTACACGCGCGGTACGGCATCCTGAGGGAATTCGACGTCa       <  1:2341111/93‑1 (MQ=255)

                                                                                   tctgAACCTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCATCGTACACGCGCGGTACGGCATCCTGAGGGAATTCGACGTCaa      <  2:1392984/92‑1 (MQ=255)

                                                                                   gctgAACCTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCATCGTACACGCGCGGTACGGCATCCTGAGGGAATTCGACGTCaa      >  2:349178/1‑93 (MQ=255)

                                                                                    ctgAACCTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCATCGTACACGCGCGGTACGGCATCCTGAGGGAATTCGACGTCAAc     <  1:983049/93‑1 (MQ=255)

                                                                                    ctgAACCTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCATCGTACACGCGCGGTACGGCATCCTGAGGGAATTCGACGTCAAc     <  1:2891423/93‑1 (MQ=255)

                                                                                    ctgAACCTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCATCGTACACGCGCGGTACGGCATCCTGAGGGAATTCGACGTCAAc     <  2:1990502/93‑1 (MQ=255)

                                                                                     tgAACCTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCATCGTACACGCGCGGTACGGCATCCTGAGGGAATTCGACGTCAACt    >  1:4389915/1‑93 (MQ=255)

                                                                                       aaCCTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCATCGTACACGCGCGGTACGGCATCCTGAGGGAATTCGACGTCAACTAc  <  1:316273/93‑1 (MQ=255)

                                                                                       aaCCTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCATCGTACACGCGCGGTACGGCATCCTGAGGGAATTCGACGTCAACTAc  <  2:2887894/93‑1 (MQ=255)

GAGGTCTTTTACATCCAGACCGCGACGCAGACCGTCGGAGGAACCCATTGCGATGGTACGTACGATACCGCCGCCGAGCTGCTGCTGAACTTCCAGCACCAGACGCTCATTACCATTTTGCACCTCAAGAGCATCGTACACGCGCGGTACGGCATCCTGAGGGAATTCGACGTCAACTAC  >  NC_000913/3914788‑3914967

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