BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	332,019	T→C	V142A (GTG→GCG)	yahA →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	332,019	0	T	C	99.7%	1174.2 / ‑5.1	367	V142A (GTG→GCG)	yahA	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (161/205); total (161/206)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.22e-01

Only 100 of 386 total aligned reads displayed.

CCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAA  >  NC_000913/331929‑332108

ccATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGgcgc                                                                                         <  1:6932361/93‑1 (MQ=255)

 cATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCt                                                                                        <  1:6161801/93‑1 (MQ=255)

  aTCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCtt                                                                                       >  1:9567201/1‑93 (MQ=255)

   tCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGCGGCGCTTg                                                                                      >  2:9362165/1‑93 (MQ=255)

   tCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTg                                                                                      >  1:9812946/1‑93 (MQ=255)

   tCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTg                                                                                      <  1:7961193/93‑1 (MQ=255)

     aGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTc                                                                                    <  1:327780/93‑1 (MQ=255)

      gTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGttc                                                                                   >  1:570506/1‑91 (MQ=255)

      gTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTcc                                                                                   <  2:9752909/93‑1 (MQ=255)

       tCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTcct                                                                                  >  1:4709258/1‑92 (MQ=255)

        cTGGCGTTGGAAAACCACGAATTCAATCCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGc                                                                                 <  1:3783690/93‑1 (MQ=255)

        cTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGc                                                                                 <  2:275890/93‑1 (MQ=255)

         tGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCt                                                                                >  2:7503207/1‑93 (MQ=255)

         tGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCt                                                                                <  1:4443091/93‑1 (MQ=255)

           gCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTgg                                                                              <  2:927442/93‑1 (MQ=255)

            cGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTggg                                                                             <  2:794879/93‑1 (MQ=255)

             gTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGa                                                                            <  1:9270177/93‑1 (MQ=255)

              ttGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGaa                                                                           <  2:3110013/93‑1 (MQ=255)

              ttGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGaa                                                                           >  1:7623761/1‑93 (MQ=255)

              ttGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGaa                                                                           <  1:5995251/93‑1 (MQ=255)

              gtGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGaa                                                                           <  1:7522384/92‑1 (MQ=255)

                ggAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACa                                                                         <  1:4043030/93‑1 (MQ=255)

                ggAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACa                                                                         <  2:8913749/93‑1 (MQ=255)

                       cACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACaa                                                                  <  2:7409638/93‑1 (MQ=255)

                        aCGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACaaa                                                                 <  2:8572862/93‑1 (MQ=255)

                        aCGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGT‑CGCTGGGAACATCCACAAAc                                                                <  1:3683755/93‑1 (MQ=255)

                         cGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAAc                                                                <  1:8454543/93‑1 (MQ=255)

                         cGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAAc                                                                <  1:3584590/93‑1 (MQ=255)

                          gAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACg                                                               <  1:6652319/93‑1 (MQ=255)

                           aaTTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTCGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCCCAAACgg                                                              >  1:888002/1‑93 (MQ=255)

                             ttCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGa                                                            >  1:2437384/1‑93 (MQ=255)

                               cAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAAt                                                          <  1:2691475/93‑1 (MQ=255)

                                 aaCCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTa                                                        <  2:2229249/93‑1 (MQ=255)

                                 aaCCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTa                                                        >  1:4150013/1‑93 (MQ=255)

                                   ccGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATc                                                      >  2:3332749/1‑93 (MQ=255)

                                   ccGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATc                                                      <  1:1789227/93‑1 (MQ=255)

                                    cGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAAGATCCACAAACGGGAATTAtac                                                     >  2:5292239/1‑91 (MQ=255)

                                    cGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATcc                                                     <  2:4624987/93‑1 (MQ=255)

                                     gTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATccc                                                    >  1:2257740/1‑93 (MQ=255)

                                      tGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATcccc                                                   >  2:6311874/1‑92 (MQ=255)

                                       gcaTCCAACCGGGTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCAc                                                  <  2:8548981/91‑1 (MQ=255)

                                        gATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCAcc                                                 <  2:1496635/93‑1 (MQ=255)

                                          tCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCgg                                               >  1:9985403/1‑93 (MQ=255)

                                            cAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCgggt                                             >  2:9092795/1‑91 (MQ=255)

                                            cAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGAt                                             >  1:9455756/1‑93 (MQ=255)

                                              aCCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCa                                           <  2:2420592/93‑1 (MQ=255)

                                              aCCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCa                                           <  1:1162234/93‑1 (MQ=255)

                                              aCCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCa                                           >  2:706891/1‑93 (MQ=255)

                                               ccGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTTCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAg                                          >  1:3040909/1‑93 (MQ=255)

                                                 ggTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGtt                                        <  2:6624487/93‑1 (MQ=255)

                                                  gTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTTCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGttt                                       <  1:2129445/93‑1 (MQ=255)

                                                  gTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGttt                                       >  2:4018363/1‑93 (MQ=255)

                                                  gTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGttt                                       >  1:1315619/1‑93 (MQ=255)

                                                   ttttCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTa                                      <  2:5921419/93‑1 (MQ=255)

                                                    tttCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTAt                                     <  2:5217598/93‑1 (MQ=255)

                                                       cTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTcc                                  <  2:7975729/93‑1 (MQ=255)

                                                       cTGCGCGCAGACTGGCGTACTGACGGGCTGTGACGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTcc                                  <  2:8328790/93‑1 (MQ=255)

                                                        tGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCct                                 >  1:5819583/1‑93 (MQ=255)

                                                        tGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCct                                 <  1:1070733/93‑1 (MQ=255)

                                                         gcgcgcAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCctc                                >  1:5331326/1‑93 (MQ=255)

                                                         gcgcgcAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCctc                                <  1:8693782/93‑1 (MQ=255)

                                                          cgcgcAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCctct                               <  2:4542966/93‑1 (MQ=255)

                                                          cgcgcAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCctct                               <  2:3900369/93‑1 (MQ=255)

                                                            cgcAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCGCCGGATCAGTTTATTCCTCTgg                             <  2:6277293/93‑1 (MQ=255)

                                                             gcAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGc                            >  2:2278646/1‑93 (MQ=255)

                                                             gcAGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGc                            <  2:4386134/93‑1 (MQ=255)

                                                               aGACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCgg                          <  1:454662/93‑1 (MQ=255)

                                                                gACTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGa                         >  2:7218651/1‑93 (MQ=255)

                                                                 aCTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAg                        <  1:90633/93‑1 (MQ=255)

                                                                  cTGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGt                       >  1:3887469/1‑93 (MQ=255)

                                                                   tGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTc                      <  2:5727712/93‑1 (MQ=255)

                                                                   tGGCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTc                      >  2:8077777/1‑93 (MQ=255)

                                                                     gCGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCAt                    >  2:501274/1‑93 (MQ=255)

                                                                      cGTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATc                   <  1:9323203/93‑1 (MQ=255)

                                                                       gTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATcc                  <  2:651794/93‑1 (MQ=255)

                                                                       gTACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATcc                  <  1:1903242/93‑1 (MQ=255)

                                                                        tACTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCg                 >  1:8589295/1‑93 (MQ=255)

                                                                         aCTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCgg                >  1:8350269/1‑93 (MQ=255)

                                                                          cTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGt               <  2:1592278/93‑1 (MQ=255)

                                                                          cTGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGt               >  2:873956/1‑93 (MQ=255)

                                                                           tGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTc              >  1:6203344/1‑93 (MQ=255)

                                                                           tGACGGGCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTc              <  1:5005600/93‑1 (MQ=255)

                                                                               gggCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTAt          >  2:5395113/1‑93 (MQ=255)

                                                                               gggCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTAt          <  1:1595008/93‑1 (MQ=255)

                                                                               gggCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTAt          >  1:4232995/1‑93 (MQ=255)

                                                                                 gCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCATCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTg        <  2:1832074/93‑1 (MQ=255)

                                                                                 gCTGTGAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTg        >  2:2808213/1‑93 (MQ=255)

                                                                                   tgtgAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTc      <  2:556669/93‑1 (MQ=255)

                                                                                   tgtgAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTc      >  1:7465007/1‑93 (MQ=255)

                                                                                    gtgAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCa     <  2:8195577/93‑1 (MQ=255)

                                                                                    gtgAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCa     >  1:7377695/1‑93 (MQ=255)

                                                                                    gtgAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCa     >  2:1655923/1‑93 (MQ=255)

                                                                                     tgAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCAt    <  2:3567151/93‑1 (MQ=255)

                                                                                     tgAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCAt    <  2:9853921/93‑1 (MQ=255)

                                                                                      gAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATa   >  1:8103774/1‑93 (MQ=255)

                                                                                      gAGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATa   >  1:6359103/1‑93 (MQ=255)

                                                                                       aGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATaa  <  2:256680/93‑1 (MQ=255)

                                                                                       aGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATaa  <  2:2992291/93‑1 (MQ=255)

                                                                                       aGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATaa  <  2:6768264/93‑1 (MQ=255)

                                                                                       aGGCGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATaa  >  2:9497912/1‑93 (MQ=255)

CCATCAGTCTGGCGTTGGAAAACCACGAATTCAAACCGTGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAA  >  NC_000913/331929‑332108

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