BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,350,598	C→T	M138I (ATG→ATA)	glpQ ←	glycerophosphoryl diester phosphodiesterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,350,598	0	C	T	100.0%	1089.1 / NA	261	M138I (ATG→ATA)	glpQ	glycerophosphoryl diester phosphodiesterase
Reads supporting (aligned to +/- strand): ref base C (0/0); new base T (130/131); total (130/131)

Only 100 of 261 total aligned reads displayed.

TCTGGATAAATACCGATATTTTTCCCGGTAGAGTGATTTAACCCCTGAACAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCCATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTCAACGACTTAATTTCATCCAGCGTGAAAT  >  NC_000913/2350500‑2350695

 cTGGATAAATACCGATATTTTTCCCGGTAGAGTGATTTAACCCCTGAACAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCTAtt                                                                                                >  1:3048603/1‑101 (MQ=255)

  tGGATAAATACCGATATTTTTCCCGGTAGAGTGATTTAACCCCTGAACAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCTATTg                                                                                               >  2:3048605/1‑101 (MQ=255)

  tGGATAAATACCGATATTTTTCCCGGTAGAGTGATTTAACCCCTGAACAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCTATTg                                                                                               <  1:3048222/101‑1 (MQ=255)

    gATAAATACCGATATTTTTCCCGGTAGAGTGATTTAACCCCTGAACAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCTATTggg                                                                                             >  1:3048608/1‑101 (MQ=255)

     aTAAATACCGATATTTTTCCCGGTAGAGTGATTTAACCCCTGAACAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCTATTGGGa                                                                                            >  2:3048610/1‑101 (MQ=255)

       aaaTACCGATATTTTTCCCGGTAGAGTGATTTAACCCCTGAACAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCTATTGGGaaa                                                                                          <  2:3048079/101‑1 (MQ=255)

       aaaTACCGATATTTTTCCCGGTAGAGTGATTTAACCCCTGAACAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCTATTGGGaaa                                                                                          <  1:3048120/101‑1 (MQ=255)

        aaTACCGATATTTTTCCCGGTAGAGTGATTTAACCCCTGAACAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCTATTGGGAAAc                                                                                         >  1:3048614/1‑101 (MQ=255)

           aCCGATATTTTTCCCGGTAGAGTGATTTAACCCCTGAACAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCTATTGGGAAACGcc                                                                                      >  1:3048617/1‑101 (MQ=255)

            ccGATATTTTTCCCGGTAGAGTGATTTAACCCCTGAACAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCTATTGGGAAACGccc                                                                                     <  2:3048201/101‑1 (MQ=255)

             cGATATTTTTCCCGGTAGAGTGATTTAACCCCTGAACAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCTATTGGGAAACGcccc                                                                                    <  1:3048184/101‑1 (MQ=255)

              gATATTTTTCCCGGTAGAGTGATTTAACCCCTGAACAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCTATTGGGAAACGCCCCg                                                                                   >  2:3048621/1‑101 (MQ=255)

              gATATTTTTCCCGGTAGAGTGATTTAACCCCTGAACAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCTATTGGGAAACGCCCCg                                                                                   <  2:3048424/101‑1 (MQ=255)

                 atTTTTCCCGGTAGAGTGATTTAACCCCTGAACAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCTATTGGGAAACGCCCCGGAt                                                                                <  2:3047999/101‑1 (MQ=255)

                 atTTTTCCCGGTAGAGTGATTTAACCCCTGAACAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCTATTGGGAAACGCCCCGGAt                                                                                >  1:3048623/1‑101 (MQ=255)

                 atTTTTCCCGGTAGAGTGATTTAACCCCTGAACAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCTATTGGGAAACGCCCCGGAt                                                                                >  1:3048626/1‑101 (MQ=255)

                  tttttCCCGGTAGAGTGATTTAACCCCTGAACAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCTATTGGGAAACGCCCCGGATa                                                                               <  1:3048202/101‑1 (MQ=255)

                   ttttCCCGGTAGAGTGATTTAACCCCTGAACAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCTATTGGGAAACGCCCCGGATaa                                                                              <  2:3048254/101‑1 (MQ=255)

                   ttttCCCGGTAGAGTGATTTAACCCCTGAACAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCTATTGGGAAACGCCCCGGATaa                                                                              <  2:3048575/101‑1 (MQ=255)

                   ttttCCCGGTAGAGTGATTTAACCCCTGAACAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCTATTGGGAAACGCCCCGGATaa                                                                              >  2:3048628/1‑101 (MQ=255)

                   ttttCCCGGTAGAGTGATTTAACCCCTGAACAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCTATTGGGAAACGCCCCGGATaa                                                                              <  1:3048292/101‑1 (MQ=255)

                     ttCCCGGTAGAGTGATTTAACCCCTGAACAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCTATTGGGAAACGCCCCGGATAAGt                                                                            <  1:3048059/101‑1 (MQ=255)

                     ttCCCGGTAGAGTGATTTAACCCCTGAACAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCTATTGGGAAACGCCCCGGATAAGt                                                                            >  1:3048631/1‑101 (MQ=255)

                       cccGGTAGAGTGATTTAACCCCTGAACAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCTATTggg                                                                                             <  1:3048081/82‑1 (MQ=255)

                       cccGGTAGAGTGATTTAACCCCTGAACAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCTATTGGGAAACGCCCCGGATAAGTCt                                                                          >  2:3048633/1‑101 (MQ=255)

                       cccGGTAGAGTGATTTAACCCCTGAACAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCTATTGGGAAACGCCCCGGATAAGTCt                                                                          >  1:3048636/1‑101 (MQ=255)

                        ccGGTAGAGTGATTTAACCCCTGAACAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCTATTGGGAAACGCCCCGGATAAGTCTg                                                                         <  2:3048297/101‑1 (MQ=255)

                         cGGTAGAGTGATTTAACCCCTGAACAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCTATTg                                                                                               >  2:3048638/1‑78 (MQ=255)

                             agagTGATTTAACCCCTGAACAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCTATTGGGAAACGCCCCGGATAAGTCTGCACtt                                                                    <  1:3048255/101‑1 (MQ=255)

                              gagTGATTTAACCCCTGAACAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCTATTGGGAAACGCCCCGGATAAGTCTGCACttt                                                                   <  1:3048447/101‑1 (MQ=255)

                              gagTGATTTAACCCCTGAACAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCTATTGGGAAACGCCCCGGATAAGTCTGCACttt                                                                   >  1:3048641/1‑101 (MQ=255)

                               agTGATTTAACCCCTGAACAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCTATTGGGAAACGCCCCGGATAAGTCTGCACtttt                                                                  <  1:3048068/101‑1 (MQ=255)

                                gTGATTTAACCCCTGAACAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCTATTGGGAAACGCCCCGGATAAGTCTGCACttttt                                                                 <  1:3048279/101‑1 (MQ=255)

                                 tGATTTAACCCCTGAACAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCTATTGGGAAACGCCCCGGATAAGTCTGCACtttttt                                                                >  2:3048644/1‑101 (MQ=255)

                                   aTTTAACCCCTGAACAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCTATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTAc                                                              >  2:3048646/1‑101 (MQ=255)

                                        aCCCCTGAACAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCTATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTACCGttt                                                         >  1:3048649/1‑101 (MQ=255)

                                        aCCCCTGAACAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCTATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTACCGttt                                                         <  2:3048106/101‑1 (MQ=255)

                                        aCCCCTGAACAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCTATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTACCGttt                                                         <  2:3048175/101‑1 (MQ=255)

                                         ccccTGAACAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCTATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTACCGtttt                                                        >  1:3048651/1‑101 (MQ=255)

                                          cccTGAACAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCTATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTACCGTTTTc                                                       >  2:3048653/1‑101 (MQ=255)

                                           ccTGAACAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCTATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTACCGTTTTCa                                                      >  1:3048658/1‑101 (MQ=255)

                                           ccTGAACAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCTATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTACCGTTTTCa                                                      >  2:3048656/1‑101 (MQ=255)

                                               aaCAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCTATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTACCGTTTTCAatat                                                  >  2:3048662/1‑101 (MQ=255)

                                               aaCAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCTATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTACCGTTTTCAatat                                                  <  1:3048242/101‑1 (MQ=255)

                                                aCAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCTATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTACCGTTTTCAATATc                                                 <  2:3048123/101‑1 (MQ=255)

                                                 cAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCTATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTACCGTTTTCAATATCg                                                <  2:3048383/101‑1 (MQ=255)

                                                 cAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCTATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTACCGTTTTCAATATCg                                                <  2:3048265/101‑1 (MQ=255)

                                                   aaTTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCTATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTACCGTTTTCAATATCGaa                                              >  2:3048666/1‑101 (MQ=255)

                                                    aTTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCTATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTACCGTTTTCAATATCGaaa                                             <  2:3048451/101‑1 (MQ=255)

                                                    aTTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCTATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTACCGTTTTCAATATCGaaa                                             <  2:3048405/101‑1 (MQ=255)

                                                    aTTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCTATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTACCGTTTTCAATATCGaaa                                             >  2:3048668/1‑101 (MQ=255)

                                                    aTTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCTATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTACCGTTTTCAATATCGaaa                                             >  1:3048669/1‑101 (MQ=255)

                                                     ttCAATCTCTTCTTCAAAGGTGAGCACCCGGAAGTCGGACTTACCTAtt                                                                                                >  2:3048672/1‑49 (MQ=38)

                                                       cAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCTATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTACCGTTTTCAATATCGAAACCt                                          <  1:3048211/101‑1 (MQ=255)

                                                       cAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCTATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTACCGTTTTCAATATCGAAACCt                                          >  1:3048674/1‑101 (MQ=255)

                                                        aaTCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCTATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTACCGTTTTCAATATCGAAACCtt                                         <  1:3048238/101‑1 (MQ=255)

                                                        aaTCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCTATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTACCGTTTTCAATATCGAAACCt                                          <  2:3048221/100‑1 (MQ=255)

                                                         aTCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCTATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTACCGTTTTCAATATCGAAACCTTc                                        >  1:3048677/1‑101 (MQ=255)

                                                           ctctTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCTATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTACCGTTTTCAATATCGAAACCTTCgg                                      >  2:3048679/1‑101 (MQ=255)

                                                            tcttctTCAAAGGTGTGCACCCGGAAGTCGGACTTACCTATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTACCGTTTTCAATATCGAAACCTTCGGt                                     <  1:3048261/101‑1 (MQ=255)

                                                            tcttctTCAAAGGTGTGCACCCGGAAGTCGGACTTACCTATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTACCGTTTTCAATATCGAAACCTTCGGt                                     <  1:3048296/101‑1 (MQ=255)

                                                             cttcttCAAAGGTGTGCACCCGGAAGTCGGACTTACCTATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTACCGTTTTCAATATCGAAACCTTCGGTa                                    >  2:3048682/1‑101 (MQ=255)

                                                              ttcttcAAAGGTGTGCACCCGGAAGTCGGACTTACCTATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTACCGTTTTCAATATCGAAACCTTCGGTaa                                   <  2:3048272/101‑1 (MQ=255)

                                                               tcttcAAAGGTGTGCACCCGGAAGTCGGACTTACCTATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTACCGTTTTCAATATCGAAACCTTCGGTaaa                                  <  1:3048355/101‑1 (MQ=255)

                                                                cttcAAAGGTGTGCACCCGGAAGTCGGACTTACCTATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAAt                                 >  1:3048686/1‑101 (MQ=255)

                                                                 ttcAAAGGTGTGCACCCGGAAGTCGGACTTACCTATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAAtt                                >  1:3048689/1‑101 (MQ=255)

                                                                 ttcAAAGGTGTGCACCCGGAAGTCGGACTTACCTATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAAtt                                >  2:3048690/1‑101 (MQ=255)

                                                                 ttcAAAGGTGTGCACCCGGAAGTCGGACTTACCTATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAAtt                                <  2:3048336/101‑1 (MQ=255)

                                                                 ttcAAAGGTGTGCACCCGGAAGTCGGACTTACCTATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAAtt                                <  2:3048311/101‑1 (MQ=255)

                                                                  tcAAAGGTGTGCACCCGGAAGTCGGACTTACCTATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAAttt                               <  1:3048486/101‑1 (MQ=255)

                                                                   cAAAGGTGTGCACCCGGAAGTCGGACTTACCTATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTc                              >  2:3048692/1‑101 (MQ=255)

                                                                    aaaGGTGTGCACCCGGAAGTCGGACTTACCTATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTCa                             <  2:3048133/101‑1 (MQ=255)

                                                                     aaGGTGTGCACCCGGAAGTCGGACTTACCTATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTCaa                            >  2:3048697/1‑101 (MQ=255)

                                                                      aGGTGTGCACGCGGAAGTCGGACTTACCTATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTACCGTTTTc                                                       <  1:3048337/73‑1 (MQ=255)

                                                                      aGGTGTGCACCCGGAAGTCGGACTTACCTATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTCAAc                           >  1:3048698/1‑101 (MQ=255)

                                                                         tgtgCACCCGGAAGTCGGACTTACCTATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTACCGTTTTCAATATCGAAACCtt                                         >  2:3048701/1‑84 (MQ=255)

                                                                          gtgCACCCGGAAGTCGGACTTACCTATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTCAAc                           >  1:3048702/1‑97 (MQ=255)

                                                                            gCACCCGGAAGTCGGACTTACCTATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTCAACGACTTa                     <  2:3048345/101‑1 (MQ=255)

                                                                              aCCCGGAAGTCGGACTTACCTATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTCAACGACTTAAt                   <  2:3048362/101‑1 (MQ=255)

                                                                              aCCCGGAAGTCGGACTTACCTATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTCAACGACTTAAt                   >  2:3048706/1‑101 (MQ=255)

                                                                               cccGGAAGTCGGACTTACCTATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTCAACGACTTAAtt                  <  1:3048341/101‑1 (MQ=255)

                                                                                ccGGAAGTCGGACTTACCTATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTCAACGACTTAAttt                 <  1:3048318/101‑1 (MQ=255)

                                                                                ccGGAAGTCGGACTTACCTATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTCAACGACTTAAttt                 >  1:3048710/1‑101 (MQ=255)

                                                                                  ggAAGTCGGACTTACCTATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTCAACGACTTAATTTCa               >  2:3048712/1‑101 (MQ=255)

                                                                                    aaGTCGGACTTACCTATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTCAACGACTTAATTTCATc             >  1:3048714/1‑101 (MQ=255)

                                                                                      gTCGGACTTACCTATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTCAACGACTTAATTTCATCCa           >  2:3048717/1‑101 (MQ=255)

                                                                                       tCGGACTTACCTATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTCAACGACTTAAt                   <  1:3048369/92‑1 (MQ=255)

                                                                                       tCGGACTTACCTATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTCAACGACTTAATTTCATCCAg          <  1:3047512/101‑1 (MQ=255)

                                                                                          gACTTACCTATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTCAACGACTTAATTTCATCCAg          >  2:3048723/1‑98 (MQ=255)

                                                                                          gACTTACCTATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTCAACGACTTAATTTCATCCAGCGt       <  2:3048613/101‑1 (MQ=255)

                                                                                          gACTTACCTATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTCAACGACTTAATTTCATCCAGCGt       >  1:3048720/1‑101 (MQ=255)

                                                                                           aCTTACCTATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTCAACGACTTAATTTCATCCAGCGTg      <  2:3048471/101‑1 (MQ=255)

                                                                                           aCTTACCTATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTCAACGACTTAATTTCATCCAGCGTg      <  2:3048389/101‑1 (MQ=255)

                                                                                            cTTACCTATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTCAACGACTTAATTTCATCCAGCGTGa     >  1:3048727/1‑101 (MQ=255)

                                                                                            cTTACCTATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTCAACGACTTAATTTCATCCAGCGTGa     >  1:3048725/1‑101 (MQ=255)

                                                                                            cTTACCTATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTCAACGACTTAATTTCATCCAGCGTGa     <  1:3048505/101‑1 (MQ=255)

                                                                                            cTTACCTATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTCAACGACTTAATTTCATCCAGCGTGa     <  1:3048399/101‑1 (MQ=255)

                                                                                             ttACCTATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTCAACGACTTAATTTCATCCAGCGTGaa    >  1:3048730/1‑101 (MQ=255)

                                                                                              tACCTATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTCAACGACTTAATTTCATCCAGCGTGaaa   <  2:3048359/101‑1 (MQ=255)

                                                                                               aCCTATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTCAACGACTTAATTTCATCCAGCGTGAAAt  <  1:3048178/101‑1 (MQ=255)

TCTGGATAAATACCGATATTTTTCCCGGTAGAGTGATTTAACCCCTGAACAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCCATTGGGAAACGCCCCGGATAAGTCTGCACTTTTTTACCGTTTTCAATATCGAAACCTTCGGTAAATTTCAACGACTTAATTTCATCCAGCGTGAAAT  >  NC_000913/2350500‑2350695

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