BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,608,034	A→G	V224A (GTG→GCG)	uxaB ←	altronate oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,608,034	0	A	G	83.0%	604.6 / 109.7	270	V224A (GTG→GCG)	uxaB	altronate oxidoreductase
Reads supporting (aligned to +/- strand): ref base A (21/25); new base G (114/110); total (135/135)
Fisher's exact test for biased strand distribution p-value = 6.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 272 total aligned reads displayed.

GATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCATTGAATAAATGCTTCTGGC  >  NC_000913/1607944‑1608126

  tttCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTc                                                                                          >  1:4844199/1‑93 (MQ=255)

   ttCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCa                                                                                         <  1:966459/93‑1 (MQ=255)

    tCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCAt                                                                                        <  2:5664001/93‑1 (MQ=255)

    tCGGGCCCTGAATCACAAACAAGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCAt                                                                                        <  1:4334843/93‑1 (MQ=255)

      ggTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATcg                                                                                      <  1:6691684/93‑1 (MQ=255)

        tcgcTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATcgcg                                                                                    >  2:4736891/4‑93 (MQ=255)

        tCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATcgcg                                                                                    >  1:2587651/1‑93 (MQ=255)

         cccTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATcgcgc                                                                                   <  2:6028969/93‑1 (MQ=255)

          ccTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATcgcgcg                                                                                  <  2:6093120/93‑1 (MQ=255)

           cTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCgg                                                                                 >  2:639632/1‑93 (MQ=255)

           cTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCgg                                                                                 >  2:5081417/1‑93 (MQ=255)

                tCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAAc                                                                            <  2:3775706/93‑1 (MQ=255)

                tCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAAc                                                                            <  1:4419815/93‑1 (MQ=255)

                  acaAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCg                                                                          >  1:3484818/1‑93 (MQ=255)

                   caAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCgg                                                                         >  1:5293806/1‑93 (MQ=255)

                       cAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAAc                                                                     >  1:1659903/1‑93 (MQ=255)

                        aGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACg                                                                    >  2:3188081/1‑93 (MQ=255)

                         ggTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGa                                                                   <  2:5526977/93‑1 (MQ=255)

                         ggTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGa                                                                   <  1:5068244/93‑1 (MQ=255)

                         ggTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGa                                                                   <  1:3639555/93‑1 (MQ=255)

                         ggTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGa                                                                   >  2:7432302/1‑93 (MQ=255)

                         ggTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGa                                                                   >  2:3218841/1‑93 (MQ=255)

                           tAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATa                                                                 <  2:339697/93‑1 (MQ=255)

                           tAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATa                                                                 >  2:7063644/1‑93 (MQ=255)

                            aaaaGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATAc                                                                >  1:116929/1‑93 (MQ=255)

                             aaaGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACg                                                               <  2:3857288/93‑1 (MQ=255)

                              aaGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACgg                                                              <  1:1796949/93‑1 (MQ=255)

                              aaGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACgg                                                              <  1:1510478/93‑1 (MQ=255)

                               aGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGt                                                             <  2:2495086/93‑1 (MQ=255)

                                gtgtTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTc                                                            >  2:6823921/1‑93 (MQ=255)

                                  gtTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTcca                                                          >  1:6189098/1‑93 (MQ=255)

                                  gtTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTcca                                                          >  2:1292804/1‑93 (MQ=255)

                                  gtTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTcca                                                          >  1:5430600/1‑93 (MQ=255)

                                   ttCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTccac                                                         <  1:4249632/93‑1 (MQ=255)

                                   ttCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTccac                                                         >  1:274950/1‑93 (MQ=255)

                                   ttCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTccac                                                         <  2:5142960/93‑1 (MQ=255)

                                    tCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTccacc                                                        >  2:4927663/1‑93 (MQ=255)

                                     cAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTccacca                                                       <  1:343930/93‑1 (MQ=255)

                                     cAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTccacca                                                       <  2:2798251/93‑1 (MQ=255)

                                      aGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAg                                                      <  1:3816884/93‑1 (MQ=255)

                                       gCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATTCGGTCCACCAGc                                                     <  2:6627639/93‑1 (MQ=255)

                                          gtgtCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTa                                                  <  2:6791088/93‑1 (MQ=255)

                                          gtgtCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTa                                                  >  2:4374206/1‑93 (MQ=255)

                                           tgtCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAg                                                 >  1:1045816/1‑93 (MQ=255)

                                            gtCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGa                                                <  1:1739120/93‑1 (MQ=255)

                                             tCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGaa                                               >  1:5386867/1‑93 (MQ=255)

                                             tCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGaa                                               <  2:5768952/93‑1 (MQ=255)

                                             tCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGaa                                               >  1:3945965/1‑93 (MQ=255)

                                                agaAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAg                                            <  1:1922083/93‑1 (MQ=255)

                                                 gaAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGa                                           <  1:4517487/93‑1 (MQ=255)

                                                 gaAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGa                                           >  1:4471521/1‑93 (MQ=255)

                                                  aaaaCCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGaa                                          <  1:4042085/93‑1 (MQ=255)

                                                  aaaaCCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGaa                                          <  1:328429/93‑1 (MQ=255)

                                                   aaaCCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAg                                         >  2:2913925/1‑93 (MQ=255)

                                                    aaCCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGc                                        <  1:2029092/93‑1 (MQ=255)

                                                    aaCCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGc                                        <  2:627177/93‑1 (MQ=255)

                                                    aaCCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGc                                        <  2:3670748/93‑1 (MQ=255)

                                                      ccATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTg                                      >  1:7178605/1‑93 (MQ=255)

                                                      ccATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTg                                      <  2:1538210/93‑1 (MQ=255)

                                                      ccATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTg                                      >  2:7312701/1‑93 (MQ=255)

                                                         tCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTa                                   <  1:2299191/93‑1 (MQ=255)

                                                          cGTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAg                                  <  1:5937120/93‑1 (MQ=255)

                                                           gTGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGc                                 <  2:172987/93‑1 (MQ=255)

                                                            tGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGcc                                <  1:4171821/93‑1 (MQ=255)

                                                            tGATAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGcc                                >  1:5606508/1‑93 (MQ=255)

                                                            tGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGcc                                <  1:6836633/93‑1 (MQ=255)

                                                            tGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGcc                                <  1:7002276/93‑1 (MQ=255)

                                                               tAACCCAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGa                             >  1:3596571/1‑93 (MQ=255)

                                                                  cccAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATcc                          <  1:6777181/93‑1 (MQ=255)

                                                                   ccAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCa                         <  2:4731875/93‑1 (MQ=255)

                                                                    cAGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAg                        >  2:3461145/1‑93 (MQ=255)

                                                                     aGTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGc                       <  1:283478/93‑1 (MQ=255)

                                                                      gTTCTTCTTCCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGcc                      >  1:4664410/1‑93 (MQ=255)

                                                                      gTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGcc                      >  1:7379681/1‑93 (MQ=255)

                                                                       ttcttcttcCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCa                     >  2:6680711/1‑93 (MQ=255)

                                                                         cttcttcCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCAtt                   <  2:6868273/93‑1 (MQ=255)

                                                                         cttcttcCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCAtt                   <  1:4955985/93‑1 (MQ=255)

                                                                         cttcttcCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCAtt                   <  1:3177731/93‑1 (MQ=255)

                                                                          ttcttcCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCATTg                  >  2:5427941/1‑93 (MQ=255)

                                                                          ttcttcCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCATTg                  <  2:261572/93‑1 (MQ=255)

                                                                          ttcttcCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCATTg                  >  1:3223867/1‑93 (MQ=255)

                                                                           tcttcCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCATTGa                 >  1:2927356/1‑93 (MQ=255)

                                                                            cttcCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCATTGaa                <  2:6255161/93‑1 (MQ=255)

                                                                            cttcCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCATTGaa                <  1:5789277/93‑1 (MQ=255)

                                                                             ttcCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCATTGAAt               <  1:4649098/93‑1 (MQ=255)

                                                                             ttcCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCATTGAAt               >  1:3008968/1‑93 (MQ=255)

                                                                              tcCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCATTGAATa              >  1:6463477/1‑93 (MQ=255)

                                                                              tcCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCATTGAATa              <  2:4557685/93‑1 (MQ=255)

                                                                              tcCAGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCATTGAATa              >  2:1751165/1‑93 (MQ=255)

                                                                                 aGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCATTGAATAAAt           <  1:803255/93‑1 (MQ=255)

                                                                                 aGTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCATTGAATAAAt           >  1:814905/1‑93 (MQ=255)

                                                                                  gTTTAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCATTGAATAAATg          >  2:178607/1‑93 (MQ=255)

                                                                                   tttAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCATTGAATAAATGc         <  1:566942/93‑1 (MQ=255)

                                                                                     tAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCATTGAATAAATGCtt       <  2:4259965/93‑1 (MQ=255)

                                                                                     tAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCATTGAATAAATGCtt       <  2:416631/93‑1 (MQ=255)

                                                                                     tAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCATTGAATAAATGCtt       >  2:1398312/1‑93 (MQ=255)

                                                                                     tAGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCATTGAATAAATGCtt       >  2:2241217/1‑93 (MQ=255)

                                                                                      aGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCATTGAATAAATGCTTc      <  1:2386568/93‑1 (MQ=255)

                                                                                      aGCCGCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCATTGAATAAATGCTTc      <  1:6019815/93‑1 (MQ=255)

                                                                                          aCTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCATTGAATAAATGCTTCTGGc  >  1:2099517/1‑93 (MQ=255)

GATTTCGGTCCCTGAATCACAAACAGGTAAAAGTGTTCAGCGGTGTCGAGAAAACCATCGTGATAACCCAGTTCTTCTTCCAGTTTAGCCACTTCATCGCGCGGATAACCGGTAACGATACGGTCCACCAGCGTAGAACAGAAGCTGTTAGCCTGATCCAGCCATTGAATAAATGCTTCTGGC  >  NC_000913/1607944‑1608126

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