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.2%	381.7 / ‑3.8	125	coding (1792/2004 nt)	tktB	transketolase 2 isozyme
Reads supporting (aligned to +/- strand): ref base T (0/0); major base C (62/62); minor base . (1/0); total (63/62)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.28e-01

Only 100 of 130 total aligned reads displayed.

GAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCA  >  NC_000913/2579357‑2579536

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

  gaAATTAGCAGGAGAAGGTCGCAATGTTCGCCTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGgt                                                                                       <  2:921699/93‑1 (MQ=255)

  gaAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGCCCGATATTTTCGACGCCCAGGATGAGGAATATTGGGAGCCggg                                                                                       >  2:670122/1‑92 (MQ=255)

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

     aTTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGtt                                                                                    >  2:4175261/1‑93 (MQ=255)

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

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

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

               gaAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAAc                                                                          >  1:1084942/1‑93 (MQ=255)

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

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

                aaGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACg                                                                         >  1:3940242/1‑93 (MQ=255)

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

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

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

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

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

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

                       cAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGc                                                                  <  2:93988/93‑1 (MQ=255)

                       cAATGTACGCGTAGTTTCCACGCCCTCGACAGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGc                                                                  <  1:2260795/93‑1 (MQ=255)

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

                         aTGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTc                                                                >  2:1604297/1‑93 (MQ=255)

                          tGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGATTCGGTGTTGCCTTCTATCGTTGCGGCTcg                                                               >  2:2219204/1‑93 (MQ=255)

                          tGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTcg                                                               >  1:2143500/1‑93 (MQ=255)

                          tGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTcg                                                               >  1:1677987/1‑93 (MQ=255)

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

                             aCGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTcgcgg                                                            >  2:43967/1‑92 (MQ=255)

                             aCGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTcgccg                                                            >  2:2569525/1‑91 (MQ=255)

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

                               gcgTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTgg                                                          <  2:3431170/93‑1 (MQ=255)

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

                                 gTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCg                                                        <  1:2618129/93‑1 (MQ=255)

                                 gTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCg                                                        <  1:2618130/93‑1 (MQ=255)

                                  tAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGGGTGgggg                                                       >  2:1138986/1‑93 (MQ=255)

                                  tAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCgg                                                       <  1:833633/93‑1 (MQ=255)

                                  tAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCgg                                                       >  2:2463068/1‑93 (MQ=255)

                                  tAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCgg                                                       >  2:4265337/1‑93 (MQ=255)

                                   aGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCggg                                                      >  1:3370445/1‑92 (MQ=255)

                                   aGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCCCGTGGCggg                                                      >  1:3472591/1‑92 (MQ=255)

                                    gTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTg                                                     >  1:778922/1‑93 (MQ=255)

                                     tttCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTgg                                                    <  1:3101609/93‑1 (MQ=255)

                                         ccTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGc                                                <  2:1552822/93‑1 (MQ=255)

                                         ccTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGc                                                <  2:2871284/93‑1 (MQ=255)

                                         ccTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGGGGCGGTGGGAGc                                                >  1:2889514/1‑93 (MQ=255)

                                          cTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGGAGCa                                               >  2:4356529/1‑93 (MQ=255)

                                           tGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTGTAACGTTGCGGCTCGCGTGGCGGGGGGAGCAg                                              >  2:3668986/1‑93 (MQ=255)

                                           tGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGGGGAAGCAg                                              >  2:3436960/1‑93 (MQ=255)

                                              ccTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTa                                           <  1:3045607/93‑1 (MQ=255)

                                               cTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGGAGCAGGTAt                                          >  2:4276452/1‑93 (MQ=255)

                                               cTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTAt                                          <  1:4348200/93‑1 (MQ=255)

                                               cTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTAt                                          <  1:2034311/93‑1 (MQ=255)

                                               cTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTAt                                          <  1:1512200/93‑1 (MQ=255)

                                               cTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGGGGGAGCAGGTAt                                          >  2:305522/1‑93 (MQ=255)

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

                                                tCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGGGGCGGCGGGAGCCGGGAtt                                         >  1:969909/1‑93 (MQ=255)

                                                  gACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGGGGGAGCCGGTAtttc                                       >  1:726703/1‑91 (MQ=255)

                                                  gACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGGGGGGGGGGGAGGAGGGAtttg                                       >  1:3249535/1‑91 (MQ=255)

                                                   aCCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGcc                                      >  1:1245845/1‑93 (MQ=255)

                                                    ccGATATTTTGGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGAGGGGGGGGAAGAAGGGATTGctg                                     >  2:2133099/1‑91 (MQ=255)

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

                                                    ccGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGGGGGAGCAGGTATTGCCg                                     >  2:1522072/1‑93 (MQ=255)

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

                                                      gATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGAt                                   <  2:4646106/93‑1 (MQ=255)

                                                      gATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCCTGGCGGGGGAAGCCGGGATTGCCggt                                   >  1:3107571/1‑91 (MQ=255)

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

                                                       atatTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGGATTGCCGAtt                                  >  1:4542095/1‑93 (MQ=255)

                                                       atatTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGGGGGAGCAGGGATTGCCGGtt                                  >  1:1203969/1‑93 (MQ=255)

                                                        tatTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGGAGCAGGTATTGGCGATTa                                 >  1:3508491/1‑93 (MQ=255)

                                                        tatTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGGGGGAGCAGGTAATGCCGATTa                                 >  1:4497973/1‑93 (MQ=255)

                                                         atTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTAc                                <  2:906752/93‑1 (MQ=255)

                                                          ttttCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGGAGGCGGGATTGCCGATTTCt                               >  2:3459249/1‑93 (MQ=255)

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

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

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

                                                              cGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGGGGGCGCAGGGATTGTCGATTACTGGTa                           >  2:4406253/1‑93 (MQ=255)

                                                                aCGCCCAGGATGAGCAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACa                         <  2:1380359/93‑1 (MQ=255)

                                                                 cGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACaa                        <  1:4613079/93‑1 (MQ=255)

                                                                 cGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACaa                        <  1:4173449/93‑1 (MQ=255)

                                                                     cAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTAt                    >  2:1532639/1‑93 (MQ=255)

                                                                     cAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTAt                    <  2:391527/93‑1 (MQ=255)

                                                                      aGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATg                   <  1:636797/93‑1 (MQ=255)

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

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

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

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

                                                                          tGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTgg               <  1:4219183/93‑1 (MQ=255)

                                                                          tGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTgg               >  1:4037663/1‑93 (MQ=255)

                                                                           gAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGGGGTGGGAGCAGGGATTGGCGATTACTGGTACAAGTATGTTGGt              >  1:4243711/1‑93 (MQ=255)

                                                                            aGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGGAGCAGGTATTGCCGATTACTGGTACCAGTATGTTGGTc             >  2:3000976/1‑93 (MQ=255)

                                                                            aGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGGATTGCCGCTTACTGGTACAAGTATGTTgggc             >  1:277084/1‑91 (MQ=38)

                                                                             ggAATATCGGGAG‑CGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGGGGGAGCAGGGATTGCCGATTACTGGGACAAGTATGGTGGGCTg           >  2:4118254/1‑93 (MQ=255)

                                                                                atatCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGGGGGAGCAGGGATTGCCGAGTACTGGTACAAGTATGGTGGTCTGaa         >  1:2322451/1‑93 (MQ=38)

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

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

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

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

                                                                                    cGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGGAGCAGGGATTGCCGATTACTGGTACCAGTATGTTGGTCTGAAAggg     >  2:2406171/1‑93 (MQ=255)

                                                                                    cGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAggg     <  1:811161/93‑1 (MQ=255)

                                                                                     gggAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAgggg    <  1:1628072/93‑1 (MQ=255)

                                                                                       gAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGGAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCa  >  2:1805553/1‑93 (MQ=255)

GAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCA  >  NC_000913/2579357‑2579536

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