BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,828,320	T→C	D472G (GAC→GGC)	astD ←	aldehyde dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,828,320	0	T	C	100.0%	1119.3 / NA	279	D472G (GAC→GGC)	astD	aldehyde dehydrogenase
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (136/143); total (136/143)

Only 100 of 283 total aligned reads displayed.

CCTACCAGCCCGTCGAAATTGACTTCCCAGGCGTTCATCGCACCACCTCATCGGAAAAATCCAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGTCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCGGCATACCAGGCGCTGGGGCGATGGTTACCGGAAGCACCAATGCCGCCGAATGGCGC  >  NC_000913/1828222‑1828417

 cTACCAGCCCGTCGAAATTGACTTCCCAGGCGTTCATCGCACCACCTCATCGGAAAAATCCAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGCCCg                                                                                                >  2:3063821/1‑101 (MQ=255)

   aCCAGCCCGTCGAAATTGACTTCCCAGGCGTTCATCGCACCACCTCATCGGAAAAATCCAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGCCCGAc                                                                                              >  1:3063824/1‑101 (MQ=255)

    ccAGCCCGTCGAAATTGACTTCCCAGGCGTTCATCGCACCACCTCATCGGAAAAATCCAGCCCGGGGTTAAGCGTGGCGGGCAATGTTACCGAGCCCGACt                                                                                             <  1:3063363/101‑1 (MQ=255)

    ccAGCCCGTCGAAATTGACTTCCCAGGCGTTCATCGCACCACCTCATCGGAAAAATCCAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGCCCGACt                                                                                             >  1:3063827/1‑101 (MQ=255)

     cAGCCCGTCGAAATTGACTTCCCAGGCGTTCATCGCACCACCTCATCGGAAAAATCCAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGCCCGACTc                                                                                            >  2:3063831/1‑101 (MQ=255)

     cAGCCCGTCGAAATTGACTTCCCAGGCGTTCATCGCACCACCTCATCGGAAAAATCCAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGCCCGACTc                                                                                            >  2:3063829/1‑101 (MQ=255)

      aGCCCGTCGAAATTGACTTCCCAGGCGTTCATCGCACCACCTCATCGGAAAAATCCAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGCCCGACTcc                                                                                           <  2:3063255/101‑1 (MQ=255)

       gCCCGTCGAAATTGACTTCCCAGGCGTTCATCGCACCACCTCATCGGAAAAATCCAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGCCCGACTCCa                                                                                          >  2:3063835/1‑101 (MQ=255)

        cccGTCGAAATTGACTTCCCAGGCGTTCATCGCACCACCTCATCGGAAAAATCCAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGCCCGACTCCAg                                                                                         <  2:3063019/101‑1 (MQ=255)

           gTCGAAATTGACTTCCCAGGCGTTCATCGCACCACCTCATCGGAAAAATCCAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGCCCGAc                                                                                              <  2:3063052/93‑1 (MQ=255)

            tCGAAATTGACTTCCCAGGCGTTCATCGCACCACCTCATCGGAAAAATCCAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGCCCGACTCCAGGCTc                                                                                     <  2:3063202/101‑1 (MQ=255)

            tCGAAATTGACTTCCCAGGCGTTCATCGCACCACCTCATCGGAAAAATCCAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGCCCGACTCCAGGCTc                                                                                     <  1:3063207/101‑1 (MQ=255)

             cGAAATTGACTTCCCAGGCGTTCATCGCACCACCTCATCGGAAAAATCCAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGCCCGACTCCAGGCTCg                                                                                    >  2:3063840/1‑101 (MQ=255)

              gAAATTGACTTCCCAGGCGTTCATCGCACCACCTCATCGGAAAAATCCAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGCCCGACTCCAGGCTCGc                                                                                   >  1:3063842/1‑101 (MQ=255)

               aaaTTGACTTCCCAGGCGTTCATCGCACCACCTCATCGGAAAAATCCAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGCCCGACTCCAGGCTCGcc                                                                                  >  2:3063846/1‑101 (MQ=255)

                aaTTGACTTCCCAGGCGTTCATCGCACCACCTCATCGGAAAAATCCAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGCCCGACTCCAGGCTCGCCa                                                                                 <  2:3063359/101‑1 (MQ=255)

                  ttGACTTCCCAGGCGTTCATCGCACCACCTCATCGGAAAAATCCAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGCCCGACTCCAGGCTCGCCATc                                                                               <  1:3063506/101‑1 (MQ=255)

                  ttGACTTCCCAGGCGTTCATCGCACCACCTCATCGGAAAAATCCAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGCCCGACTCCAGGCTCGCCATc                                                                               <  1:3063293/101‑1 (MQ=255)

                  ttGACTTCCCAGGCGTTCATCGCACCAACTCATCGGAAAAATCCAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGCCCGACTCCAGGCTCGCCATc                                                                               <  1:3063237/101‑1 (MQ=255)

                   tGACTTCCCAGGCGTTCATCGCACCACCTCATCGGAAAAATCCAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGCCCGACTCCAGGCTCGCCATCg                                                                              >  1:3063851/1‑101 (MQ=255)

                    gACTTCCCAGGCGTTCATCGCACCACCTCATCGGAAAAATCCAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGCCCGACTCCAGGCTc                                                                                     >  2:3063853/1‑93 (MQ=255)

                    gACTTCCCAGGCGTTCATCGCACCACCTCATCGGAAAAATCCAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGCCCGACTCCAGGCTCGCCATCgg                                                                             <  1:3063013/101‑1 (MQ=255)

                    gACTTCCCAGGCGTTCATCGCACCACCTCATCGGAAAAATCCAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGCCCGACTCCAGGCTCGCCATCgg                                                                             <  2:3063501/101‑1 (MQ=255)

                      cTTCCCAGGCGTTCATCGCACCACCTCATCGGAAAAATCCAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGCCCGACTCCAGGCTCGCCATCGGcc                                                                           <  2:3063398/101‑1 (MQ=255)

                       ttCCCAGGCGTTCATCGCACCACCTCATCGGAAAAATCCAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGCCCGACTCCAGGCTCGCCa                                                                                 >  1:3063857/1‑94 (MQ=255)

                       ttCCCAGGCGTTCATCGCACCACCTCATCGGAAAAATCCAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGCCCGACTCCAGGCTCGCCATCGGCCa                                                                          <  2:3063309/101‑1 (MQ=255)

                        tCCCAGGCGTTCATCGCACCACCTCATCGGAAAAATCCAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGCCCGACTCCAGGCTCGCCATCGGCCAt                                                                         <  2:3063651/101‑1 (MQ=255)

                           cAGGCGTTCATCGCACCACCTCATCGGAAAAATCCAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGCCCGACTCCAGGCTCGCCATCGGCCATgcg                                                                      >  1:3063862/1‑101 (MQ=255)

                            aGGCGTTCATCGCACCACCTCATCGGAAAAATCCAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGCCCGACTCCAGGCTCGCCATCGGcc                                                                           >  2:3063868/1‑95 (MQ=255)

                            aGGCGTTCATCGCACCACCTCATCGGAAAAATCCAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGCCCGACTCCAGGCTCGCCATCGGCCATgcgc                                                                     >  1:3063867/1‑101 (MQ=255)

                             ggCGTTCATCGCACCACCTCATCGGAAAAATCCAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGCCCGACTCCAGGCTCGcc                                                                                  >  2:3063870/1‑87 (MQ=255)

                              gCGTTCATCGCACCACCTCATCGGAAAAATCCAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGCCCGACTCCAGGCTCGCCATCGGCCATGCGCAg                                                                   <  1:3063308/101‑1 (MQ=255)

                               cGTTCATCGCACCACCTCATCGGAAAAATCCAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGCCCGACTCCAGGCTCGCCATCGGCCATgc                                                                       >  2:3063873/1‑96 (MQ=255)

                                gTTCATCGCACCACCTCATCGGAAAAATCCAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGCCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTa                                                                 <  2:3063323/101‑1 (MQ=255)

                                gTTCATCGCACCACCTCATCGGAAAAATCCAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGCCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTa                                                                 >  2:3063878/1‑101 (MQ=255)

                                 ttCATCGCACCACCTCATCGGAAAAATCCAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGCCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTaa                                                                >  1:3063883/1‑101 (MQ=255)

                                  tCATCGCACCACCTCATCGGAAAAATCCAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGCCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAAt                                                               >  2:3063887/1‑101 (MQ=255)

                                   cATCGCACCACCTCATCGGAAAAATCCAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGCCCGa                                                                                               >  1:3063889/1‑68 (MQ=255)

                                   cATCGCACCACCTCATCGGAAAAATCCAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGCCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATc                                                              <  2:3063433/101‑1 (MQ=255)

                                    aTCGCACCACCTCATCGGAAAAATCCAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGCCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCt                                                             <  1:3063583/101‑1 (MQ=255)

                                       gCACCACCTCATCGGAAAAATCCAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGCCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCg                                                          >  2:3063899/1‑101 (MQ=255)

                                       gCACCACCTCATCGGAAAAATCCAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGCCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCg                                                          >  1:3063896/1‑101 (MQ=255)

                                       gCACCACCTCATCGGAAAAATCCAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGCCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCg                                                          <  2:3063578/101‑1 (MQ=255)

                                        caccacCTCATCGGAAAAATCCAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGCCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCgg                                                         >  1:3063901/1‑101 (MQ=255)

                                          ccaccTCATCGGAAAAATCCAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGCCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCGGCa                                                       >  2:3063907/1‑101 (MQ=255)

                                           caccTCATCGGAAAAATCCAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGCCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAAt                                                               >  2:3063909/1‑92 (MQ=255)

                                             ccTCATCGGAAAAATCCAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGCCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATc                                                              <  2:3063518/91‑1 (MQ=255)

                                               tCATCGGAAAAATCCAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGCCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCGGCATACCa                                                  <  2:3063489/101‑1 (MQ=255)

                                               tCATCGGAAAAATCCAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGCCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCGGCATACCa                                                  >  1:3063912/1‑101 (MQ=255)

                                                 aTCGGAAAAATCCAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGCCCGACTCCAGGCTCGcc                                                                                  >  1:3063917/1‑67 (MQ=255)

                                                 aTCGGAAAAATCCAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGCCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCGGc                                                        <  1:3063660/93‑1 (MQ=255)

                                                 aTCGGAAAAATCCAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGCCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCGGCATACCAgg                                                <  2:3063412/101‑1 (MQ=255)

                                                 aTCGGAAAAATCCAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGCCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCGGCATACCAgg                                                <  2:3063540/101‑1 (MQ=255)

                                                 aTCGGAAAAATCCAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGCCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCGGCATACCAgg                                                >  2:3063915/1‑101 (MQ=255)

                                                 aTCGGAAAAATCCAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGCCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCGGCATACCAgg                                                <  2:3063368/101‑1 (MQ=255)

                                                  tCGGAAAAATCCAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGCCCGACTCCAGGCTCGCCATCGGTCATGCGCAGTAATCTGCGGCATACCAGgc                                               <  1:3063614/101‑1 (MQ=255)

                                                   cGGAAAAATCCAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGCCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCGGCATACCAGgcg                                              <  2:3063451/101‑1 (MQ=255)

                                                    ggAAAAATCCAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGCCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCGGCATACCa                                                  <  1:3063563/96‑1 (MQ=255)

                                                     gAAAAATCCAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGCCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCGGCATACCAGGCGCt                                            >  1:3063922/1‑101 (MQ=255)

                                                     gAAAAATCCAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGCCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCGGCATACCAGGCGCt                                            >  2:3063923/1‑101 (MQ=255)

                                                       aaaaTCCAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGCCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCGGCATACCAGGCGCTgg                                          <  2:3063473/101‑1 (MQ=255)

                                                       aaaaTCCAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGCCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCGGCATACCAGGCGCTgg                                          >  1:3063926/1‑101 (MQ=255)

                                                        aaaTCCAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGCCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCGGCATACCAGGCGCTggg                                         <  1:3063188/101‑1 (MQ=255)

                                                         aaTCCAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGCCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCGGCATACCAGGCGCTgggg                                        >  1:3063931/1‑101 (MQ=255)

                                                          aTCCAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGCCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCGGCATACCAGGCGCTGGGGc                                       <  1:3063523/101‑1 (MQ=255)

                                                          aTCCAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGCCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCGGCATACCAGGCGCTGGGGc                                       >  2:3063933/1‑101 (MQ=255)

                                                             cAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGCCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCGGCATACCAGGCGCTGGGGCGAt                                    <  1:3063555/101‑1 (MQ=255)

                                                             cAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGCCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCGGCATACCAGGCGCTGGGGCGAt                                    <  2:3063441/101‑1 (MQ=255)

                                                             cAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGCCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCGGCATACCAGGCGCTGGGGCGAt                                    >  2:3063937/1‑101 (MQ=255)

                                                              aGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGCCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCGGCATACCAGGCGCTGGGGCGATg                                   <  1:3063160/101‑1 (MQ=255)

                                                                cccGGGGTTAAGCGTGGCGGGCAATGTTAACGAGCCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCGGCATACCAGGCGCTGGGGCGATGGt                                 >  2:3063941/1‑101 (MQ=255)

                                                                 ccGGGGTTAAGCGTGGCGGGCAATGTTAACGAGCCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCGGCATACCAGGCGCTGGGGCGATGGtt                                >  1:3063943/1‑101 (MQ=255)

                                                                   ggggTTAAGCGTGGCGGGCAATGTTAACGAGCCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCGGCATACCAGGCGCTGGGGCGATGGTTAc                              >  1:3063946/1‑101 (MQ=255)

                                                                    gggTTAAGCGTGGCGGGCAATGTTAACGAGCCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCGGCATACCAGGCGCTGGGGCGATGGTTAcc                             <  2:3063461/101‑1 (MQ=255)

                                                                    gggTTAAGCGTGGCGGGCAATGTTAACGAGCCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCGGCATACCAGGCGCTGGGGCGATGGTTAcc                             <  1:3063457/101‑1 (MQ=255)

                                                                      gTTAAGCGTGGCGGGCAATGTTAACGAGCCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCGGCATACCAGGCGCTGGGGCGATGGTTACCgg                           >  1:3063949/1‑101 (MQ=255)

                                                                      gTTAAGCGTGGCGGGCAATGTTAACGAGCCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCGGCATACCAGGCGCTGGGGCGATGGTTACCgg                           <  1:3063406/101‑1 (MQ=255)

                                                                      gTTAAGCGTGGCGGGCAATGTTAACGAGCCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCGGCATACCAGGCGCTGGGGCGATGGTTACCgg                           <  2:3063556/101‑1 (MQ=255)

                                                                      gTTAAGCGTGGCGGGCAATGTTAACGAGCCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCGGCATACCAGGCGCTGGGGCGATGGTTACCgg                           <  2:3063736/101‑1 (MQ=255)

                                                                       ttAAGCGTGGCGGGCAATGTTAACGAGCCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCGGCATACCAGGCGCTGGGGCGATgg                                  >  2:3063950/1‑93 (MQ=255)

                                                                       ttAAGCGTGGCGGGCAATGTTAACGAGCCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCGGCATACCAGGCGCTGGGGCGATGGTTACCGGa                          <  2:3063394/101‑1 (MQ=255)

                                                                        tAAGCGTGGCGGGCAATGTTAACGAGCCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCGGCATACCAGGCGCTGGGGCGATGGTTACCGGaa                         <  1:3063538/101‑1 (MQ=255)

                                                                        tAAGCGTGGCGGGCAATGTTAACGAGCCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCGGCATACCAGGCGCTGGGGCGATGGTTACCGGaa                         <  2:3063382/101‑1 (MQ=255)

                                                                          aGCGTGGCGGGCAATGTTAACGAGCCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCGGCATACCAGGCGCTGGGGCGATGGTTACCGGAAGc                       <  2:3063268/101‑1 (MQ=255)

                                                                            cGTGGCGGGCAATGTTAACGAGCCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCGGCATACCAGGCGCt                                            <  2:3063339/78‑1 (MQ=255)

                                                                            cGTGGCGGGCAATGTTAACGAGCCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCGGCATACCAGGCGCTgggg                                        >  2:3063954/1‑82 (MQ=255)

                                                                             gTGGCGGGCAATGTTAACGAGCCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCGGCATACCAGGCGCTGGGGCGATGGTTAcc                             >  1:3063955/1‑92 (MQ=255)

                                                                              tGGCGGGCAATGTTAACGAGCCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCg                                                          >  2:3063958/1‑62 (MQ=255)

                                                                                 cGGGCAATGTTAACGAGCCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCGGCATACCAGGCGCTGGGGCGATGGTTACCGGAAGCACCAATg                <  1:3063466/101‑1 (MQ=255)

                                                                                    gCAATGTTAACGAGCCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCGGCATACCAGGCGCTGGGGCGATGGTTACCGGAAGCACCAATgccg             <  2:3063625/101‑1 (MQ=255)

                                                                                     cAATGTTAACGAGCCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCGGCATACCAGGCGCTGGGGCGATGGTTACCGGAAGCACCAATgccgc            >  1:3063963/1‑101 (MQ=255)

                                                                                      aaTGTTAACGAGCCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCGGCATACCAGGCGCTGGGGCGATGGTTACCGGAAGCACCAATgccgcc           <  1:3063481/101‑1 (MQ=255)

                                                                                      aaTGTTAACGAGCCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCGGCATACCAGGCGCTGGGGCGATGGTTACCGGAAGCACCAATgccgcc           <  2:3063111/101‑1 (MQ=255)

                                                                                       aTGTTAACGAGCCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCGGCATACCAGGCGCTGGGGCGATGGTTACCGGAAGCACCAATGccgccg          >  1:3063966/1‑101 (MQ=255)

                                                                                        tGTTAACGAGCCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCGGCATACCAGGCGCTGGGGCGATgg                                  >  2:3063968/1‑76 (MQ=255)

                                                                                         gTTAACGAGCCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCGGCATACCAGGCGCTGGGGCGATGGTTACCGGAAGCACCAATGCCGCCGaa        <  1:3063415/101‑1 (MQ=255)

                                                                                           tAACGAGCCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCGGCATACCAGGCGCTGGGGCGATGGTTACCGGAAGCACCAATGCCGCCGAATg      <  1:3063335/101‑1 (MQ=255)

                                                                                           tAACGAGCCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCGGCATACCAGGCGCTGGGGCGATGGTTACCGGAAGCACCAATGCCGCCGAATg      >  1:3063972/1‑101 (MQ=255)

                                                                                             aCGAGCCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCGGCATAc                                                    <  1:3063431/53‑1 (MQ=255)

                                                                                              cGAGCCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCGGCATACCAGGCGCTGGGGCGATGGTTACCGGAAGCACCAATGCCGCCGAATGgcg   >  2:3063974/1‑101 (MQ=255)

CCTACCAGCCCGTCGAAATTGACTTCCCAGGCGTTCATCGCACCACCTCATCGGAAAAATCCAGCCCGGGGTTAAGCGTGGCGGGCAATGTTAACGAGTCCGACTCCAGGCTCGCCATCGGCCATGCGCAGTAATCTGCGGCATACCAGGCGCTGGGGCGATGGTTACCGGAAGCACCAATGCCGCCGAATGGCGC  >  NC_000913/1828222‑1828417

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