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	99.2%	372.3 / ‑5.9	133	noncoding (71/87 nt)	leuP	tRNA‑Leu
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (93/39); total (94/39)
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.96e-01

Only 100 of 155 total aligned reads displayed.

TGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑‑GG‑GGGGGGACT‑TGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAA  >  NC_000913/4603699‑4603874

tGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑‑GG‑GGGGGGGCT‑TGaa                                                                                     >  2:551380/1‑93 (MQ=255)

  aaGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GG‑GGGGGGGCT‑TGAAc                                                                                    >  2:2402959/1‑93 (MQ=255)

   aGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GG‑GGGGGGGCT‑TGAAcc                                                                                   <  1:2841457/93‑1 (MQ=37)

     cTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GG‑GGGGGGGCT‑TGAAcccc                                                                                 <  1:670066/93‑1 (MQ=38)

      tAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GG‑GGGGGGGCT‑TCAAccccc                                                                                >  2:648009/1‑93 (MQ=38)

       aGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GG‑GGGGGGGCT‑TGAACCCCCa                                                                               <  2:695475/93‑1 (MQ=38)

          gcgTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GG‑GGGGGGGCT‑TTAAccccccctc                                                                            >  2:4329815/1‑89 (MQ=38)

          gcgTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GG‑GGGGGGGCT‑TGAACCCCCACgc                                                                            >  1:1036462/1‑92 (MQ=39)

          gcgTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GG‑GGGGGGGCT‑TGAACCCCCACGt                                                                            >  1:3314533/1‑93 (MQ=255)

           cgTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GG‑GGGGGGGCT‑TGAACCCCCACGTc                                                                           <  1:1840375/93‑1 (MQ=39)

           cgTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GG‑GGGGGGGCT‑TGAACCCCCACGTc                                                                           <  1:3559278/93‑1 (MQ=39)

            gTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GG‑GGGGGGGCT‑TGAATCCCCACGTcc                                                                          >  2:3963345/1‑93 (MQ=38)

            gTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GG‑GGGGGGGCT‑TGAACCCCCACGTcc                                                                          <  1:1778288/93‑1 (MQ=39)

            gTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GG‑GGGGGGGCT‑TGAACCCCCACGTcc                                                                          >  2:3089501/1‑93 (MQ=255)

            gTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GG‑GGGGGGGCT‑TGAACCCCCACGTcc                                                                          >  1:97255/1‑93 (MQ=255)

            gTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GG‑GGGGGGGCT‑TGAACCCCCACGTcc                                                                          >  1:983041/1‑93 (MQ=255)

             tCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GG‑GGGGGGGCT‑TGAACCCCCACGTCCg                                                                         >  1:4167522/1‑93 (MQ=255)

             tCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGCTGCGA‑‑GGTGCGGGGGGT‑TGAACCCCCACGTCCg                                                                         >  1:363312/1‑93 (MQ=38)

                cccAACTCCGCCCCCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTGGTTTGGTGCGAG‑GG‑GGGGGGGCT‑TGACCCCCCACGTCCGTaa                                                                      <  2:1284496/92‑1 (MQ=37)

                aCCAATTCCGCCCCCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GG‑GGGGGGGCT‑TGAACCCCCACGTCCGTaa                                                                      <  1:1882379/93‑1 (MQ=38)

                  cAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GG‑GGGGGGGCT‑TGAACCCCCACGTCCGTAAGa                                                                    >  2:2162199/1‑93 (MQ=255)

                  cAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑‑GG‑GGGGGGGCT‑TGAACCCCCACGTCCGTAAGaa                                                                   >  2:1644148/1‑93 (MQ=255)

                   aaTTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GG‑GGGGGGGCT‑TGAACCCCCACGTCCGTAAGaa                                                                   >  1:234276/1‑93 (MQ=255)

                   aaTTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑‑GG‑GGGGGGGCT‑TGAACCCCCACGTCCGTAAGaaa                                                                  >  2:2774573/1‑92 (MQ=255)

                   aaTTCCGCCACCGTCGCACATTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GG‑GGGGGGGCT‑TGAACCCCCACGTCCGTAAGaa                                                                   <  2:4600366/93‑1 (MQ=37)

                     ttCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GG‑GGGGGGGCT‑TGAACCCCCACGTCCGTAAGAaca                                                                 >  1:2232826/1‑93 (MQ=255)

                        cGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GG‑GGGGGGGCT‑TGAACCCCCACGTCCGTAAGAACACTa                                                              >  1:886467/1‑93 (MQ=255)

                         cccACCTTCGCACAGTCATCTTCCTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GG‑GGGGGGGCT‑TGAACCCCCACGTCCGTAAGAACACTaa                                                             <  1:2341994/92‑1 (MQ=37)

                           cACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GG‑GGGGGGGCT‑TGAACCCCCACGTCCGTAAGAACACTAaca                                                           >  2:1651138/1‑93 (MQ=255)

                           cACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GG‑GGGGGGGCT‑TGAACCCCCACGTCCGTAAGAACACTAaca                                                           >  1:1537915/1‑93 (MQ=255)

                            aCCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑‑GG‑GGGGGGGCT‑TGAACCCCCACGTCCGTAAGAACACTAACAcc                                                         >  1:318407/1‑93 (MQ=255)

                              cTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GG‑GGGGGGGCG‑TGAACCCCCACGGCCGTAAGAAGACTAACACCt                                                        >  2:2875128/1‑93 (MQ=25)

                               ttCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GG‑GGGGGGGCT‑TGAACCCCCACGTCCGTAAGAACACTAACACCTg                                                       >  2:2654721/1‑93 (MQ=255)

                                 cGCACAGTCTTTTTTCTTTTTTTGATATCGCCTTGTTTGGTGCGAG‑GG‑GGGGGGGCT‑TGACCCCCCACGTCCGTAAGAACACTAACACCTGaa                                                     <  2:4030810/93‑1 (MQ=18)

                                  gCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GG‑GGGGGGGCT‑TGAACCCCCACGTCCGTAAGAACACTAACACCTGAAg                                                    >  2:4403130/1‑93 (MQ=255)

                                  gCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑‑GG‑GGGGGGGCT‑TGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGc                                                   >  1:848618/1‑93 (MQ=255)

                                    acaGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAC‑GG‑GGGGGGGGTGTGAACCCCCACCTCCGTAAGAACACTAACACCTaaatc                                                   >  1:2516554/1‑91 (MQ=255)

                                     caGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGG‑GGGGGGGCT‑TGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCt                                                  <  1:1209157/93‑1 (MQ=25)

                                     caGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GG‑GGGGGGGCT‑TGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTa                                                 >  1:2884865/1‑93 (MQ=255)

                                     caGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GG‑GGGGGGGAT‑TGAACCCCCACGGCCGTAAGAACACTAACACCTGAAGATa                                                 >  2:1326549/1‑93 (MQ=21)

                                      aGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GG‑GGGGGGGCT‑TGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAg                                                <  2:719635/93‑1 (MQ=255)

                                      aGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GG‑GGGGGGGCT‑TGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAg                                                >  1:4332418/1‑93 (MQ=255)

                                       gTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GG‑GGGGGGGCT‑TGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAgg                                               >  1:512311/1‑92 (MQ=21)

                                       gTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GG‑GGGGGGGCT‑TGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAgc                                               >  2:3176312/1‑93 (MQ=255)

                                        tCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GG‑GGGGGGGCT‑TGAACCCCCACGTCCGGAAGAACACTAACACCTGAAGCTAgcg                                              >  1:4676766/1‑93 (MQ=21)

                                         cATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GG‑GGGGGGGCT‑TGAACCCCCACGTCCGTAAGAAAACTAACACCTGAAGCTaaccg                                             >  1:3837748/1‑89 (MQ=21)

                                          aTCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GG‑GGGGGGGCT‑TGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGcgcg                                            >  1:3352644/1‑93 (MQ=255)

                                          aTCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GG‑GGGGGGGCT‑TGAACCCCCACGTCCGTAAGAACACTAACACCTGAAACTAGcgcg                                            >  2:1537465/1‑93 (MQ=21)

                                           tCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GG‑GGGGGGGCT‑TGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGt                                           >  1:4180217/1‑93 (MQ=255)

                                           tCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GG‑GGGGGGGCT‑TGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGt                                           >  1:3977857/1‑93 (MQ=255)

                                           tCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GG‑GGGGGGGCT‑TGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGt                                           >  1:2385403/1‑93 (MQ=255)

                                            cTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GG‑GGGGGGGCT‑TTAAACCCCACGGCCGCTAAAAAACTTAAAACTGGAGCTAGCGCGTc                                          >  1:1128775/1‑93 (MQ=16)

                                            cTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GG‑GGGGGGGCT‑TGAGCCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTc                                          >  2:826572/1‑93 (MQ=18)

                                             ttACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GG‑GGGGGGGCT‑TGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCt                                         >  1:1434715/1‑93 (MQ=255)

                                             ttACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GG‑GGGGGGGCT‑TGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCt                                         >  2:701578/1‑93 (MQ=255)

                                             ttACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GG‑GGGGGGGCT‑TGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCt                                         >  2:64306/1‑93 (MQ=255)

                                             ttACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GG‑GGGGGGGCT‑TGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCt                                         >  2:2151808/1‑93 (MQ=255)

                                              ttcTTTTTTTGTTATGACGTCGTTCGGTTGGAG‑GG‑GGGGGGGCG‑TGAACCCCCGCGTCCGTAAGAACACTAGCGCCTGAAGCTAGCGCGTCTa                                        <  2:708338/91‑1 (MQ=9)

                                               aCTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GG‑GGGGGGGCT‑TGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTAc                                       >  2:3474269/1‑93 (MQ=255)

                                               aCTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GG‑GGGGGGGCT‑TGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTAc                                       >  2:651661/1‑93 (MQ=255)

                                               aCTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GG‑GGGGGGGCT‑TGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTAc                                       >  2:4523006/1‑93 (MQ=255)

                                                 tttttttGATATCGCCTCGTTTGGTGCGAG‑GG‑GGGGGGGCT‑TGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCa                                     >  2:970469/1‑93 (MQ=255)

                                                 tttttttGATATCGCCTCGTTTGGTGCGAG‑GG‑GGGGGGGCT‑TGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCa                                     >  1:4049608/1‑93 (MQ=255)

                                                      ttGATATCGCCTCGTTTGGTGCGA‑‑GG‑GGGGGGGCT‑TGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCg                               >  2:466404/1‑93 (MQ=255)

                                                       tGATATCGCCTCGTTTGGTGCAA‑‑GG‑GGGGGGACT‑TGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGc                              >  1:794469/1‑93 (MQ=255)

                                                       tGAAATCGCCTCGTTTGGGGCGAG‑GG‑GGGGGGGCT‑TGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCg                               <  2:3674605/93‑1 (MQ=17)

                                                        gATATCGCCTCGTTTGGTGCGAG‑GG‑GGGGGGGCT‑TGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGc                              >  2:1830036/1‑93 (MQ=255)

                                                        gATATCGCCTCGTTTGGTGCGAG‑GG‑GGGGGGGCT‑TGAACCCCCACGCCCGTAAGAACACTAACACCTGAAGCTAGGGCGTCTAACAATTCCGc                              >  2:4674526/1‑93 (MQ=17)

                                                          tatCGCCTCGTTTGGTGCGAGGGG‑GGGGGGGCT‑TGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGcc                             >  1:2273749/1‑93 (MQ=14)

                                                           atCGCCTCGTTTGGTGCGAG‑GG‑GGGGGGGCT‑TGAACCCCCACGTCCGTAAGAACACTAAAACCTAAAGCTAGCGCGCCTACCAATTCCCCCAc                           >  1:1174133/1‑93 (MQ=14)

                                                            tCGCCTCGTTTGGTGCGAG‑GG‑GGGGGGGCT‑TGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCAct                          >  2:3814161/1‑92 (MQ=14)

                                                            tCGCCTCGTTTGGTGCGA‑‑GG‑GGGGGGGCT‑TGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCt                         >  1:2287090/1‑93 (MQ=255)

                                                             cGCCTGGTTGGGTGGGGG‑GG‑GGGGGGGCT‑TGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCt                         <  2:1734414/93‑1 (MQ=11)

                                                             cGCCTCGTTTGGTGCGAG‑GG‑GGGGGGGCT‑TGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCt                         <  1:4398499/93‑1 (MQ=255)

                                                              gCCTCGTTTGGTGCGA‑‑GG‑GGGGGGGCT‑TGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTc                       >  2:2326096/1‑93 (MQ=255)

                                                                cTCGTTTGGTGCGAGGGG‑GGGGGGGCT‑TGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTc                       >  1:435062/1‑93 (MQ=14)

                                                                cTCGTTTGGTGCGAG‑GG‑GGGGGGGCT‑TGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTcc                      >  1:2493300/1‑92 (MQ=14)

                                                                cTCGTTTGGTGCGAG‑GG‑GGGGGGGCT‑TGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCg                      >  1:826919/1‑93 (MQ=255)

                                                                cTCGTTTGGTGCGAG‑GG‑GGGGGGGCT‑TGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCg                      <  2:1579505/93‑1 (MQ=255)

                                                                 tCGTTTGGTGCGAG‑GG‑GGGGGGGCT‑TGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCCCCACCTTCGc                     >  2:1275274/1‑93 (MQ=11)

                                                                   gTTTGGTGCGAG‑GG‑GGGGGGGCT‑TGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGcca                   >  2:1986628/1‑92 (MQ=11)

                                                                    tttGGTGCGAG‑GG‑GGGGGGGCT‑TGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCa                  >  1:3139005/1‑93 (MQ=255)

                                                                    tttGGTGCGAG‑GG‑GGGGGGGCT‑TGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCCTCGCCCa                  >  1:1325805/1‑93 (MQ=11)

                                                                        gTGCGAG‑GG‑GGGGGGGCT‑TGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGc              <  1:3496064/93‑1 (MQ=255)

                                                                        gTGCGA‑‑GG‑GGGGGGGCT‑TGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCg             <  2:1092833/93‑1 (MQ=255)

                                                                         tGCGAG‑GG‑GGGGGGGCT‑TGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCg             >  1:115063/1‑93 (MQ=255)

                                                                          gCGA‑‑GG‑GGGGGGGCT‑TGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAg           >  1:3614943/1‑93 (MQ=255)

                                                                            gggggg‑ggggggGCT‑TGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAg           <  2:4465026/89‑1 (MQ=255)

                                                                             gaggg‑ggggggGCT‑TGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGc          <  2:1346306/90‑1 (MQ=255)

                                                                             gaggg‑ggggggGCT‑TGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGc          >  2:2109710/4‑93 (MQ=255)

                                                                               ggg‑ggggggGCT‑TGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCaa        <  2:115077/92‑1 (MQ=255)

                                                                                gg‑gggtggGCT‑TGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAat       <  2:2830068/93‑1 (MQ=255)

                                                                                 g‑gggggggCT‑TGATCCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAata      <  1:529761/93‑1 (MQ=255)

                                                                                 g‑gggggggCT‑TGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAata      <  1:218795/93‑1 (MQ=255)

                                                                                   ggggggGCT‑TGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAatat     <  2:1551469/93‑1 (MQ=255)

                                                                                    ggggggCT‑TGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAatata    >  1:1476757/1‑92 (MQ=255)

                                                                                    ggggggCT‑TGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATc    <  2:2044792/93‑1 (MQ=255)

                                                                                      ggggCT‑TGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCaa  <  2:4457765/93‑1 (MQ=255)

                                                                                      ggggCT‑TGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCaa  <  1:4685281/93‑1 (MQ=255)

                                                                                      ggggCT‑TGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCaa  <  1:1076363/93‑1 (MQ=255)

TGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑‑GG‑GGGGGGACT‑TGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAA  >  NC_000913/4603699‑4603874

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