BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,032,068	A→C	Q303P (CAG→CCG)	pepQ →	proline dipeptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,032,068	0	A	C	97.4%	803.9 / 2.0	231	Q303P (CAG→CCG)	pepQ	proline dipeptidase
Reads supporting (aligned to +/- strand): ref base A (3/3); new base C (112/112); total (115/116)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 233 total aligned reads displayed.

CGCGCAGCTGGTGAAGGACGTTAATGATGAACAACTGGCGCTGATCGCCACCATGAAAGCAGGCGTCAGCTATGTGGATTACCACATCCAGTTCCATCAGCGCATCGCCAAACTGCTGCGTAAACATCAAATCATCACCGATATGAGTGAAGAAGCGATGGTCGAAAACGATCTTACCGGGCCGTTTATGCCGCAT  >  NC_000913/4031970‑4032165

  cgcAGCTGGTGAAGGACGTTAATGATGAACAACTGGCGCTGATCGCCACCATGAAAGCAGGCGTCAGCTATGTGGATTACCACATCCAGTTCCATCcgcg                                                                                                >  2:1439994/1‑100 (MQ=255)

   gcAGCTGGTGAAGGACGTTAATGATGAACAACTGGCGCTGATCGCCACCATGAAAGCAGGCGTCAGCTATGTGGATTACCACATCCAGTTCCATCCgcgc                                                                                               >  1:1792758/1‑100 (MQ=255)

    cAGCTGGTGAAGGACGTTAATGATGAACAACTGGCGCTGATCGCCACCATGAAAGCAGGCGTCAGCTATGTGGATTACCACATCCAGTTCCATCCGCGCa                                                                                              >  2:3737222/1‑100 (MQ=255)

      gCTGGTGAAGGACGTTAATGATGAACAACTGGCGCTGATCGCCACCATGAAAGCAGGCGTCAGCTATGTGGGTTCCCGTATCCAGTTCCATGCGCGCATc                                                                                            >  2:3168157/1‑100 (MQ=255)

      gCTGGTGAAGGACGTTAATGATGAACAACTGGCGCTGATCGCCACCATGAAAGCAGGCGTCAGCTATGTGGATTACCACATCCAGTTCCATCCGCGCATc                                                                                            <  2:4627161/100‑1 (MQ=255)

       cTGGTGAAGGACGTTAATGATGAACAACTGGCGCTGATCGCCACCATGAAAGCAGGCGTCAGCTATGTGGATTACCACATCCAGTTCCATCCGCGCATCg                                                                                           >  2:3443059/1‑100 (MQ=255)

         ggTGAAGGACGTTAATGATGAACAACTGGCGCTGATCGCCACCATGAAAGCAGGCGTCAGCTATGTGGATTACCACATCCAGTTCCATCCGCGCATCGcg                                                                                         >  2:4288025/1‑99 (MQ=255)

          ntGAAGGACGTTAATGATGAACAACTGGCGCTGATCGCCACCATGAAAGCAGGCGTCAGCTATGTGGATTACCACATCCAGTTCCATCCGCGCATCGCCa                                                                                        >  1:4448365/2‑100 (MQ=255)

            gAAGGACGTTAATGATGAACAACTGGCGCTGATCGCCACCATGAAAGCAGGCGTCAGCTATGTGGATTACCACATCCAGTTCCATCCGCGCATCGCCaaa                                                                                      >  1:1073281/1‑100 (MQ=255)

             aaGGACGTTAATGATGAACAACTGGCGCTGATCGCCACCATGAAAGCAGGCGTCAGCTATGTGGATTACCACATCCAGTTCCATCCGCGCATCGCCAAAc                                                                                     <  1:888011/100‑1 (MQ=255)

             aaGGACGTTAATGATGAACAACTGGCGCTGATCGCCACCATGAAAGCAGGCGTCAGCTATGTGGATTACCACATCCAGTTCCATCCGCGCATCGCCAAAc                                                                                     <  2:23163/100‑1 (MQ=255)

             aaGGACGTTAATGATGAACAACTGGCGCTGATCGCCACCATGAAAGCAGGCGTCAGCTATGTGGATTACCACATCCAGTTCCATCCGCGCATCGCCAAAc                                                                                     <  2:2663207/100‑1 (MQ=255)

               ggACGTTAATGATGAACAACTGGCGCTGATCGCCACCATGAAAGCAGGCGTCAGCTATGTGGATTACCACATCCAGTTCCATCCGCGCATCGCCAAActg                                                                                   >  2:4703495/1‑100 (MQ=255)

                 aCGTTAATGATGAACAACTGGCGCTGATCGCCACCATGAAAGCAGGCGTCAGCTATGTGGATTACCACATCCAGTTCCATCCGCGCATCGCCAAActgct                                                                                 >  2:2188102/1‑100 (MQ=255)

                    ttAATGATGAACAACTGGCGCTGATCGCCACCATGAAAGCAGGCGTCAGCTATGTGGATTACCACATCCAGTTCCATCCGCGCATCGCCAAACTGCTGCg                                                                              >  1:3685484/1‑100 (MQ=255)

                     tAATGATGAACAACTGGCGCTGATCGCCACCATGAAAGCAGGCGTCAGCTATGTGGATTACCACATCCAGTTCCATCCGCGCATCGCCAAACTGCTGCGt                                                                             >  1:1511876/1‑100 (MQ=255)

                        tgatgaACAACTGGCGCTGATCGCCACCATGAAAGCAGGCGTCAGCTATGTGGATTACCACATCCAGTTCCATCCGCGCATCGCCAAACTGCTGCGTaaa                                                                          >  1:4633645/1‑100 (MQ=255)

                        agatgaACAACTGGCGCTGATCGCCACCATGAAAGCAGGCGTCAGCTATGTGGATTACCACATCCAGTTCCATCCGCGCATCGCCAAACTGCTGCGTaaa                                                                          <  2:3547019/99‑1 (MQ=255)

                          atgaACCACTGGCGCTGATCGCCACCATGAAAGCAGGCGTCAGCTATGTAGATTACCACATCCAGTTCCATCCGCGCATCGCCAAACTGCTGCCTAAACa                                                                        >  2:5183001/1‑100 (MQ=255)

                          atgaACAACTGGCGCTGATCGCCACCATTAAAGCAGGCGTCAGCTATGTGGATTACCACATCCAGTTCCATCCGCGCATCGCCAAACTGCTGCGTAAACa                                                                        <  2:4864773/100‑1 (MQ=255)

                          atgaACAACTGGCGCTGATCGCCACCATGAAAGCAGGCGTCAGCTATGTGGATTACCACATCCAGTTCCATCCGCGCATCGCCAAACTGCTGCGTAAACa                                                                        <  2:4680646/100‑1 (MQ=255)

                          atgaACAACTGGCGCTGATCGCCACCATGAAAGCAGGCGTCAGCTATGTGGATTACCACATCCAGTTCCATCCGCGCATCGCCAAACTGCTGCGTAAACa                                                                        >  2:973419/1‑100 (MQ=255)

                           tgaACAACTGGCGCTGATCGCCACCATGAAAGCAGGCGTCAGCTATGTGGATTACCACATCCAGTTCCATCCGCGCATCGCCAAACTGCTGCGTAAACAt                                                                       <  1:5515461/100‑1 (MQ=255)

                           tgaACAACTGGCGCTGATCGCCACCATGAAAGCAGGCGTCAGCTATGTGGATTACCACATCCAGTTCCATCCGCGCATCGCCAAACTGCTGCGTAAACAt                                                                       >  1:1676087/1‑100 (MQ=255)

                            gaACAACTGGCGCTGATCGCCACCATGAAAGCAGGCGTCAGCTATGTGGATTACCACATCCAGTTCCATCCGCGCATCGCCAAACTGCTGCGTAAACATc                                                                      <  2:2723944/100‑1 (MQ=255)

                               caacTGGCGCTGATCGCCACCATGAAAGCAGGCGTCAGCTATGTGGATTACCACATCCAGTTCCATCCGCGCATCGCCAAACTGCTGCGTAAACATCaaa                                                                   <  1:4768748/100‑1 (MQ=255)

                               caacTGGCGCTGATCGCCACCATGAAAGCAGGCGTCAGCTATGTGGATTACCACATCCAGTTCCATCCGCGCATCGCCAAACTGCTGCGTAAACATCaaa                                                                   <  2:2275271/100‑1 (MQ=255)

                                 acTGGCGCTGATCGCCACCATGAAAGCAGGCGTCAGCTATGTGGATTACCACATCCAGTTCCATCCGCGCATCGCCAAACTGCTGCGTAAACATCAAatc                                                                 >  2:1959243/1‑100 (MQ=255)

                                  cTGGCGCTGATCGCCACCATGAAAGCAGGCGTCAGCTATGTGGATTACCACATCCAGTTCCATCCGCGCATCGCCAAACTGCTGCGTAAACATCAAatca                                                                <  1:2435039/100‑1 (MQ=255)

                                  cTGGCGCTGATCGCCACCATGAAAGCAGGCGTCAGCTATGTGGATTACCACATCCAGTTCCATCCGCGCATCGCCAAACTGCTGCGTAAACATCAAatca                                                                >  1:1478725/1‑100 (MQ=255)

                                    ggCGCTGATCGCCACCATGAAAGCAGGCGTCAGCTATGTGGATTACCACATCCAGTTCCATCCGCGCATCGCCAAACTGCTGCGTAAACATCAAatcatc                                                              >  2:2499122/1‑100 (MQ=255)

                                    ggCGCTGATCGCCACCATGAAAGCAGGCGTCAGCTATGTGGATTACCACATCCAGTTCCATCCGCGCATCGCCAAACTGCTGCGTAAACATCAAatcatc                                                              >  1:2800073/1‑100 (MQ=255)

                                    ggCGCTGATCGCCACCATGAAAGCAGGCGTCAGCTATGTGGATTACCACATCCAGTTCCATCCGCGCATCGCCAAACTGCTGCGTAAACATCAAatcatc                                                              <  1:2588233/100‑1 (MQ=255)

                                      cgcTGATCGCCACCATGAAAGCAGGCGTCAGCTATGTGGATTACCACATCCAGTTCCATCCGCGCATCGCCAAACTGCTGCGTAAACATCAAATCATCAc                                                            <  2:2802641/100‑1 (MQ=255)

                                      cgcTGATCGCCACCATGAAAGCAGGCGTCAGCTATGTGGATTACCACATCCAGTTCCATCCGCGCATCACCAAACTGCTGCGTAAACATCAAATCATCAc                                                            >  2:3761157/1‑100 (MQ=255)

                                       gcTGATCGCCACCATGAAAGCAGGCGTCAGCTATGTGGATTACCACATCCAGTTCCATCCGCGCATCGCCAAACTGCTGCGTAAACATCAAATCATCAcn                                                           <  1:1859781/100‑2 (MQ=255)

                                       gcTGATCGCCACCATGAAAGCAGGCGTCAGCTATGTGGATTACCACATCCAGTTCCATCCGCGCATCGCCAAACTGCTGCGTAAACATCAAATCATCAcc                                                           >  1:1242210/1‑100 (MQ=255)

                                          gATCGCCACCATGAAAGCAGGCGTCAGCTATGTGGATTACCACATCCAGTTCCATCCGCGCATCGCCAAACTGCTGCGTAAACATCAAATCATCACCGat                                                        <  2:852791/100‑1 (MQ=255)

                                          gATCGCCACCATGAAAGCAGGCGTCAGCTATGTGGATTACCACATCCAGTTCCATCCGCGCATCGCCAAACTGCTGCGTAAACATCAAATCATCACCGat                                                        >  2:3610571/1‑100 (MQ=255)

                                              gCCACCATGAAAGCAGGCGTCAGCTATGTGGATTACCACATCCAGTTCCATCCGCGCATCGCCAAACTGCTGCGTAAACATCAAATCATCACCGCTATGa                                                    >  1:5347446/1‑100 (MQ=255)

                                              gCCACCATGAAAGCAGGCGTCAGCTATGTGGATTACCACATCCAGTTCCATCCGCGCATCGCCAAACTGCTGCGTAAACATCAAATCATCACCGATATGa                                                    >  1:3939260/1‑100 (MQ=255)

                                               ccaccaTGAAAGCAGGCGTCAGCTATGTGGATTACCACATCCAGTTCCATCCGCGCATCGCCAAACTGCTGCGTAAACATCAAATCATCACCGATATGAg                                                   <  2:1931519/100‑1 (MQ=255)

                                               ccaccaTGAAAGCAGGCGTCAGCTATGTGGATTACCACATCCAGTTCCATCCGCGCATCGCCAAACTGCTGCGTAAACATCAAATCATCACCGATATGAg                                                   <  2:3690690/100‑1 (MQ=255)

                                               ccaccaTGAAAGCAGGCGTCAGCTATGTGGATTACCACATCCAGTTCCATCCGCGCATCGCCAAACTGCTGCGTAAACATCAAATCATCACCGATATGAg                                                   <  2:4794852/100‑1 (MQ=255)

                                               ccaccaTGAAAGCAGGCGGCAGCTATGTGGATTACCACATCCAGTTCCATCCGCGCATCGCCAAACTGCTGCGTAAACATCAAATCATCACCGATATGAg                                                   <  1:628324/100‑1 (MQ=255)

                                                caccaTGAAAGCAGGCGTCAGCTATGTGGATTACCACATCCAGTTCCATCCGCGCATCGCCAAACTGCTGCGTAAACATCAAATCATCACCGATATGAGt                                                  >  2:3976638/1‑100 (MQ=255)

                                                caccaTGAAAGCAGGCGTCAGCTATGTGGATTACCACATCCAGTTCCATCCGCGCATCGCCAAACTGCTGCGTAAACATCAAATCATCACCGATATGAGt                                                  >  2:2934590/1‑100 (MQ=255)

                                                  ccaTGAAAGCAGGCGTCAGCTATGTGGATTACGACATCAAGTTCCAGCCGCGCATCGCCAAACTGCTGCGTAAACATCAACTCATCACCGATATGAGTga                                                >  1:4940884/1‑100 (MQ=255)

                                                  ccaTGAAAGCAGGCGTCAGCTATGTGGATTACCACATCCAGTTCCATCCGCGCATCGCCAAACTGCTGCGTAAACATCAAATCATCACCGATATGAGTga                                                >  1:5144142/1‑100 (MQ=255)

                                                   caTGAAAGCAGGCGTCAGCTATGTGGATTACCACATCCAGTTCCATCCGCGCATCGCCAAACTGCTGCGTAAACATCAAATCATCACCGATATGAGTgaa                                               >  2:1443974/1‑100 (MQ=255)

                                                     tGAAAGCAGGCGTCAGCTATGTGGATTACCACATCCAGTTCCATCCGCGCATCGCCAAACTGCTGCGTAAACATCAAATCATCACCGATATGAGTgaaga                                             <  1:3629689/100‑1 (MQ=255)

                                                     tGAAAGCAGGCGTCAGCTATGTGGATTACCACATCCAGTTCCATCCGCGCATCGCCAAACTGCTGCGTAAACATCAAATCATCACCGATATCAGTgaaga                                             <  2:1459753/100‑1 (MQ=255)

                                                      gAAAGCAGGCGTCAGCTATGTGGATTACCACATCCAGTTCCATCCGCGCATCGCCAAACTGCTGCGTAAACATCAAATCATCACCGATATGAGTgaagaa                                            >  1:1349347/1‑100 (MQ=255)

                                                       aaaGCAGGCGTCAGCTATGTGGATTACCACATCCAGTTCCATCCGCGCATCGCCAAACTGCTGCGTAAACATCAAATCATCACCGATATGAGTGaagaag                                           <  1:848669/100‑1 (MQ=255)

                                                        aaGCAGGCGTCAGCTATGTGGATTACCACATCCAGTTCCATCCGCGCATCGCCAAACTGCTGCGTAAACATCAAATCATCACCGATATGAGTGAAGAAGc                                          >  1:835900/1‑100 (MQ=255)

                                                          gCAGGCGTCAGCTATGTGGATTACCACATCCAGTTCCATCCGCGCATCGCCAAACTGCTGCGTAAACATCAAATCATCACCGATATGAGTGAAGAAGCGa                                        >  1:474014/1‑100 (MQ=255)

                                                            aGGCGTCAGCTATGTGGATTACCACATCCAGTTCCATCCGCGCATCGCCAAACTGCTGCGTAAACATCAAATCATCACCGATATGAGTGAAGAAGCGATg                                      <  1:4503525/100‑1 (MQ=255)

                                                            aGGCGTCAGCTATGTGGATTACCACATCCAGTTCCATCCGCGCATCGCCAAACTGCTGCGTAAACATCAAATCATCACCGATATGAGTGAAGAAGCGATg                                      <  1:5064052/100‑1 (MQ=255)

                                                            aGGCGTCAGCTATGTGGATTACCACATCCAGTTCCATCCGCGCATCGCCAAACTGCTGCGTAAACATCAAATCATCACCGATATGAGTGAAGAAGCGATg                                      <  1:986594/100‑1 (MQ=255)

                                                            aGGCGTCAGCTATGTGGATTACCACAGCCAGTTCCATCCGCGCATCGCCAGACTGCTGCGTAAACATCAAATCATCACCGATATGAGTGAATAAGCGATg                                      <  2:5275143/100‑1 (MQ=255)

                                                              gCGTCAGCTATGTGGATTACCACATCCAGTTCCATCCGCGCATCGCCAAACTGCTGCGTAAACATCAAATCATCACCGATATGAGTGAAGAAGCGATGGt                                    <  1:3900432/100‑1 (MQ=255)

                                                               cGTCAGCTATGTGGATTACCACATCCAGTTCCATCCGCGCATCGCCAAACTGCTGCGTAAACATCAAATCATCACCGATATGAGTGAAGAAGCGATGGTc                                   <  1:5301956/100‑1 (MQ=255)

                                                                gTCAGCTATGTGGATTACCACATCCAGTTCCATCCGCGCATCGCCAAACTGCTGCGTAAACATCAAATCATCACCGATATGAGTGAAGAAGCGATGGTCg                                  <  2:2838444/100‑1 (MQ=255)

                                                                gTCAGCTATGTGGATTACCACATCCAGTTCCATCCGCGCATCGACAAACTGCTGCGTAAACATCAAATCATCACCGATATGAGTGAAGAAGCGATGGTCg                                  >  2:1049809/1‑100 (MQ=255)

                                                                 tCAGCTATGTGGATTACCACATCCAGTTCCATCCGCGCATCGCCAAACTGCTGCGTAAACATCAAATCATCACCGATATGAGTGAAGAAGCGATGGTCGa                                 <  2:3278811/100‑1 (MQ=255)

                                                                  cAGCTATGTGGATTACCACATCCAGTTCCATCCGCGCATCGCCAAACTGCTGCGTAAACATCAAATCATCACCGATATGAGTGAAGAAGCGATGGTCGaa                                >  1:944654/1‑100 (MQ=255)

                                                                  cAGCTATGTGGATTACCACATCCAGTTCCATCCGCGCATCGCCAAACTGCTGCGTAAACATCAAATCATCACCGACATGAGAGAAGAAGCGATGGCCGaa                                >  2:4252429/1‑100 (MQ=255)

                                                                   aGCTATGTGGATTACCACATCCAGTTCCATCCGCGCATCGCCAAACTGCTGCGTAAACATCAAATCATCACCGATATGAGTGAAGAAGCGATGGTCGaaa                               >  1:2133973/1‑100 (MQ=255)

                                                                    gCTATGTGGATTACCACATCCAGTTCCATCCGCGCATCGCCAAACTGCTGCGTAAACATCAAATCATCACCGATATGAGTGAAGAAGCGATGGTCGaaaa                              <  2:405222/100‑1 (MQ=255)

                                                                     cTATGTGGATTACCACATCCAGTTCCATCCGCGCATCGCCAAACTGCTGCGTAAACATCAAATCATCACCGATATGAGTGAAGAAGCGATGGTCGAAAAc                             <  1:1970121/100‑1 (MQ=255)

                                                                        tgtgGATTACCACATCCAGTTCCATCCGCGCATCGCCAAACTGCTGCGTAAACATCAAATCATCACCGATATGAGTGAAGAAGCGATGGTCGAAAACGAt                          >  1:2201928/1‑100 (MQ=255)

                                                                        tgtgGATTACCACATCCAGTTCCATCCGCGCATCGCCAAACTGCTGCGTAAACATCAAATCATCACCGATATGAGTGAAGAAGCGATGGTCGAAAACGAt                          <  2:3018223/100‑1 (MQ=255)

                                                                        tgtgGATTACCACATCCAGTTCCATCCGCGCATCGCCAAACTGCTGCGTAAACATCAAATCATCACCGATATGAGTGAAGAAGCGATGGTCGAAAACGAt                          <  1:2294520/100‑1 (MQ=255)

                                                                         gtgGATTACCACATCCATTTCCATCCGGGTATCGCCAAACTGCTGCGCAAACATCAAATCATCACCGATATGAGTGATGAATCGATGGTCGAAAACgcgc                         <  2:916231/100‑4 (MQ=255)

                                                                         gtgGATTACCACATCCAGTTCCATCCGCGCATCGCCAAACTGCTGCGTAAACATCAAATCATCACCGATATGAGTGAAGAAGCGATGGTCGAAAACGATc                         <  1:3430028/100‑1 (MQ=255)

                                                                         gtgGATTACCACATCCAGTTCCATCCGCGCATCGCCAAACTGCTGCGTAAACATCAAATCATCACCGATATGAGTGAAGAAGCGATGGTCGAAAACGATc                         <  1:2513724/100‑1 (MQ=255)

                                                                         gtgGATTACCACATCCAGTTCCATCCGCGCATAGCCAAACTGCTGCGTAAACATCAAATCATCACCGATATGAGTGAAGAAGCGATGGTCGAAAACGATc                         <  2:5498238/100‑1 (MQ=255)

                                                                                aCCACATCCAGTTCCATCCGCGCATCGCCAAACTGCTGCGTAAACATCAAATCATCACCGATATGAGTGAAGAAGCGATGGTCGAAAACGATCTTACCgg                  <  2:4406095/100‑1 (MQ=255)

                                                                                aCCACATCCAGTTCCATCCGCGCATCGCCAAACTGCTGCGTAAACATCAAATCATCACCGATATGAGTGAAGAAGCGATGGTCGAAAACGATCTTACCgg                  >  2:1639180/1‑100 (MQ=255)

                                                                                aCCACATCCAGTTCCATCCGCGCATCGCCAAACTGCTGCGTAAACATCAAATCATCACCGATATGAGTGAAGAAGCGATGGTCGAAAACGATCTTACCgg                  >  1:2963542/1‑100 (MQ=255)

                                                                                aCCACATCCAGTTCCATCCGCGCATCGCCAAACTGCTGCGTAAACATCAAATCATCACCGATATGAGTGAAGAAGCGATGGTCGAAAACGATCTTACCgg                  >  2:1126619/1‑100 (MQ=255)

                                                                                aCCACATCCAGTTCCATCCGCGCATCGCCAAACTGCTGCGTAAACATCAAATCATCACCGATATGAGTGAAGAAGCGATGGTCGAAAACGATCTTACCgg                  >  1:754749/1‑100 (MQ=255)

                                                                                aCCACATCCAGTTCCATCCGCGCATCGCCAAACTGCTGCGTAAACATCAAATCATCACCGATATGAGTGAAGAAGCGATGGTCGAAAACGATCTTACCgg                  <  1:4570989/100‑1 (MQ=255)

                                                                                aCCACATCCAGTTCCATCCGCGCATCGCCAAACTGCTGCGTAAACATCAAATCATCACCGATATGAGTGAAGAAGCGATGGTCGAAAACGATCTTACCgg                  >  2:4875573/1‑100 (MQ=255)

                                                                                 ccACATCCAGTTCCATCCGCGCATCGCCAAACTGCTGCGTAAACATCAAATCATCACCGATATGAGTGAAGAAGCGATGGTCGAAAACGATCTTACCggg                 >  2:4118244/1‑100 (MQ=255)

                                                                                  cacaTCCAGTTCCATCCGCGCATCGCCAAACTGCTGCGTAAACATCAAATCATCACCGATATGAGTGAAGAAGCGATGGTCGAAAACGATCTTACCGGGc                >  2:1992404/1‑100 (MQ=255)

                                                                                   acaTCCAGTTCCATCCGCGCATCGCCAAACTGCTGCGTAAACATCAAATCATCACCGATATGAGTGAAGAAGCGATGGTCGAAAACGATCTTACCGGGcc               <  1:4327898/100‑1 (MQ=255)

                                                                                    caTCCAGTTCCATCCGCGCATCGCCAAACTGCTGCGTAAACATCAAATCATCACCGATATGAGTGAAGAAGCGATGGTCGAAAACGATCTTACCGGGCCg              >  1:451766/1‑100 (MQ=255)

                                                                                        cAGTTCCATCCGCGCATCGCCAAACTGCTGCGTAAACATCAAATCATCACCGATATGAGTGAAGAAGCGATGGTCGAAAACGATCTTACCGGGCCGTTTa          <  2:1238898/100‑1 (MQ=255)

                                                                                        cAGTTCCATCCGCGCATCGCCAAACTGCTGCGTAAACATCAAATCATCACCGATATGAGTGAAGAAGCGATGGTCGAAAACGATCTTACCGGGCCGTTTa          >  2:1763052/1‑100 (MQ=255)

                                                                                         aGTTCCATCCGCGCATCGCCAAACTGCTGCGTAAATATCAAATCATCACCGATATGAGTGAAGAAGCGATGGTCGAAAACGATCTTACCGGGCCGTTTAt         <  2:5424916/100‑1 (MQ=255)

                                                                                         aGTTCCATCCGCGCATCGCCAAACTGCTGCGTAAACATCAAATCATCACCGATATGAGTGAAGAAGCGATGGTCGAAAACGATCTTACCGGGCCGTTTAt         <  2:1331521/100‑1 (MQ=255)

                                                                                         aGTTCCATCCGCGCATCGCCAAACTGCTGCGTAAACATCAAATCATCACCGATATGAGTGAAGAAGCGATGGTCGAAAACGATCTTACCGGGCCGTTTAt         <  1:1159461/100‑1 (MQ=255)

                                                                                          gTTCCATCCGCGCATCGCCAAACTGCTGCGTAAACATCAAATCATCACCGATATGAGTGAAGAAGCGATGGTCGAAAACGATCTTACCGGGCCGTTTATg        <  2:1656291/100‑1 (MQ=255)

                                                                                          gTTCCATCCGCGCATCGCCAAACTGCTGCGTAAACATCAAATCATCACCGATATGAGTGAAGAAGCGATGGTCGAAAACGATCTTACCGGGCCGTTTATg        >  2:598897/1‑100 (MQ=255)

                                                                                          gTTCCATCCGCGCATCGCCAAACTGCTGCGTAAACATCAAATCATCACCGATATGAGTGAAGAAGCGATGGTCGAAAACGATCTTACCGGGCCGTTTATg        >  1:4086938/1‑100 (MQ=255)

                                                                                          gTTCCATCCGCGCATCGCCAAACTGCTGCGTAAACATCAAATCATCACCGATATGAGTGAAGAAGCGATGGTCGAAAACGATCTTACCGGGCCGTTTATg        >  1:2937543/1‑100 (MQ=255)

                                                                                           ttCCATCCGCGCATCGCCAAACTGCTGCGTAAACATCAAATCATCACCGATATGAGTGAAGAAGCGATGGTCGAAAACGATCTTACCGGGCCGTTTATGc       <  2:5073359/100‑1 (MQ=255)

                                                                                              cATCCGCGCATCGCCAAACTGCTGCGTAAACATCAAATCATCACCGATATGAGTGAAGAAGCGATGGTCGAAAACGATCTTACCGGGCCGTTTATGCCGc    >  1:4672865/1‑100 (MQ=255)

                                                                                               aTCCGCGCATAGCCAAGCTGCTGCGTAAACATCAAATCATCACCAATATGAGTGAAGAAGCGATGGTCGAAAACGATCTTACCGGGCCGTTTATGCCGCa   >  1:4878248/1‑100 (MQ=255)

CGCGCAGCTGGTGAAGGACGTTAATGATGAACAACTGGCGCTGATCGCCACCATGAAAGCAGGCGTCAGCTATGTGGATTACCACATCCAGTTCCATCAGCGCATCGCCAAACTGCTGCGTAAACATCAAATCATCACCGATATGAGTGAAGAAGCGATGGTCGAAAACGATCTTACCGGGCCGTTTATGCCGCAT  >  NC_000913/4031970‑4032165

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