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.6%	1024.7 / ‑5.0	273	D79D (GAC→GAT)	ybaB	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (0/0); major base T (133/139); minor base G (0/1); total (133/140)
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.97e-01

Only 100 of 283 total aligned reads displayed.

CACAACTGCCGTCGCGTAGAGATCGACCCGAGCCTGCTGGAAGACGACA‑AAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGACGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTT  >  NC_000913/493447‑493625

cacaACTGCCGTCGCGTAGAGATCGACCCGAGCCTGCTGGAAGACGACA‑AAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATgca                                                                                        >  2:4908423/1‑93 (MQ=255)

 acaACTGCCGTCGCGTAGAGATCGACCCGAGCCTGCTGGAAGACGACA‑AAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATgcag                                                                                       >  1:678728/1‑93 (MQ=255)

 acaACTGCCGTCGCGTAGAGATCGACCCGAGCCTGCTGGAAGACGACA‑AAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATgcag                                                                                       >  2:3856133/1‑93 (MQ=255)

  caACTGCCGTCGCGTAGAGATCGACCCGAGCCTGCTGGAAGACGACA‑AAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATgcagc                                                                                      >  2:4065875/1‑93 (MQ=255)

  caACTGCCGTCGCGTAGAGATCGACCCGAGCCTGCTGGAAGACGACA‑AAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATgcagc                                                                                      >  1:870010/1‑93 (MQ=255)

   agcTGCCGTCGCGTAGAGATCGACCCGAGCCTGCTGGAAGACGACA‑AAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATgcagca                                                                                     <  2:2018016/91‑1 (MQ=255)

   aaCTGCCGTCGCGTAGAGATCGACCCGAGCCTGCTGGAAGACGACA‑AAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATgcagct                                                                                     >  2:4464076/1‑92 (MQ=255)

   aaCTGCCGTCGCGTAGAGATCGACCCGAGCCTGCTGGAAGACGACA‑AAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATgcagca                                                                                     <  1:1006806/93‑1 (MQ=255)

   aaCTGCCGTCGCGTAGAGATCGACCCGAGCCTGCTGGAAGACGACA‑AAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATgcagca                                                                                     >  2:2720618/1‑93 (MQ=255)

    agtcccGTGGAGTGGAGATCGACCCGAGCCTGCTGGAAGACGACA‑AGGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCAc                                                                                    <  2:5091508/89‑1 (MQ=255)

     cTGCCGTCGCGTAGAGATCGACCCGAGCCTGCTGGAAGACGACA‑AAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCAcg                                                                                   >  2:5968347/1‑93 (MQ=255)

      tGCCGTCGCGTAGAGATCGACCCGAGCCTGCTGGAAGACGACA‑AAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCAcgt                                                                                  <  2:2273861/93‑1 (MQ=255)

       gCCGTCGCGTAGAGATCGACCCGAGCCTGCTGGAAGACGACA‑AAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCAcgtc                                                                                 <  1:1812249/93‑1 (MQ=255)

       gCCGTCGCGTAGAGATCGACCCGAGCCTGCTGGAAGACGACA‑AAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCAcgtc                                                                                 <  1:352882/93‑1 (MQ=255)

        ccGTCGCGTAGAGATCGACCCGAGCCTGCTGGAAGACGACA‑AAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCAcgtcg                                                                                <  1:5628889/93‑1 (MQ=255)

          gTCGCGTAGAGATCGACCCGAGCCTGCTGGAAGGCGACA‑AAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGCCGTa                                                                              >  2:1278078/1‑93 (MQ=255)

            cgcgTAGAGATCGACCCGAGCCTGCTGGAAGACGACA‑AAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTAtt                                                                            <  2:325890/93‑1 (MQ=255)

            cgcgTAGAGATCGACCCGAGCCTGCTGGAAGACGACA‑AAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTAtt                                                                            <  1:5116700/93‑1 (MQ=255)

            cgcgTAGAGATCGACCCGAGCCTGCTGGAAGACGACA‑AAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTAtt                                                                            <  1:3692476/93‑1 (MQ=255)

                tAGAGATCGACCCGAGCCTGCTGGAAGACGACA‑AAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTgaag                                                                        >  1:256035/1‑93 (MQ=255)

                tAGAGATCGACCCGAGCCTGCTGGAAGACGACA‑AAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTgaag                                                                        >  1:824640/1‑93 (MQ=255)

                  gagaTCGACCCGAGCCTGCTGGAAGACGACA‑AAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTgaagaa                                                                      >  1:1072626/1‑93 (MQ=255)

                   agaTCGACCCGAGCCTGCTGGAAGACGACA‑AAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGaaa                                                                     <  1:1460452/93‑1 (MQ=255)

                   agaTCGACCCGAGCCTGCTGGAAGACGACA‑AAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGaaa                                                                     <  2:4386900/93‑1 (MQ=255)

                         aCCCGAGCCTGCTGGAAGACTACA‑AAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCaga                                                               >  1:5410310/1‑93 (MQ=255)

                          cccGAGCCTGCTGGAAGACGACA‑AAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTTAAGAAACGCagaa                                                              >  2:4296771/1‑93 (MQ=255)

                           ccGAGCCTGCTGGAAGACGACA‑AAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCagaaa                                                             <  1:1783514/93‑1 (MQ=255)

                            cGAGCCTGCTGGAAGACGACA‑AAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCagaaag                                                            <  2:2200453/93‑1 (MQ=255)

                             gAGCCTGCTGGAAGACGACA‑AAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCagaaaga                                                           >  2:626397/1‑93 (MQ=255)

                              aGCCTGCTGGAAGACGACA‑AAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCagaaagaa                                                          >  1:1224017/1‑93 (MQ=255)

                              aGCCTGCTGGAAGACGACA‑AAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCagaaagaa                                                          <  1:4836158/93‑1 (MQ=255)

                               gCCTGCTGGAAGACGACA‑AAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAgaaagaaa                                                         <  1:2059306/93‑1 (MQ=255)

                                ccTGCTGGAAGACGACA‑AAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGaaaa                                                        <  1:6194515/93‑1 (MQ=255)

                                ccTGCTGGAAGACGACA‑AAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGaaaa                                                        <  2:4971040/93‑1 (MQ=255)

                                  tgctgGAAGACGACA‑AAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGaaaaaa                                                      <  1:4982387/93‑1 (MQ=255)

                                   gctgGAAGACGACA‑AAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAAt                                                     <  2:4562999/93‑1 (MQ=255)

                                   gctgGAAGACGACA‑AAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAAt                                                     <  2:4777116/93‑1 (MQ=255)

                                    gcagAAGACGGCATACGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAAt                                                     <  1:2921534/90‑1 (MQ=255)

                                    ctgGAAGACGACA‑AAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATg                                                    <  1:36750/93‑1 (MQ=255)

                                     tgGAAGACGACA‑AAGAGATGCTGGAAGACCTGGTGGCTTCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATgg                                                   >  2:978341/1‑93 (MQ=255)

                                     tgGAAGACGACA‑AAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAAttg                                                   >  2:5423369/1‑91 (MQ=255)

                                     tgGAAGACGACA‑AAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATgg                                                   >  1:4724618/1‑93 (MQ=255)

                                     tgGAAGACGACA‑AAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATgg                                                   >  2:5242173/1‑93 (MQ=255)

                                     tgGAAGACGACA‑AAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATgg                                                   >  2:679419/1‑93 (MQ=255)

                                      ggAAGACGACA‑AAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGc                                                  <  1:1885525/93‑1 (MQ=255)

                                      ggAAGACGACA‑AAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGc                                                  <  1:2161153/93‑1 (MQ=255)

                                      ggAAGACGACA‑AAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGc                                                  <  1:3419503/93‑1 (MQ=255)

                                       gAAGACGACA‑AAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGcc                                                 >  2:3034253/1‑93 (MQ=255)

                                       gAAGACGACA‑AAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGcc                                                 <  1:5540397/93‑1 (MQ=255)

                                       gAAGACGACA‑AAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGcc                                                 <  1:2877434/93‑1 (MQ=255)

                                       gAAGACGACA‑AAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGcc                                                 <  1:986585/93‑1 (MQ=255)

                                         aGACGACA‑AAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCctc                                               >  2:5868436/1‑93 (MQ=255)

                                          gacgacA‑AAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCctct                                              <  2:239265/93‑1 (MQ=255)

                                           acgacA‑AAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTg                                             <  2:3448898/93‑1 (MQ=255)

                                           acgacA‑AAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTg                                             <  2:3777097/93‑1 (MQ=255)

                                             gacA‑AAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTa                                           <  1:4897696/93‑1 (MQ=255)

                                             gacA‑AAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTa                                           >  1:3056079/1‑93 (MQ=255)

                                              acA‑AAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTAt                                          <  1:4783489/93‑1 (MQ=255)

                                               cA‑AAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTAtc                                         >  1:4106839/1‑93 (MQ=255)

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

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

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

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

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

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

                                                         gCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAAt                                <  1:6120475/93‑1 (MQ=255)

                                                         gCTGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAAt                                <  2:4167393/93‑1 (MQ=255)

                                                           tGGAAGACCTGGTGGCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGTAATGc                              >  1:4609002/1‑93 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

                                                                     ggtggCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAAGGGCCTCTGTATCCTCCGGAATGCAGCTgccgcc                    <  2:830478/93‑1 (MQ=255)

                                                                       tggCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTg                  <  1:5868922/93‑1 (MQ=255)

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

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

                                                                        ggCTGCAGCATTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTgg                 >  2:5808759/1‑93 (MQ=255)

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

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

                                                                          cTGCAGCATTCACCGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTGTCCTCCGGAATGCAGCTGCCGCCTGGCt               <  1:3363511/93‑1 (MQ=255)

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

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

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

                                                                             cagcaTTCAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTa            >  2:2887560/1‑93 (MQ=255)

                                                                              agcaTTCAACGATGCAGCAAGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTaa           >  2:6148254/1‑93 (MQ=255)

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

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

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

                                                                                    cAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGcc     <  2:2572931/93‑1 (MQ=255)

                                                                                    cAACGATGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGcc     <  2:3919750/93‑1 (MQ=255)

                                                                                      acgatgCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGt   <  1:1543899/93‑1 (MQ=255)

                                                                                       cgatgCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCTTGGCTTTAAGATGCCGtt  >  1:4257320/1‑93 (MQ=255)

CACAACTGCCGTCGCGTAGAGATCGACCCGAGCCTGCTGGAAGACGACA‑AAGAGATGCTGGAAGACCTGGTGGCTGCAGCATTCAACGACGCAGCACGTCGTATTGAAGAAACGCAGAAAGAAAAAATGGCCTCTGTATCCTCCGGAATGCAGCTGCCGCCTGGCTTTAAGATGCCGTT  >  NC_000913/493447‑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: `‑`