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	99.6%	718.4 / ‑4.9	224	T93A (ACC→GCC)	sbcB	exonuclease I
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (98/125); total (99/125)
Fisher's exact test for biased strand distribution p-value = 4.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.92e-01

Only 100 of 246 total aligned reads displayed.

GCCGTATTAATTACCGGTATTACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTACCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTCACACGCAACATTTTTTATCGTAATTTCTACGAT  >  NC_000913/2080964‑2081143

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

gCCGTATTAATTACCGGTATTACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTtgcc                                                                                         >  2:3306200/1‑93 (MQ=255)

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

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

     attaattaCCGGTATTACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTtgccgtacc                                                                                    <  1:6715086/93‑1 (MQ=255)

       taattaCCGGTATTACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGa                                                                                  <  1:6467873/93‑1 (MQ=255)

        aattaCCGGTATTACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGaa                                                                                 <  1:4859910/93‑1 (MQ=255)

        aattaCCGGTATTACACCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGaa                                                                                 <  2:6651504/93‑1 (MQ=255)

         attaCCGGTATTACCCCGCAGGAAGCACGGGCGTAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAg                                                                                <  2:4880119/93‑1 (MQ=255)

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

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

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

            aCCGGTATTACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGAcc                                                                             <  2:6595652/93‑1 (MQ=255)

            aCCGGTATTACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGAcc                                                                             <  2:5792880/93‑1 (MQ=255)

             ccGGTATTACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCt                                                                            >  1:6263300/1‑93 (MQ=255)

              cGGTATTACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTg                                                                           <  1:5967126/93‑1 (MQ=255)

                  aTTACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTAtt                                                                       <  1:1671288/93‑1 (MQ=255)

                  aTTACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTAtt                                                                       >  1:3653178/1‑93 (MQ=255)

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

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

                   ttACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACAGAAGACCTGTATTc                                                                      >  1:1132760/1‑93 (MQ=255)

                    tACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCt                                                                     >  2:3007630/1‑93 (MQ=255)

                     aCCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTg                                                                    <  2:1085853/93‑1 (MQ=255)

                      ccccGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTgg                                                                   <  1:324607/93‑1 (MQ=255)

                       cccGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTggg                                                                  >  2:3622350/1‑93 (MQ=255)

                       cccGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATCCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTggg                                                                  <  1:1382494/93‑1 (MQ=255)

                         cGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCt                                                                >  2:4889567/1‑93 (MQ=255)

                         cGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCt                                                                <  2:484231/93‑1 (MQ=255)

                          gCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTa                                                               >  1:4930910/1‑93 (MQ=255)

                            aGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTaca                                                             >  2:1256218/1‑93 (MQ=255)

                             ggAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTGACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTacaa                                                            >  2:3810199/1‑93 (MQ=255)

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

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

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

                               aaGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTCGGCTacaaca                                                          >  1:1289680/1‑93 (MQ=255)

                                 gCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAAt                                                        <  2:1443629/93‑1 (MQ=255)

                                 gCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAAt                                                        <  1:5358627/93‑1 (MQ=255)

                                 gCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAAt                                                        >  2:2805715/1‑93 (MQ=255)

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

                                    cGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAAtgtg                                                     <  2:3183776/93‑1 (MQ=255)

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

                                      ggCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGACGACCTGTATTCTGGGCTACAACAATGTGCg                                                   >  2:5171017/1‑93 (MQ=255)

                                      ggCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCg                                                   >  2:3371028/1‑93 (MQ=255)

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

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

                                        cGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCAGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGtt                                                 <  1:801332/93‑1 (MQ=255)

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

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

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

                                          aaaGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTc                                               <  1:2378009/93‑1 (MQ=255)

                                          aaaGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTc                                               <  1:4694604/93‑1 (MQ=255)

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

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

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

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

                                                    aCGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAg                                     >  2:3248235/1‑93 (MQ=255)

                                                    aCGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAg                                     <  2:2654656/93‑1 (MQ=255)

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

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

                                                       aaGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGTCGAAGTca                                  <  1:6526896/93‑1 (MQ=255)

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

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

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

                                                         gCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTcaca                                <  1:3916928/93‑1 (MQ=255)

                                                              gTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTCACACGCaa                           >  1:5609382/1‑93 (MQ=255)

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

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

                                                                  gccgccCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTCACACGCAACAtt                       <  2:5149145/93‑1 (MQ=255)

                                                                    cgccCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTCACACGCAACAtttt                     <  1:6116913/93‑1 (MQ=255)

                                                                      cccGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTCACACGCAACAtttttt                   <  1:568202/93‑1 (MQ=255)

                                                                       ccGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTCACACGCAACATTTTTTa                  >  2:1006893/1‑93 (MQ=255)

                                                                       ccGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTCACACGCAACATTTTTTa                  >  1:2916620/1‑93 (MQ=255)

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

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

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

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

                                                                          tATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTCACACGCAACATTTTTTATCg               <  2:1393962/93‑1 (MQ=255)

                                                                          tATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTCACACGCAACATTTTTTATCg               >  1:221827/1‑93 (MQ=255)

                                                                          tATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTCACACGCAACATTTTTTATCg               <  1:2634620/93‑1 (MQ=255)

                                                                          tATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTCACACGCAACATTTTTTATCg               <  1:5522703/93‑1 (MQ=255)

                                                                           aTTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTCACACGCAACATTTTTTATCGt              >  2:3555075/1‑93 (MQ=255)

                                                                           aTTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTCACACGCAACATTTTTTATCGt              <  1:2786704/93‑1 (MQ=255)

                                                                             tCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTCACACGCAACATTTTTTATCGTaa            <  2:2000267/93‑1 (MQ=255)

                                                                              cACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTCACACGCAACATTTTTTATCGTAAt           >  2:93913/1‑93 (MQ=255)

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

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

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

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

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

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

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

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

                                                                                  cGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTCACACGCAACATTTTTTATCGTAATTTCt       >  2:2174557/1‑93 (MQ=255)

                                                                                   gCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTCACACGCAACATTTTTTATCGTAATTTCTa      >  2:2414165/1‑93 (MQ=255)

                                                                                    cTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTCACACGCAACATTTTTTATCGTAATTTCTAc     <  1:6017670/93‑1 (MQ=255)

                                                                                    cTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTCACACGCAACATTTTTTATCGTAATTTCTAc     <  2:6438738/93‑1 (MQ=255)

                                                                                     tttttGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTCACACGCAACATTTTTTATCGTAATTTCTACg    >  2:6181083/1‑93 (MQ=255)

                                                                                     tttttGCCGTACCGAAGACATGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTCACACGCAACATTTTTTATCGTAATTTCTACg    <  2:3996179/93‑1 (MQ=255)

                                                                                      ttttGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTCACACGCAACATTTTTTATCGTAATTTCTACGa   <  1:1509238/93‑1 (MQ=255)

                                                                                       tttGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTCACACGCAACATTTTTTATCGTAATTTCTACGAt  <  1:6687799/93‑1 (MQ=255)

GCCGTATTAATTACCGGTATTACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTACCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTCACACGCAACATTTTTTATCGTAATTTCTACGAT  >  NC_000913/2080964‑2081143

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