BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,740,514	T→C	R51R (CGT→CGC)	yiaK →	putative dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,740,514	0	T	C	97.8%	578.8 / ‑1.8	181	R51R (CGT→CGC)	yiaK	putative dehydrogenase
Reads supporting (aligned to +/- strand): ref base T (4/0); new base C (74/103); total (78/103)
Fisher's exact test for biased strand distribution p-value = 3.30e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 190 total aligned reads displayed.

TGACAGCGAAACGGCTGACGCCTGTGCAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGTTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCAATTGAACAGTGGGACGC  >  NC_000913/3740424‑3740603

tGACAGCGAAACGGCTGACGCCTGTGCAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGctt                                                                                         >  2:1086125/1‑93 (MQ=255)

tGACAGCGAAACGGCTGACGCCTGTGCAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGctt                                                                                         <  1:3170998/93‑1 (MQ=255)

 gACAGCGAAACGGCTGACGCCTGTGCAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCttc                                                                                        <  2:447359/93‑1 (MQ=255)

  aCAGCGAAACGGCTGACGCCTGTGCAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGcttat                                                                                       >  2:3534053/1‑91 (MQ=255)

  aCAGCGAAACGGCTGACGCCTGTGCAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCttca                                                                                       <  1:2043546/93‑1 (MQ=255)

       gAAACGGCTGACGCCTGTGCAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTcaa                                                                                  >  2:2511123/1‑93 (MQ=255)

       gAAACGGCTGACGCCTGTGCAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTcaa                                                                                  <  1:514681/93‑1 (MQ=255)

       gAAACGGCTGACGCCTGTGCAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTcaa                                                                                  <  1:482972/93‑1 (MQ=255)

       gAAACGGCTGACGCCTGTGCAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTcaa                                                                                  >  2:4362486/1‑93 (MQ=255)

         aaCGGCTGACGCCTGTGCAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGGTTCCCTCGCTTCATTcaaca                                                                                >  2:3836545/1‑93 (MQ=255)

          aCGGCTGACGCCTGTGCAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTcaacaa                                                                               <  2:4178753/93‑1 (MQ=255)

           cGGCTGACGCCTGTGCAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCaacaac                                                                              <  2:564775/93‑1 (MQ=255)

            ggCTGACGCCTGTGCAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACt                                                                             >  2:1684523/1‑93 (MQ=255)

            ggCTGACGCCTGTGCAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACt                                                                             <  1:640959/93‑1 (MQ=255)

            ggCTGACGCCTGTGCAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACt                                                                             <  1:4090611/93‑1 (MQ=255)

             gCTGACGCCTGTGCAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTg                                                                            >  2:4594188/1‑93 (MQ=255)

              gtGACGCCTGTGCAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTgg                                                                           <  1:2802519/92‑1 (MQ=255)

               tGACGCCTGTGCAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGa                                                                          >  1:2132539/1‑93 (MQ=255)

               tGACGCCTGTGCAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGa                                                                          >  1:1106253/1‑93 (MQ=255)

               tGACGCCTGTGCAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGa                                                                          >  2:2726029/1‑93 (MQ=255)

                gACGCCTGTGCAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGaa                                                                         >  2:2755229/1‑93 (MQ=255)

                          cAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACTACTGGAAAACGGCGata                                                               >  2:1610601/1‑93 (MQ=255)

                           agagATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGatat                                                              >  1:1502118/1‑93 (MQ=255)

                           agagATGTTAGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGatat                                                              <  1:3822496/93‑1 (MQ=255)

                            gagaTGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATc                                                             >  2:2858826/1‑93 (MQ=255)

                            gagaTGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATc                                                             <  2:3422915/93‑1 (MQ=255)

                             agaTGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCa                                                            <  1:4148118/93‑1 (MQ=255)

                             agaTGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCa                                                            <  1:1786550/93‑1 (MQ=255)

                              gaTGTTCGCCCGCACCCCCGAATCCGGCGTTTATGCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCAt                                                           <  1:626907/93‑1 (MQ=255)

                              gaTGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCAt                                                           >  2:1799361/1‑93 (MQ=255)

                              gaTGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCAt                                                           <  1:559106/93‑1 (MQ=255)

                               aTGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCAtt                                                          <  1:771031/93‑1 (MQ=255)

                                tGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTc                                                         <  2:47188/93‑1 (MQ=255)

                                tGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTc                                                         <  2:3298708/93‑1 (MQ=255)

                                 gTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTcc                                                        <  1:3855978/93‑1 (MQ=255)

                                 gTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTcc                                                        <  1:4453255/93‑1 (MQ=255)

                                  ttCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCt                                                       <  1:1141536/93‑1 (MQ=255)

                                     gCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGAt                                                    <  2:2090350/93‑1 (MQ=255)

                                     ccccGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGAt                                                    <  2:1024641/92‑1 (MQ=255)

                                      cccGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACACCTGGAAAACGGCGATATCATTCCTGATg                                                   <  2:2953918/93‑1 (MQ=255)

                                      cccGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATg                                                   <  1:1163105/93‑1 (MQ=255)

                                        cGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGcc                                                 >  1:3150200/1‑93 (MQ=255)

                                        cGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGcc                                                 <  2:1303245/93‑1 (MQ=255)

                                          caccacCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGccca                                               >  1:2423738/1‑93 (MQ=255)

                                           accaccGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGcccaa                                              >  2:677200/1‑93 (MQ=255)

                                             caccGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGcccaacc                                            <  1:4664290/93‑1 (MQ=255)

                                             caccGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGcccaacc                                            <  2:4490417/93‑1 (MQ=255)

                                             caccGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGcccaacc                                            >  1:2826325/1‑93 (MQ=255)

                                             cacaGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGcccaacc                                            >  2:4409563/1‑93 (MQ=255)

                                                 gAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCaaa                                        >  2:1834/1‑93 (MQ=255)

                                                 gAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCaaa                                        >  2:719191/1‑93 (MQ=255)

                                                 gAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCaaa                                        <  1:1015632/93‑1 (MQ=255)

                                                      cGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGtat                                   <  2:114735/93‑1 (MQ=255)

                                                       ggCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGtata                                  >  1:1631973/1‑93 (MQ=255)

                                                       ggCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGtata                                  >  1:1835181/1‑93 (MQ=255)

                                                       ggCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGtata                                  >  1:291483/1‑93 (MQ=255)

                                                        gCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCCAACGTATaa                                 >  2:394215/1‑93 (MQ=255)

                                                        gCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATaa                                 <  1:3455526/93‑1 (MQ=255)

                                                        gCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATaa                                 >  2:4098179/1‑93 (MQ=255)

                                                          gTTTATTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAAcc                               <  2:2967553/93‑1 (MQ=255)

                                                            ttATTCTCACGGCGTTAATCGTTTCCCTCGTTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCac                             >  1:4103447/1‑92 (MQ=255)

                                                              aTTCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACCACTGGAAAGCGGCGATATCATTCCTGATGTCCAACCCAACCGTATATCCAGcc                           >  2:3745945/1‑93 (MQ=255)

                                                               ttCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCt                          <  1:180732/93‑1 (MQ=255)

                                                               ttCTCACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCt                          >  2:3646532/1‑93 (MQ=255)

                                                                tctcACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTc                         <  2:3095692/93‑1 (MQ=255)

                                                                tctcACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTc                         <  1:1049977/93‑1 (MQ=255)

                                                                tctcACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTc                         <  1:4362486/93‑1 (MQ=255)

                                                                 ctcACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCg                        <  2:3357819/93‑1 (MQ=255)

                                                                  tcACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCgg                       >  1:4232165/1‑93 (MQ=255)

                                                                   cACGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGgc                      <  1:3856915/93‑1 (MQ=255)

                                                                    aCGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGgcg                     >  1:2950828/1‑93 (MQ=255)

                                                                     cGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGgcgc                    <  2:2248011/93‑1 (MQ=255)

                                                                     cGGCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGgcgc                    >  2:3582839/1‑93 (MQ=255)

                                                                      ggCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCa                   <  1:1201532/93‑1 (MQ=255)

                                                                       gCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCaa                  <  1:430775/93‑1 (MQ=255)

                                                                       gCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCaa                  <  1:3230946/93‑1 (MQ=255)

                                                                       gCGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCaa                  >  1:1289916/1‑93 (MQ=255)

                                                                        cGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCAAt                 >  2:770514/1‑93 (MQ=255)

                                                                        cGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCAAt                 <  1:980092/93‑1 (MQ=255)

                                                                        cGTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCAAt                 <  1:980093/93‑1 (MQ=255)

                                                                         gTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCAAtt                >  1:3534277/1‑93 (MQ=255)

                                                                         gTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCAAtt                <  2:39175/93‑1 (MQ=255)

                                                                         gTTAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCAAtt                <  2:4033619/93‑1 (MQ=255)

                                                                           tAATCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCAATTGa              <  1:4117641/93‑1 (MQ=255)

                                                                            aaTCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCAATTGaa             <  2:2982845/93‑1 (MQ=255)

                                                                             aTCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCAATTGAAc            <  2:3726348/93‑1 (MQ=255)

                                                                             aTCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCAATTGAAc            <  2:1173647/93‑1 (MQ=255)

                                                                              tCGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCAATTGAACa           >  1:4488557/1‑93 (MQ=255)

                                                                               cGTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCAATTGAACAg          <  1:4377789/93‑1 (MQ=255)

                                                                                gTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCAATTGAACAGt         <  1:4551272/93‑1 (MQ=255)

                                                                                gTTTCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCAATTGAACAGt         <  2:2207439/93‑1 (MQ=255)

                                                                                 tttCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCAATTGAACAGTg        >  2:1463012/1‑93 (MQ=255)

                                                                                 tttCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCAATTGAACAGTg        <  1:3702407/93‑1 (MQ=255)

                                                                                  ttCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCAATTGAACAGTgg       <  1:4250747/93‑1 (MQ=255)

                                                                                   tCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCAATTGAACAGTggg      >  1:4274367/1‑93 (MQ=255)

                                                                                   tCCCTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCAATTGAACAGTggg      <  1:2726197/93‑1 (MQ=255)

                                                                                    cccTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCAATTGAACAGTGGGa     <  2:1956888/93‑1 (MQ=255)

                                                                                    cccTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCAATTGAACAGTGGGa     <  2:826842/93‑1 (MQ=255)

                                                                                     ccTCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCAATTGAACAGTGGGAc    >  1:3021879/1‑93 (MQ=255)

                                                                                       tCGCTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCAATTGAACAGTGGGACGc  <  2:468821/93‑1 (MQ=255)

TGACAGCGAAACGGCTGACGCCTGTGCAGAGATGTTCGCCCGCACCACCGAATCCGGCGTTTATTCTCACGGCGTTAATCGTTTCCCTCGTTTCATTCAACAACTGGAAAACGGCGATATCATTCCTGATGCCCAACCCAAACGTATAACCAGCCTCGGCGCAATTGAACAGTGGGACGC  >  NC_000913/3740424‑3740603

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