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.7%	928.4 / ‑5.3	290	T93A (ACC→GCC)	sbcB	exonuclease I
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (148/141); total (149/141)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.03e-01

Only 100 of 334 total aligned reads displayed.

GCCGTATTAATTACCGGTATTACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTACCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTCACACGCAACATTTTTTATCGTAATTTCTACGAT  >  NC_000913/2080964‑2081143

gCCGTATTAATTACCGGTATTACCCCGCAGGAAGCATGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTtgcc                                                                                         <  1:6416560/93‑1 (MQ=255)

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

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

  cGTATTAATTCCCGGTATTACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTtgccgt                                                                                       <  2:614270/93‑1 (MQ=255)

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

   gTATTAATTTCCGGTATTTCCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTtgccgta                                                                                      <  2:8271937/93‑1 (MQ=255)

   gTATTAATTACCGGTATTCCCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTtgccgta                                                                                      <  1:4773767/93‑1 (MQ=255)

   gTATTAATTACCGGTATTACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTtgccgta                                                                                      >  2:9398/1‑93 (MQ=255)

   gTATTAATTACCGGTATTACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTtgccgta                                                                                      >  1:459162/1‑93 (MQ=255)

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

    tATTAATTACCGGTATTACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTtgccgtac                                                                                     <  1:3752640/93‑1 (MQ=255)

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

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

     attaattaCCGGTATTACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTtaccgtacc                                                                                    >  2:6590802/1‑93 (MQ=255)

      ttaattaCCGGTATTACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTtgccgtaccg                                                                                   >  2:6519052/1‑93 (MQ=255)

      ttaattaCCGGTATTACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTtgccgtaccg                                                                                   >  2:2319723/1‑93 (MQ=255)

       taattaCCGGTATTACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGa                                                                                  >  2:3453967/1‑93 (MQ=255)

       taattaCCGGTATTACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGa                                                                                  >  1:5113681/1‑93 (MQ=255)

       taattaCCGGTATTACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCAGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGa                                                                                  >  1:2463806/1‑93 (MQ=255)

            aCCGGTATTACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGAcc                                                                             >  2:3215498/1‑93 (MQ=255)

             ccGGTATTACCCCCCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCCCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCt                                                                            <  2:2160098/93‑1 (MQ=255)

                 tATTACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCTCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTAt                                                                        <  2:7669897/93‑1 (MQ=255)

                 tATTACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTAt                                                                        >  2:6463244/1‑93 (MQ=255)

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

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

                    tACCCCGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCt                                                                     <  1:9820047/93‑1 (MQ=255)

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

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

                       cccGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTggg                                                                  <  2:8657871/93‑1 (MQ=255)

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

                        ccGCAGGAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGc                                                                 >  2:384196/1‑93 (MQ=255)

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

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

                              gAAGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGTCTacaac                                                           >  2:9475758/1‑93 (MQ=255)

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

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

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

                                aGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTAcaacaa                                                         <  1:6891616/93‑1 (MQ=255)

                                aGCACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTAcaacaa                                                         <  1:2380997/93‑1 (MQ=255)

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

                                  cACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAAtg                                                       <  1:8388947/93‑1 (MQ=255)

                                  cACGGGCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAAtg                                                       >  2:4851080/1‑93 (MQ=255)

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

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

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

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

                                      ggCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCg                                                   <  1:80322/93‑1 (MQ=255)

                                       gCGAAAGGAGAAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGt                                                  <  1:901796/93‑1 (MQ=255)

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

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

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

                                              gagaAAACGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTccaca                                           >  1:5670114/1‑92 (MQ=255)

                                               agaAAACGAAGCCGCGTTTGCCGCCCGTATTCCCTCGTTTTTTGCCGTACCGAAGACCTGTATTCTGAGCTACAACAATGTGCGTTTcgacga                                          <  2:4129944/93‑1 (MQ=255)

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

                                                 aaaaCGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTcgacgacg                                        >  2:3060220/1‑93 (MQ=255)

                                                 aaaaCGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTcgacgacg                                        >  2:1175312/1‑93 (MQ=255)

                                                  aaaCGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTcgacgacga                                       <  2:5469591/93‑1 (MQ=255)

                                                   aaCGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGaa                                      >  2:4798421/1‑93 (MQ=255)

                                                   aaCGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGaa                                      >  1:2662196/1‑93 (MQ=255)

                                                    aCGAAGCCGCGTTTGCCGCCCGTATTCCCTCGTTTTTTCCCTTACCG‑AGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGt                                    <  1:4844773/93‑1 (MQ=255)

                                                    aCGAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAATACCTGTCTTCTGGGCTACAACAATGTGCGGTTCGATGACGaaa                                     >  1:2734070/1‑92 (MQ=255)

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

                                                      gAAGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTc                                   >  1:4081770/1‑93 (MQ=255)

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

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

                                                        aGCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTcac                                 <  1:1894233/93‑1 (MQ=255)

                                                         gCCGCGTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTcaca                                >  2:2658156/1‑93 (MQ=255)

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

                                                           cgcgTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTCacacc                              >  2:2069053/1‑92 (MQ=255)

                                                            gcgTTTGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTCACACGc                             >  1:2074664/1‑93 (MQ=255)

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

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

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

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

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

                                                                ttGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTCACACGCAACa                         <  1:6065634/93‑1 (MQ=255)

                                                                ttGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTCACACGCAACa                         <  1:7004472/93‑1 (MQ=255)

                                                                 tGCCGCCCGTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTCACACGCAACAt                        <  1:7723729/93‑1 (MQ=255)

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

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

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

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

                                                                         gTATTCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTCACACGCAACATTTTTTATc                <  2:312313/93‑1 (MQ=255)

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

                                                                           aTTCACTCGATTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGCCACACGCAAAATTTTTTATCGt              >  2:5613842/1‑93 (MQ=255)

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

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

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

                                                                            ttCACTCGCTTTTTGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGGTTCGACGACGAAGTCACACCCAACATTTTTTATCGTa             >  1:8432218/1‑93 (MQ=255)

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

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

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

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

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

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

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

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

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

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

                                                                                      ttttGCCGTACCGAAGACCTGTATTCTGGGCTACAACAATGTGCGTTTCGACGACGAAGTCACACGCAACATTTTTTATCGTAATTTCTACGa   <  1:1045891/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` `< 28 ≤` `ATCG/ATCG` `< 33 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`