BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,755,127	T→C	V469A (GTC→GCC)	pykF →	pyruvate kinase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,755,127	0	T	C	82.6%	479.7 / 99.6	224	V469A (GTC→GCC)	pykF	pyruvate kinase I
Reads supporting (aligned to +/- strand): ref base T (13/26); new base C (70/115); total (83/141)
Fisher's exact test for biased strand distribution p-value = 7.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

Only 100 of 226 total aligned reads displayed.

AGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCAATAAAA  >  NC_000913/1755037‑1755217

agagCGGTCTGGCGCACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGccc                                                                                          <  2:2710857/93‑1 (MQ=255)

 gagCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGCCCt                                                                                         >  2:4221728/1‑93 (MQ=255)

   gCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGCCCTGt                                                                                       <  1:3537405/93‑1 (MQ=255)

    cGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTa                                                                                      <  2:6018973/93‑1 (MQ=255)

    cGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTa                                                                                      <  1:2921508/93‑1 (MQ=255)

     ggTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGCCCTGTaa                                                                                     <  1:363349/93‑1 (MQ=255)

      gTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAat                                                                                    >  2:4470485/1‑93 (MQ=255)

      gTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGCCCTGTAat                                                                                    >  2:3256478/1‑93 (MQ=255)

       tCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGCCCTGTAata                                                                                   <  2:6139189/93‑1 (MQ=255)

        cTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGCCCTGTAatat                                                                                  <  1:4815127/93‑1 (MQ=255)

         tGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATAtt                                                                                 <  1:1527157/93‑1 (MQ=255)

         tGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGCCCTGTAATAtt                                                                                 >  2:4945604/1‑93 (MQ=255)

         tGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGCCCTGTAATAtt                                                                                 <  2:4093291/93‑1 (MQ=255)

         tGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGCCCTGTAATAtt                                                                                 <  2:6902126/93‑1 (MQ=255)

            cacacaAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGCCCTGTAATATTGCt                                                                              >  2:3172154/1‑93 (MQ=255)

                caAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGCCCTGTAATATTGCTTTtg                                                                          >  2:2532532/1‑93 (MQ=255)

                caAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGCCCTGTAATATTGCTTTtg                                                                          <  1:6895396/93‑1 (MQ=255)

                 aaaGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGCCCTGTAATATTGCTTTtgt                                                                         >  2:191221/1‑93 (MQ=255)

                  aaGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGCCCTGTAATATTGCTTTtgtg                                                                        <  2:1765741/93‑1 (MQ=255)

                  aaGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGCCCTGTAATATTGCTTTtgtg                                                                        <  2:1463572/93‑1 (MQ=255)

                     gTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGaat                                                                     <  1:1270884/93‑1 (MQ=255)

                      tGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGCCCTGTAATATTGCTTTTGTGaatt                                                                    >  1:1231636/1‑93 (MQ=255)

                       gtcGTTGTAGTTATGGTTTCTGGTGCATTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGCCCTGTAATATTGCTTTTGTGaatta                                                                   <  1:130518/91‑1 (MQ=255)

                       gACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGCCCTGTAATATTGCTTTTGTGaatta                                                                   >  1:2602018/1‑93 (MQ=255)

                        aCGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGaattaa                                                                  <  2:6334624/93‑1 (MQ=255)

                        aCGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGCCCTGTAATATTGCTTTTGTGaattaa                                                                  >  2:6354161/1‑93 (MQ=255)

                           ttGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAAttt                                                               >  1:3314443/1‑93 (MQ=255)

                            tGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTg                                                              <  1:6567760/93‑1 (MQ=255)

                              tAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGCCCTGTAATATTGCTTTTGTGAATTAATTTGta                                                            <  2:5603674/93‑1 (MQ=255)

                               aGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGCCCTGTAATATTGCTTTTGTGAATTAATTTGtat                                                           <  2:3481274/93‑1 (MQ=255)

                                gTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGtata                                                          <  1:4734435/93‑1 (MQ=255)

                                gTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGCCCTGTAATATTGCTTTTGTGAATTAATTTGtata                                                          <  2:2885311/93‑1 (MQ=255)

                                gTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGCCCTGTAATATTGCTTTTGTGAATTAATTTGtata                                                          >  2:1233709/1‑93 (MQ=255)

                                gTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGCCCTGTAATATTGCTTTTGTGAATTAATTTGtata                                                          <  1:7429145/93‑1 (MQ=255)

                                 ttATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTatat                                                         >  2:401225/1‑93 (MQ=255)

                                  tATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGCCCTGTAATATTGCTTTTGTGAATTAATTTGTATATc                                                        <  2:2783281/93‑1 (MQ=255)

                                   aTGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGCCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCg                                                       <  2:4605721/93‑1 (MQ=255)

                                    tGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGCCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGa                                                      >  2:117828/1‑93 (MQ=255)

                                      gTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGCCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAg                                                    <  1:4223231/93‑1 (MQ=255)

                                       tttCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGCCCTGTAATATTGCTTTTGTGAATTAAGTTGTATATCGAAgc                                                   <  2:5182853/93‑1 (MQ=255)

                                        ttCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGCCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAgcg                                                  <  1:663340/93‑1 (MQ=255)

                                         tCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGCCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAgcgc                                                 >  2:5841453/1‑93 (MQ=255)

                                           tGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGCCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGccc                                               <  2:3833545/93‑1 (MQ=255)

                                           tGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGCCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGccc                                               <  1:5138524/93‑1 (MQ=255)

                                           tGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGCCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGccc                                               <  2:4726547/93‑1 (MQ=255)

                                           tGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGCCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGccc                                               <  1:2179618/93‑1 (MQ=255)

                                           tGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGCCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGccc                                               <  1:7395567/93‑1 (MQ=255)

                                            ggTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGCCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCt                                              <  2:103167/93‑1 (MQ=255)

                                            ggTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGCCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCt                                              <  2:1486272/93‑1 (MQ=255)

                                            ggTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGCCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCt                                              <  1:1766562/93‑1 (MQ=255)

                                              tGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGa                                            <  2:723903/93‑1 (MQ=255)

                                               gCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATCGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGAt                                           <  1:4014351/93‑1 (MQ=255)

                                               gCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGCCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGAt                                           <  1:281760/93‑1 (MQ=255)

                                                   tGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGCCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATAGgc                                       >  2:795057/1‑93 (MQ=255)

                                                      tACCGAGCGGCACTACTAACACCGCATCTGTTCACGCCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCt                                    >  2:6726312/1‑93 (MQ=255)

                                                        ccGAGCGGCACTACTAACACCGCATCTGTTCACGCCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCttt                                  >  2:3775917/1‑93 (MQ=255)

                                                         cGAGCGGCACTACTAACACCGCATCTGTTCACGCCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCtttt                                 >  2:2445273/1‑93 (MQ=255)

                                                         cGAGCGGCACTACTAACACCGCATCTGTTCACGCCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCtttt                                 <  1:2068711/93‑1 (MQ=255)

                                                         cGAGCGGCACTACTAACACCGCATCTGTTCACGCCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCtttt                                 <  1:6024691/93‑1 (MQ=255)

                                                          gAGCGGCACTACTAACACCGCATCTGTTCACGCCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCttttt                                >  1:6737583/1‑93 (MQ=255)

                                                          gAGCGGCACTACTAACACCGCATCTGTTCACGCCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCttttt                                <  2:2046460/93‑1 (MQ=255)

                                                           aGCGGCACTACTAACACCGCATCTGTTCACGCCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCtttttt                               <  1:5842693/93‑1 (MQ=255)

                                                           aGCGGCACTACTAACACCGCATCTGTTCACGCCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCtttttt                               <  1:7235359/93‑1 (MQ=255)

                                                             cGGCACTACTAACACCGCATCTGTTCACGCCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTttta                             <  1:4120980/93‑1 (MQ=255)

                                                              ggCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTtttat                            >  1:3414104/1‑93 (MQ=255)

                                                              ggCACTACTAACACCGCATCTGTTCACGCCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCtttttttgt                            <  2:704409/93‑3 (MQ=255)

                                                               gCACTACTAACACCGCATCTGTTCACGCCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTtttatt                           <  2:7212027/93‑1 (MQ=255)

                                                                cACTACTAACACCGCATCTGTTCACGCCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTtttattt                          <  2:4768076/93‑1 (MQ=255)

                                                                  ctactaACACCGCATCTGTTCACGCCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTaa                        >  1:6692971/1‑93 (MQ=255)

                                                                   tactaACACCGCATCTGTTCACGCCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAAt                       >  1:7030623/1‑93 (MQ=255)

                                                                   tactaACACCGCATCTGTTCACGCCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAAt                       >  1:4860796/1‑93 (MQ=255)

                                                                   tactaACACCGCATCTGTTCACGCCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAAt                       <  2:4657296/93‑1 (MQ=255)

                                                                   tactaACACCGCATCTGTTCACGCCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAAt                       >  2:5282263/1‑93 (MQ=255)

                                                                   tactaACACCGCATCTGTTCACGCCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAAt                       >  1:7070251/1‑93 (MQ=255)

                                                                    actaACACCGCATCTGTTCACGCCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATc                      <  2:1327249/93‑1 (MQ=255)

                                                                     ctaACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCg                     >  2:751345/1‑93 (MQ=255)

                                                                     ctaACACCGCATCTGTTCACGCCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCg                     >  1:7473545/1‑93 (MQ=255)

                                                                      taACACCGCATCTGTTCACGCCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGa                    >  2:846080/1‑93 (MQ=255)

                                                                      taACACCGCATCTGTTCACGCCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGa                    <  1:1635563/93‑1 (MQ=255)

                                                                      taACACCGCATCTGTTCACGCCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGa                    <  2:6631565/93‑1 (MQ=255)

                                                                      taACACCGCATCTGTTCACGCCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGa                    <  2:5775105/93‑1 (MQ=255)

                                                                      taACACCGCATCTGTTCACGCCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGACCTGATGGGCGCTTTTTTTATTTAATCGa                    >  1:1175062/1‑93 (MQ=255)

                                                                       aaCACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGAt                   <  2:5529109/93‑1 (MQ=255)

                                                                        acacCGCATCTGTTCACGCCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATa                  >  2:3033326/1‑93 (MQ=255)

                                                                        acacCGCATCTGTTCACGCCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATa                  <  1:6777787/93‑1 (MQ=255)

                                                                         cacCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATaa                 <  1:3904172/93‑1 (MQ=255)

                                                                         cacCGCATCTGTTCACGCCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATaa                 <  1:6316329/93‑1 (MQ=255)

                                                                            cGCATCTGTTCACGCCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCa              <  1:5581673/93‑1 (MQ=255)

                                                                            cGCATCTGTTCACGCCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCa              <  1:548266/93‑1 (MQ=255)

                                                                              cATCTGTTCACGCCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGa            >  1:181843/1‑93 (MQ=255)

                                                                                  tGTTCACGCCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCa        >  1:2391452/1‑93 (MQ=255)

                                                                                  tGTTCACGCCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCa        >  1:2216188/1‑93 (MQ=255)

                                                                                   gTTCACGCCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCaa       <  1:5734841/93‑1 (MQ=255)

                                                                                   gTTCACGCCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCaa       <  1:3676825/93‑1 (MQ=255)

                                                                                   gTTCACGCCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCaa       <  2:4926063/93‑1 (MQ=255)

                                                                                   gTTCACGCCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCaa       >  2:2155838/1‑93 (MQ=255)

                                                                                    ttCACGCCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCAAt      >  2:419045/1‑93 (MQ=255)

                                                                                    ttCACGCCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCAAt      <  1:7275728/93‑1 (MQ=255)

                                                                                    ttCACGCCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCAAt      <  2:1593084/93‑1 (MQ=255)

                                                                                     tCACGCCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCAATa     >  2:3592072/1‑93 (MQ=255)

AGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTATTTAATCGATAACCAGAAGCAATAAAA  >  NC_000913/1755037‑1755217

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