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.3%	1034.9 / ‑4.9	289	D79D (GAC→GAT)	ybaB	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (2/0); new base T (128/159); total (130/159)
Fisher's exact test for biased strand distribution p-value = 2.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.62e-01

Only 100 of 296 total aligned reads displayed.

ACACAACTGCCGTCGCGTAGAGATCGACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGACGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCG  >  NC_000913/493446‑493623

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

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

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

  acaACTGCCGTCGCGTAGAGATCGACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGACgcag                                                                                     >  2:3513828/1‑93 (MQ=255)

  acaACTGCCGTCGCGGAGAGGTCGACCCGAGCCTGCTGGAAGGCGCCAAAGAGATGCTGGAAGACCTGGTGGCTGCAGTCTTCAACGATgcac                                                                                     >  2:7770413/1‑92 (MQ=255)

   caACTGCCGTCGCGTAGAGATCGACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATgcagc                                                                                    >  1:2549498/1‑93 (MQ=255)

    aaCTGCCGTCGCGTAGAGATCGACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATgcagca                                                                                   <  2:3215769/93‑1 (MQ=255)

    aaCTGCCGTCGCGTAGAGATCGACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATgcagca                                                                                   <  1:3910045/93‑1 (MQ=255)

     aCTGCCGTCGCGTAGAGATCGACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCAc                                                                                  >  1:6453295/1‑93 (MQ=255)

       tGCCGTCGCGTAGAGATCGACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGACTGCAGCATTCAACGATGCAGCAcgt                                                                                <  2:5661760/93‑1 (MQ=255)

       cgCCGTCGCGTAGAGATCGACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCAcgt                                                                                <  1:6163294/92‑1 (MQ=255)

           gTCGCGTAGAGATCGACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGGAGCATTCAACGATGCAGCACGTCGTa                                                                            <  1:936125/93‑1 (MQ=255)

           gTCGCGTAGAGATCGACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTa                                                                            >  1:5682665/1‑93 (MQ=255)

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

            tCGCGTAGAGATCGACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTAt                                                                           >  2:3861606/1‑93 (MQ=255)

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

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

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

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

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

              gcgTAGAGATCGACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTg                                                                         <  1:403504/93‑1 (MQ=255)

               cgTAGAGATCGACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTga                                                                        >  1:4161587/1‑93 (MQ=255)

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

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

                 tAGAGATCGACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTgaaa                                                                      >  2:2263643/1‑92 (MQ=255)

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

                     gaTCGACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAAc                                                                  <  1:7104823/93‑1 (MQ=255)

                         gACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCag                                                              >  1:501795/1‑93 (MQ=255)

                          aCCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCaga                                                             <  2:2426190/93‑1 (MQ=255)

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

                           cccGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCagaa                                                            >  2:7806435/1‑93 (MQ=255)

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

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

                                gCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAgaaagaaa                                                       <  2:6146702/93‑1 (MQ=255)

                                 ccTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGaaaa                                                      <  2:3490142/93‑1 (MQ=255)

                                 ccTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGaaaa                                                      >  1:7362701/1‑93 (MQ=255)

                                  ctgctgGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGaaaaa                                                     <  1:1593573/93‑1 (MQ=255)

                                  ctgctgGAAGACGACAAAGAGATGCTGGAAGACCTGGGGGATGCAGCAATCAACGATGCAGCAAGTCGTGTTGTAGAAACGCAGAAAGaaaaa                                                     >  2:847138/1‑93 (MQ=255)

                                   tgctgGAAGCCGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTAGTATTGAAGAAACGCAGAAAGaaaaaa                                                    >  2:6518476/1‑93 (MQ=255)

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

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

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

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

                                   tgctgGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGaaaaaa                                                    >  2:2983541/1‑93 (MQ=255)

                                    gctgGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCATCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAAt                                                   <  2:3706018/93‑1 (MQ=255)

                                     gtgGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATg                                                  <  1:2818745/92‑1 (MQ=255)

                                       tgAAGACGACAAAGATATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGc                                                <  1:3096494/92‑1 (MQ=255)

                                       ggAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGc                                                >  2:4715788/1‑93 (MQ=255)

                                         aaGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCct                                              <  2:2525710/93‑1 (MQ=255)

                                         aaGACGAAAGAGAGATGCTGGAGGACCTGGTAGCTGCAGGATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAAGGGCct                                              <  1:1068607/93‑1 (MQ=255)

                                          aGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCctc                                             <  1:1019480/93‑1 (MQ=255)

                                          aGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCctc                                             <  1:1474038/93‑1 (MQ=255)

                                           gacgacAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCctct                                            <  2:6235996/93‑1 (MQ=255)

                                           gacgacAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCctct                                            <  2:1810011/93‑1 (MQ=255)

                                             cgacAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGt                                          >  2:6648749/1‑93 (MQ=255)

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

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

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

                                               acAAAGAGATGCTGGAAGACCAGGTGGCTGCAGCCTTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTAt                                        >  2:2794602/1‑93 (MQ=255)

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

                                                 aaaGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTAtcc                                      <  1:3542027/93‑1 (MQ=255)

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

                                                  aaGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTAtcct                                     <  2:1266093/93‑1 (MQ=255)

                                                   agagATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTAtcctc                                    <  1:5550494/93‑1 (MQ=255)

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

                                                    gagaTGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTAtcctcc                                   <  1:2631615/93‑1 (MQ=255)

                                                    gagaTGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTAtcctcc                                   <  2:6075538/93‑1 (MQ=255)

                                                      gaTGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCgg                                 >  1:5230915/1‑93 (MQ=255)

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

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

                                                           tGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGc                            <  1:840743/93‑1 (MQ=255)

                                                           tGGAAGACCTGGTGGCTGCAGCAGTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGc                            <  2:3889157/93‑1 (MQ=255)

                                                              aaGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGc                         <  2:5933353/93‑1 (MQ=255)

                                                              aaGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGc                         <  2:5146782/93‑1 (MQ=255)

                                                              aaGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGc                         <  2:2678119/93‑1 (MQ=255)

                                                                 aCCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTgc                      >  2:5330324/1‑93 (MQ=255)

                                                                  ccTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAACAATGGCCTCTGTATCCTCCGGAATGCAGCTgcc                     >  1:7650834/1‑93 (MQ=255)

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

                                                                    tggtggCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTgccgc                   >  2:6908837/1‑93 (MQ=255)

                                                                    tggtggCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTgccgc                   >  1:59084/1‑93 (MQ=255)

                                                                    tggtggCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTgccgc                   >  2:7503955/1‑93 (MQ=255)

                                                                    tggtggCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTgccgc                   >  1:4730548/1‑93 (MQ=255)

                                                                    tggtggCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTgccgc                   >  1:5129191/1‑93 (MQ=255)

                                                                     ggtggCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTgccgcc                  >  1:2926420/1‑93 (MQ=255)

                                                                     ggtggCTGCAGCATCCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTgccgcc                  >  2:1644695/1‑93 (MQ=255)

                                                                      gtggCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCt                 >  2:2854205/1‑93 (MQ=255)

                                                                         gCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGc              >  2:4732157/1‑93 (MQ=255)

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

                                                                           tGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCtt            >  2:3808330/1‑93 (MQ=255)

                                                                           tGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATACTCCGGAATGCAGCTGCCGCCCGGCtt            >  2:4127456/1‑93 (MQ=255)

                                                                               gcaTTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAg        >  2:2091762/1‑93 (MQ=255)

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

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

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

                                                                                caTTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGa       >  1:6075818/1‑93 (MQ=255)

                                                                                caTTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGa       >  2:1985856/1‑93 (MQ=255)

                                                                                 aTTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGAt      >  2:575921/1‑93 (MQ=255)

                                                                                   tCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGc    >  1:4514849/1‑93 (MQ=255)

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

                                                                                     aaCGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCg  >  2:2891072/1‑93 (MQ=255)

ACACAACTGCCGTCGCGTAGAGATCGACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGACGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCG  >  NC_000913/493446‑493623

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: `‑`