BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,603,663	(G)_8→9	noncoding (72/88 nt)	leuV ←	tRNA‑Leu

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,603,655	1	.	G	94.3%	558.5 / 25.4	176	noncoding (80/88 nt)	leuV	tRNA‑Leu
Reads supporting (aligned to +/- strand): ref base . (6/4); new base G (107/59); total (113/63)
Fisher's exact test for biased strand distribution p-value = 7.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.32e-01

Only 100 of 216 total aligned reads displayed.

TTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GA‑‑GGGGGGG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCAT  >  NC_000913/4603563‑4603742

ttATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑Ga                                                                                             >  1:4833347/1‑93 (MQ=255)

 tATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GA‑‑g                                                                                          >  2:1572501/1‑93 (MQ=255)

 tATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GA‑‑g                                                                                          >  2:1754543/1‑93 (MQ=255)

  aTTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GA‑‑gg                                                                                         <  2:6015955/93‑1 (MQ=255)

   ttGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GA‑‑ggg                                                                                        <  2:4977888/93‑1 (MQ=255)

   ttGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GA‑‑ggg                                                                                        <  1:5015292/93‑1 (MQ=255)

   ttGCTGGCATCATGCTGATGATAATGAATTAGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GA‑‑ggg                                                                                        <  2:6463631/93‑1 (MQ=255)

     gCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GA‑‑ggggg                                                                                      <  1:1341747/93‑1 (MQ=255)

      cTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GA‑‑gggggg                                                                                     >  1:63326/1‑93 (MQ=255)

        ggCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAAAAAGGAAGGAAAGATAATTGGTGC‑GA‑‑ggggggg‑g                                                                                  <  2:5750071/93‑1 (MQ=255)

           atcatGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GA‑‑ggggggg‑ggcc                                                                               >  2:2746247/1‑90 (MQ=255)

            tcatGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GA‑‑ggggggg‑ggact                                                                              >  2:723067/1‑89 (MQ=255)

            tcatGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GA‑‑ggggggg‑ggact                                                                              <  2:509558/93‑5 (MQ=255)

            tcatGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GA‑‑ggggggg‑ggact                                                                              <  2:3018883/93‑5 (MQ=255)

            tcatGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GA‑‑ggggggg‑ggact                                                                              >  2:2533563/1‑89 (MQ=255)

              atGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GA‑‑GGGGGGG‑G‑ACTTGa                                                                           >  1:4412201/1‑93 (MQ=255)

               tGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GAG‑GGGGGGG‑G‑ACTTGa                                                                           >  2:5306639/1‑93 (MQ=255)

                 cTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GA‑‑GGGGGGG‑G‑ACTTGAAcc                                                                        >  2:1511676/1‑93 (MQ=255)

                   gatgatAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAGAGGAAGTAAAGATAATTGGTGC‑GAG‑GGGGGGG‑G‑ACTTGAAccc                                                                       >  2:4289401/1‑93 (MQ=255)

                   gatgatAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GAG‑GGGGGGG‑G‑ACTTGAAccc                                                                       >  1:4358114/1‑93 (MQ=255)

                   gatgatAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GAG‑GGGGGGG‑G‑ACTTGAAccc                                                                       >  1:2146086/1‑93 (MQ=255)

                     tgatAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GAG‑GGGGGGG‑G‑ACTTGAAccccc                                                                     >  1:2398528/1‑93 (MQ=255)

                      gatAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GAG‑GGGGGGG‑G‑ACTTGAACCCCCa                                                                    <  2:4309279/93‑1 (MQ=255)

                      gatAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GAG‑GGGGGGG‑G‑ACTTGAACCCCCa                                                                    <  2:758151/93‑1 (MQ=255)

                      gatAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GAG‑GGGGGGG‑G‑ACTTGAACCCCCa                                                                    >  2:5940899/1‑93 (MQ=255)

                      gatAATGAATGTGAATGAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GAG‑GGGGGGG‑G‑ACTTGAACCCCCa                                                                    <  2:1975164/93‑1 (MQ=39)

                       atAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GAG‑GGGGGGG‑G‑ACTTGAACCCCCAc                                                                   >  1:3434601/1‑93 (MQ=255)

                         aaTGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GAG‑GGGGGGG‑G‑ACTTGAACCCCCACGt                                                                 <  2:2649484/93‑1 (MQ=255)

                         aaTGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GAG‑GGGGGGG‑G‑ACTTGAACCCCCACGt                                                                 >  2:3962974/1‑93 (MQ=255)

                          aTGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GAG‑GGGGGGG‑G‑GCTTGAACCCCCACGTc                                                                >  1:5726839/1‑93 (MQ=255)

                          aTGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GAG‑GGGGGGG‑G‑ACTTGAACCCCCACGTc                                                                >  1:3631432/1‑93 (MQ=255)

                          aTGAATTTAAAACAGGAATAAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGAAAATTGGTGG‑GAG‑GGGGGGG‑G‑ACTTGAACCCCCACGTc                                                                <  2:1362212/93‑1 (MQ=37)

                           tGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GAG‑GGGGGGG‑G‑ACTTGAACCCCCACGTcc                                                               >  1:3318119/1‑93 (MQ=255)

                            gAATTTGAATCAGGGATAAAAAATCAAGGAAGAAAAAAGAAAGGAAGAAAAGATAATTGGTGC‑GAG‑GGGGGGG‑G‑ACTTGAACCCCCACGTCCg                                                              <  1:1521119/93‑1 (MQ=37)

                             aaTTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GAG‑GGGGGGG‑G‑ACTTGAACCCCCACGTCCga                                                             >  1:2720389/1‑92 (MQ=255)

                             aaTTTGAAAAAGGGATAAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGGGC‑GAG‑GGGGGGG‑G‑ACTTGAACCCCCACGTCCGt                                                             <  2:5821605/93‑1 (MQ=37)

                                  gAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GAG‑GGGGGGG‑G‑ACTTGAACCCCCACGTCCGTAAGGa                                                        >  2:3234025/1‑93 (MQ=255)

                                   aaTCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GAG‑GGGGGGG‑G‑ACTTGAACCCCCACGTCCGTAAGGac                                                       <  1:6373330/93‑1 (MQ=255)

                                    aTCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GAG‑GGGGGGG‑G‑ACTTGAACCCCCACGTCCGTAAGGaca                                                      >  1:58508/1‑93 (MQ=255)

                                    aTCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GAG‑GGGGGGG‑G‑ACTTGAACCCCCACGTCCGTAACGaaa                                                      >  1:1086173/1‑91 (MQ=37)

                                     tCAGGGATAAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GAG‑GGGGGGG‑G‑ACTTGAACCCCCACGTCCGTAAGGacac                                                     <  1:707747/93‑1 (MQ=255)

                                      aaGGGATGAAAAATAAAGGAAGAAACAAGAAAGGAAGTAAAGATAATGGGTGC‑GAG‑GGGGGGG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACt                                                    <  1:1795442/92‑1 (MQ=25)

                                       aGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GAG‑GGGGGGG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTa                                                   <  2:764915/93‑1 (MQ=255)

                                       aGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GAG‑GGGGGGG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTa                                                   >  2:6342593/1‑93 (MQ=255)

                                        gagATAAAAAAAAAAGAAAAAAACAAAAAAGAAAGAAAAGAAAATTGGTGC‑GA‑‑GGGGGGG‑GTACTTGACCCCCCACGTCCGTAAGGACACTaa                                                  <  2:3615485/91‑1 (MQ=12)

                                          gATGAAAAAACAAAGAAAAAACAAAAAAAGAAAGAAAAGTAATTGGGGCGGG‑‑GGGGGGG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAaca                                                <  2:2142628/93‑1 (MQ=9)

                                           aTGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GAG‑GGGGGGG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAacac                                               >  1:288190/1‑93 (MQ=255)

                                            agAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGGGC‑GA‑‑GGGGGGGAG‑ATTTGACCCCCCACGTCCGTAAGGACACTAACAcc                                              <  1:2913678/92‑1 (MQ=25)

                                             gAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GAG‑GGGGGGG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCt                                             <  2:1807695/93‑1 (MQ=255)

                                             gAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GAG‑GGGGGGG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCt                                             >  1:4773536/1‑93 (MQ=255)

                                               aaaaTCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GAG‑GGGGGGG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGa                                           >  1:6259374/1‑93 (MQ=255)

                                               aaaaTCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GAG‑GGGGGGG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGa                                           >  1:2235758/1‑93 (MQ=255)

                                               aaaaTCAAGGAAGAAAAAAGAAAGGAAGTAAAGATAATTGGTGC‑GAG‑GGGGGGG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGa                                           <  2:4826942/93‑1 (MQ=255)

                                                 aaTCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GAG‑GGGGGGG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAg                                         <  2:6273681/93‑1 (MQ=255)

                                                  aTCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GAG‑GGGGGGG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGc                                        >  1:389031/1‑93 (MQ=255)

                                                  aTCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GAG‑GGGGGGG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGc                                        >  2:2257866/1‑93 (MQ=255)

                                                  aTCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GAG‑GGGGGGG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGc                                        >  1:3018994/1‑93 (MQ=255)

                                                   tCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GAG‑GGGGGGG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCt                                       >  1:5704298/1‑93 (MQ=255)

                                                      aGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GAG‑GGGGGGG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAgc                                    >  1:1905289/1‑93 (MQ=255)

                                                        gaagaaACAAGAAAGGAAGTAAAGATAATTGGTGC‑GAG‑GGGGGGG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAgcgc                                  <  1:4302625/93‑1 (MQ=255)

                                                        gaagaaACAAGAAAGGAAGTAAAGATAATTGGTGC‑GAG‑GGGGGGG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAgcgc                                  <  2:696558/93‑1 (MQ=255)

                                                        gaagaaACAAGAAAGGAAGTAAAGATAATTGGTGC‑GAG‑GGGGGGG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAgcgc                                  <  1:5398441/93‑1 (MQ=255)

                                                           gaaACAAGAAAGGAAGTAAAGATAATTGGTGC‑GAG‑GGGGGGG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTc                               >  2:6412312/1‑93 (MQ=255)

                                                              aCAAGAAAGGAAGTAAAGATAATTGGTGC‑GAG‑GGGGGGG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTAc                            >  1:5583648/1‑93 (MQ=255)

                                                               cAAGAAAGGAAGTAAAGATAATTGGTGC‑GAG‑GGGGGGG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTAcc                           >  2:5482839/1‑93 (MQ=255)

                                                                 aGAAAGGAAGTAAAGATAATTGGTGC‑GAG‑GGGGGGG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCaa                         >  2:902591/1‑93 (MQ=255)

                                                                 aGAAAGGAAGTAAAGATAATTGGTGC‑GAG‑GGGGGGG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCaa                         >  2:4888963/1‑93 (MQ=255)

                                                                   aaaGGAAGTAAAGATATTGGGGGC‑GAG‑GGGGGGG‑G‑ACTTGACCCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAAtt                       <  2:3871146/93‑1 (MQ=14)

                                                                   aaaGGAAGTAAAGATAATTGGTGC‑GAGGGGGGGGG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAAt                        >  2:5398455/1‑93 (MQ=17)

                                                                    aaGGAAGTAAAGATAATTGGTGC‑GAG‑GGGGGGG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTc                      <  2:4076412/93‑1 (MQ=255)

                                                                     aGGAAGTAAAGATAATTGGTGC‑GAG‑GGGGGGG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTcc                     <  1:794915/93‑1 (MQ=255)

                                                                     aGGAAGTAAAGATAATTGGTGC‑GAG‑GGGGGGG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTcc                     >  2:3212116/1‑93 (MQ=255)

                                                                      ggAAGTAAAGATAATTGGTGC‑GAG‑GGGGGGG‑G‑TTTTGAACCCCCACGTCCGTAAGGACACTAACACTTCAAGCTAGCACTGCAACAAATTCCg                    >  2:2772244/1‑93 (MQ=17)

                                                                      ggAAGTAAAGATAATTGGTGC‑GAG‑GGGGGGG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCg                    >  2:5802632/1‑93 (MQ=255)

                                                                        aaGTAAAGAAAGTGGGTGC‑GAG‑GGGGGGG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGcc                  <  2:238830/93‑1 (MQ=11)

                                                                         aGTAAAGATATTTGGGGC‑GA‑‑GGGGGGG‑GAACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCa                 <  2:380502/93‑1 (MQ=14)

                                                                         aGTAAAGATAATTGGTGC‑GAG‑GGGGGGG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCa                 >  1:3614514/1‑93 (MQ=255)

                                                                         aGTAAAGATAATTGGTGC‑GAG‑GGGGGGG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCa                 >  1:377582/1‑93 (MQ=255)

                                                                         aGTAAAGATAATTGGTGC‑GA‑‑GGGGGGG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCAc                >  1:3232352/1‑93 (MQ=255)

                                                                          gTAAAGATAATTGGTGC‑GAG‑GGGGGGG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCAc                >  1:1249089/1‑93 (MQ=255)

                                                                          atAAAGATAATTGGTGC‑GAG‑GGGGGGG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCAc                <  2:5178218/92‑1 (MQ=255)

                                                                            aaaGATAATTGGTGC‑GAG‑GGGGGGG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCt              <  1:858802/93‑1 (MQ=255)

                                                                            aaaGATAATTGGTGC‑GAG‑GGGGGGG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCt              <  2:5896209/93‑1 (MQ=255)

                                                                             aggATAATTGGTGC‑GAG‑GGGGGGG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCtt             <  2:6368625/91‑1 (MQ=255)

                                                                             aaaaTAATTGGGGC‑GAG‑GGGGGGG‑G‑ACTTGACCCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCtt             <  2:3680573/90‑1 (MQ=11)

                                                                             aaGATAATTGGTGC‑GAG‑GGGGGGG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCtt             <  1:2014069/93‑1 (MQ=255)

                                                                             aaGATAATTGGTGC‑GAG‑GGGGGGG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCtt             >  2:2576457/1‑93 (MQ=255)

                                                                             aaGATAATTGGTGC‑GAG‑GGGGGGG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCtt             >  2:3288462/1‑93 (MQ=255)

                                                                             aaGATAATTGGTGC‑GAG‑GGGGGGG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCtt             <  1:1662303/93‑1 (MQ=255)

                                                                             aaGATAATTGGTGC‑GA‑‑GGGGGGG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTc            >  1:6398355/1‑93 (MQ=38)

                                                                              aGATAATTGGTGC‑GAG‑GGGGGGG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTc            >  2:1191502/1‑93 (MQ=255)

                                                                              aGATAATGGGTGC‑GAG‑GGGGGGG‑G‑ACTTGACCCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGCCTACCAATTCCGCCACCTTc            <  1:3652765/93‑1 (MQ=11)

                                                                              aGATAAGTGGTGC‑GAG‑GGGGGGG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTc            <  2:2009735/93‑1 (MQ=255)

                                                                               gATAATTGGTGC‑GAG‑GGGGGGG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGCTTACCAATTCCGCAACCttgc           >  1:5300364/1‑91 (MQ=11)

                                                                                  aaTTGGTGC‑GAG‑GGGGGGG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGcac        >  2:2169227/1‑93 (MQ=16)

                                                                                   aTTGGTGC‑GAG‑GGGGGGG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGcaca       >  1:1013259/1‑93 (MQ=16)

                                                                                    ttGGTGC‑GAG‑GGGGGGG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGcacac      >  1:5265141/1‑92 (MQ=12)

                                                                                    ttGGTGC‑GAG‑GGGGGGG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAg      <  1:5341970/93‑1 (MQ=16)

                                                                                     tgttGC‑GAG‑GGGGGGG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGt     <  2:6358257/90‑1 (MQ=255)

                                                                                        tGC‑GAG‑GGGGGGG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCAt  >  2:1865032/1‑93 (MQ=255)

TTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGC‑GA‑‑GGGGGGG‑G‑ACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCAT  >  NC_000913/4603563‑4603742

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