BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,032,327	A→G	Q44R (CAG→CGG)	hemG →	protoporphyrin oxidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,032,327	0	A	G	95.2%	431.4 / 8.8	147	Q44R (CAG→CGG)	hemG	protoporphyrin oxidase
Reads supporting (aligned to +/- strand): ref base A (3/4); new base G (76/64); total (79/68)
Fisher's exact test for biased strand distribution p-value = 7.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.76e-01

Only 100 of 149 total aligned reads displayed.

CGCGCGAGATTGCCTCCTACCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACAGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGCGTTCCAGGAATTTGTCAAAAAACATG  >  NC_000913/4032237‑4032416

cgcgcgAGATTGCCTCCTACCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGt                                                                                         <  1:882645/93‑1 (MQ=255)

   gcgAGATTGCCTCCTACCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTggg                                                                                      >  1:1039458/1‑93 (MQ=255)

   gcgAGATTGCCTCCTACCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTggg                                                                                      <  2:4312729/93‑1 (MQ=255)

     gagaTTGCCTCCTACCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGaa                                                                                    <  1:4319325/93‑1 (MQ=255)

     gagaTTGCCTCCTACCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGaa                                                                                    >  2:2143017/1‑93 (MQ=255)

          tGCCTCCTACCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTa                                                                               <  1:3467665/93‑1 (MQ=255)

          tGCCTCCTACCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTa                                                                               <  2:189547/93‑1 (MQ=255)

             ctcctACCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGTATTGAAGAACCACGGTGGGAAAACTATGa                                                                            <  1:263090/93‑1 (MQ=255)

             ctcctACCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGa                                                                            >  1:3800891/1‑93 (MQ=255)

             ctcctACCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGa                                                                            <  1:280201/93‑1 (MQ=255)

             ctcctACCTGGCTTCGGAACTGAAAGAACTGCCGGTCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGa                                                                            >  2:3748233/1‑93 (MQ=255)

              tcctACCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGAc                                                                           >  2:3946212/1‑93 (MQ=255)

               cctACCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGAcc                                                                          >  1:4055213/1‑93 (MQ=255)

                ctACCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCg                                                                         >  2:3551223/1‑93 (MQ=255)

                  aCCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCgtg                                                                       <  2:4540106/93‑1 (MQ=255)

                  aCCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCgtg                                                                       >  2:2601134/1‑93 (MQ=255)

                    cTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGtgtg                                                                     >  1:3700084/1‑93 (MQ=255)

                    cTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGtgtg                                                                     >  1:381017/1‑93 (MQ=255)

                    cTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGtgtg                                                                     >  2:2581882/1‑93 (MQ=255)

                       gCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTc                                                                  <  2:3763241/93‑1 (MQ=255)

                       gCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTc                                                                  >  1:183274/1‑93 (MQ=255)

                       gCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTc                                                                  <  1:1892785/93‑1 (MQ=255)

                       gCTTCGGAACTGAAAGAACTGGCGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTc                                                                  <  1:1892784/93‑1 (MQ=255)

                         ttCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCAt                                                                >  2:1371418/1‑93 (MQ=255)

                           cGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTg                                                              >  1:573744/1‑93 (MQ=255)

                            ggAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCACTGAAGAACCACGGTGGGAAAACTATGACCGTGGGGTCATTgg                                                             >  1:4197050/1‑93 (MQ=255)

                             gAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGt                                                            <  1:2565823/93‑1 (MQ=255)

                             gAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACAGTGGGAAAACTATGACCGTGTGGTCATTGGt                                                            <  1:3298125/93‑1 (MQ=255)

                               aCTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGc                                                          >  1:3671682/1‑93 (MQ=255)

                               aCTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGc                                                          <  2:4732581/93‑1 (MQ=255)

                                cTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCt                                                         <  1:438093/93‑1 (MQ=255)

                                cTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCt                                                         <  2:4693546/93‑1 (MQ=255)

                                 tGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCtt                                                        >  1:100891/1‑93 (MQ=255)

                                  gAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTc                                                       >  1:4334033/1‑93 (MQ=255)

                                   aaaGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCt                                                      <  1:762924/93‑1 (MQ=255)

                                   aaaGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCt                                                      <  1:3764431/93‑1 (MQ=255)

                                     aGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTAt                                                    <  1:2307529/93‑1 (MQ=255)

                                      gAACTGGGGATCCAGGCGGATGGCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTTTAtt                                                   <  2:1156712/93‑1 (MQ=255)

                                         cTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGc                                                <  2:163728/93‑1 (MQ=255)

                                          tGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCt                                               >  2:2641834/1‑93 (MQ=255)

                                           ggggATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTa                                              <  1:1464063/93‑1 (MQ=255)

                                           ggggATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACAGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTa                                              >  1:4278274/1‑93 (MQ=255)

                                            gggATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTTTGACCGTGTGGTCATTGGTGCTTCTATTCGCTAt                                             <  1:2606089/93‑1 (MQ=255)

                                            gggATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGCCCGTGTGGTCATTGGTGCTTCTATTCGCTAt                                             <  2:4719606/93‑1 (MQ=255)

                                             ggATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATg                                            <  2:4582700/93‑1 (MQ=255)

                                             ggATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATg                                            >  1:3047352/1‑93 (MQ=255)

                                               aTCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGt                                          <  1:1444969/93‑1 (MQ=255)

                                                   aGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACt                                      >  2:3001551/1‑93 (MQ=255)

                                                    ggCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTa                                     >  2:4017275/1‑93 (MQ=255)

                                                    ggCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTa                                     >  1:979621/1‑93 (MQ=255)

                                                      cGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTAcc                                   <  2:2622491/93‑1 (MQ=255)

                                                       ggATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCa                                  >  2:1729824/1‑93 (MQ=255)

                                                       ggATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCa                                  <  2:1604097/93‑1 (MQ=255)

                                                       ggATGTCGCCAATGTGCACCGCATTGAAGAACCACAGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCa                                  <  2:4013823/93‑1 (MQ=255)

                                                         aTGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTTTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCAtt                                <  2:2270325/93‑1 (MQ=255)

                                                         aTGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCAtt                                >  1:2205221/1‑93 (MQ=255)

                                                          tGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTc                               <  2:2829265/93‑1 (MQ=255)

                                                            tCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAg                             <  2:4399374/93‑1 (MQ=255)

                                                            tCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAg                             >  1:780554/1‑93 (MQ=255)

                                                            tCGCCAATGTGCACCGCATTGAAGAACCACAGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAg                             >  1:1557112/1‑93 (MQ=255)

                                                              gCCAATGTGCACCGCATTGAAGCACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAgtg                           >  2:661559/1‑91 (MQ=255)

                                                              gCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAgtg                           <  1:4327627/93‑3 (MQ=255)

                                                               ccAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAgtgt                          <  1:1194688/93‑1 (MQ=255)

                                                               ccAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAgtgt                          <  1:1072553/93‑1 (MQ=255)

                                                                 aaTGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTc                        >  1:3163477/1‑93 (MQ=255)

                                                                 aaTGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTc                        >  2:1558292/1‑93 (MQ=255)

                                                                 aaTGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTc                        >  2:3326790/1‑93 (MQ=255)

                                                                  aTGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTcc                       >  2:3829863/1‑93 (MQ=255)

                                                                   tgtgCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCa                      >  1:1752156/1‑93 (MQ=255)

                                                                   tgtgCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCa                      <  1:2626679/93‑1 (MQ=255)

                                                                   tgtgCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCa                      >  1:2255552/1‑93 (MQ=255)

                                                                   tgtgCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCa                      <  2:502986/93‑1 (MQ=255)

                                                                   tgtgCACCGCATTGAAGAACCACAGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGCGTTCCa                      >  2:3825127/1‑93 (MQ=255)

                                                                       cACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGaa                  <  2:724041/93‑1 (MQ=255)

                                                                        aCCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAAt                 <  1:665440/93‑1 (MQ=255)

                                                                        aCCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAAt                 >  1:3939494/1‑93 (MQ=255)

                                                                         ccGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGACtt                >  1:269571/1‑93 (MQ=255)

                                                                         ccGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAAtt                >  1:2531359/1‑93 (MQ=255)

                                                                          cccATTAAAGAACCACGGTGGGAAAACTATGACCGAGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTATCATTCCGTGTTGCAGGAAttt               >  2:3855866/3‑93 (MQ=255)

                                                                          cGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAAttt               >  1:4187684/1‑93 (MQ=255)

                                                                           gCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTg              >  2:3365922/1‑93 (MQ=255)

                                                                           gCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTg              <  1:1269509/93‑1 (MQ=255)

                                                                           gCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTg              >  1:2889574/1‑93 (MQ=255)

                                                                           gCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTg              >  2:2363407/1‑93 (MQ=255)

                                                                             aTTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTc            >  2:4730426/1‑93 (MQ=255)

                                                                             aTTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTc            >  2:1106569/1‑93 (MQ=255)

                                                                              ttGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCa           >  1:1503589/1‑93 (MQ=255)

                                                                              ttGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCa           >  2:2547797/1‑93 (MQ=255)

                                                                              ttGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCa           >  1:463125/1‑93 (MQ=255)

                                                                                gaagaaCCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCaaa         >  2:2566042/1‑93 (MQ=255)

                                                                                gaagaaCCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCaaa         <  2:3455483/93‑1 (MQ=255)

                                                                                gaagaaCCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCaaa         >  1:923618/1‑93 (MQ=255)

                                                                                  agaaCCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCaaaaa       <  2:64596/93‑1 (MQ=255)

                                                                                   gaaCCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCaaaaaa      <  1:2376565/93‑1 (MQ=255)

                                                                                    aaCCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAAc     >  1:419039/1‑93 (MQ=255)

                                                                                    aaCCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAAc     <  2:78518/93‑1 (MQ=255)

                                                                                     aCCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTCCCATTCAGTGTTCCAGGAATTTGTCAAAAAACa    >  1:4590540/1‑93 (MQ=255)

                                                                                     aCCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACa    >  2:466705/1‑93 (MQ=255)

                                                                                     aCCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACa    >  1:2540796/1‑93 (MQ=255)

                                                                                      ccACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACAt   <  1:755094/93‑1 (MQ=255)

CGCGCGAGATTGCCTCCTACCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACAGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGCGTTCCAGGAATTTGTCAAAAAACATG  >  NC_000913/4032237‑4032416

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