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.4%	1148.1 / ‑4.7	313	D79D (GAC→GAT)	ybaB	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (0/2); new base T (144/167); total (144/169)
Fisher's exact test for biased strand distribution p-value = 5.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

Only 100 of 328 total aligned reads displayed.

ACACAACTGCCGTCGCGTAGAGATCGACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGACGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTT  >  NC_000913/493446‑493625

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

  acaACTGCCGTCGCGTAGAGATCGACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATgcag                                                                                       >  2:2433730/1‑93 (MQ=255)

  acaACTGCCGTCGCGTAGAGATCGACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATgcag                                                                                       <  1:4695943/93‑1 (MQ=255)

  acaACTGCCGTCGCGTAGAGATCGACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATgcag                                                                                       <  2:4383993/93‑1 (MQ=255)

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

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

       tGCCGTCGCGTAGAGATCGACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCAcgt                                                                                  <  2:3828709/93‑1 (MQ=255)

       tGCCGTCGCGTAGAGATCGACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCAcgt                                                                                  <  1:126197/93‑1 (MQ=255)

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

          cGTCGCGTAGAGATCGACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCAcgtcgt                                                                               <  1:358167/93‑1 (MQ=255)

           gtggcgTAGAGATCGACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTa                                                                              <  1:1451276/90‑1 (MQ=255)

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

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

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

               cgTAGAGATCGCCCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTga                                                                          <  1:9838907/93‑1 (MQ=255)

               cgTAGAGATCGACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTga                                                                          <  1:6995539/93‑1 (MQ=255)

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

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

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

                  agagATCGGCCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTgaaga                                                                       <  2:5550993/93‑1 (MQ=255)

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

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

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

                      aTCGACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACg                                                                   >  2:4549341/1‑93 (MQ=255)

                        cGACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCa                                                                 <  1:8353050/93‑1 (MQ=255)

                        cGACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCa                                                                 <  1:1859569/93‑1 (MQ=255)

                          aCCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAGCGATGCAGCACGTCGTATTGAAGAAACGCaga                                                               <  2:8363481/93‑1 (MQ=255)

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

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

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

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

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

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

                                     ctgGAAGACGACAAAGAGATGCTGGAAGCCCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATg                                                    <  1:5186145/93‑1 (MQ=255)

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

                                      tgGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATgg                                                   >  2:3957998/1‑93 (MQ=255)

                                       ggAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGc                                                  <  2:1482881/93‑1 (MQ=255)

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

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

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

                                         aaGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCct                                                >  1:242287/1‑93 (MQ=255)

                                         aaGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCct                                                >  2:1144543/1‑93 (MQ=255)

                                         aaGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCct                                                >  1:950615/1‑93 (MQ=255)

                                         aaGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGACGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCct                                                <  1:2623178/93‑1 (MQ=255)

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

                                           gacgacAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCctct                                              >  2:7660537/1‑93 (MQ=255)

                                            acgacAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTg                                             <  2:360709/93‑1 (MQ=255)

                                            acgacAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTg                                             <  1:3812266/93‑1 (MQ=255)

                                             cgacAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGt                                            <  1:6333135/93‑1 (MQ=255)

                                              gacAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCGGAAAGAAAAAATGGCCTCTGTa                                           >  2:3062208/1‑93 (MQ=255)

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

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

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

                                                aaaaGAGAGGCTAGAAGGGCTGGTGGATGCAGCATTCAACCATGCAGCACGTCGTATTGAAGAAACGGAGAAAGAAAGAATGGCCTCTGTagc                                         <  1:1479127/92‑3 (MQ=255)

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

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

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

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

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

                                                        tGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGaa                                 <  1:7040669/93‑1 (MQ=255)

                                                        tGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGaa                                 <  1:48700/93‑1 (MQ=255)

                                                         gCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTCTTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAAt                                >  1:8186312/1‑93 (MQ=255)

                                                         gCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAAt                                >  2:2089607/1‑93 (MQ=255)

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

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

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

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

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

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

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

                                                               aGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCt                          >  1:7567951/1‑93 (MQ=255)

                                                                 aCCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTgc                        <  2:3481076/93‑1 (MQ=255)

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

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

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

                                                                    tggtggCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTggcgc                     >  2:255823/1‑93 (MQ=255)

                                                                    tggtggCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTgccgc                     <  2:9052950/93‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

                                                                                caTTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGACCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGa         <  2:8927761/93‑1 (MQ=255)

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

                                                                                  ttCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATg       >  1:5126510/1‑93 (MQ=255)

                                                                                  ttCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATg       <  2:8463725/93‑1 (MQ=255)

                                                                                  ttCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATg       >  2:537819/1‑93 (MQ=255)

                                                                                  ttCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATg       >  2:7403000/1‑93 (MQ=255)

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

                                                                                      acgatgCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGt   >  1:7363424/1‑93 (MQ=255)

                                                                                      acgatgCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGt   >  1:10021096/1‑93 (MQ=255)

                                                                                      acgatgCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCGGGCTTTAAGATGCCGt   <  2:8007124/93‑1 (MQ=255)

                                                                                       cgatgCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGtt  >  1:1328651/1‑93 (MQ=255)

                                                                                       cgatgCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGtt  <  1:6777511/93‑1 (MQ=255)

                                                                                       cgatgCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGtt  >  2:5991031/1‑93 (MQ=255)

ACACAACTGCCGTCGCGTAGAGATCGACCCGAGCCTGCTGGAAGACGACAAAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGACGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTT  >  NC_000913/493446‑493625

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