BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,196,394	A→G	Q11Q (CAA→CAG)	nfi →	endonuclease V

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,196,394	0	A	G	88.8%	807.1 / 78.0	278	Q11Q (CAA→CAG)	nfi	endonuclease V
Reads supporting (aligned to +/- strand): ref base A (19/12); new base G (126/121); total (145/133)
Fisher's exact test for biased strand distribution p-value = 3.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.96e-01

Only 100 of 287 total aligned reads displayed.

TGCTGGCGTGTTCGTGGAGGCAGTGCATCGACTGTCTGAACAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGA  >  NC_000913/4196304‑4196483

tGCTGGCGTGTTCGTGGAGGCAGTGCATCGACTGTCTGAACAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTcagca                                                                                         >  2:3763039/1‑93 (MQ=255)

 gCTGGCGTGTTCGTGGAGGCAGTGCATCGACTGTCTGAACAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTcagcaa                                                                                        >  1:5503146/1‑93 (MQ=255)

   tGGCGTGTTCGTGGAGGCAGTGCATCGACTGTCTGAACAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAGCAAAt                                                                                      <  1:3066963/93‑1 (MQ=255)

   tGGCGTGTTCGTGGAGGCAGTGCATCGACTGTCTGAACAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAGCAAAt                                                                                      >  1:2359039/1‑93 (MQ=255)

   tGGCGTGTTCGTGGAGGCAGTGCATCGACTGTCTGAACAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAAt                                                                                      <  1:1104312/93‑1 (MQ=255)

    ggCGTGTTCGTGGAGGCAGTGCATCGACTGTCTGAACAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAGCAAATc                                                                                     >  2:2575571/1‑93 (MQ=255)

    ggCGTGTTCGTGGAGGCAGTGCATCGACTGTCTGAACAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAGCAAATc                                                                                     <  1:5121872/93‑1 (MQ=255)

     gCGTGTTCGTGGAGGCAGTGCATCGACTGTCTGAACAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAGCAAATCg                                                                                    >  2:706082/1‑93 (MQ=255)

        tgtTCGTGGAGGCAGTGCATCGACTGTCTGAACAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAATCGAAc                                                                                 >  2:6538306/1‑93 (MQ=255)

         gtTCGTGGAGGCAGTGCATCGACTGTCTGAACAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAGCAAATCGAACt                                                                                <  2:5258549/93‑1 (MQ=255)

         gtTCGTGGAGGCAGTGCATCGACTGTCTGAACAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAGCAAATCGAACt                                                                                >  2:5695839/1‑93 (MQ=255)

          ttCGTGGAGGCAGTGCATCGACTGTCTGAACAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAGCAAATCGAACTg                                                                               >  1:4077881/1‑93 (MQ=255)

            cGTGGAGGCAGTGCATCGACTGTCTGAACAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGc                                                                             <  1:3770500/93‑1 (MQ=255)

             gTGGAGGCAGTGCATCGACTGTCTGAACAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAGCAAATCGAACTGGct                                                                            >  1:2189507/1‑93 (MQ=255)

             gTGGAGGCAGTGCATCGACTGTCTGAACAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAGCAAATCGAACTGGct                                                                            <  1:2533084/93‑1 (MQ=255)

              tGGAGGCAGTGCATCGACTGTCTGAACAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAGCAAATCGAACTGGctt                                                                           <  2:4946294/93‑1 (MQ=255)

                gAGGCAGTGCATCGACTGTCTGAACAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAGCAAATCGAACTGGcttct                                                                         >  1:1517430/1‑93 (MQ=255)

                gAGGCAGTGCATCGACTGTCTGAACAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAGCAAATCGAACTGGcttct                                                                         >  1:3868070/1‑93 (MQ=255)

                gAGGCAGTGCATCGACTGTCTGAACAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGcttct                                                                         >  1:91805/1‑93 (MQ=255)

                 aGGCAGTGCATCGACTGTCTGAACAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAGCAAATCGAACTGGcttctt                                                                        <  2:14498/93‑1 (MQ=255)

                 aGGCAGTGCATCGACTGTCTGAACAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAGCAAATCGAACTGGcttctt                                                                        >  2:2159787/1‑93 (MQ=255)

                  ggCAGTGCATCGACTGTCTGAACAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAGCAAATCGAACTGGCttcttc                                                                       <  1:6022403/93‑1 (MQ=255)

                  ggCAGTGCATCGACTGTCTGAACAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAGCAAATCGAACTGGCttcttc                                                                       <  2:3351827/93‑1 (MQ=255)

                  ggCAGTGCATCGACTGTCTGAACAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAGCAAATCGAACTGGCttcttc                                                                       <  1:1047135/93‑1 (MQ=255)

                  ggCAGTGCATCGACTGTCTGAACAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAGCAAATCGAACTGGCttcttc                                                                       >  1:4243885/1‑93 (MQ=255)

                     aGTGCATCGACTGTCTGAACAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAGCAAATCGAACTGGCTTCTTCtgt                                                                    <  2:1200186/93‑1 (MQ=255)

                      gTGCATCGACTGTCTGAACAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAGCAAATCGAACTGGCTTCTTCtgtg                                                                   <  2:632183/93‑1 (MQ=255)

                       tGCATCGACTGTCTGAACAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCGGCAAATCGAACTGGCTTCTTCTGTGa                                                                  <  2:5514557/93‑1 (MQ=255)

                        gCATCGACTGTCTGAACAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAGCAAATCGAACTGGCTTCTTCTGTGAt                                                                 >  1:6263970/1‑93 (MQ=255)

                         cATCGACTGTCTGAACAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAGCAAATCGAACTGGCTTCTTCTGTGATc                                                                <  1:3335818/93‑1 (MQ=255)

                         cATCGACTGTCTGAACAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAGCAAATCGAACTGGCTTCTTCTGTGATc                                                                >  1:2263645/1‑93 (MQ=255)

                         cATCGACTGTCTGAACAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAGCAAATCGAACTGGCTTCTTCTGTGATc                                                                <  2:2985980/93‑1 (MQ=255)

                         cATCGACTGTCTGAACAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAGCAAATCGAACTGGCTTCTTCTGTGATc                                                                <  2:3849445/93‑1 (MQ=255)

                           tCGACTGTCTGAACAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAGCAAATCGAACTGGCTTCTTCTGTGATCcg                                                              >  2:594018/1‑93 (MQ=255)

                           tCGACTGTCTGAACAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAGCAAATCGAACTGGCTTCTTCTGTGATCcg                                                              <  2:2953194/93‑1 (MQ=255)

                            cGACTGTCTGAACAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAGCAAATCGAACTGGCTTCTTCTGTGATCcgc                                                             <  2:2729093/93‑1 (MQ=255)

                             gACTGTCTGAACAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAGCAAATCGAACTGGCTTCTTCTGTGATCcgcg                                                            >  2:1935818/1‑93 (MQ=255)

                              aCTGTCTGAACAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAGCAAATCGAACTGGCTTCTTCTGTGATCCGCGa                                                           >  2:4431546/1‑93 (MQ=255)

                              aCTGTCTGAACAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAGCAAATCGAACTGGCTTCTTCTGTGATCCGCGa                                                           >  1:3009636/1‑93 (MQ=255)

                               cTGTCTGAACAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAGCAAATCGAACTGGCTTCTTCTGTGATCCGCGAg                                                          <  1:1244893/93‑1 (MQ=255)

                                tGTCTGAACAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAGCAAATCGAACTNGCTTCTTCTGTGATCCGCGNgg                                                         >  2:4689036/1‑93 (MQ=255)

                                 gTCTGAACAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAGCAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGa                                                        >  2:2478161/1‑93 (MQ=255)

                                 gTCTGAACAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGa                                                        >  2:2877128/1‑93 (MQ=255)

                                  tCTGAACAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAGCAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGAt                                                       <  1:2597326/93‑1 (MQ=255)

                                  tCTGAACAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAGCAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGAt                                                       <  1:762282/93‑1 (MQ=255)

                                    tGAACAGTCTCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAGCAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCg                                                     >  1:6467249/1‑93 (MQ=255)

                                      aaCAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGAc                                                   <  2:17830/93‑1 (MQ=255)

                                       aCAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAGCAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACt                                                  <  2:47612/93‑1 (MQ=255)

                                        cAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAGCAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTc                                                 <  2:3949283/93‑1 (MQ=255)

                                        cAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAGCAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTc                                                 <  2:1389042/93‑1 (MQ=255)

                                          gTATCACCGCTAAGGAGTGATTATGGATCTCTCGTCATTACGCGCTCAGCAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGGTCGCCTCGa                                               >  2:4511785/1‑93 (MQ=255)

                                            aTCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAGCAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATa                                             >  2:3462830/1‑93 (MQ=255)

                                            aTCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAGCAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATa                                             >  1:1448936/1‑93 (MQ=255)

                                             tCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAGCAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATaa                                            >  2:3639614/1‑93 (MQ=255)

                                              cACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAGCAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATaaa                                           <  2:4846829/93‑1 (MQ=255)

                                              cACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAGCAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATaaa                                           >  2:4979787/1‑93 (MQ=255)

                                              cACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAGCAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATaaa                                           <  2:419563/93‑1 (MQ=255)

                                                ccGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGa                                         >  1:5195576/1‑93 (MQ=255)

                                                 cGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAGCAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGAt                                        >  1:3414414/1‑93 (MQ=255)

                                                 cGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGAt                                        <  2:978715/93‑1 (MQ=255)

                                                  gCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAGCAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATc                                       <  1:1765668/93‑1 (MQ=255)

                                                  gCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATc                                       <  1:843049/93‑1 (MQ=255)

                                                   cTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAGCAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATcc                                      >  2:4784605/1‑93 (MQ=255)

                                                   cTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAGCAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATcc                                      <  2:3248808/93‑1 (MQ=255)

                                                     aaGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAGCAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCAc                                    >  1:3564800/1‑93 (MQ=255)

                                                      aGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAGCAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCAcc                                   >  1:3688091/1‑93 (MQ=255)

                                                       ggAGTGATTATGGATCTCGCGTCATTACGCGCTCAGCAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCg                                  <  2:3069674/93‑1 (MQ=255)

                                                        gAGTGATTATGGATCTCGCGTCATTACGCGCTCAGCAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCgg                                 <  2:5167485/93‑1 (MQ=255)

                                                           ttaTTATGGATCTCGCGTCATTACGCGCTCAGCAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATc                              <  1:4947586/91‑1 (MQ=255)

                                                            gATTATGGATCTCGCGTCATTACGCGCTCAGCAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCCCCGGATCt                             <  1:669270/93‑1 (MQ=255)

                                                              ttATGGATCTCGCGTCATTACGCGCTCAGCAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGa                           >  2:6636604/1‑93 (MQ=255)

                                                               tacggATCTCGCGTCATTACGCGCTCAGCAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGAt                          <  1:4635788/90‑1 (MQ=255)

                                                               tATGGATTTCGCGTCATTACGCGCTCAGCAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGAt                          <  1:2396793/93‑1 (MQ=255)

                                                               tATGGATCTCGCGTCATTACGCGCTCAGCAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGAt                          <  2:3581135/93‑1 (MQ=255)

                                                               tATGGATCTCGCGTCATTACGCGCTCAGCAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGAt                          <  1:4557073/93‑1 (MQ=255)

                                                               tATGGATCTCGCGTCATTACGCGCTCAGCAAATCGAACTGGCTTCTGCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGAt                          >  1:5363603/1‑93 (MQ=255)

                                                                    aTCTCGCGTCATTACGCGCTCAGCAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCg                     <  1:5395541/93‑1 (MQ=255)

                                                                     tctcGCGTCATTACGCGCTCAGCAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCgg                    >  1:273807/1‑93 (MQ=255)

                                                                      ctcGCGTCATTACGCGCTCAGCAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGa                   <  2:5817146/93‑1 (MQ=255)

                                                                      ctcGCGTCATTACGCGCTCAGCAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGa                   <  1:1931028/93‑1 (MQ=255)

                                                                       tcGCGTCATTACGCGCTCAGCAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAg                  <  1:6184455/93‑1 (MQ=255)

                                                                       tcGCGTCATTACGCGCTCAGCAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAg                  <  2:1589122/93‑1 (MQ=255)

                                                                       tcGCGTCATTACGCGCTCAGCAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAg                  >  2:2351536/1‑93 (MQ=255)

                                                                        cgcgTCATTACGCGCTCAGCAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGc                 >  2:60748/1‑93 (MQ=255)

                                                                         gcgTCATTACGCGCTCAGCAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGcc                <  1:2023963/93‑1 (MQ=255)

                                                                         gcgTCATTACGCGCTCAGCAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGcc                <  2:5137351/93‑1 (MQ=255)

                                                                          cgTCATTACGCGCTCAGCAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCg               <  2:225573/93‑1 (MQ=255)

                                                                           gTCATTACGCGCTCAGCAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGa              <  1:187590/93‑1 (MQ=255)

                                                                            tCATTACGCGCTCAGCAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGAt             <  1:2503300/93‑1 (MQ=255)

                                                                            tCATTACGCGCTCAGCAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGAt             >  1:5726185/1‑93 (MQ=255)

                                                                            tCATTACGCGCTCAGCAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGAt             >  2:3309742/1‑93 (MQ=255)

                                                                             cATTACGCGCTCAGCAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATg            >  1:3728526/1‑93 (MQ=255)

                                                                               ttACGCGCTCAGCAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTc          >  1:4158471/1‑93 (MQ=255)

                                                                               ttACGCGCTCAGCAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTc          >  1:1159566/1‑93 (MQ=255)

                                                                                  cgcgCTCAGCAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAGAGATCTACCGGATCTGATCGCCGGAGCCGATGTCggg       >  1:2899789/1‑93 (MQ=255)

                                                                                   gcgcTCAGCAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGt      >  1:6585447/1‑93 (MQ=255)

                                                                                      cTCAGCAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTg   >  1:1883558/1‑93 (MQ=255)

                                                                                      cTCAGCAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTg   >  2:2223133/1‑93 (MQ=255)

                                                                                       tCAGCAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGa  >  2:6457311/1‑93 (MQ=255)

                                                                                       tCAACAAATCGAACTGGCTTCTTCTGTTATCCGCGAGGATCGACTCGATAAAGATCCACCGGTTCTGATCGCCGGAGCCGATGTCGGGTTTGa  <  1:4404328/93‑1 (MQ=255)

TGCTGGCGTGTTCGTGGAGGCAGTGCATCGACTGTCTGAACAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGA  >  NC_000913/4196304‑4196483

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