BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,972,641	G→A	R145H (CGC→CAC)	wecD →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,972,641	0	G	A	99.4%	1166.6 / ‑4.5	327	R145H (CGC→CAC)	wecD	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (1/1); new base A (135/190); total (136/191)
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.79e-01

Only 100 of 332 total aligned reads displayed.

AATGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATCGTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGAATGTAG  >  NC_000913/3972553‑3972724

aaTGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGt                                                                                 >  1:5138336/1‑93 (MQ=255)

aaTGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGt                                                                                 >  2:1741553/1‑93 (MQ=255)

aaTGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGt                                                                                 >  2:2684884/1‑93 (MQ=255)

    cGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTaaaa                                                                             <  1:3268686/93‑1 (MQ=255)

     gACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAaac                                                                            <  1:3894576/93‑1 (MQ=255)

     gACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAaac                                                                            <  2:311421/93‑1 (MQ=255)

       cAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAaacaa                                                                          >  1:4892190/1‑93 (MQ=255)

        aGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAaacaac                                                                         <  1:1254043/93‑1 (MQ=255)

        aGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAaacaac                                                                         >  1:2833105/1‑93 (MQ=255)

         taTGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACt                                                                        >  2:293241/2‑93 (MQ=255)

         gATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACt                                                                        >  1:445074/1‑93 (MQ=255)

         gATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACt                                                                        <  2:7336844/93‑1 (MQ=255)

          aTGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACtt                                                                       <  1:6168017/93‑1 (MQ=255)

          aTGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACtt                                                                       >  2:6435908/1‑93 (MQ=255)

           tGCGCGAATTGGCCTGCTGGCTGGACGCGGTTCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACttt                                                                      <  1:2536275/93‑1 (MQ=255)

           tGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTGTGCTAACGGTAAAACAGCtta                                                                      <  2:7027850/93‑2 (MQ=255)

           tGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACttt                                                                      <  2:232566/93‑1 (MQ=255)

           tGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACttt                                                                      <  1:6625073/93‑1 (MQ=255)

           tGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACttt                                                                      <  2:5587729/93‑1 (MQ=255)

           tGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACttt                                                                      >  2:1691192/1‑93 (MQ=255)

            gcgcGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTg                                                                     <  2:4037819/93‑1 (MQ=255)

            gcgcGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTg                                                                     <  2:4924519/93‑1 (MQ=255)

              gcgAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCg                                                                   <  2:725346/93‑1 (MQ=255)

              gcgAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCg                                                                   <  2:1878063/93‑1 (MQ=255)

              gcgAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCg                                                                   <  1:1165992/93‑1 (MQ=255)

                gAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCggg                                                                 >  2:3917643/1‑93 (MQ=255)

                 aaTTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCgggg                                                                >  2:4193743/1‑92 (MQ=255)

                 aaTTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGt                                                                <  1:6566948/93‑1 (MQ=255)

                 aaTTGGCCTGCTGGCTGGACGCGGTGCAGGTGCCGAGCTGATGCAAACAGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTATGCgggc                                                                >  2:7861340/1‑92 (MQ=255)

                   ttGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTgg                                                              <  1:7078763/93‑1 (MQ=255)

                      gCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGa                                                           <  1:731874/93‑1 (MQ=255)

                         tgctgGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGAccc                                                        >  1:7504369/1‑93 (MQ=255)

                         tgctgGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGAGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGAccc                                                        >  1:4365977/1‑93 (MQ=255)

                            tggctggACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCaaa                                                     <  1:4777918/93‑1 (MQ=255)

                             ggctggACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAAt                                                    <  1:5643412/93‑1 (MQ=255)

                             ggctggACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAAt                                                    <  2:3483247/93‑1 (MQ=255)

                              gctggACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATg                                                   >  1:1930908/1‑93 (MQ=255)

                              gctggACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATg                                                   >  2:2412520/1‑93 (MQ=255)

                               ctggACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATgg                                                  <  2:6683167/93‑1 (MQ=255)

                                 ggACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGc                                                <  2:7527434/93‑1 (MQ=255)

                                   aCGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCaa                                              <  1:405021/93‑1 (MQ=255)

                                   aCGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCaa                                              >  1:4509569/1‑93 (MQ=255)

                                       ggTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACAcc                                          <  2:1326064/93‑1 (MQ=255)

                                       ggTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACAcc                                          <  2:5321516/93‑1 (MQ=255)

                                       ggTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACAcc                                          <  1:3360573/93‑1 (MQ=255)

                                       ggTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACAcc                                          >  2:6270116/1‑93 (MQ=255)

                                        gTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACAccg                                         >  2:859190/1‑93 (MQ=255)

                                         tGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACAccgc                                        <  2:6065452/93‑1 (MQ=255)

                                         tGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACAccgc                                        >  2:479468/1‑93 (MQ=255)

                                             ggTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCcgcg                                    >  1:6269419/1‑93 (MQ=255)

                                              gTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCgcgc                                   <  2:2037890/93‑1 (MQ=255)

                                               tGCGGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCt                                  <  1:2175395/93‑1 (MQ=255)

                                                gctgagctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCtt                                 <  1:825347/93‑1 (MQ=255)

                                                 ctgagctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTa                                <  1:5878360/93‑1 (MQ=255)

                                                 ctgagctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTa                                <  2:6957012/93‑1 (MQ=255)

                                                 ctgagctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTa                                >  2:3642974/1‑93 (MQ=255)

                                                  tgagctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTaa                               <  1:3119683/93‑1 (MQ=255)

                                                   gagctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTaaa                              >  2:3295175/1‑93 (MQ=255)

                                                     gctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACg                            <  2:5784044/93‑1 (MQ=255)

                                                      ctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAATGa                           <  1:4090702/93‑1 (MQ=255)

                                                      ctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGa                           >  1:6673486/1‑93 (MQ=255)

                                                       tgaTGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGat                          >  2:6069442/1‑93 (MQ=255)

                                                        gaTGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGata                         >  1:5365215/1‑93 (MQ=255)

                                                         aTGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGatac                        >  1:7809967/1‑93 (MQ=255)

                                                          tGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGataca                       <  2:2579100/93‑1 (MQ=255)

                                                             aaaCGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGatacatac                    >  2:5970804/1‑93 (MQ=255)

                                                             aaaCGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGatacatac                    >  2:4583233/1‑93 (MQ=255)

                                                              aaCGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGAtacataca                   >  1:5063889/1‑93 (MQ=255)

                                                              aaCGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGAtacataca                   <  2:1422078/93‑1 (MQ=255)

                                                              aaCGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGAtacataca                   >  1:167897/1‑93 (MQ=255)

                                                               aCGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACaa                  >  2:7876549/1‑93 (MQ=255)

                                                                cGGGGCTAACCTGCGCGTATGCTCACGGTAAACCAACTTTGCGGTTGGCGACCCAAATGGGCAGCACCGCCGCGCTTAAACGATACATGCaaa                 <  2:3209791/93‑1 (MQ=255)

                                                                cGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACaaa                 <  2:3406845/93‑1 (MQ=255)

                                                                cGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACaaa                 <  1:2328554/93‑1 (MQ=255)

                                                                 ggCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAg                <  1:6949283/93‑1 (MQ=255)

                                                                 ggCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAg                <  1:2715937/93‑1 (MQ=255)

                                                                  gcgcTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAgt               <  2:5482845/93‑1 (MQ=255)

                                                                  gcgcTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAgt               <  1:1442268/93‑1 (MQ=255)

                                                                   cgcTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAgtg              <  1:4658514/93‑1 (MQ=255)

                                                                    gcTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAgtgg             >  2:4608181/1‑93 (MQ=255)

                                                                    gcTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAgtgg             >  2:5138298/1‑93 (MQ=255)

                                                                    gcTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAgtgg             >  2:1209459/1‑93 (MQ=255)

                                                                    gcTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAgtgg             <  2:4409456/93‑1 (MQ=255)

                                                                     cTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAgtggt            <  1:913863/93‑1 (MQ=255)

                                                                     cTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAgtggt            <  2:7652685/93‑1 (MQ=255)

                                                                     cTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAgtggt            >  1:7576441/1‑93 (MQ=255)

                                                                     cTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAgtggt            >  1:6772886/1‑93 (MQ=255)

                                                                     cTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAgtggt            <  2:6821444/93‑1 (MQ=255)

                                                                      tAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAgtggtg           <  2:2870671/93‑1 (MQ=255)

                                                                      tAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAgtggtg           <  2:7822691/93‑1 (MQ=255)

                                                                       aaaCTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGc          <  2:3828653/93‑1 (MQ=255)

                                                                       aaaCTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGc          <  2:7429618/93‑1 (MQ=255)

                                                                          cTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAAGGATACATACAAAGTGGTGCGaa       <  2:4980300/93‑1 (MQ=255)

                                                                          cTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGaa       >  1:5204063/1‑93 (MQ=255)

                                                                          cTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGaa       <  1:2886161/93‑1 (MQ=255)

                                                                           tGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGAAt      <  1:6483223/93‑1 (MQ=255)

                                                                            gggCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGAATg     <  2:4326178/93‑1 (MQ=255)

                                                                            gggCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGAATg     <  1:5952718/93‑1 (MQ=255)

                                                                            gggCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGAATg     <  1:2081804/93‑1 (MQ=255)

                                                                             ggCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGAATGt    <  1:2786036/93‑1 (MQ=255)

AATGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATCGTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGAATGTAG  >  NC_000913/3972553‑3972724

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: `‑`