BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,415,130	C→T	P201S (CCT→TCT)	yjfP →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,415,130	0	C	T	99.5%	753.0 / NA	209	P201S (CCT→TCT)	yjfP	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (0/1); new base T (94/114); total (94/115)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 230 total aligned reads displayed.

CCACTGGCAGAGTGGGAAGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGCCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAGACGGGACGGGATAAACTGCTAACCTGTTCATGGCAGCCAGGCGTGCGT  >  NC_000913/4415040‑4415219

ccACTGGCAGAGTGGGAAGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGtgtct                                                                                         >  1:2551022/1‑93 (MQ=255)

  aCTGGCAGAGTGGGAAGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGtgtctgc                                                                                       <  1:3747151/93‑1 (MQ=255)

   cTGGCAGAGTGGGAAGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGtgtctgcc                                                                                      <  2:4020889/93‑1 (MQ=255)

   cTGGCAGAGTGGAAAGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGtgtctgcc                                                                                      <  2:3490800/93‑1 (MQ=255)

     ggCAGAGTGGGAAGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGTCTGCcgc                                                                                    >  1:4590363/1‑92 (MQ=255)

         gagTGGGAAGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGTCTGCCGACGaa                                                                                >  2:144616/1‑93 (MQ=255)

         gagTGGGAAGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGTCTGCCGACGaa                                                                                >  2:1558183/1‑93 (MQ=255)

         gagTGGGAAGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGTCTGCCGACGaa                                                                                >  1:3695309/1‑93 (MQ=255)

          agTGGGAAGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGTCTGCCGACGAAt                                                                               >  1:3434818/1‑93 (MQ=255)

           gTGGGAAGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGTCTGCCGACGAATc                                                                              <  1:1560194/93‑1 (MQ=255)

            tGGGAAGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGTCTGCCGACGAATCa                                                                             >  2:3111485/1‑93 (MQ=255)

            tGGGAAGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGTCTGCCGACGAATCa                                                                             >  1:2916105/1‑93 (MQ=255)

            tGGGAAGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGTCTGCCGACGAATCa                                                                             >  1:511196/1‑93 (MQ=255)

              ggAAGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGTCTGCCGACGAATCACt                                                                           >  2:2855446/1‑93 (MQ=255)

               gAAGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGTCTGCCGACGAATCACTa                                                                          <  2:899176/93‑1 (MQ=255)

               gAAGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGTCTGCCGACGAATCACTa                                                                          <  2:3291799/93‑1 (MQ=255)

               gAAGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGTCTGCCGACGAATCACTa                                                                          <  1:2112284/93‑1 (MQ=255)

               gAAGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGTCTGCCGACGAATCACTa                                                                          <  2:266647/93‑1 (MQ=255)

                aaGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGTCTGCCGACGAATCACTAc                                                                         <  2:2458067/93‑1 (MQ=255)

                aaGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGTCTGCCGACGAATCACTAc                                                                         <  1:33375/93‑1 (MQ=255)

                 aGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGTCTGCCGACGAATCACTACg                                                                        >  1:4163285/1‑93 (MQ=255)

                  gCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGTCTGCCGACGAATCACTACGt                                                                       <  2:4634035/93‑1 (MQ=255)

                  acGACAAACCACCTGCAACAGCTTAGTGACAGACTTCTACTGCTGTGGCATGGCTTCGACGATGACGTTGTGTCTGGCGACGAATCACTGCGt                                                                       <  2:4360227/92‑1 (MQ=255)

                   cGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGTCTGCCGACGAATCACTACGtt                                                                      <  1:2989069/93‑1 (MQ=255)

                   cGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGTCTGCCGACGAATCACTACGtt                                                                      <  1:3248987/93‑1 (MQ=255)

                   cGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGTCTGCCGACGAATCACTACGtt                                                                      <  2:2073385/93‑1 (MQ=255)

                      cAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGTCTGCCGACGAATCACTACGTTTgc                                                                   <  2:3848759/93‑1 (MQ=255)

                       aaaCCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGTCTGCCGACGAATCACTACGTTTgca                                                                  >  2:4137844/1‑93 (MQ=255)

                        aaCCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGTCTGCCGACGAATCACTACGTTTgcag                                                                 >  1:3970961/1‑93 (MQ=255)

                         accaccTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGTCTGCCGACGAATCACTACGTTTgcagc                                                                >  2:3710472/1‑93 (MQ=255)

                          ccaccTGGAACAACTTAGTGACAGACCTCTACTGCTGTTGCATGGCCTCGACGATGACGTTTTGCCTGCCGACGTATCACTACGTTTgcagca                                                               <  2:3182838/93‑1 (MQ=255)

                          ccaccTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGTCTGCCGACGAATCACTACGTTTgcagca                                                               >  2:2597863/1‑93 (MQ=255)

                          ccaccTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGTCTGCCGACGAATCACTACGTTTgcagca                                                               <  2:3590825/93‑1 (MQ=255)

                          ccaccTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGTCTGCCGACGAATCACTACGTTTgcagca                                                               <  2:2744732/93‑1 (MQ=255)

                           caccTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGTCTGCCGACGAATCACTACGTTTGcagcag                                                              <  1:849365/93‑1 (MQ=255)

                            cccTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGTCTGCCGACGAATCACTACGTTTGCAGCAgg                                                             <  1:1332202/92‑1 (MQ=255)

                            cccTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGTCTGCCGACGAATCACTACGTTTGCAGCAgg                                                             <  2:2565700/92‑1 (MQ=255)

                              cTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGTCTGCCGACGAATCACTACGTTTGCAGCAGGcc                                                           <  2:1787291/93‑1 (MQ=255)

                                ggAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGTCTGCCGACGAATCACTACGTTTGCAGCAGGCCtt                                                         >  1:1070534/1‑93 (MQ=255)

                                 gAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGTCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTa                                                        <  2:78795/93‑1 (MQ=255)

                                  aacaacTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGTCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTaa                                                       <  1:506631/93‑1 (MQ=255)

                                   acaacTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGTCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAg                                                      <  1:1593289/93‑1 (MQ=255)

                                    caacTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGTCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGc                                                     >  2:554735/1‑93 (MQ=255)

                                     aacTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGTCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCg                                                    >  2:4473392/1‑93 (MQ=255)

                                     aacTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGTCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCg                                                    >  2:4286991/1‑93 (MQ=255)

                                      acTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGTCTGCCGACGAATTACTACGTTTGCAGCAGGCCTTAAGCga                                                   >  1:1535504/1‑93 (MQ=255)

                                      acTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGTCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAgaga                                                   >  2:3450380/1‑93 (MQ=255)

                                      acTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGTCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCga                                                   <  1:1500593/93‑1 (MQ=255)

                                      acTTAGTGACAGACCTCTACTGCTGTGGCATGGCCCCGACGATGACGTTGTGTCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCga                                                   >  1:1712924/1‑93 (MQ=255)

                                       cTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGTCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCgag                                                  >  2:697450/1‑93 (MQ=255)

                                        ttAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGTCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCgaga                                                 <  2:2247905/93‑1 (MQ=255)

                                         tAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGTCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAgac                                                <  1:2375629/93‑1 (MQ=255)

                                          aGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGTCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAgacg                                               <  1:3527934/93‑1 (MQ=255)

                                           gTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGTCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAgacgg                                              <  1:1452316/93‑1 (MQ=255)

                                             gACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGTCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAgacggga                                            <  1:3261736/93‑1 (MQ=255)

                                             gACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGTCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAgacggga                                            <  1:3891380/93‑1 (MQ=255)

                                               cAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGTCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAgacgggacg                                          >  1:383968/1‑93 (MQ=255)

                                                aGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGTCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAgacgggacgg                                         <  2:735591/93‑1 (MQ=255)

                                                 gACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGTCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAGacgggacggg                                        >  1:3188400/1‑93 (MQ=255)

                                                 gACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGTCTGCCGACGAATCACTACGTTTGAAGCAGGCCTTAAGCGAGacgggacggg                                        >  1:2433535/1‑93 (MQ=255)

                                                   ccTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGTCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAGACGGGACGGGAt                                      <  2:4442032/93‑1 (MQ=255)

                                                   ccTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGTCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAGACGGGACGGGAt                                      <  1:4339960/93‑1 (MQ=255)

                                                    ctctACTGCTGTGGCATGGCCTCGACGATGACGTTGTGTCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAGACGGGACGGGATa                                     >  2:1842803/1‑93 (MQ=255)

                                                    ctctACTGCTGTGGCATGGCCTCGACGATGACGTTGTGTCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAGACGGGACGGGATa                                     >  2:2341324/1‑93 (MQ=255)

                                                     tctACTGCTGTGGCATGGCCTCGACGATGACGTTGTGTCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAGACGGGACGGGATaa                                    <  2:1678793/93‑1 (MQ=255)

                                                     cctACTGCTGTGGCATGGCCTCGACGATGACGTTGTGTCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAGACGGGACGGGATaa                                    <  2:228777/92‑1 (MQ=255)

                                                      ctACTGCTGTGGCATGGCCTCGACGATGACGTTGTGTCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAGACGGGACGGGATaaa                                   <  1:953548/93‑1 (MQ=255)

                                                      ctACTGCTGTGGCATGGCCTCGACGATGACGTTGTGTCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAGACGGGACGGGATaaa                                   <  1:3623444/93‑1 (MQ=255)

                                                       tACTGCTGTGGCATGGCCTCGACGATGACGTTGTGTCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAGACGGGACGGGATAAAc                                  >  2:1832706/1‑93 (MQ=255)

                                                       tACTGCTGTGGCATGGCCTCGACGATGACGTTGTGTCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAGACGGGACGGGATAAAc                                  >  2:2087435/1‑93 (MQ=255)

                                                       tACTGCTGTGGCATGGCCTCGACGATGACGTTGTGTCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAGACGGGACGGGATAAAc                                  <  2:138583/93‑1 (MQ=255)

                                                        aCTGCTGTGGCATGGCCTCGACGATGACGTTGTGTCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAGACGGGACGGGATAAACt                                 <  1:2308468/93‑1 (MQ=255)

                                                        aCTGCTGTGGCATGGCCTCGACGATGACGTTGTGTCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAGACGGGACGGGATAAACt                                 <  1:3622371/93‑1 (MQ=255)

                                                           gctgTGGCATGGCCTCGACGATGACGTTGTGTCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAGACGGGACGGGATAAACTGCt                              <  2:1259925/93‑1 (MQ=255)

                                                             tgtgGCATGGCCTCGACGATGACGTTGTGTCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAGACGGGACGGGATAAACTGCTaa                            >  1:723059/1‑93 (MQ=255)

                                                               tgGCATGGCCTCGACGATGACGTTGTGTCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAGACGGGACGGGATAAACTGCTAAcc                          <  2:3664854/93‑1 (MQ=255)

                                                               tgGCATGGCCTCGACGATGACGTTGTGTCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAGACGGGACGGGATAAACTGCTAAcc                          <  2:1915224/93‑1 (MQ=255)

                                                                 gCATGGCCTCGACGATGACGTTGTGTCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAGACGGGACGGGATAAACTGCTAACCTg                        <  1:1868682/93‑1 (MQ=255)

                                                                  gaTGGCCTCGACGATGACGTTGTGTCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAGACGGGACGGGATAAACTGCTAACCTGt                       <  1:1224213/92‑1 (MQ=255)

                                                                   aTGGCCTCGACGATGACGTTGTGTCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAGACGGGACGGGATAAACTGCTAACCTGtt                      >  2:925545/1‑93 (MQ=255)

                                                                    tGGCCTCGACGATGACGTTGTGTCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAGACGGGACGGGATAAACTGCTAACCTGTTc                     >  2:2631903/1‑93 (MQ=255)

                                                                    tGGCCTCGACGATGACGTTGTGTCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAGACGGGACGGGATAAACTGCTAACCTGTTc                     >  2:1603453/1‑93 (MQ=255)

                                                                     ggCCTCGACGATGACGTTGTGTCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAGACGGGACGGGATAAACTGCTAACCTGTTCa                    >  2:3907260/1‑93 (MQ=255)

                                                                      gCCTCGACGATGACGTTGTGTCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAGACGGGACGGGATAAACTGCTAACCTGTTCAt                   <  2:723685/93‑1 (MQ=255)

                                                                       ccTCGACGATGACGTTGTGTCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAGACGGGACGGGATAAACTGCTAACCTGTTCATg                  <  2:2467369/93‑1 (MQ=255)

                                                                       ccTCGACGATGACGTTGTGTCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAGACGGGACGGGATAAACTGCTAACCTGTTCATg                  >  2:1175213/1‑93 (MQ=255)

                                                                       ccTCGACGATGACGTTGTGTCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAGACGGGACGGGATAAACTGCTAACCTGTTCATg                  <  1:4467464/93‑1 (MQ=255)

                                                                          cgacgaTGACGTTGTGTCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAGACGGGACGGGATAAACTGCTAACCTGTTCATGGCa               >  2:2712198/1‑93 (MQ=255)

                                                                              gaTGACGTTGTGTCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAGACGGGACGGGATAAACTGCTAACCTGTTCATGGCAGCCa           >  2:1125806/1‑93 (MQ=255)

                                                                              gaTGACGTTGTGTCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAGACGGGACGGGATAAACTGCTAACCTGTTCATGGCAGCCa           >  1:4269851/1‑93 (MQ=255)

                                                                               aTGACGTTGTGTCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAGACGGGACGGGATAAACTGCTAACCTGTTCATGGCAGCCAg          >  1:2498505/1‑93 (MQ=255)

                                                                                tGACGTTGTGTCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAGACGGGACGGGATAAACTGCTAACCTGTTCATGGCAGCCAgg         >  1:2684059/1‑93 (MQ=255)

                                                                                tGACGTTGTGTCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAGACGGGACGGGATAAACTGCTAACCTGTTCATGGCAGCCAgg         >  2:4648798/1‑93 (MQ=255)

                                                                                 gACGTTGTGTCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAGACGGGACGGGATAAACTGCTAACCTGTTCATGGCAGCCAGgc        <  1:3160269/93‑1 (MQ=255)

                                                                                  aCGTTGTGTCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAGACGGGACGGGATAAACTGCTAACCTGTTCATGGCAGCCAGgcg       <  2:4586458/93‑1 (MQ=255)

                                                                                  aCGTTGTGTCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAGACGGGACGGGATAAACTGCTAACCTGTTCATGGCAGCCAGgcg       >  1:4105953/1‑93 (MQ=255)

                                                                                    gTTGTGTCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAGACGGGACGGGATAAACTGCTAACCTGTTCATGGCAGCCAGgcgtg     >  1:4177852/1‑93 (MQ=255)

                                                                                     ttGTGTCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAGACCGGACGGGATAAACTGCTAACCTGTTCATGGCAGCCAGgcgtgc    >  1:4611623/1‑93 (MQ=255)

                                                                                       gtgtCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAGACGGGACGGGATAAACTGCTAACCTGTTCATGGCAGCCAGgcgtgcgt  >  2:1314399/1‑93 (MQ=255)

                                                                                       gtgtCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAGACGGGACGGGATAAACTGCTAACCTGTTCATGGCAGCCAGgcgggcgt  >  1:1274915/1‑93 (MQ=255)

CCACTGGCAGAGTGGGAAGCGACAAACCACCTGGAACAACTTAGTGACAGACCTCTACTGCTGTGGCATGGCCTCGACGATGACGTTGTGCCTGCCGACGAATCACTACGTTTGCAGCAGGCCTTAAGCGAGACGGGACGGGATAAACTGCTAACCTGTTCATGGCAGCCAGGCGTGCGT  >  NC_000913/4415040‑4415219

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