BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	932,454	A→G	Q3R (CAG→CGG)	ftsK →	cell division protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	932,454	0	A	G	99.7%	926.0 / ‑5.2	284	Q3R (CAG→CGG)	ftsK	cell division protein
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (118/165); total (118/166)
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.09e-01

Only 100 of 297 total aligned reads displayed.

TACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGCTGATCCTTATTGTCCTGT  >  NC_000913/932364‑932543

tACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCggg                                                                                         >  2:7385164/1‑93 (MQ=255)

 aCAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGa                                                                                        <  2:7615210/93‑1 (MQ=255)

  cAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGaa                                                                                       >  1:6598152/1‑93 (MQ=255)

   aGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAAt                                                                                      >  1:6969778/1‑93 (MQ=255)

   aGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAAt                                                                                      <  1:7549634/93‑1 (MQ=255)

    gCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATa                                                                                     <  1:4579464/93‑1 (MQ=255)

    gCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATa                                                                                     >  1:4456458/1‑93 (MQ=255)

     cAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATAc                                                                                    <  2:333007/93‑1 (MQ=255)

      aaCAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACa                                                                                   <  1:8192675/93‑1 (MQ=255)

      aaCAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACa                                                                                   >  2:4424682/1‑93 (MQ=255)

      aaCAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACa                                                                                   >  1:3424159/1‑93 (MQ=255)

       aCAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACAt                                                                                  >  1:3529508/1‑93 (MQ=255)

       aCAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACAt                                                                                  <  2:5885678/93‑1 (MQ=255)

       aCAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACAt                                                                                  >  2:6675358/1‑93 (MQ=255)

       aCAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACAt                                                                                  >  2:1530155/1‑93 (MQ=255)

         aGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTg                                                                                >  1:9977474/1‑93 (MQ=255)

         aGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTg                                                                                <  1:5248768/93‑1 (MQ=255)

         aGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTg                                                                                <  1:1763057/93‑1 (MQ=255)

          gTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGa                                                                               <  2:4033819/93‑1 (MQ=255)

           tACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGaa                                                                              <  2:5527571/93‑1 (MQ=255)

           tACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGaa                                                                              >  1:1083684/1‑93 (MQ=255)

             cTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGa                                                                            <  2:7043283/93‑1 (MQ=255)

             cTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGa                                                                            <  1:6528734/93‑1 (MQ=255)

              tGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGAc                                                                           <  2:846563/93‑1 (MQ=255)

                   tAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAaaga                                                                      <  2:6779695/93‑1 (MQ=255)

                    aaCCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAaagaa                                                                     <  2:8682700/93‑1 (MQ=255)

                    aaCCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAaagaa                                                                     <  2:9324963/93‑1 (MQ=255)

                      tcgggTACGGTTGTCCGGTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACAATGTAGACAAAGAAGt                                                                   >  2:8981967/4‑93 (MQ=255)

                       cTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTc                                                                  >  1:1247168/1‑93 (MQ=255)

                       cTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTc                                                                  >  2:4732112/1‑93 (MQ=255)

                        tGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTca                                                                 <  2:5291691/93‑1 (MQ=255)

                        tGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTca                                                                 <  2:396918/93‑1 (MQ=255)

                         ggTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTcac                                                                >  1:604918/1‑93 (MQ=255)

                         ggTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTcac                                                                >  1:7903964/1‑93 (MQ=255)

                         ggTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTcac                                                                <  1:8715105/93‑1 (MQ=255)

                          gTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTcaca                                                               <  2:6942784/93‑1 (MQ=255)

                          gTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTcaca                                                               <  1:3841197/93‑1 (MQ=255)

                           tACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACAt                                                              <  2:8326753/93‑1 (MQ=255)

                            aCTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACAtt                                                             <  1:667455/93‑1 (MQ=255)

                             cTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTg                                                            <  2:4197495/93‑1 (MQ=255)

                             cTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTg                                                            >  1:5345491/1‑93 (MQ=255)

                              tgttgtCCGTTTTTGCATCGGGCAGGAAACGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGCATACATTGAAGACACAGAAGTCACATTGa                                                           >  2:9806684/1‑93 (MQ=255)

                              tgttgtCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGa                                                           >  1:5170051/1‑93 (MQ=255)

                                 tgtCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACaa                                                        >  1:5993064/1‑93 (MQ=255)

                                  gtCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACaaa                                                       >  1:4668258/1‑93 (MQ=255)

                                    tcGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGt                                                     <  1:9798272/92‑1 (MQ=255)

                                    ccGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGt                                                     <  2:5743329/93‑1 (MQ=255)

                                     cGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGtt                                                    <  2:2336422/93‑1 (MQ=255)

                                      gTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTa                                                   <  2:9263928/93‑1 (MQ=255)

                                      gTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTa                                                   <  1:6791479/93‑1 (MQ=255)

                                      gTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTa                                                   <  1:7417857/93‑1 (MQ=255)

                                      gTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTa                                                   <  2:9105679/93‑1 (MQ=255)

                                       ttttAGCATGGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTaa                                                  <  2:3529508/93‑1 (MQ=255)

                                       ttttAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTaa                                                  >  1:6231489/1‑93 (MQ=255)

                                       ttttAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTaa                                                  <  2:5969824/93‑1 (MQ=255)

                                       ttttAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTaa                                                  <  1:8100015/93‑1 (MQ=255)

                                         ttATCATCGGGCAGCAAAAGCCTGTACCCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGt                                                <  2:3407955/93‑1 (MQ=255)

                                               tCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGgc                                          >  2:1789149/1‑93 (MQ=255)

                                               tCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGgc                                          >  1:9654070/1‑93 (MQ=255)

                                               tCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGgc                                          <  2:2990739/93‑1 (MQ=255)

                                                cGGGCAGGAAAAGCCTGTATCCTGGAGAGCCTTTCTTGAGCCGTGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGgcc                                         <  1:1443435/93‑1 (MQ=255)

                                                cGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGgcc                                         <  1:3297071/93‑1 (MQ=255)

                                                  ggCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGgccgc                                       >  1:5786851/1‑93 (MQ=255)

                                                   gCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGgccgcc                                      <  2:5004983/93‑1 (MQ=255)

                                                    cAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGgccgccg                                     >  2:6617202/1‑93 (MQ=255)

                                                      ggAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGgccgccgcc                                   >  1:9450816/1‑93 (MQ=255)

                                                      ggAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGgccgccgcc                                   >  2:3649516/1‑93 (MQ=255)

                                                      ggAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGgccgccgcc                                   >  1:7594025/1‑93 (MQ=255)

                                                      ggAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGgccgccgcc                                   <  2:9576566/93‑1 (MQ=255)

                                                       gAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCt                                  <  2:4400078/93‑1 (MQ=255)

                                                       gAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCt                                  >  1:41731/1‑93 (MQ=255)

                                                        aaaaGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCtt                                 <  2:6071106/93‑1 (MQ=255)

                                                         aaaTCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTc                                <  1:2252932/93‑1 (MQ=255)

                                                          aaGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCt                               <  2:7514876/93‑1 (MQ=255)

                                                           aGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTg                              >  1:8297097/1‑93 (MQ=255)

                                                           aGCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTg                              <  2:8052044/93‑1 (MQ=255)

                                                            gCCTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTgg                             <  1:8056076/93‑1 (MQ=255)

                                                              cTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAg                           >  2:4691142/1‑93 (MQ=255)

                                                              cTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAg                           <  2:3509639/93‑1 (MQ=255)

                                                              cTGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAg                           >  2:9375996/1‑93 (MQ=255)

                                                               tGTAACCTGGAGAGCCTTTCTTGGGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAgg                          >  1:8497018/1‑93 (MQ=255)

                                                               tGTAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAgg                          >  1:9053366/1‑93 (MQ=255)

                                                                 tAACCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCg                        <  1:7244104/93‑1 (MQ=255)

                                                                  aaCCTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTCGCGGCCGCCGCCTTCTGGAGGCGt                       <  2:8901442/93‑1 (MQ=255)

                                                                    ccTGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTg                     <  2:1416738/93‑1 (MQ=255)

                                                                      tGGAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGCt                   >  1:134605/1‑93 (MQ=255)

                                                                       ggAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGCTg                  >  1:4933681/1‑93 (MQ=255)

                                                                       ggAGAGCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGCTg                  <  2:653859/93‑1 (MQ=255)

                                                                           agCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGCTGATcc              <  1:2950548/93‑1 (MQ=255)

                                                                            gCCTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGCTGATCCt             >  2:9154649/1‑93 (MQ=255)

                                                                            gCCTTTCTTGAGCCGGGAATACATGGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGCTGATCCt             <  1:4411698/93‑1 (MQ=255)

                                                                              cTTTCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGCTGATCCTTa           >  1:9520522/1‑93 (MQ=255)

                                                                               tttCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGCTGATCCTTAt          >  2:686725/1‑93 (MQ=255)

                                                                               tttCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGCTGATCCTTAt          <  1:6181093/93‑1 (MQ=255)

                                                                               tttCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGCTGATCCTTAt          >  1:9323027/1‑93 (MQ=255)

                                                                               tttCTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGCTGATCCTTAt          >  1:9168688/1‑93 (MQ=255)

                                                                                  cTTGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGCTGATCCTTATTGt       >  1:3094076/1‑93 (MQ=255)

                                                                                   ttGAGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGCTGATCCTTATTGTc      >  1:6813194/1‑93 (MQ=255)

                                                                                      aGCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGCTGATCCTTATTGTCCTg   <  1:7360967/93‑1 (MQ=255)

                                                                                       gCCGGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGCTGATCCTTATTGTCCTGt  <  1:275834/93‑1 (MQ=255)

TACAGCAACAGTACTGGGGTAACCTGGTACTGTTGTCCGTTTTAGCATCGGGCAGGAAAAGCCTGTAACCTGGAGAGCCTTTCTTGAGCCAGGAATACATTGAAGACAAAGAAGTCACATTGACAAAGTTAAGTAGCGGCCGCCGCCTTCTGGAGGCGTTGCTGATCCTTATTGTCCTGT  >  NC_000913/932364‑932543

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