BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,603,785	A→G	noncoding (71/87 nt)	leuP ←	tRNA‑Leu

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,603,785	0	A	G	96.8%	257.1 / ‑1.5	97	noncoding (71/87 nt)	leuP	tRNA‑Leu
Reads supporting (aligned to +/- strand): ref base A (2/1); new base G (57/33); total (60/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

Only 100 of 112 total aligned reads displayed.

ACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GA‑‑GGGGGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAAC  >  NC_000913/4603696‑4603875

        cTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GA‑‑GGGGGGG‑GGCTTGAAccccc                                                                                 <  1:3098499/93‑1 (MQ=255)

         tAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAGGGGGGGGG‑GGCTTGAAcccc                                                                                  >  1:2590302/1‑93 (MQ=38)

         tAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GGGGGGG‑GGCTTGAAccccc                                                                                 >  1:3568740/1‑93 (MQ=255)

          accgccTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GGGGGGG‑GGCTTGAACCCCCa                                                                                <  1:3676543/91‑1 (MQ=38)

          aGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GGGGGGG‑GGCTTGAAcccccc                                                                                >  2:2569818/1‑92 (MQ=38)

          aGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GGGGGGG‑GGCTTGAACCCCCa                                                                                <  2:3829325/93‑1 (MQ=38)

           gcgcGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAGGGGGGGGG‑GGCTTGAAcccccc                                                                                >  2:3780433/1‑92 (MQ=38)

           gcgcGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GGGGGGG‑GGCTTGAAccccccc                                                                               >  2:4389187/1‑91 (MQ=38)

            cgcgTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GGGGGGG‑GGCTTGAACCCCCACg                                                                              >  1:1256454/1‑93 (MQ=255)

             gcgTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GGGGGGG‑GGCTTGAAccccccggc                                                                             >  2:3673879/1‑89 (MQ=38)

             gcgTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GGGGGGG‑GGCTTGAACCCCCCCgc                                                                             >  2:1758630/1‑92 (MQ=38)

              cgTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GGGGGGG‑GGCTTGAACCCCCACGTc                                                                            >  1:2739044/1‑93 (MQ=255)

               gTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GA‑‑GGGGGGG‑GGCTTGAACCCCCACGTCCg                                                                          >  1:299559/1‑93 (MQ=255)

               gTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GA‑‑GGGGGGG‑GGCTTGAACCCCCACGTCCg                                                                          <  1:571210/93‑1 (MQ=255)

                  tACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GGGGGGG‑GGCTTGAACCCCCACGTCCGTa                                                                        >  1:2843826/1‑93 (MQ=255)

                  tACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GGGGGGG‑GGCTTGAACCCCCACGTCCGTa                                                                        >  1:4355969/1‑93 (MQ=255)

                  tACAAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GGGGGGG‑GGCTTGAACCCCCACGTCCGTa                                                                        <  2:2931224/93‑1 (MQ=38)

                     cAATTCCGCCACCTTCGCACAGTCATCTTCCTTTTTTTCATATCGCCTCGTTTGGTGC‑GA‑‑GGGGGGG‑G‑CTTGAACCCCCACGTCCGTAAGAac                                                                   <  2:197084/93‑1 (MQ=37)

                      aaTTCCGCCACCTGCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GGGGGGG‑GGCTTGAACCCCCACGTCCGTAAGaa                                                                    <  2:1654088/93‑1 (MQ=37)

                       aTTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GGGGGGG‑GGGTTGAACCCCCCCGGCCCTAAGaaa                                                                   >  2:3855958/1‑92 (MQ=255)

                        ttCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GGGGGGG‑GGGTTGAACCCCCACCGCCCGAAGAaca                                                                  >  2:2473454/1‑93 (MQ=255)

                        ttCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GGGGGGG‑GGCTTGAACCCCCACGTCCGTAAGAaca                                                                  <  1:2425974/93‑1 (MQ=38)

                              cACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GGGGGGG‑GGCTTGAACCCCCACGTCCGTAAGAACACTAaca                                                            >  1:4440171/1‑93 (MQ=255)

                              cACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GGGGGGG‑GGCTTGAACCCCCACGTCCGTAAGAACACTAaca                                                            >  1:515131/1‑93 (MQ=255)

                              cACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GGGGGGG‑GGCTTGAACCCCCACGTCCGTAAGAACACTAaca                                                            >  2:672505/1‑93 (MQ=255)

                               aCCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GGGGGGG‑GGCTTGAACCCCCACGTCCGTAAGAACAATAacac                                                           >  2:831195/1‑93 (MQ=37)

                               aCCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GA‑‑GGGGGGG‑GACTTGACCCCCCACGTCCGTAAGAACACTAACAcc                                                          <  1:3789699/93‑1 (MQ=255)

                                ccTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GGGGGGG‑GGCTTGAACCCCCACGTCCGTAAGAACACTAACAcc                                                          >  1:2501110/1‑93 (MQ=255)

                                 cTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GGGGGGG‑GGCTTGAACCCCCACGTCCGTAAGAACACTAACACCt                                                         >  2:4456592/1‑93 (MQ=255)

                                  ttCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GGGGGGG‑GGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTg                                                        >  1:3403817/1‑93 (MQ=21)

                                  ttCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GGGGGGG‑GGCTTGAACCCCCACGTCCGCACGAATACTAACACCTg                                                        >  2:3284953/1‑93 (MQ=25)

                                  ttCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GA‑‑GGGGGGG‑GGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGa                                                       <  2:2810135/93‑1 (MQ=255)

                                  ctCGCACAGTCATCTTACTTTTGTTGATATCGCCTCGTTTGGTGC‑GA‑‑GGGGGGG‑G‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGaa                                                      <  1:3029175/92‑1 (MQ=255)

                                   tCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GGGGGGG‑GGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGa                                                       >  1:1482726/1‑93 (MQ=255)

                                        caGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GA‑‑GGGGGGG‑GGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAg                                                 <  2:1690844/93‑1 (MQ=255)

                                        caGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GA‑‑GGGGGGG‑GGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAg                                                 <  1:1831355/93‑1 (MQ=255)

                                         aGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GGGGGGG‑GGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAg                                                 >  2:3461533/1‑93 (MQ=255)

                                            cATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GGGGGGG‑GGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAgcgc                                              <  1:4032663/93‑1 (MQ=255)

                                            cATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GGGGGGG‑GGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAgcgc                                              >  2:1616656/1‑93 (MQ=255)

                                              tCTTCCTTTTTTTGATATCCCCTTGTTTGGTGC‑GA‑‑GGGGGGG‑GGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTc                                           <  2:602757/93‑1 (MQ=17)

                                              tCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GA‑‑GGGGGGG‑GGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTc                                           >  2:535779/1‑93 (MQ=255)

                                               cTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GGGGGGG‑GGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTc                                           >  1:2283771/1‑93 (MQ=255)

                                               cTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GGGGGGG‑GGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTc                                           >  1:2543237/1‑93 (MQ=255)

                                               cTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GGGGGGG‑GGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTc                                           <  1:2537180/93‑1 (MQ=255)

                                                ttACTTTTTTTGATATCGCCTCGTTTGGTGC‑GAGGGGGGGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTc                                           >  1:3582801/1‑93 (MQ=18)

                                                 tACTTTTTTTGATATCGCCTCGTTTGGTGC‑GA‑‑GGGGGGG‑GGCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTAc                                        >  1:1832867/1‑93 (MQ=255)

                                                 tACTTTTTTTGATATCGCCTCGTTTGGTGC‑GA‑‑GGGGGGG‑GGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTAc                                        <  1:327789/93‑1 (MQ=255)

                                                  aCTTTTTTTGATATCGCCTCGTTTGGTGC‑GAG‑GGGGGGG‑GGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTAc                                        >  1:4640917/1‑93 (MQ=255)

                                                  aCTTTTTTTGATATCGCCTCGTTTGGTGC‑GA‑‑GGGGGGG‑AGCTTGCACCCCCACGCCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTAcc                                       >  2:2243228/1‑93 (MQ=18)

                                                   cTTTTTTTGATATCGCCTCGTTTGGTGC‑GA‑‑GGGGGGG‑GGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCa                                      <  2:596122/93‑1 (MQ=255)

                                                    tttttttGATATCGCCTCGTTTGGTGC‑GAG‑GGGGGGG‑GGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCa                                      >  1:4068278/1‑93 (MQ=255)

                                                    tttttttGATATCGCCTCGTTTGGTGC‑GAG‑GGGGGGG‑GGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCa                                      >  1:1302093/1‑93 (MQ=255)

                                                    tttttttGATATCGCCTCGTTTGGTGC‑GA‑‑GGGGGGG‑GGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCaa                                     >  2:1064065/1‑93 (MQ=255)

                                                     ttttttGATATCGCCTCGTTTGGTGC‑GAG‑GGGGGGG‑GGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCaa                                     >  2:2709980/1‑93 (MQ=255)

                                                      tttttGATATCGCCTCGTTTGGTGC‑GA‑‑GGGGGGG‑GGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAAtt                                   <  2:3071549/93‑1 (MQ=255)

                                                      tttttGATAACGCCTCGTTTGGTGC‑GAG‑GGGGGGG‑GGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAAt                                    <  1:3270912/93‑1 (MQ=17)

                                                        tttGATATCGCCTCGTTTGGTGC‑GAG‑GGGGGGG‑GGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTc                                  <  1:1828931/93‑1 (MQ=255)

                                                        tttGATATCGCCTCGTTTGGTGC‑GAG‑GGGGGGG‑GGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTc                                  >  1:2293817/1‑93 (MQ=255)

                                                         ttGATATCGCCTCGTTTGGTGC‑GAG‑GGGGGGG‑GGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTcc                                 <  1:574775/93‑1 (MQ=255)

                                                           gATATCGCCTCGTTTGGTGC‑GAG‑GGGGGGG‑GGCTTGAACCCCCCCGCCCGTAAAACAATTAAAACCGAAACCTACGGGGCCTACCAATTcccc                               >  2:3760651/1‑91 (MQ=12)

                                                             tatCGCCTCGTTTGGTGC‑GAGGGGGGGGG‑GTCTTGACCCCCCACGTCCGTAAGACCACTAAAACCTGAAGCTAGCGCGTCTACCAATTCCGcc                              >  1:630903/1‑93 (MQ=14)

                                                              atCGCCTCGTTTGGTGC‑GAG‑GGGGGGG‑GGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCAc                            >  2:3557263/1‑93 (MQ=255)

                                                              atCGCCTCGTTTGGTGC‑GAG‑GGGGGGG‑GGCTTGAACCCCCACGTCCGTAAGAAAACTAACACCTGAGGTTAGCGGGTCTACCAATTACGCCAc                            >  1:1893547/1‑93 (MQ=14)

                                                               tCGCCTCGTTTGGTGC‑GAG‑GGGGGGG‑GGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCAcc                           >  1:4660078/1‑93 (MQ=255)

                                                               tCGCCTCGTTTGGTGC‑GAG‑GGGGGGG‑GGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCAcc                           >  2:1218173/1‑93 (MQ=255)

                                                               tCGCCTCGTTTGGTGC‑GAG‑GGGGGGG‑GGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAACTAGCGCCTCTACCAATTCCCCCAcc                           >  1:2612108/1‑93 (MQ=14)

                                                               tCGCCTCGTTTGGTGC‑GAG‑GGGGGGG‑GGCTTGAACCCCCACGTCCGTAAGAAAACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCAcc                           >  2:4651057/1‑93 (MQ=14)

                                                                cGCCTTGTTTGGTGC‑GA‑‑GGGGGGG‑GGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCtt                         <  2:4710870/93‑1 (MQ=255)

                                                                cGCCTCGTTTGGTGC‑GAG‑GGGGGGG‑GGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCt                          <  2:4528659/93‑1 (MQ=255)

                                                                 gCCTCGTTTGGTGC‑GA‑‑GGGGGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTc                        >  2:4137916/1‑93 (MQ=40)

                                                                  ccTCGTTTGGTGC‑GAG‑GGGGGGG‑GGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTc                        >  1:1708487/1‑93 (MQ=255)

                                                                  ccTCGTTTGGTGC‑GAG‑GGGGGGG‑GGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTc                        >  2:2239071/1‑93 (MQ=255)

                                                                  ccTCGTTTGGTGC‑GA‑‑GGGGGGG‑GGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCg                       >  2:346246/1‑93 (MQ=255)

                                                                   cTCGTTTGGTGC‑GAGGGGGGGGG‑GGCTTGAACCCCCCCGTCCGTAAAAACACAAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCCCCTTc                        >  2:3582200/1‑93 (MQ=14)

                                                                   cTCGTTTGGTGC‑GAG‑GGGGGGG‑GGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCg                       <  1:611614/93‑1 (MQ=255)

                                                                   cTCGTTTGGTGC‑GAG‑GGGGGGG‑GGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCg                       <  1:2464417/93‑1 (MQ=255)

                                                                      gTTTGGTGC‑GAG‑GGGGGGG‑GGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGccc                    >  2:1827168/1‑93 (MQ=255)

                                                                      gTTTGGTGC‑GAG‑GGGGGGG‑GGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGccc                    <  2:719018/93‑1 (MQ=255)

                                                                      gTTTGGTGC‑GAG‑GGGGGGG‑GGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGcat                    >  1:4061083/1‑91 (MQ=11)

                                                                       tttGGTGCGGG‑‑GGGGGGG‑GGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCa                   <  2:2717891/93‑1 (MQ=11)

                                                                       tttGGTGC‑GAG‑GGGGGGG‑GGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGcata                   >  1:450074/1‑90 (MQ=1)

                                                                       tttGGTGC‑GAG‑GGGGGGG‑GGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCa                   <  1:780146/93‑1 (MQ=255)

                                                                              c‑GA‑‑GGGGGGGTGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAg            <  2:2529680/93‑1 (MQ=255)

                                                                              c‑GA‑‑GGGGGGGTGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAg            <  2:2529681/93‑1 (MQ=255)

                                                                              c‑GA‑‑GGGGGGG‑GGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGc           >  2:808297/1‑93 (MQ=255)

                                                                                ggggggggggg‑gCCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAg            <  2:2063281/89‑1 (MQ=255)

                                                                                 gagggggggg‑gGCTTGAACCCCCACGTCAGTAAGAACACTAACTCCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGc           <  1:342505/90‑1 (MQ=14)

                                                                                 a‑‑GGGGGGG‑GGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCaa         <  2:2936154/93‑1 (MQ=255)

                                                                                  ggggggggg‑gGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCa          >  2:1943320/3‑93 (MQ=255)

                                                                                  agggggggg‑gGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCa          >  2:178169/3‑93 (MQ=255)

                                                                                  agggggggg‑gGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGCGCGAgaa          >  2:4729402/3‑91 (MQ=14)

                                                                                    ggggggg‑gGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCact        >  2:2628115/1‑91 (MQ=255)

                                                                                     gggggg‑ggGTGGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAata       <  1:2683819/93‑1 (MQ=255)

                                                                                     gggggg‑ggCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAata       <  2:673137/93‑1 (MQ=255)

                                                                                      ggggg‑gGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAatat      >  2:3268625/1‑93 (MQ=255)

                                                                                      ggggg‑gGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAatat      <  1:756368/93‑1 (MQ=255)

                                                                                      ggggg‑gGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAatat      >  1:2699043/1‑93 (MQ=255)

                                                                                        ggg‑gGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCa    <  1:1065894/93‑1 (MQ=255)

                                                                                         gg‑ggCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCaa   >  1:2021056/1‑93 (MQ=255)

                                                                                           gcaCTTGACCCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAAc  <  2:4018457/91‑1 (MQ=255)

ACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑GA‑‑GGGGGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAAC  >  NC_000913/4603696‑4603875

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