BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	618,600	A→G	H375R (CAC→CGC)	fepE →	FepE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	618,600	0	A	G	99.7%	1129.1 / NA	327	H375R (CAC→CGC)	fepE	FepE
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (144/182); total (144/183)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 343 total aligned reads displayed.

CCGCGTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGC  >  NC_000913/618510‑618689

ccGCGTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCt                                                                                         >  1:6215643/1‑93 (MQ=255)

ccGCGTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCt                                                                                         <  2:2481108/93‑1 (MQ=255)

 cgcgTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCtt                                                                                        <  1:5588898/93‑1 (MQ=255)

  gcgTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTa                                                                                       >  1:3757780/1‑93 (MQ=255)

  gcgTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTa                                                                                       >  1:6758193/1‑93 (MQ=255)

  gcgTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTa                                                                                       >  1:14146/1‑93 (MQ=255)

  gcgTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTa                                                                                       <  2:8106064/93‑1 (MQ=255)

   cgTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAg                                                                                      <  1:6023626/93‑1 (MQ=255)

    gTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGt                                                                                     <  2:7217311/93‑1 (MQ=255)

    gTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGt                                                                                     <  2:4940701/93‑1 (MQ=255)

     ttGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGtt                                                                                    >  2:5847636/1‑93 (MQ=255)

     ttGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGtt                                                                                    >  1:4456021/1‑93 (MQ=255)

     ttGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGtt                                                                                    >  1:6840594/1‑93 (MQ=255)

      tGATCGGCGGGATGGTGGCTTGTGGTGGCGCGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGttt                                                                                   <  2:770132/93‑1 (MQ=255)

      tGATCGGCGGGATGGTGGCGTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTGAGttt                                                                                   <  1:6980758/93‑1 (MQ=255)

       gATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTa                                                                                  <  1:453005/93‑1 (MQ=255)

       gATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTa                                                                                  >  1:742137/1‑93 (MQ=255)

        aTCGGCGGGATGGTGGCTTGTGGTGGCGGGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGCCCGCTGAGTTTaa                                                                                 >  1:6471046/1‑93 (MQ=255)

         tCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAAc                                                                                <  2:2402842/93‑1 (MQ=255)

             cggGATTGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTAtt                                                                            <  2:3327058/93‑1 (MQ=255)

             cggGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTAtt                                                                            <  2:3738107/93‑1 (MQ=255)

             cggGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTAtt                                                                            <  2:4003688/93‑1 (MQ=255)

              gggATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTAttt                                                                           <  1:3569342/93‑1 (MQ=255)

              gggATGGTGGCTTGTGGTGGCGTGCTGTTGCACTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTAttt                                                                           <  2:7746684/93‑1 (MQ=255)

               ggATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTg                                                                          >  1:6293632/1‑93 (MQ=255)

               ggATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTg                                                                          <  2:8354876/93‑1 (MQ=255)

                  tggtggCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCgg                                                                       >  1:1053525/1‑93 (MQ=255)

                  tggtggCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCgg                                                                       <  1:5250859/93‑1 (MQ=255)

                   ggtggCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGt                                                                      <  1:5165020/93‑1 (MQ=255)

                    gtggCTTGTGGTGGCGTGCTGTTGCGCTATGCTATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGtt                                                                     >  1:8005643/1‑93 (MQ=255)

                     tggCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTg                                                                    <  1:8076962/93‑1 (MQ=255)

                     tggCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTg                                                                    <  1:1896917/93‑1 (MQ=255)

                      ggCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGa                                                                   >  2:255029/1‑93 (MQ=255)

                      ggCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGa                                                                   >  2:161944/1‑93 (MQ=255)

                      ggCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGa                                                                   >  2:1286427/1‑93 (MQ=255)

                      ggCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGa                                                                   >  1:1016781/1‑93 (MQ=255)

                      ggCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGa                                                                   >  2:4333077/1‑93 (MQ=255)

                       gCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGaa                                                                  <  2:2051971/93‑1 (MQ=255)

                       gCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGaa                                                                  <  1:942154/93‑1 (MQ=255)

                        cTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCCGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAg                                                                 >  1:3642747/1‑93 (MQ=255)

                        cTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAg                                                                 <  1:7536044/93‑1 (MQ=255)

                         ttGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAgg                                                                >  2:4149421/1‑93 (MQ=255)

                         ttGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAgg                                                                <  2:3447116/93‑1 (MQ=255)

                            tggtggCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGc                                                             <  1:4830335/93‑1 (MQ=255)

                              gtggCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCg                                                           <  2:7541941/93‑1 (MQ=255)

                              gtggCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCg                                                           >  2:3671057/1‑93 (MQ=255)

                                ggCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGtt                                                         <  1:6559274/93‑1 (MQ=255)

                                 gCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTc                                                        >  1:5351056/1‑93 (MQ=255)

                                  ggTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCg                                                       <  1:2261998/92‑1 (MQ=255)

                                   gTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGt                                                      <  1:4708807/93‑1 (MQ=255)

                                        gTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAg                                                 >  2:868224/1‑93 (MQ=255)

                                         ttGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGc                                                >  2:1207473/1‑93 (MQ=255)

                                              cTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGcac                                           >  1:394227/1‑93 (MQ=255)

                                              cTATGCGATGGCATCAAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGcac                                           <  1:3431277/93‑1 (MQ=255)

                                                 tGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGcaccac                                        <  2:6967792/93‑1 (MQ=255)

                                                    gATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAg                                     <  2:3871542/93‑1 (MQ=255)

                                                     aTGGCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAgc                                    >  2:7383464/1‑93 (MQ=255)

                                                       ggCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAgcgg                                  >  1:5436374/1‑93 (MQ=255)

                                                       ggCATCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAgcgg                                  <  2:1413593/93‑1 (MQ=255)

                                                           tCCAGAAAACAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGtt                              >  2:5033911/1‑93 (MQ=255)

                                                                  aaCAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCaa                       <  2:3066520/93‑1 (MQ=255)

                                                                   aCAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAAc                      <  2:4640540/93‑1 (MQ=255)

                                                                   aCAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCAAGCGGCACCACAAGCGGCGTTCCGGCAAc                      <  2:5888732/93‑1 (MQ=255)

                                                                   aCAGGAGGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAAc                      <  1:8384419/93‑1 (MQ=255)

                                                                    cAGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACt                     <  1:1152248/93‑1 (MQ=255)

                                                                     aGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTg                    >  1:5654724/1‑93 (MQ=255)

                                                                     aGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTg                    >  1:873571/1‑93 (MQ=255)

                                                                     aGGATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCAAGCGGCACCACAAGCGGCGTTCCGGCAACTg                    <  1:4173229/93‑1 (MQ=255)

                                                                      ggATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTgg                   <  2:672191/93‑1 (MQ=255)

                                                                      ggATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTgg                   <  2:5287636/93‑1 (MQ=255)

                                                                       gATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGa                  <  1:1379878/93‑1 (MQ=255)

                                                                       gATGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGa                  <  2:4397922/93‑1 (MQ=255)

                                                                        gtGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGAt                 <  1:6628174/92‑1 (MQ=255)

                                                                        aTGCCATGATGTCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCTTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGAt                 <  1:4448797/93‑1 (MQ=255)

                                                                        aTGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGAt                 <  2:6563788/93‑1 (MQ=255)

                                                                        aTGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGAt                 <  2:7141709/93‑1 (MQ=255)

                                                                         tGCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATc                >  2:2176627/1‑93 (MQ=255)

                                                                          gCCATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCg               <  1:2977051/93‑1 (MQ=255)

                                                                            cATGATGGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTc             <  2:794958/93‑1 (MQ=255)

                                                                              tgatgGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCaa           <  1:7266464/93‑1 (MQ=255)

                                                                              tgatgGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCaa           >  2:6856896/1‑93 (MQ=255)

                                                                              tgatgGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCaa           <  2:1098181/93‑1 (MQ=255)

                                                                               gatgGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAAt          <  1:2510455/93‑1 (MQ=255)

                                                                                atgGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATg         >  1:5012349/1‑93 (MQ=255)

                                                                                atgGCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATg         <  1:7959171/93‑1 (MQ=255)

                                                                                  ggCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGAt       <  1:3049713/93‑1 (MQ=255)

                                                                                  ggCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGAt       <  1:2089201/93‑1 (MQ=255)

                                                                                   gCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATc      >  1:5882132/1‑93 (MQ=255)

                                                                                   gCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATc      >  1:1682598/1‑93 (MQ=255)

                                                                                   gCAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATc      >  1:3332775/1‑93 (MQ=255)

                                                                                    cAGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCa     <  2:8102404/93‑1 (MQ=255)

                                                                                     aGACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCAt    <  1:1574815/93‑1 (MQ=255)

                                                                                     aGACCGCTTAGTTTAACTATTTGCGGTGGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCAt    <  2:6462090/93‑1 (MQ=255)

                                                                                      gACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATg   >  2:5699805/1‑93 (MQ=255)

                                                                                      gACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATg   >  2:5645766/1‑93 (MQ=255)

                                                                                      gACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATg   >  1:7055972/1‑93 (MQ=255)

                                                                                      gACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATg   <  1:779685/93‑1 (MQ=255)

                                                                                      gACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATg   >  2:4099867/1‑93 (MQ=255)

                                                                                      gACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATg   <  1:2826304/93‑1 (MQ=255)

                                                                                      gACCGCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATg   >  2:2978213/1‑93 (MQ=255)

CCGCGTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGC  >  NC_000913/618510‑618689

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