BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	493,536	C→T	D79D (GAC→GAT)	ybaB →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	493,536	0	C	T	99.5%	690.4 / ‑5.7	190	D79D (GAC→GAT)	ybaB	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (0/1); new base T (90/99); total (90/100)
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.41e-01

Only 100 of 197 total aligned reads displayed.

ACACAACTGCCGTCGCGTAGAGATCGACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGACGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTT  >  NC_000913/493446‑493625

acacaACTGCCGTCGCGTAGAGATCGACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATgc                                                                                         >  1:4627068/1‑93 (MQ=255)

acacaACTGCCGTCGCGTAGAGATCGACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATgc                                                                                         <  1:4557446/93‑1 (MQ=255)

acacaACTGCCGTCGCGTAGAGATCGACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATgc                                                                                         >  1:3544255/1‑93 (MQ=255)

acacaACTGCCGTCGCGTAGAGATCGACCAGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATgc                                                                                         <  1:581623/93‑1 (MQ=255)

 cacaACTTCCGTCGCGTAGAGATCGACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATgca                                                                                        <  2:3103821/93‑1 (MQ=255)

 cacaACTGCCGTCGCGTAGAGATCGACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATgca                                                                                        >  2:1682854/1‑93 (MQ=255)

 cacaACTGCCGTCGCGTAGAGATCGACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATgca                                                                                        >  2:764734/1‑93 (MQ=255)

 cacaACTGCCGTCGCGTAGAGATCGACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATgca                                                                                        <  1:2577826/93‑1 (MQ=255)

 cacaACTCCCGTCGCGTAGAGATCGCCCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATgca                                                                                        <  2:3263620/93‑1 (MQ=255)

  acaACTGCCGTCGCGTAGAGATCGGCCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATgcag                                                                                       <  1:2894230/93‑1 (MQ=255)

    aaCTGCCGTCGCGTAGAGATCGCCCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATgcagca                                                                                     <  1:2229778/93‑1 (MQ=255)

     aCTGCCGTCGCGTAGAGATCGACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCAc                                                                                    <  2:2314802/93‑1 (MQ=255)

     aCTGCCGTCGCGTAGAGATCGACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCAc                                                                                    <  2:193461/93‑1 (MQ=255)

       tGCCGTCGCGTAGAGATCGACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGCGGCTGCAGCATTCAACGATGCAGCAcgt                                                                                  <  1:3118727/93‑1 (MQ=255)

        gCCGTCGCGTAGAGATCGACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCAcgtc                                                                                 <  2:3900420/93‑1 (MQ=255)

        gCCGTCGCGTAGAGATCGACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCAcgtc                                                                                 <  1:326200/93‑1 (MQ=255)

         ccGGCGCGAAGAGATCAACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCAcgtcg                                                                                <  2:490959/93‑1 (MQ=255)

           gTCGCGTAGAGATCGACCCGAGCCTGCTGGAAGACGACCAAGAGATGCTGGAAGACCTGGTGGCCGCAGCATTCAACGATGCAGCAcgtcgct                                                                              >  2:2010159/1‑91 (MQ=255)

            tCGCGTAGAGATCGACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTAt                                                                             <  1:1376415/93‑1 (MQ=255)

             cgcgTAGAGATCGACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTAtt                                                                            <  2:560373/93‑1 (MQ=255)

              gcgTAGAGATCGACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTg                                                                           >  1:1097042/1‑93 (MQ=255)

               cgTAGCGATCGACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTga                                                                          <  2:3704400/93‑1 (MQ=255)

                gTAGAGATCGACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTgaa                                                                         <  2:1856355/93‑1 (MQ=255)

                gTAGAGATCGACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTgaa                                                                         >  1:3029286/1‑93 (MQ=255)

                 tAGAGATCGACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTgaag                                                                        >  2:3556849/1‑93 (MQ=255)

                 tAGAGATCGACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGCGGCTGCAGCATTCAACGATGCAGCACGTGGTATTgaaa                                                                        >  1:3419888/1‑92 (MQ=255)

                  agagATCGACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTgaaga                                                                       <  1:2414317/93‑1 (MQ=255)

                    agaTCGACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGaaa                                                                     >  1:2706858/1‑93 (MQ=255)

                    agaTCGACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGaaa                                                                     >  1:3124275/1‑93 (MQ=255)

                    agaTCGACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGaaa                                                                     <  1:4653413/93‑1 (MQ=255)

                    agaTCGACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGaaa                                                                     <  2:4400244/93‑1 (MQ=255)

                     gaTCGACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAAc                                                                    >  1:626453/1‑93 (MQ=255)

                     gaTCGACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAAc                                                                    >  1:1835745/1‑93 (MQ=255)

                       tCGACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGc                                                                  <  1:2900185/93‑1 (MQ=255)

                       tCGACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGc                                                                  <  1:4339938/93‑1 (MQ=255)

                          aCCCGAGCCTGCTGGAAGACTACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCaga                                                               >  2:3140624/1‑93 (MQ=255)

                           cccGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCagaa                                                              <  1:2004824/93‑1 (MQ=255)

                            ccGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCagaaa                                                             <  2:406710/93‑1 (MQ=255)

                            ccGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCagaaa                                                             >  1:4219224/1‑93 (MQ=255)

                            ccGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCagaaa                                                             <  1:1876716/93‑1 (MQ=255)

                             cGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCagaaag                                                            <  1:2980487/93‑1 (MQ=255)

                              gAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCagaaaga                                                           <  1:1275440/93‑1 (MQ=255)

                              gAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCagaaaga                                                           <  2:4369011/93‑1 (MQ=255)

                               aGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCagaaagaa                                                          <  2:3511810/93‑1 (MQ=255)

                               aGCCGGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCagaaagaa                                                          <  2:3683751/93‑1 (MQ=255)

                                gCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAgaaagaaa                                                         >  2:2500408/1‑93 (MQ=255)

                                   tgctgGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGaaaaaa                                                      <  2:573400/93‑1 (MQ=255)

                                   tgctgGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGaaaaaa                                                      <  2:1104462/93‑1 (MQ=255)

                                      tgGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATgg                                                   <  2:2959998/93‑1 (MQ=255)

                                        gAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGcc                                                 <  1:4168809/93‑1 (MQ=255)

                                          aGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCctc                                               >  1:4464488/1‑93 (MQ=255)

                                          aGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCctc                                               >  1:3790263/1‑93 (MQ=255)

                                           aacgacAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCctct                                              <  1:229430/92‑1 (MQ=255)

                                            acgacAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTg                                             >  1:2182768/1‑93 (MQ=255)

                                            acgacAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTg                                             >  1:3506926/1‑93 (MQ=255)

                                              gacAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTa                                           >  2:631340/1‑93 (MQ=255)

                                              gacAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTa                                           <  2:1820116/93‑1 (MQ=255)

                                              gacAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTa                                           >  1:1591645/1‑93 (MQ=255)

                                               acAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTAt                                          >  1:3596740/1‑93 (MQ=255)

                                                cAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTAtc                                         <  2:18782/93‑1 (MQ=255)

                                                 aaaGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTAtcc                                        >  2:201614/1‑93 (MQ=255)

                                                  aaGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTAtcct                                       >  2:2220830/1‑93 (MQ=255)

                                                   agagATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTAtcctc                                      >  1:1652732/1‑93 (MQ=255)

                                                    gagaTGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTAtcctcc                                     >  2:2608756/1‑93 (MQ=255)

                                                     agaTGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCg                                    <  1:898030/93‑1 (MQ=255)

                                                       aTGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGa                                  >  2:1526065/1‑93 (MQ=255)

                                                        tGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGaa                                 >  2:1045552/1‑93 (MQ=255)

                                                           tGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGc                              >  2:1547664/1‑93 (MQ=255)

                                                            ggAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCa                             <  2:954219/93‑1 (MQ=255)

                                                            ggAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCa                             >  2:46416/1‑93 (MQ=255)

                                                            ggAAGACCTGGTGGCTGCAGAATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCa                             <  2:777791/93‑1 (MQ=255)

                                                             gAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAg                            >  2:129323/1‑93 (MQ=255)

                                                               aGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCt                          <  1:489673/93‑1 (MQ=255)

                                                                 gccgggtggCTGCAGCATACAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTgc                        <  1:1961944/89‑1 (MQ=255)

                                                                  ccTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTgcc                       <  1:2718138/93‑1 (MQ=255)

                                                                   cTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTgccg                      <  2:161697/93‑1 (MQ=255)

                                                                   cTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTgccg                      >  1:1141588/1‑93 (MQ=255)

                                                                    tggtggCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTGTGTATCCTCCGGAATGCAGCTgccgc                     >  2:3793759/1‑93 (MQ=255)

                                                                     ggtggCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTgccgcc                    <  1:2801178/93‑1 (MQ=255)

                                                                     ggtgcCTGCAGCATACAGCGATGCAGCGCTTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTgccccc                    <  2:2073364/93‑1 (MQ=255)

                                                                      gtggCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGGAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCt                   <  1:4494565/93‑1 (MQ=255)

                                                                      gtggCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGATACGCATTAAGAAAATATGGTCTTTGTATCCTCCGGAATGCAGCTGcctcct                   >  1:4683738/1‑93 (MQ=255)

                                                                      gtggCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCt                   <  1:1708029/93‑1 (MQ=255)

                                                                       tggCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTg                  >  1:411372/1‑93 (MQ=255)

                                                                        ggCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTgg                 <  2:2827757/93‑1 (MQ=255)

                                                                         gCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGc                <  2:1343843/93‑1 (MQ=255)

                                                                          cTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCt               >  1:42587/1‑93 (MQ=255)

                                                                          cTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCt               <  2:1409944/93‑1 (MQ=255)

                                                                          cTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCt               <  1:294569/93‑1 (MQ=255)

                                                                           tGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCtt              <  1:2533147/93‑1 (MQ=255)

                                                                            gcagcaTTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCttt             >  1:386791/1‑93 (MQ=255)

                                                                              agcaTTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTaa           <  1:3161277/93‑1 (MQ=255)

                                                                              agcaTTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTaa           <  2:937942/93‑1 (MQ=255)

                                                                              agcaTTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTaa           <  1:2691629/93‑1 (MQ=255)

                                                                              agcaTTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTaa           >  2:2911725/1‑93 (MQ=255)

                                                                                caTTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGa         <  1:2332174/93‑1 (MQ=255)

                                                                                 aTTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGAt        <  2:3989069/93‑1 (MQ=255)

                                                                                   tCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGc      <  1:3586403/93‑1 (MQ=255)

                                                                                   tCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGc      <  1:816627/93‑1 (MQ=255)

                                                                                      acgacgCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGt   >  2:2706457/1‑93 (MQ=255)

ACACAACTGCCGTCGCGTAGAGATCGACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGACGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTT  >  NC_000913/493446‑493625

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