BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,579,447	T→C	S598P (TCG→CCG)	tktB →	transketolase 2 isozyme

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,579,447	0	T	C	99.5%	596.9 / NA	190	S598P (TCG→CCG)	tktB	transketolase 2 isozyme
Reads supporting (aligned to +/- strand): ref base T (0/0); major base C (91/98); minor base A (1/0); total (92/98)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 195 total aligned reads displayed.

GAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGC  >  NC_000913/2579357‑2579535

gagaAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCg                                                                                        <  1:121682/93‑1 (MQ=255)

  gaAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGgt                                                                                      <  2:3919188/93‑1 (MQ=255)

  gaAATGAGCAGGAGAAGGTCGCAATGTACGCTTAGTTTCCCTGCCCTCGACCGATATTTTCGCCGCCCAGGATGAGGAATATCGGGAGCCGgt                                                                                      <  2:3182913/93‑1 (MQ=255)

  gaAATGAGCAGGAGAAGGTCGAAATGTACGCGTAGTTTCCCTGCCCGCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGgt                                                                                      <  2:4535115/93‑1 (MQ=255)

   aaaTTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGgtg                                                                                     <  1:5511442/93‑1 (MQ=255)

    aaTTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGgtgt                                                                                    >  1:3550371/1‑93 (MQ=255)

    aaTTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGgggt                                                                                    >  1:1828451/1‑93 (MQ=255)

      ttAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTg                                                                                  <  1:5322626/93‑1 (MQ=255)

        aGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGcc                                                                                >  1:6189712/1‑93 (MQ=255)

           aGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTc                                                                             <  1:5862962/93‑1 (MQ=255)

            ggAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCt                                                                            >  2:2944096/1‑93 (MQ=255)

            ggAGAAGGTCGCAATGTACGCGTAGTTCCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCt                                                                            >  2:5812837/1‑93 (MQ=255)

               gaAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAAc                                                                         <  2:4834033/93‑1 (MQ=255)

                aaGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACg                                                                        <  1:1028829/93‑1 (MQ=255)

                 aGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGt                                                                       >  2:1425625/1‑93 (MQ=255)

                  ggTCGGAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGtt                                                                      <  1:4675992/93‑1 (MQ=255)

                  ggTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGtt                                                                      >  1:22779/1‑93 (MQ=255)

                   gTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTg                                                                     <  2:3868273/93‑1 (MQ=255)

                     cGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCg                                                                   >  1:974620/1‑93 (MQ=255)

                      gCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCgg                                                                  >  1:251965/1‑93 (MQ=255)

                      gCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCgg                                                                  <  2:5063386/93‑1 (MQ=255)

                       cAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGc                                                                 >  2:4153323/1‑93 (MQ=255)

                       cAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGc                                                                 >  2:3544267/1‑93 (MQ=255)

                        aaTGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCt                                                                <  1:6248700/93‑1 (MQ=255)

                        aaTGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCt                                                                <  2:2838188/93‑1 (MQ=255)

                         aTGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTc                                                               <  2:5187141/93‑1 (MQ=255)

                            tACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTcgcg                                                            <  1:4115909/93‑1 (MQ=255)

                            tACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTcgcg                                                            >  1:3702149/1‑93 (MQ=255)

                            tACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTcgcg                                                            <  1:1386410/93‑1 (MQ=255)

                              cgcgTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTg                                                          >  1:5341341/1‑93 (MQ=255)

                              cgcgTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTg                                                          <  2:2181411/93‑1 (MQ=255)

                               gcgTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTgg                                                         >  1:5190532/1‑93 (MQ=255)

                                cgTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGc                                                        <  2:4738071/93‑1 (MQ=255)

                                cgTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGc                                                        <  1:394255/93‑1 (MQ=255)

                                     tttCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTgg                                                   >  2:586533/1‑93 (MQ=255)

                                      tccccTGCCCTCGACCGATATTTTCTACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGa                                                  <  2:4269616/91‑1 (MQ=255)

                                       tCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGaa                                                 <  1:2655854/93‑1 (MQ=255)

                                       tCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGaa                                                 <  2:5492644/93‑1 (MQ=255)

                                        cccTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAg                                                <  1:4611796/93‑1 (MQ=255)

                                         ccTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGGGGGAgg                                               >  2:1872371/1‑92 (MQ=255)

                                         ccTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGGGGGAGc                                               >  2:208148/1‑93 (MQ=255)

                                           tGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGGGGCGGGGGGAGCAg                                             >  2:2551907/1‑93 (MQ=255)

                                            gCCCTCGACCGATATTTTCGACGGCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAgg                                            <  2:4689862/93‑1 (MQ=255)

                                               cTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCCCGTGGCGGGGGGAGAAgggtt                                         >  1:1617662/1‑90 (MQ=255)

                                                tCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTAtt                                        <  2:5398759/93‑1 (MQ=255)

                                                tCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTAtt                                        <  1:1117751/93‑1 (MQ=255)

                                                 cGACTGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTg                                       <  1:4507465/93‑1 (MQ=255)

                                                  gACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGc                                      <  1:5963482/93‑1 (MQ=255)

                                                   aCCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGcc                                     <  1:2144224/93‑1 (MQ=255)

                                                   aCCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGGGGCGGTGGGAGCAGGGATTGcc                                     >  2:3347501/1‑93 (MQ=255)

                                                   aCCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGACTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAACAGGTATTGcg                                     >  1:4254773/1‑92 (MQ=255)

                                                    ccGATTTTTTCTACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCg                                    <  2:4474243/93‑1 (MQ=255)

                                                    ccGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCg                                    <  1:3935379/93‑1 (MQ=255)

                                                    ccGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGGGGGAGCAGGGATTGCCg                                    >  2:129860/1‑93 (MQ=255)

                                                    ccGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCCCTCGGTGGTCGACGGGGGGATTGCCg                                    >  2:5587859/1‑93 (MQ=255)

                                                     cGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGa                                   <  2:144605/93‑1 (MQ=255)

                                                     cGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGGATTGCCGa                                   >  1:4638003/1‑93 (MQ=255)

                                                     cGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGGGGGAGCAGGGATTGCCGa                                   >  2:3761214/1‑93 (MQ=255)

                                                     cGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGGGGAAGCAGGTATTGCCGa                                   >  2:752535/1‑93 (MQ=255)

                                                       atatTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGAtt                                 <  1:3082347/93‑1 (MQ=255)

                                                        tatTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTa                                <  2:3039460/93‑1 (MQ=255)

                                                        tatTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTa                                <  1:1472778/93‑1 (MQ=255)

                                                         atTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGGGGGGGGGGGAGGAGGGGTTGCCGGtttt                               >  2:4502405/1‑91 (MQ=255)

                                                           tttCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTg                             >  2:2747294/1‑93 (MQ=255)

                                                            ttCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTgg                            <  2:4804074/93‑1 (MQ=255)

                                                             tCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGt                           <  1:5142668/93‑1 (MQ=255)

                                                             tCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGt                           <  2:2288254/93‑1 (MQ=255)

                                                               gACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTAc                         <  2:4973119/93‑1 (MQ=255)

                                                                  gCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCCTGGCGGTGGAAGCAGGTGTTGCCGCTTACTGGTACAAg                      >  1:4817379/1‑93 (MQ=255)

                                                                     cAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGGAGCAGGTATTGCCGATTACTGGTACAAGTAt                   >  2:17017/1‑93 (MQ=255)

                                                                     cAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGGAGCAGGGATTGCCGATTACTGGTACAAGTAt                   >  1:4999268/1‑93 (MQ=255)

                                                                      aGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGGAGCAGGGATTGCCGATTACTGGGACAAATATg                  >  2:2659400/1‑93 (MQ=255)

                                                                      aGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATg                  >  1:5889701/1‑93 (MQ=255)

                                                                       ggATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGt                 <  2:1637798/93‑1 (MQ=255)

                                                                       ggATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGGATTGCCGATTACTGGTACAAGTATGt                 >  1:3037503/1‑93 (MQ=255)

                                                                       ggATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGGATTGCCGATTACTGGTACAAGTATGt                 >  1:5691698/1‑93 (MQ=255)

                                                                        gATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGtt                <  2:2818710/93‑1 (MQ=255)

                                                                        gATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGCGGCGGTGGGAGCAGGGATTGCCGCTTACTGGGCCCAGTATGtt                >  2:2497370/1‑93 (MQ=255)

                                                                         aTGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTg               <  2:5356651/93‑1 (MQ=255)

                                                                         aTGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGGGGAAGCAAGTATTGCCGATTACTGGGACAAGTATGGTg               >  2:6286539/1‑93 (MQ=255)

                                                                          tGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGGAGCAGGTATTGCCGATTACTGGGACAAGTATGTTgg              >  2:361841/1‑93 (MQ=255)

                                                                           gAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGCATTGCCGGTTACTGGTACCAGGATGTTTGt             >  2:4635415/1‑93 (MQ=255)

                                                                            aGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTc            <  1:4159880/93‑1 (MQ=255)

                                                                             ggAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCt           <  1:5163025/93‑1 (MQ=255)

                                                                             ggAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCt           <  1:473285/93‑1 (MQ=255)

                                                                              gAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGGAGCAGGGATTGCCGATTACTGGGACACGTATGTTGGTCTg          >  2:527375/1‑93 (MQ=255)

                                                                              gAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTg          <  1:3851325/93‑1 (MQ=255)

                                                                               aaTATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGGATTGCCGATTACTGGTACAAGTATGTTGGTCTGa         >  1:2311815/1‑93 (MQ=255)

                                                                                atatCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGGAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGaa        >  2:1931679/1‑93 (MQ=255)

                                                                                atatCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGaa        <  1:5443535/93‑1 (MQ=255)

                                                                                 tatCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGaaa       <  2:243749/93‑1 (MQ=255)

                                                                                 tatCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGaaa       <  2:4252942/93‑1 (MQ=255)

                                                                                 tatCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGaaa       >  2:3475278/1‑93 (MQ=255)

                                                                                  atCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAg      <  1:1031944/93‑1 (MQ=255)

                                                                                  atCGGGAGCCGGTGTGGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAg      <  1:1066555/93‑1 (MQ=255)

                                                                                   tCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAgg     <  1:6057548/93‑1 (MQ=255)

                                                                                   tCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGGGGGAGCAGGTATTGCCGATTACTTGTACAAGTATGTTGGTCTGAAAgg     >  2:3277211/1‑93 (MQ=38)

                                                                                      ggAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGGAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGc  >  2:3539704/1‑93 (MQ=255)

                                                                                      ggAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGc  >  1:6094806/1‑93 (MQ=255)

                                                                                      ggAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGACGGTGGAAGCATGTATTGCCGATTACTGGTACAAGTATGTAGGTCTGAAAGGGGc  <  1:3195938/93‑1 (MQ=38)

GAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGC  >  NC_000913/2579357‑2579535

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