BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,081,054	A→G	T93A (ACC→GCC)	sbcB →	exonuclease I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,081,054	0	A	G	98.7%	497.5 / ‑2.9	158	T93A (ACC→GCC)	sbcB	exonuclease I
Reads supporting (aligned to +/- strand): ref base A (0/2); new base G (79/77); total (79/79)
Fisher's exact test for biased strand distribution p-value = 4.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.12e-01

Only 100 of 178 total aligned reads displayed.

GCCGTATTAATTACCGGTATTACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTACCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTCACACGCAACATTTTTTATCGTAATTTCTACGAT  >  NC_000913/2080964‑2081143

gCCGTATTAATTACCGGTATTACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTtgcc                                                                                         <  1:4170781/93‑1 (MQ=255)

  cGTATTAATTACCGGTATTACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTtgccgt                                                                                       <  2:1626274/93‑1 (MQ=255)

    tATTAATTACCGGTATTACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTtgccgtac                                                                                     >  2:2329932/1‑93 (MQ=255)

     attaattaCCGGTATTACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTtgccgtacc                                                                                    >  2:3867660/1‑93 (MQ=255)

      ttaattaCCGGTATTACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTtgccgtaccg                                                                                   <  2:2831977/93‑1 (MQ=255)

         attaCCGGTATTACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAg                                                                                >  1:3902852/1‑93 (MQ=255)

          ttaCCGGTATTACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGa                                                                               <  2:4316047/93‑1 (MQ=255)

                gTATTACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGtc                                                                         >  1:1796935/1‑92 (MQ=255)

                gTATTACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTa                                                                         >  2:1397712/1‑93 (MQ=255)

                gTATTACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTa                                                                         >  2:1429413/1‑93 (MQ=255)

                gTATTACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTa                                                                         >  1:1796806/1‑93 (MQ=255)

                 tATTACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTAt                                                                        <  2:711892/93‑1 (MQ=255)

                   ttACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTc                                                                      <  1:823082/93‑1 (MQ=255)

                   ttACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTc                                                                      >  1:2088787/1‑93 (MQ=255)

                   ttACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTc                                                                      <  1:274833/93‑1 (MQ=255)

                     aCCCCGCAAGAAGCCCGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCCCTCGCTTTTTGCCGTACCGAAGACCTGTATTCTg                                                                    <  1:2275852/93‑1 (MQ=255)

                      ccccGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTgg                                                                   >  2:3864322/1‑93 (MQ=255)

                        ccGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGc                                                                 <  1:1059816/93‑1 (MQ=255)

                        ccGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGc                                                                 <  1:1474503/93‑1 (MQ=255)

                          gCAGGAAGCACGGGCGAAAGGAGAAAACGAAGGCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTa                                                               >  2:3941995/1‑93 (MQ=255)

                          gCAGGAAGCACGGGCGAAAGGAGAAAAAGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTa                                                               >  2:4614829/1‑93 (MQ=255)

                           cAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTac                                                              <  1:4625816/93‑1 (MQ=255)

                           cAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTac                                                              >  1:3742193/1‑93 (MQ=255)

                             ggAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTacaa                                                            >  2:3479146/1‑93 (MQ=255)

                             ggAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTacaa                                                            >  1:1224273/1‑93 (MQ=255)

                             ggAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTacaa                                                            >  1:4006968/1‑93 (MQ=255)

                              gAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTacaac                                                           >  2:1687434/1‑93 (MQ=255)

                              gAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTacaac                                                           >  1:20315/1‑93 (MQ=255)

                                aGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTacaacac                                                         >  2:2536224/1‑92 (MQ=255)

                                 gCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTACCGTACCGAAGACCTGTATTCTGGGCTACAACAAt                                                        <  2:4041661/93‑1 (MQ=255)

                                   aCGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAAtgt                                                      <  1:3418836/93‑1 (MQ=255)

                                    cGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAAtgtg                                                     >  2:1081572/1‑93 (MQ=255)

                                     gggCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGc                                                    <  1:1490340/93‑1 (MQ=255)

                                        cGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGtt                                                 <  1:524534/93‑1 (MQ=255)

                                         gAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGttt                                                <  2:4194098/93‑1 (MQ=255)

                                         gAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGttt                                                >  1:3042572/1‑93 (MQ=255)

                                          aaaGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTc                                               >  1:669619/1‑93 (MQ=255)

                                           aaGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTcg                                              >  2:396204/1‑93 (MQ=255)

                                            aGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTcga                                             <  1:109375/93‑1 (MQ=255)

                                            aGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTcga                                             >  1:4130004/1‑93 (MQ=255)

                                            aGGAGAAAACGAAGCCGCGTTTGACGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTcga                                             <  2:2602983/93‑1 (MQ=255)

                                             ggAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTcgac                                            <  2:4645439/93‑1 (MQ=255)

                                             ggAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTcgac                                            >  2:3809167/1‑93 (MQ=255)

                                              gagaAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTcgacg                                           <  2:486254/93‑1 (MQ=255)

                                              gagaAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTcgacg                                           >  1:3988717/1‑93 (MQ=255)

                                               agaAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTcgacga                                          >  1:1703630/1‑93 (MQ=255)

                                                gaAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTcgacgac                                         <  2:3259814/93‑1 (MQ=255)

                                                gaAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTcgacgac                                         >  2:3552210/1‑93 (MQ=255)

                                                  aaaCGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTcgacgacga                                       >  1:3168572/1‑93 (MQ=255)

                                                   agcGAATCCGGGTTTGCCGCCTGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGaa                                      <  2:2971930/91‑1 (MQ=255)

                                                   aaCGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGaa                                      <  1:4428814/93‑1 (MQ=255)

                                                   aaCGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGaa                                      <  2:3357630/93‑1 (MQ=255)

                                                     cGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGt                                    <  2:1262966/93‑1 (MQ=255)

                                                      gAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTc                                   <  2:1085566/93‑1 (MQ=255)

                                                      gAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTc                                   <  1:3074313/93‑1 (MQ=255)

                                                       aaGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTca                                  >  2:3425518/1‑93 (MQ=255)

                                                        aGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTcac                                 >  1:3867288/1‑93 (MQ=255)

                                                          ccGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTCacac                               <  2:552801/93‑1 (MQ=255)

                                                            gcgTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTCACACGc                             <  2:4444069/93‑1 (MQ=255)

                                                              gTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTCACACGCaa                           <  1:1997417/93‑1 (MQ=255)

                                                               tttGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTCACACGCAAc                          <  1:2095727/93‑1 (MQ=255)

                                                               tttGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTCACACGCAAc                          >  2:2181966/1‑93 (MQ=255)

                                                                  gccgccCGTATGCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTCACACGCAACAtt                       <  2:2292365/93‑1 (MQ=255)

                                                                   ccgccCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTCACACGCAACAttt                      >  1:3614112/1‑93 (MQ=255)

                                                                     gccCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTCACACGCAACAttttt                    >  1:1636149/1‑93 (MQ=255)

                                                                     gccCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTCACACCTCACAttttt                    >  2:1913632/1‑93 (MQ=255)

                                                                      cccGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTCACACGCAACAtttttt                   >  2:3510006/1‑93 (MQ=255)

                                                                       ccGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTCACACGCAACATTTTTTa                  <  2:933796/93‑1 (MQ=255)

                                                                       ccGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTCACACGCAACATTTTTTa                  <  1:557997/93‑1 (MQ=255)

                                                                        cGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTCACACGCAACATTTTTTAt                 >  1:3101948/1‑93 (MQ=255)

                                                                        cGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTCACACGCAACATTTTTTAt                 <  1:2291500/93‑1 (MQ=255)

                                                                        cGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTCACACGCAACATTTTTTAt                 <  2:923380/93‑1 (MQ=255)

                                                                         ggaTTCGCTCGCTTTTTGCCGTCACGAAGACCCGGATTCTTGGCTACAACAATGTGCGTTTCGACGACGAAGTCACCCGCAACATTTTTTATc                <  1:4480375/91‑1 (MQ=255)

                                                                         gTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTCACACGCAACATTTTTTATc                >  1:1153633/1‑93 (MQ=255)

                                                                         gTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTCACACGCAACATTTTTTATc                >  2:3100828/1‑93 (MQ=255)

                                                                          tATTCACTCACTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTCACACGCAACATTTTTTATCg               <  2:2099086/93‑1 (MQ=255)

                                                                           aTTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGACGTCACACGCAACATTTTTTATCGt              >  1:2652419/1‑93 (MQ=255)

                                                                            ttCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTCACACGCAACATTTTTTATCGTa             >  1:1567206/1‑93 (MQ=255)

                                                                            ttCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTCACACGCAACATTTTTTATCGTa             >  1:3478547/1‑93 (MQ=255)

                                                                              cACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGaa                                      >  2:891058/1‑66 (MQ=255)

                                                                              cACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGaa                                      <  1:891058/66‑1 (MQ=255)

                                                                              cACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTCACACGCAACATTTTTTATCGTAAt           <  1:331861/93‑1 (MQ=255)

                                                                              cACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTCACAAGCAACATTTTTTATCGTAAt           >  2:1839182/1‑93 (MQ=255)

                                                                               ccTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTCACACGCAACATTTTTTATCGTAAtt          <  1:642966/92‑1 (MQ=255)

                                                                               aCTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTCACACGCAACATTTTTTATCGTAAtt          <  2:2283659/93‑1 (MQ=255)

                                                                               aCTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTCACACGCAACATTTTTTATCGTAAtt          <  1:3006615/93‑1 (MQ=255)

                                                                               aCTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTCACACGCAACATTTTTTATCGTAAtt          >  2:3312977/1‑93 (MQ=255)

                                                                                ctagCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTCTCGACGACCAAGTCACACGCAACAGTTTTTATCGTAAttt         >  2:413172/4‑93 (MQ=255)

                                                                                cTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTCACACGCAACATTTTTTATCGTAAttt         >  2:1769010/1‑93 (MQ=255)

                                                                                cTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTCACACGCAACATTTTTTATCGTAAttt         >  2:1961313/1‑93 (MQ=255)

                                                                                 tCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTCACACGCAACATTTTTTATCGTAATTTc        >  1:2794186/1‑93 (MQ=255)

                                                                                 tCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTCACACGCAACATTTTTTATCGTAATTTc        >  2:1721234/1‑93 (MQ=255)

                                                                                 tCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTCACACGCAACATTTTTTATCGTAATTTc        <  2:704995/93‑1 (MQ=255)

                                                                                  cGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTCACACGCAACATTTTTTATCGTAAtttat       >  2:3673517/1‑91 (MQ=255)

                                                                                  cGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTCACACGCAACATTTTTTATCGTAATTTCt       <  1:3902469/93‑1 (MQ=255)

                                                                                   gCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTCACACGCAACATTTTTTATCGTAATTTCTa      <  1:789732/93‑1 (MQ=255)

                                                                                   gCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTCACACGCAACATTTTTTATCGTAATTTCTa      <  1:4173973/93‑1 (MQ=255)

                                                                                   gCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTCACACGCAACATTTTTTATCGTAATTTCTa      <  1:4312110/93‑1 (MQ=255)

                                                                                   gCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTCACACGCAACATTTTTTATCGTAATTTCTa      <  2:672359/93‑1 (MQ=255)

                                                                                       tttGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTCACACGCAACATTTTTTATCGTAATTTCTACGAt  >  1:62322/1‑93 (MQ=255)

GCCGTATTAATTACCGGTATTACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTACCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTCACACGCAACATTTTTTATCGTAATTTCTACGAT  >  NC_000913/2080964‑2081143

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