BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,525,330	T→C	T408T (ACT→ACC)	yrfF →	putative dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,525,330	0	T	C	99.6%	876.7 / NA	271	T408T (ACT→ACC)	yrfF	putative dehydrogenase
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (129/141); total (130/141)
Fisher's exact test for biased strand distribution p-value = 4.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 275 total aligned reads displayed.

GCCACCAGCGTAAAACAACTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACTTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCGCTCATTGCCATT  >  NC_000913/3525241‑3525419

gCCACCAGCGTAAAACAACTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTcc                                                                                        <  1:256340/93‑1 (MQ=255)

gCCACCAGCGTAAAACAACTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTcc                                                                                        >  2:1334901/1‑93 (MQ=255)

 ccaccaGCGTAAAACAACTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCg                                                                                       <  1:1065438/93‑1 (MQ=255)

  caccaGCGTAAAACAACTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCgg                                                                                      >  1:5466538/1‑93 (MQ=255)

  caccaGCGTAAAACAACTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCgg                                                                                      >  2:3130523/1‑93 (MQ=255)

  caccaGCGTAAAACAACTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCgg                                                                                      <  2:1698881/93‑1 (MQ=255)

   accaGCGTAAAACAACTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGa                                                                                     >  1:3732418/1‑93 (MQ=255)

    ccaGCGTAAAACAACTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGaa                                                                                    >  2:4207118/1‑93 (MQ=255)

    ccaGCGTAAAACAACTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGaa                                                                                    >  1:2214073/1‑93 (MQ=255)

      aGCGTAAAACAACTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAAcc                                                                                  <  1:901922/93‑1 (MQ=255)

      aGCGTAAAACAACTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAAcc                                                                                  <  1:2644108/93‑1 (MQ=255)

      aGCATAAAACAACTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAAcc                                                                                  <  1:4162590/93‑1 (MQ=255)

       gCGTAAAACAACTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCt                                                                                 <  1:1563268/93‑1 (MQ=255)

         gTAAAACAACTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTgg                                                                               <  1:5179281/93‑1 (MQ=255)

           aaaaCAACTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAg                                                                             <  2:1922698/93‑1 (MQ=255)

           aaaaCAACTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAg                                                                             <  1:3947413/93‑1 (MQ=255)

            aaaCAACTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAgc                                                                            <  1:4192061/93‑1 (MQ=255)

             aacaacTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAgcg                                                                           <  1:5268196/93‑1 (MQ=255)

              acaacTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAgcgc                                                                          >  2:2723733/1‑93 (MQ=255)

               caacTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGcgcg                                                                         <  1:5677554/93‑1 (MQ=255)

               caacTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGcgcg                                                                         <  2:3924594/93‑1 (MQ=255)

                aacTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGa                                                                        >  1:6440297/1‑93 (MQ=255)

                aacTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGa                                                                        >  1:856399/1‑93 (MQ=255)

                aacTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGa                                                                        >  1:5408866/1‑93 (MQ=255)

                aacTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGa                                                                        >  2:415808/1‑93 (MQ=255)

                 acTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGaa                                                                       >  1:1541511/1‑93 (MQ=255)

                 acTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGaa                                                                       <  2:37226/93‑1 (MQ=255)

                   tGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGaaaa                                                                     >  2:4571427/1‑93 (MQ=255)

                    ggCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAAc                                                                    <  1:3683314/93‑1 (MQ=255)

                     gCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAAcc                                                                   >  1:5907927/1‑93 (MQ=255)

                      cTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCa                                                                  <  1:771874/93‑1 (MQ=255)

                        gACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAAt                                                                >  1:679188/1‑93 (MQ=255)

                         aCGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAAtt                                                               >  1:5839520/1‑93 (MQ=255)

                          cGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTc                                                              >  1:4699327/1‑93 (MQ=255)

                           ccTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGATCCTGGAGCGCGAAAACCAATTCa                                                             <  1:3287385/92‑1 (MQ=255)

                            cTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCAc                                                            <  2:1197628/93‑1 (MQ=255)

                            cTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCAc                                                            >  2:3446725/1‑93 (MQ=255)

                            cTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCAc                                                            >  2:6086163/1‑93 (MQ=255)

                             tGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACTTCCGGAACCTGGAGCGCGAAAACCAATTCAct                                                           >  1:3439123/1‑92 (MQ=255)

                             tGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCAcc                                                           <  1:2302941/93‑1 (MQ=255)

                              ggCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCt                                                          <  1:1036318/93‑1 (MQ=255)

                               gCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCtt                                                         <  1:5811179/93‑1 (MQ=255)

                               gCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCtt                                                         <  2:480994/93‑1 (MQ=255)

                                 gTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCtttt                                                       <  1:5367610/93‑1 (MQ=255)

                                  tGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCttttt                                                      >  1:547306/1‑93 (MQ=255)

                                   gCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCtttttt                                                     <  2:4545173/93‑1 (MQ=255)

                                    cGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTa                                                    >  1:3203614/1‑93 (MQ=255)

                                    cGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTa                                                    <  1:4842151/93‑1 (MQ=255)

                                      ggTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTAcc                                                  <  2:2143588/93‑1 (MQ=255)

                                       gTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCg                                                 <  1:2539684/93‑1 (MQ=255)

                                        tAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGt                                                <  1:3113222/93‑1 (MQ=255)

                                        tAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGt                                                >  2:5744211/1‑93 (MQ=255)

                                        tAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGt                                                <  2:2366424/93‑1 (MQ=255)

                                         aGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGtt                                               <  2:2940459/93‑1 (MQ=255)

                                          ggCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGttt                                              <  1:6153274/93‑1 (MQ=255)

                                             gATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGAc                                           >  2:2778853/1‑93 (MQ=255)

                                             gATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGAc                                           <  1:1268709/93‑1 (MQ=255)

                                             gATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGAc                                           <  2:1133266/93‑1 (MQ=255)

                                             gATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGAc                                           >  1:4900006/1‑93 (MQ=255)

                                             gATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGAc                                           >  2:6731970/1‑93 (MQ=255)

                                              aTACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACt                                          <  2:772097/93‑1 (MQ=255)

                                                aCATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGc                                        >  2:795561/1‑93 (MQ=255)

                                                  aTTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTc                                      <  1:1202372/93‑1 (MQ=255)

                                                   ttGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCg                                     >  1:165120/1‑93 (MQ=255)

                                                    tGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGc                                    <  2:5694598/93‑1 (MQ=255)

                                                    tGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCATTTTTACCGTTTGACTGCTCGc                                    <  2:5366136/93‑1 (MQ=255)

                                                       gTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGa                                 <  1:2009728/93‑1 (MQ=255)

                                                       gTATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGa                                 >  1:5249468/1‑93 (MQ=255)

                                                        tATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGat                                >  2:1731763/1‑93 (MQ=255)

                                                        tATTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGat                                >  2:6264120/1‑93 (MQ=255)

                                                         aTTAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGatc                               <  2:5485493/93‑1 (MQ=255)

                                                          ttAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGACAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGatca                              >  1:4196820/1‑93 (MQ=255)

                                                          ttAGCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGatca                              <  2:261654/93‑1 (MQ=255)

                                                             gCGGTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCt                           <  2:3507834/93‑1 (MQ=255)

                                                               ggTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTgg                         >  2:1761670/1‑93 (MQ=255)

                                                               ggTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTgg                         <  1:5058885/93‑1 (MQ=255)

                                                                gTACGGGAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGa                        >  2:5135500/1‑93 (MQ=255)

                                                                    gggAATGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAacga                    >  2:3169866/1‑93 (MQ=255)

                                                                       aaTGTGTAATATTCGAACCTCCGGCACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGc                 >  2:2059812/1‑93 (MQ=255)

                                                                       aaTGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGc                 <  1:4077912/93‑1 (MQ=255)

                                                                       aaTGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGc                 >  2:3305444/1‑93 (MQ=255)

                                                                        aTGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAGAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGcc                <  1:3515334/93‑1 (MQ=255)

                                                                        aTGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGcc                >  2:1587215/1‑93 (MQ=255)

                                                                        aTGTGTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGcc                >  2:4926889/1‑93 (MQ=255)

                                                                         tgtgTAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGccc               >  2:2512243/1‑93 (MQ=255)

                                                                             tAATATTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCGCTc           >  1:5730169/1‑93 (MQ=255)

                                                                               atatTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCGCTCAt         >  1:4993011/1‑93 (MQ=255)

                                                                               atatTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCGCTCAt         <  1:2787755/93‑1 (MQ=255)

                                                                               atatTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCGCTCAt         >  2:6488818/1‑93 (MQ=255)

                                                                                tatTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCGCTCAtt        >  1:1741997/1‑93 (MQ=255)

                                                                                 atTCGAGCCTCCGGAGCCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCGCTCATTg       <  1:3521714/93‑1 (MQ=255)

                                                                                 atTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCGCTCATTg       >  1:6622125/1‑93 (MQ=255)

                                                                                 atTCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACATTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCGCTCATTg       >  2:931534/1‑93 (MQ=255)

                                                                                   tCGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCGCTCATTGcc     >  2:1005267/1‑93 (MQ=255)

                                                                                    cGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCGCTCATTGCCa    >  1:630651/1‑93 (MQ=255)

                                                                                    cGAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCGCTCATTGCCa    >  2:6517595/1‑93 (MQ=255)

                                                                                     gAACCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCGCTCATTGCCAt   <  1:6226946/93‑1 (MQ=255)

                                                                                      aaCCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCGCTCATTGCCAtt  >  2:4837514/1‑93 (MQ=255)

                                                                                      aaCCTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCGCTCATTGCCAtt  <  1:236585/93‑1 (MQ=255)

                                                                                      aaCCTCCGGAACCTGGAGCGCGAAAACCAATCCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCGCTCATTGCCAtt  >  1:3824855/1‑93 (MQ=255)

GCCACCAGCGTAAAACAACTGGCTGACGCTGGCGTGCGGGTAGGCGATACATTGCGTATTAGCGGTACGGGAATGTGTAATATTCGAACTTCCGGAACCTGGAGCGCGAAAACCAATTCACCTTTTTTACCGTTTGACTGCTCGCAGATCATCTGGAACGACGCCCGCTCATTGCCATT  >  NC_000913/3525241‑3525419

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