BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,185,646	A→G	T194A (ACT→GCT)	pepT →	putative peptidase T

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,185,646	0	A	G	99.4%	508.5 / NA	158	T194A (ACT→GCT)	pepT	putative peptidase T
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (71/86); total (71/87)
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 165 total aligned reads displayed.

GTGATATTCGCGTCGCCTTTACCCCGGATGAAGAAGTGGGCAAAGGGGCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACACTGTTGATGGTGGTGGCGTAGGCGAACTGGAGTTTGAAAACTTCAACGCCGCGTCGGTCAATATCAAAATTGTCGGTAACAATGTTCAT  >  NC_000913/1185554‑1185735

  gATATTCGCGTCGCCTTTACCCCGGATGAAGAAGTGGGCAAAGGGGCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACGCt                                                                                         >  1:4598115/1‑93 (MQ=255)

     atTCGCGTCGCCTTTACCCCGGATGAAGAAGTGGGCAAAGGGGCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACGCTGtt                                                                                      >  1:3793012/1‑93 (MQ=255)

      ttCGCGTCGCCTTTACCCCGGATGAAGAAGTGGGCAAAGGGGCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACGCTGTTg                                                                                     >  1:1881111/1‑93 (MQ=255)

         gcgTCGCCTTTACCCCGGATGAAGAAGTGGGCAAAGGGGCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACGCTGTTGAtg                                                                                  >  2:1326519/1‑93 (MQ=255)

           gTCGCCTTTACCCCGGATGAAGAAGTGGGCAAAGGGGCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACGCTGTTGAtggt                                                                                >  1:1004919/1‑93 (MQ=255)

           gTCGCCTTTACCCCGGATGAAGAAGTGGGCAAAGGGGCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACGCTGTTGAtggt                                                                                <  2:1963861/93‑1 (MQ=255)

           gTCGCCTTTACCCCGGATGAAGAAGTGGGCAAAGGGGCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACGCTGTTGAtggt                                                                                <  1:4314187/93‑1 (MQ=255)

           gTCGCCTTTACCCCGGATGAAGAAGTGGGCAAAGGGGCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACGCTGTTGAtggt                                                                                <  1:4271270/93‑1 (MQ=255)

           gTCGCCTTTACCCCGGATGAAGAAGTGGGCAAAGGGGCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACGCTGTTGAtggt                                                                                <  1:1597168/93‑1 (MQ=255)

            tCGCCTTTACCCCGGATGAAGAAGTGGGCAAAGGGGCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACGCTGTTGAtggtg                                                                               <  1:1890804/93‑1 (MQ=255)

            tCGCCTTTACCCCGGATGAAGAAGTGGGCAAAGGGGCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACGCTGTTGAtggtg                                                                               >  1:4469194/1‑93 (MQ=255)

            gcGCCTTTCCCCCGGATGAAGAAGTGGGCAAAGGGGCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACGCTGTTGAtggtg                                                                               <  1:3435831/92‑1 (MQ=255)

             cGCCTTTACCCCGGATGAAGAAGTGGGCAAAGGGGCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACGCTGTTGAtggtgg                                                                              <  1:885262/93‑1 (MQ=255)

              gCCTTTACCCCGGATGAAGAAGTGGGCAAAGGGGCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACGCTGTTGAtggtggg                                                                             >  2:4053365/1‑92 (MQ=255)

                ctgtACCCCGGATGAAGAAGTGGGCAAAGGGGCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACGCTGTTGAtggtggtgg                                                                           <  2:3681659/90‑1 (MQ=255)

                cTTTACCCCGGATGAAGAAGTGGGCAAAGGGGCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACGCTGTTGAtggtggtgg                                                                           >  2:4050062/1‑93 (MQ=255)

                cTTTACCCCGGATGAAGAAGTGGGCAAAGGGGCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACGCTGTTGAtggtggtgg                                                                           >  1:661811/1‑93 (MQ=255)

                 tttACCCCGGATGAAGAAGTGGGCAAAGGGGCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACGCTGTTGAtggtggtggg                                                                          >  2:3951629/1‑92 (MQ=255)

                    aCCCCGGATGAAGAAGTGGGCAAAGGGGCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACGCTGTTGATGGTGGTGGCGTa                                                                       <  1:1743353/93‑1 (MQ=255)

                    aCCCCGGATGAAGAAGTGGGCAAAGGGGCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACGCTGTTGATGGTGGTGGCGTa                                                                       <  1:56161/93‑1 (MQ=255)

                     ccccGGATGAAGAAGTGGGCAAAGGGGCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACGCTGTTGATGGTGGTGGCGTAg                                                                      <  1:3473719/93‑1 (MQ=255)

                     ccccGGATGAAGAAGTGGGCAAAGGGGCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACGCTGTTGATGGTGGGGGCGTAg                                                                      >  2:1835827/1‑93 (MQ=255)

                      cccGGATGAAGAAGTGGGCAAAGGGGCGAAACAATTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACGCTGTTGATGGTGGTGGCGTAgg                                                                     >  2:3460190/1‑93 (MQ=255)

                       ccGGATGAAGAAGTGGGCAAAGGGGCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACGCTGTTGATGGTGGTGGCGTAggg                                                                    >  2:1097875/1‑92 (MQ=255)

                       ccGGATGAAGAAGTGGGCAAAGGGGCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACGCTGTTGATGGTGGTGGCGTAGGc                                                                    >  2:2665111/1‑93 (MQ=255)

                           aTGAAGAAGTGGGCAAAGGGGCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACGCTGTTGATGGTGGTGGCGTAGGCGAAc                                                                <  1:2555340/93‑1 (MQ=255)

                            tGAAGAAGTGGGCAAAGGGGCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACGCTGTTGATGGTGGTGGCGTAGGCGAACt                                                               <  1:1857479/93‑1 (MQ=255)

                            tGAAGAAGTGGGCAAAGGGGCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACGCTGTTGATGGTGGTGGCGTAGGCGAACt                                                               <  1:3626382/93‑1 (MQ=255)

                             gaagaaGTGGGCAAAGGGGCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACGCTGTTGATGGTGGTGGCGTAGGCGAACTg                                                              <  1:2285143/93‑1 (MQ=255)

                             gaagaaGTGGGCAAAGGGGCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACGCTGTTGATGGTGGTGGCGTAGGCGAACTg                                                              <  2:1451365/93‑1 (MQ=255)

                               agaagTGGGCAAAGGGGCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACGCTGTTGATGGTGGTGGCGTAGGCGAACTggg                                                            >  2:3194855/1‑92 (MQ=255)

                               agaagTGGGCAAAGGGGCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACGCTGTTGATGGTGGTGGCGTAGGCGAACTGGa                                                            <  1:281411/93‑1 (MQ=255)

                               agaagTGGGCAAAGGGGCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACGCTGTTGATGGTGGTGGCGTAGGCGAACTGGa                                                            <  1:1710734/93‑1 (MQ=255)

                                     gggCAAAGGGGCGAAACATTTTGATGTTGACTCCTTCGATGCCCGCTGGGCTTACGCTGTTGATGGTGGTGGCGTAGGCGAACTGGAGTTTGa                                                      <  1:972190/93‑1 (MQ=255)

                                     gggCAAAGGGGCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACGCTGTTGATGGTGGTGGCGTAGGCGAACTGGCGTTTGa                                                      >  1:3748785/1‑93 (MQ=255)

                                     gggCAAAGGGGCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACGCTGTTGATGGTGGTGGCGTAGGCGAACTGGAGTTTGa                                                      >  1:2485820/1‑93 (MQ=255)

                                      ggCAAAGGGGCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACGCTGTTGATGGTGGTGGCGTAGGCGAACTGGAGTTTGaa                                                     >  1:3393198/1‑93 (MQ=255)

                                      ggCAAAGGGGCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACGCTGTTGATGGTGGTGGCGGAGGCGAACTGGAGTTTGaa                                                     >  2:4720003/1‑93 (MQ=255)

                                       gCAAAGGGGCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACGCTGTTGATGGTGGTGGCGTAGGCGAACTGGAGTTTGaaa                                                    <  1:4055597/93‑1 (MQ=255)

                                        gaaaGGGGCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACGCTGTTGATGGTGGTGGCGTAGGCGAACTGGAGTTTGaaaa                                                   <  1:3521686/92‑1 (MQ=255)

                                          aaGGGGCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACGCTGTTGATGGTGGTGGCGTAGGCGAACTGGAGTTTGAAAACt                                                 >  1:3072381/1‑93 (MQ=255)

                                           aGGGGCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACGCTGTTGATGGTGGTGGCGTAGGCGAACTGGAGTTTGAAAACtt                                                >  2:1508/1‑93 (MQ=255)

                                           aGGGGCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACGCTGTTGATGGTGGTGGCGTAGGCGAACTGGAGTTTGAAAACtt                                                >  2:4263756/1‑93 (MQ=255)

                                            ggggCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACGCTGTTGATGGTGGTGGCGTAGGCGAACTGGAGTTTGAAAACTTc                                               <  2:3690480/93‑1 (MQ=255)

                                            ggggCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACGCTGTTGATGGTGGTGGCGTAGGCGAACTGGAGTTTGAAAACTTc                                               <  1:462740/93‑1 (MQ=255)

                                              ggCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACGCTGTTGATGGTGGTGGCGTAGGCGAACTGGAGTTTGAAAACTTCaa                                             >  2:2379539/1‑93 (MQ=255)

                                              ggCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACGCTGTTGATGGTGGTGGCGTAGGCGAACTGGAGTTTGAAAACTTCaa                                             <  2:3973358/93‑1 (MQ=255)

                                              ggCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACGCTGTTGATGGTGGTGGCGTAGGCGAACTGGAGTTTGAAAACTTCaa                                             <  2:304948/93‑1 (MQ=255)

                                               gCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACGCTGTTGATGGTGGTGGCGTAGGCGAACTGGAGTTTGAAAACTTCAAc                                            <  1:3041408/93‑1 (MQ=255)

                                               gCGAAACATTTTGATGTGGACGCCTTCGATGCCCGCTGGGCTTACGCTGTTGATGGTGGTGGCGTAGGCGAACTGGAGTTTGAAAACTTCAAc                                            <  1:182413/93‑1 (MQ=255)

                                                cGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACGCTGTTGATGGTGGTGGCGTAGGCGAACTGGAGTTTGAAAACTTCAAcg                                           <  1:3760027/93‑1 (MQ=255)

                                                cGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACGCTGTTGATGGTGGTGGCGTAGGCGAACTGGAGTTTGAAAACTTCAAcg                                           >  1:2153687/1‑93 (MQ=255)

                                                 gAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACGCTGTTGATGGTGGTGGCGTAGGCGAACTGGAGTTTGAAAACTTCAAcgc                                          >  1:2744651/1‑93 (MQ=255)

                                                 gAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACGCTGTTGATGGTGGTGGCGTAGGCGAACTGGAGTTTGAAAACTTCAAcgc                                          <  2:3233250/93‑1 (MQ=255)

                                                  aaaCATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACGCTGTTGATGGTGGTGGCGTAGGCGAACTGGAGTTTGAAAACTTCAAcgcc                                         >  1:979940/1‑93 (MQ=255)

                                                   aaCATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACGCTGTTGATGGTGGTGGCGTAGGCGAACTGGAGTTTGAAAACTTCAAcgccg                                        >  2:84110/1‑93 (MQ=255)

                                                    aCATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACGCTGTTGATGGTGGTGGCGTAGGCGAACTGGAGTTTGAAAACTTCAAcgccgc                                       <  1:3924420/93‑1 (MQ=255)

                                                     cATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACGCTGTTGATGGTGGTGGCGTAGGCGAACTGGAGTTTGAAAACTTCAACGCcgcg                                      <  1:77320/93‑1 (MQ=255)

                                                       ttttGATGTTGACGCCTTCGATGCCCGCTGGGCTTACGCTGTTGATGGTGGTGGCGTAGGCGAACTGGAGTTTGAAAACTTCAACGCCGCGTc                                    <  2:269903/93‑1 (MQ=255)

                                                       ttttGATGTTGACGCCTTCGATGCCCGCTGGGCTTACGCTGTTGATGGTGGTGGCGTAGGCGAACTGGAGTTTGAAAACTTCAACGCCGCGTc                                    >  1:1820118/1‑93 (MQ=255)

                                                       ttttGATGTTGACGCCTTCGATGCCCGCTGGGCTTACGCTGTTGATGGTGGTGGCGTAGGCGAACTGGAGTTTGAAAACTTCAACGCCGCGTc                                    >  2:597403/1‑93 (MQ=255)

                                                         ttGATGTTGACGCCTTCGATGCCCGCTGGGCTTACGCTGTTGATGGTGGTGGCGTAGGCGAACTGGAGTTTGAAAACTTCAACGCCGCGTCgg                                  <  2:3715950/93‑1 (MQ=255)

                                                         ttGATGTTGACGCCTTCGATGCCCGCTGGGCTTACGCTGTTGATGGTGGTGGCGTAGGCGAACTGGAGTTTGAAAACTTCAACGCCGCGTCgg                                  >  1:3012047/1‑93 (MQ=255)

                                                         ttGATGTTGACGCCTTCGATGCCCGCTGGGCTTACGCTGTTGATGGTGGTGGCGTAGGCGAACTGGAGTTTGAAAACTTCAACGCCGCGTCgg                                  <  2:4697770/93‑1 (MQ=255)

                                                           gATGTTGACGCCTTCGATGCCCGCTGGGCTTACGCTGTTGATGGTGGTGGCGTAGGCGAACTGGAGTTTGAAAACTTCAACGCCGCGTCGGTc                                <  1:336401/93‑1 (MQ=255)

                                                           gATGTTGACGCCTTCGATGCCCGCTGGGCTAACGCTGTTAATGGTGGTGACGTAGACGAACTAGGGTTTGGGTGCTTCGACGCCGAGTCGGTc                                <  2:4363826/93‑1 (MQ=255)

                                                            aTGTTGACGCCTTCGATGCCCGCTGGGCTTACGCTGTTGATGGTGGTGGCGTAGGCGAACTGGAGTTTGAAAACTTCAACGCCGCGTCGGTCa                               <  2:176284/93‑1 (MQ=255)

                                                             tGTTGACGCCTTTGATGCCCGCTGGGCTTACGCTGTTGATGGTGGTGGCGTAGGCGAACTGGAGTTTGAAAACTTCAACGCCGCGTCGGTCaa                              <  2:1462723/93‑1 (MQ=255)

                                                             tGTTGACGCCTTCGATGCCCGCTGGGCTTACGCTGTTGATGGTGGTGGCGTAGGCGAACTGGAGTTTGAAAACTTCAACGCCGCGTCGGTCaa                              <  2:3549154/93‑1 (MQ=255)

                                                             tGTTGACGCCTTCGATGCCCGCTGGGCTTACGCTGTTGATGGTGGTGGCGTAGGCGAACTGGAGTTTGAAAACTTCAACGCCGCGTCGGTCaa                              <  2:3549152/93‑1 (MQ=255)

                                                             tGTTGACGCCTTCGATGCCCGCTGGGCTTACGCTGTTGATGGTGGTGGCGTAGGCGAACTGGAGTTTGAAAACTTCAACGCCGCGTCGGTCaa                              >  1:2103766/1‑93 (MQ=255)

                                                                 gACGCCTTCGATGCCCGCTGGGCTTACGCTGTTGATGGTGGTGGCGTAGGCGAACTGGAGTTTGAAAACTTCAACGCCGCGTCGGTCAATATc                          <  2:2055853/93‑1 (MQ=255)

                                                                 gACGCCTTCGATGCCCGCTGGGCTTACGCTGTTGATGGTGGTGGCGTAGGCGAACTGGAGTTTGAAAACTTCAACGCCGCGTCGGTCAATATc                          >  2:663508/1‑93 (MQ=255)

                                                                   cGCCTTCGATGCCCGCTGGGCTTACGCTGTTGATGGTGGTGGCGTAGGCGAACTGGAGTTTGAAAACTTCAACGCCGCGTCGGTCAATATCaa                        <  2:4547548/93‑1 (MQ=255)

                                                                    gCCTTCGATGCCCGCTGGGCTTACGCTGTTGATGGTGGTGGCGTAGGCGAACTGGAGTTTGAAAACTTCAACGCCGCGTCGGTCAATATCaaa                       >  2:691661/1‑93 (MQ=255)

                                                                    gCCTTCGATGCCCGCTGGGCTTACGCTGTTGATGGTGGTGGCGTAGGCGAACTGGAGTTTGAAAACTTCAACGCCGCGTCGGTCAATATCaaa                       >  1:4027294/1‑93 (MQ=255)

                                                                      cTTCGATGCCCGCTGGGCTTACGCTGTTGATGGTGGTGGCGTAGGCGAACTGGAGTTTGAAAACTTCAACGCCGCGTCGGTCAATATCAAAAt                     <  2:956856/93‑1 (MQ=255)

                                                                      cTTCGATGCCCGCTGGGCTTACGCTGTTGATGGTGGTGGCGTAGGCGAACTGGAGTTTGAAAACTTCAACGCCGCGTCGGTCAATATCAAAAt                     >  1:1314669/1‑93 (MQ=255)

                                                                       ttCGATGCCCGCTGGGCTTACGCTGTTGATGGTGGTGGCGTAGGCGAACTGGAGTTTGAAAACTTCAACGCCGCGTCGGTCAATATCAAAAtt                    <  2:3302026/93‑1 (MQ=255)

                                                                       ttCGATGCCCGCTGGGCTTACGCTGTTGATGGTGGTGGCGTAGGCGAACTGGAGTTTGAAAACTTCAACGCCGCGTCGGTCAATATCAAAAtt                    >  1:2219523/1‑93 (MQ=255)

                                                                       ttCGATGCCCGCTGGGCTTACGCTGTTGATGGTGGTGGCGTAGGCGAACTGGAGTTTGAAAACTTCAACGCCGCGTCGGTCAATATCAAAAtt                    <  1:1295545/93‑1 (MQ=255)

                                                                       ttCGATGCCCGCTGGGCTTACGCTGTTGATGGTGGTGGCGTAGGCGAACTGGAGTTTGAAAACTTCAACGCCGCGTCGGTCAATATCAAAAtt                    <  2:3149597/93‑1 (MQ=255)

                                                                        tCGATGCCCGCTGGGCTTACGCTGTTGATGGTGGTGGCGTAGGCGAACTGGAGTTTGAAAACTTCAACGCCGCGTCGGTCACTATCAAACTTg                   >  2:854353/1‑93 (MQ=255)

                                                                          gATGCCCGCTGGGCTTACGCTGTTGATGGTGGTGGCGTAGGCGAACTGGAGTTTGAAAACTTCAACGCCGCGTCGGTCAATATCAAAATTGTc                 >  2:768751/1‑93 (MQ=255)

                                                                           aTGCCCGCTGGGCTTACGCTGTTGATGGTGGTGGCGTAGGCGAACTGGAGTTTGAAAACTTCAACGCCGCGTCGGTCAATATCAAAATTGTCg                <  1:3607325/93‑1 (MQ=255)

                                                                           aTGCCCGCTGCGCTTACGCTGTTGATGGTGGTGGCGTAGGCGAACTGGAGTTTGAAAACTTCAACGCCGCGTCGGTCAATATCAAAATTGTCg                <  1:320718/93‑1 (MQ=255)

                                                                            tGCCCGCTGGGCTTACGCTGTTGATGGTGGTGGCGTAGGCGAACTGGAGTTTGAAAACTTCAACGCCGCGTCGGTCAATATCAAAATTGTCgg               >  1:2922215/1‑93 (MQ=255)

                                                                              cccGCTGGGCTTACGCTGTTGATGGTGGTGGCGTAGGCGAACTGGAGTTTGAAAACTTCAACGCCGCGTCGGTCAATATCAAAATTGTCGGTa             >  2:1382365/1‑93 (MQ=255)

                                                                              cccGCTGGGCTTACGCTGTTGATGGTGGTGGCGTAGGCGAACTGGAGTTTGAAAACTTCAACGCCGCGTCGGTCAATATCAAAATTGTCGGTa             <  1:2449133/93‑1 (MQ=255)

                                                                               ccGCTGGGCTTACGCTGTTGATGGTGGTGGCGTAGGCGAACTGGAGTTTGAAAACTTCAACGCCGCGTCGGTCAATATCAAAATTGTCGGTaa            >  2:1722248/1‑93 (MQ=255)

                                                                                cGCTGGGCTTACGCTGTTGATGGTGGTGGCGTAGGCGAACTGGAGTTTGAAAACTTCAACGCCGCGTCGGTCAATATCAAAATTGTCGGTAAc           <  2:3882032/93‑1 (MQ=255)

                                                                                  cttggCTTACGCTGGTGATGGTGGTGGCGTAGGCGAACTGGAGTTTGAAAACTTCAACGCCGCGTCGGTCAATATCAAAATTGTCGGTAACaa         <  1:1383715/90‑1 (MQ=255)

                                                                                  cTGGGCTTACGCTGTTGATGGTGGTGGCGTAGGCGAACTGGAGTTTGAAAACTTCAACGCCGCGTCGGTCAATATCAAAATTGTCGGTAACaa         >  1:3566863/1‑93 (MQ=255)

                                                                                   tGGGCTTACGCTGTTGATGGTGGTGGCGTAGGCGAACTGGAGTTTGAAAACTTCAACGCCGCGTCGGTCAATATCAAAATTGTCGGTAACAAt        <  1:1242011/93‑1 (MQ=255)

                                                                                     ggCTTACGCTGTTGATGGTGGTGGCGTAGGCGAACTGGAGTTTGAAAACTTCAACGCCGCGTCGGTCAATATCAAAATTGTCGGTAACAATGt      >  2:2311671/1‑93 (MQ=255)

                                                                                     ggCTTACGCTGTTGATGGTGGTGGCGTAGGCGAACTGGAGTTTGAAAACTTCAACGCCGCGTCGGTCAATATCAAAATTGTCGGTAACAATGt      <  2:302032/93‑1 (MQ=255)

                                                                                     ggCTTACGCTGTTGATGGTGGTGGCGTAGGCGAACTGGAGTTTGAAAACTTCAACGCCGCGTCGGTCAATATCAAAATTGTCGGTAACAATGt      <  2:3282904/93‑1 (MQ=255)

                                                                                     ggCTTACGCTGTTGATGGTGGTGGCGTAGGCGAACTGGAGGTTGAAAACTTCAACGCCGCGTCGGTCAATATCAAAATTGTCGGTAACAATGt      >  2:4300693/1‑93 (MQ=255)

                                                                                      gCTTACGCTGTTGATGGTGGTGGCGTAGGCGAACTGGAGTTTGAAAACTTCAACGCCGCGTCGGTCAATATCAAAATTGTCGGTAACAATGtt     <  2:3388073/93‑1 (MQ=255)

                                                                                        ttACGCTGTTGATGGTGGTGGCGTAGGCGAACTGGAGTTTGAAAACTTCAACGCCGCGTCGGTCAATATCAAAATTGTCGGTAACAATGTTCa   >  1:309945/1‑93 (MQ=255)

GTGATATTCGCGTCGCCTTTACCCCGGATGAAGAAGTGGGCAAAGGGGCGAAACATTTTGATGTTGACGCCTTCGATGCCCGCTGGGCTTACACTGTTGATGGTGGTGGCGTAGGCGAACTGGAGTTTGAAAACTTCAACGCCGCGTCGGTCAATATCAAAATTGTCGGTAACAATGTTCAT  >  NC_000913/1185554‑1185735

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