BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,702,500	C→G	6.5%	H90D (CAC→GAC)	add →	adenosine deaminase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,702,500	0	C	G	6.5%	485.5 / 20.4	138	H90D (CAC→GAC)	add	adenosine deaminase
Reads supporting (aligned to +/- strand): ref base C (67/62); new base G (5/4); total (72/66)
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.15e-01

Only 100 of 142 total aligned reads displayed.

CACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGGTGTGCAGGCGAAGCTTATCGGCAT  >  NC_000913/1702357‑1702648

cacTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACgg                                                                                                                                                  >  1:754965/1‑148 (MQ=255)

  cTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGcc                                                                                                                                                <  2:892635/148‑1 (MQ=255)

  cTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGcc                                                                                                                                                <  1:362240/148‑1 (MQ=255)

    gATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTg                                                                                                                                              >  1:436972/1‑148 (MQ=255)

      ttCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCAc                                                                                                                                           <  4:823977/149‑1 (MQ=255)

        ccccACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGGCCGTGACGGGCTGCatc                                                                                                                                          <  1:727312/148‑3 (MQ=255)

          ccACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTAt                                                                                                                                        >  1:502975/1‑148 (MQ=255)

          ccACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTAt                                                                                                                                        >  2:172386/1‑148 (MQ=255)

           cACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTg                                                                                                                                              <  2:299/141‑1 (MQ=255)

           cACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATg                                                                                                                                       >  2:305935/1‑148 (MQ=255)

           cACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTCAAGATGCAGCCCGTCACGGCCTg                                                                                                                                              >  1:299/1‑141 (MQ=255)

              gTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCg                                                                                                                                    <  1:810546/148‑1 (MQ=255)

                 cAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTc                                                                                                                                     >  4:373971/1‑144 (MQ=255)

                 cAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTc                                                                                                                                     >  4:725850/1‑144 (MQ=255)

                 cAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTc                                                                                                                                     <  3:725850/144‑1 (MQ=255)

                   ggTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTg                                                                                                                               >  1:771478/1‑148 (MQ=255)

                    gTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGGCCGTGACGGGCTGCatct                                                                                                                                         >  3:890133/1‑134 (MQ=255)

                      cATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGt                                                                                                                            >  1:509970/1‑148 (MQ=255)

                       aTTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGttt                                                                                                                          <  4:700673/149‑1 (MQ=255)

                       aTTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGttt                                                                                                                          <  4:293494/149‑1 (MQ=255)

                             aaCGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCAc                                                                                                                     >  2:751166/1‑148 (MQ=255)

                              aCGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCAcc                                                                                                                    >  2:642132/1‑148 (MQ=255)

                                   cccGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACg                                                                                                                                                   <  4:552139/112‑1 (MQ=255)

                                   cccGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACg                                                                                                                                                   >  3:552139/1‑112 (MQ=255)

                                     cGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGtttt                                                                                                                         >  1:837074/1‑136 (MQ=255)

                                     cGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGtttt                                                                                                                         <  2:837074/136‑1 (MQ=255)

                                       aTCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACatg                                                                                                          <  3:238615/149‑1 (MQ=255)

                                        tCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGCTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACatg                                                                                                          >  2:110432/1‑148 (MQ=255)

                                             tGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACatggcaatg                                                                                                    <  4:874748/149‑1 (MQ=255)

                                               aGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACatggcaatggc                                                                                                  <  4:384689/149‑1 (MQ=255)

                                                gCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACatggcaatggca                                                                                                 <  4:692437/149‑1 (MQ=255)

                                                gCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACatggcaatggc                                                                                                  >  2:789129/1‑148 (MQ=255)

                                                         ccAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCt                                                                                         >  1:885956/1‑148 (MQ=255)

                                                         ccAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTg                                                                                        >  3:197079/1‑149 (MQ=255)

                                                          cAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTTAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGc                                                                                       >  4:549314/1‑149 (MQ=255)

                                                          cAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGc                                                                                       >  4:863849/1‑149 (MQ=255)

                                                          cAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGc                                                                                       >  4:427498/1‑149 (MQ=255)

                                                          cAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGc                                                                                       <  3:355990/149‑1 (MQ=255)

                                                          cAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGACTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACatgg                                                                                                         <  4:683297/131‑1 (MQ=255)

                                                          cAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGACTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACatgg                                                                                                         >  3:683297/1‑131 (MQ=255)

                                                           aaaCTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTc                                                                                                                                                      >  1:420053/1‑85 (MQ=255)

                                                           aaaCTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGc                                                                                       >  1:192324/1‑148 (MQ=255)

                                                               ttGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTa                                                                                  >  3:406429/1‑149 (MQ=255)

                                                                   cTGGCGCGCTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGTCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTTCGTTTTTCACCAGGCTACAGGGCAATGGCACATCAGCTGCCTGTAGCgg                                                                              <  4:345944/149‑1 (MQ=255)

                                                                     ggggCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGgt                                                                            <  4:724485/149‑1 (MQ=255)

                                                                       ggCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGgtg                                                                           >  2:684951/1‑148 (MQ=255)

                                                                        gCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGgtgt                                                                          <  2:462711/148‑1 (MQ=255)

                                                                          gTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGtgttgt                                                                       >  4:362485/1‑149 (MQ=255)

                                                                           ttAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTc                                                                      >  3:420251/1‑149 (MQ=255)

                                                                           ttAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTc                                                                      <  3:863849/149‑1 (MQ=255)

                                                                             aaaGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGa                                                                    >  4:296299/1‑149 (MQ=255)

                                                                                  tctcGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCg                                                                >  1:787431/1‑148 (MQ=255)

                                                                                   ctcGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTAGATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGt                                                              >  3:766735/1‑149 (MQ=255)

                                                                                   ctcGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGt                                                              >  3:673336/1‑149 (MQ=255)

                                                                                    tcGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGt                                                              <  1:245330/148‑1 (MQ=255)

                                                                                     cGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTg                                                             <  2:885956/148‑1 (MQ=255)

                                                                                     cGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGa                                                            <  4:17980/149‑1 (MQ=255)

                                                                                      gCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATc                                                                                             >  3:3557/1‑115 (MQ=255)

                                                                                      gCCTCTCTTGAGGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACAAGGCTATATGGCAGTGGCACATCAGCTGCCTGTAGCGGGTTTTGTCGAAgatcggaa                                                           >  4:82019/1‑142 (MQ=255)

                                                                                        ctctctTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAg                                                                                 >  3:836193/1‑125 (MQ=255)

                                                                                         tctctTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGa                                                        >  3:13668/1‑149 (MQ=255)

                                                                                          ctctTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCGCCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGAt                                                       <  4:717719/149‑1 (MQ=255)

                                                                                          ctctTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGa                                                        <  2:754965/148‑1 (MQ=255)

                                                                                          ctctTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGAt                                                       <  3:654965/149‑1 (MQ=255)

                                                                                          ctctTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGAGGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGa                                                        <  1:110432/148‑1 (MQ=255)

                                                                                             ttGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGt                                                                                                                                      <  3:461998/67‑1 (MQ=255)

                                                                                             ttGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGt                                                                                                                                      >  4:461998/1‑67 (MQ=255)

                                                                                             ttGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATgg                                                     <  1:789129/148‑1 (MQ=255)

                                                                                             ttGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATgg                                                     >  1:171862/1‑148 (MQ=255)

                                                                                              tGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATg                                                                                                                                       <  1:387523/65‑1 (MQ=255)

                                                                                              tGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATg                                                                                                                                       >  2:387523/1‑65 (MQ=255)

                                                                                                 tGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTAc                                                <  3:216230/149‑1 (MQ=255)

                                                                                                  gCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACg                                               >  4:157205/1‑149 (MQ=255)

                                                                                                  gCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACg                                               >  4:548772/1‑149 (MQ=255)

                                                                                                   ccTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACg                                               >  2:675721/1‑148 (MQ=255)

                                                                                                      gTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGaa                                           <  4:673336/149‑1 (MQ=255)

                                                                                                         gccgcGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGt                                        <  4:197079/149‑1 (MQ=255)

                                                                                                            gcgTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGc                                     <  4:351460/149‑1 (MQ=255)

                                                                                                             cgTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGc                                     <  1:181911/148‑1 (MQ=255)

                                                                                                                  cATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCAcc                               <  4:203062/149‑1 (MQ=255)

                                                                                                                   aTTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCAcc                               >  1:856243/1‑148 (MQ=255)

                                                                                                                     ttGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCtt                             >  2:452770/1‑148 (MQ=255)

                                                                                                                     ttGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCtt                             >  2:426682/1‑148 (MQ=255)

                                                                                                                        aaaaCATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGgt                         >  3:672309/1‑149 (MQ=255)

                                                                                                                            cATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGGtgtg                      >  1:160802/1‑148 (MQ=255)

                                                                                                                             aTTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGGTGTGc                     <  2:771478/148‑1 (MQ=255)

                                                                                                                               tGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGGTGTGCAgg                  <  3:157205/149‑1 (MQ=255)

                                                                                                                               tGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGGTGTGCAgg                  <  3:548772/149‑1 (MQ=255)

                                                                                                                                 aaGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGGTGTGCAGGc                 >  1:218188/1‑148 (MQ=255)

                                                                                                                                    aTGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGGTGTGCAGGCGaa              <  2:171862/148‑1 (MQ=255)

                                                                                                                                     tGCAGGCCGTGACGGGCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATc                                                                                             <  2:402779/68‑1 (MQ=38)

                                                                                                                                     tGCAGGCCGTGACGGGCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGGTGTGCAGGCGAAg             >  2:695544/1‑148 (MQ=255)

                                                                                                                                      gCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGGTGTGCAGGCGAAGCt           >  3:19413/1‑149 (MQ=255)

                                                                                                                                      gCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGGTGTGCAGGCGAAGCt           >  3:494036/1‑149 (MQ=255)

                                                                                                                                        aGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGGTGTGCAGGCGAAGCtt          >  2:531931/1‑148 (MQ=255)

                                                                                                                                        aGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGGTGTGCAGGCGAAGCTTa         <  4:494036/149‑1 (MQ=255)

                                                                                                                                         ggccGTGACGGGCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGGTGTGCAGGCGAAGCTTa         >  1:439814/3‑148 (MQ=255)

                                                                                                                                         gCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGaa                                           <  4:239345/114‑1 (MQ=255)

                                                                                                                                         gCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGaa                                           >  3:239345/1‑114 (MQ=255)

                                                                                                                                               cACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGGTGTGCAGGCGAAGCTTATCGGCAt  <  4:406429/149‑1 (MQ=255)

CACTGATTCCCCACGTTCAGGTCATTGCCAACGAACCCGATCTGGTGAGCTTTCTGACCAAACTTGACTGGGGCGTTAAAGTTCTCGCCTCTCTTGATGCCTGTCGCCGCGTGGCATTTGAAAACATTGAAGATGCAGCCCGTCACGGCCTGCACTATGTCGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGGTGTGCAGGCGAAGCTTATCGGCAT  >  NC_000913/1702357‑1702648

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