BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	198,418	A→G	D164G (GAC→GGC)	yaeT →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	198,418	0	A	G	98.4%	578.1 / ‑3.2	184	D164G (GAC→GGC)	yaeT	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (1/2); new base G (90/91); total (91/93)
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.98e-01

Only 100 of 186 total aligned reads displayed.

AGAAAGGTCTGGAAGACTTCTACTACAGCGTCGGTAAATATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGACCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACATTGTTGGTAACCATGCTTTCACCACCGACGAACTGA  >  NC_000913/198328‑198507

agaAAGGTCTGGAAGACTTCTACTACAGCGTCGGTAAATATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGcc                                                                                         >  1:102596/1‑93 (MQ=255)

agaAAGGTCTGGAAGACTTCTACTACAGCGTCGGTAAATATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGcc                                                                                         >  2:4554257/1‑93 (MQ=255)

 gaAAGGTCTGGAAGACTTCTACTACAGCGTCGGTAAATATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCt                                                                                        >  1:4432334/1‑93 (MQ=255)

 gaAAGGTCTGGAAGACTTCTACTACAGCGTCGGTAAATATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCt                                                                                        >  1:345812/1‑93 (MQ=255)

 gaAAGGTCTGGAAGACTTCTACTACAGCGTCGGTAAATATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCt                                                                                        >  2:502478/1‑93 (MQ=255)

  aaaGGTCTGGAAGACTTCTACTACAGCGTCGGTAAATATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTa                                                                                       >  1:195490/1‑93 (MQ=255)

    aGGTCTGGAAGACTTCTACTACAGCGTCGGTAAATATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTaaa                                                                                     >  1:902726/1‑93 (MQ=255)

     ggTCTGGAAGACTTCTACTACAGCGTCGGTAAATATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTaaaa                                                                                    >  1:2938130/1‑93 (MQ=255)

     ggTCTGGAAGACTTCTACTACAGCGTCGGTAAATATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTaaaa                                                                                    >  2:4379207/1‑93 (MQ=255)

      gTCTGGAAGACTTCTACTACAGCGTCGGTAAATATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAAc                                                                                   >  2:2069199/1‑93 (MQ=255)

      gTCTGGAAGACTTCTACTACAGCGTCGGTAAATATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAAc                                                                                   <  2:652861/93‑1 (MQ=255)

       tCTGGAAGACTTCTACTACAGCGTCGGTAAATATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACt                                                                                  <  1:3503669/93‑1 (MQ=255)

       tCTGGAAGACTTCTACTACAGCGTCGGTAAATATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACt                                                                                  >  2:3498714/1‑93 (MQ=255)

        cTGGAAGACTTCTACTACAGCGTCGGTAAATATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTg                                                                                 <  2:213345/93‑1 (MQ=255)

         tGGAAGACTTCTACTACAGCGTCGGTAAATATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTgg                                                                                >  2:1853957/1‑93 (MQ=255)

          ggAAGACTTCTACTACAGCGTCGGTAAATATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGgt                                                                               >  2:3532941/1‑93 (MQ=255)

          ggAAGACTTCTACTACAGCGTCGGTAAATATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGgt                                                                               >  2:2673417/1‑93 (MQ=255)

           gAAGACTTCTACTACAGCGTCGGTAAATATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGgtg                                                                              >  2:220347/1‑93 (MQ=255)

           gAAGACTTCTACTACAGCGTCGGTAAATATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGgtg                                                                              <  1:706733/93‑1 (MQ=255)

            gaGACGTCTACTACAGCGTCGGTAAATATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGgtgt                                                                             <  1:123610/92‑1 (MQ=255)

            aaGACTTCTACTACAGCGTCGGTAAATATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGgtgt                                                                             <  1:2695296/93‑1 (MQ=255)

            aaGACTTCTACTACAGCGTCGGTAAATATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGgtgt                                                                             <  2:3157032/93‑1 (MQ=255)

             aGACTTCTACTACAGCGTCGGTAAATATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGtt                                                                            <  2:4392058/93‑1 (MQ=255)

             aGACTTCTACTACAGCGTCGGTAAATATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGtt                                                                            >  1:312404/1‑93 (MQ=255)

             aGACTTCTACTACAGCGTCGGTAAATATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGAGtt                                                                            <  1:4377546/93‑1 (MQ=255)

                cTTCTACTACAGCGTCGGTAAATATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCa                                                                         <  1:2017951/93‑1 (MQ=255)

                 ttCTACTACAGCGTCGGTAAATATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAg                                                                        <  1:4146435/93‑1 (MQ=255)

                 ttCTACTACAGCGTCGGTAAATATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAg                                                                        <  1:3701461/93‑1 (MQ=255)

                    tactacAGCGTCGGTAAATATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGaa                                                                     >  1:1101664/1‑93 (MQ=255)

                      ctacAGCGTCGGTAAATATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAgg                                                                   <  2:3780902/93‑1 (MQ=255)

                      ctacAGCGTCGGTAAATATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGGGTTCCAGGAAgg                                                                   >  2:277406/1‑93 (MQ=255)

                       tacAGCGTCGGTAAATATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAggg                                                                  >  2:3510024/1‑92 (MQ=255)

                       tacAGCGTCGGTAAATATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGgt                                                                  >  1:2635694/1‑93 (MQ=255)

                        acAGCGTCGGTAAATATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGgtg                                                                 >  1:4699152/1‑93 (MQ=255)

                          aGCGTCGGTAAATATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGgtgtg                                                               >  1:505011/1‑93 (MQ=255)

                          aGCGTCGGTAAATATAGCGCCAGCGTAAAAGCTGTCGTGACACCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGgtgtg                                                               >  1:930592/1‑93 (MQ=255)

                           gCGTCGGTAAATATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGgtgtgt                                                              >  1:2754067/1‑93 (MQ=255)

                             gTCGGTAAATATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCa                                                            >  1:2473559/1‑93 (MQ=255)

                              tCGGTAAATATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAg                                                           >  1:4447711/1‑93 (MQ=255)

                               cGGTAAATATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGc                                                          <  2:2506660/93‑1 (MQ=255)

                               cGGTAAATATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGc                                                          >  1:1911294/1‑93 (MQ=255)

                               cGGTAAATATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGc                                                          >  1:4253173/1‑93 (MQ=255)

                               cGGTAAATATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGc                                                          <  2:1276072/93‑1 (MQ=255)

                                ggTAAATATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCt                                                         <  1:4485190/93‑1 (MQ=255)

                                  tAAATATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGa                                                       >  2:3330334/1‑93 (MQ=255)

                                    aaTATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGaaa                                                     >  2:3012927/1‑93 (MQ=255)

                                    aaTATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGaaa                                                     >  1:3667593/1‑93 (MQ=255)

                                     atatAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGaaac                                                    >  1:3952896/1‑92 (MQ=255)

                                     atatAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAAt                                                    <  2:1213916/93‑1 (MQ=255)

                                      tataGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATc                                                   >  1:4067348/1‑93 (MQ=255)

                                      tataGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATc                                                   >  1:1138998/1‑93 (MQ=255)

                                      aataGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATc                                                   <  1:1145386/92‑1 (MQ=255)

                                       ataGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATcc                                                  >  2:145505/1‑93 (MQ=255)

                                          gcgcCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGc                                               >  2:809902/1‑93 (MQ=255)

                                           cgcCAGCCTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCa                                              <  1:1105983/93‑1 (MQ=255)

                                             ccAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCaaa                                            <  1:2722460/93‑1 (MQ=255)

                                             ccAGCATAAAAGCTGTCGTGCCCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTATCAGCTGAAATCCAGCaaa                                            <  2:1731903/93‑1 (MQ=255)

                                              cAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAAt                                           <  2:2662296/93‑1 (MQ=255)

                                               aGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAAtt                                          <  1:1386381/93‑1 (MQ=255)

                                                 cGTAAAAGCTGTCGTGACGCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTaa                                        <  1:2076912/93‑1 (MQ=255)

                                                   tataaGCTGTCGTGCCCCCGCAGCCGCGCAGCCGTATTGGCCTAAAACTGGTGTTACAGGAAGGTATATCAGCTGAAATCCAGCAAATTAACa                                      <  1:4209755/90‑1 (MQ=255)

                                                   tAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACa                                      >  1:1402572/1‑93 (MQ=255)

                                                   tAAAAGCTGTCGTGACCCCGCTGACGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAACAAATTAACa                                      >  2:3840287/1‑93 (MQ=255)

                                                     aaaGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACAtt                                    >  2:1532435/1‑93 (MQ=255)

                                                     aaaGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACAtt                                    >  1:2826148/1‑93 (MQ=255)

                                                     aaaGCTGGCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTTTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACAtt                                    <  2:1876483/93‑1 (MQ=255)

                                                      aaGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACAttt                                   >  2:5698/1‑92 (MQ=255)

                                                      aaGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACAttg                                   <  2:2287532/93‑1 (MQ=255)

                                                       aGCTGTCGTGGCCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACAttgt                                  <  1:4251855/93‑1 (MQ=255)

                                                        gCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACAttgtt                                 <  1:1352362/93‑1 (MQ=255)

                                                        gCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACAttgtt                                 <  2:847681/93‑1 (MQ=255)

                                                        gCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACAttgtt                                 <  1:1352361/93‑1 (MQ=255)

                                                        gCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACAttgtt                                 <  1:3212604/93‑1 (MQ=255)

                                                        gCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACAttgtt                                 <  2:3058489/93‑1 (MQ=255)

                                                        gCTGTCGTGAACCCGCTGCCGCGCAACCGTGTTGACCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACAttgtt                                 <  2:3583755/93‑1 (MQ=255)

                                                         cTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACAttgttg                                <  2:1322460/93‑1 (MQ=255)

                                                          tGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACATTGTTgg                               <  2:2759426/93‑1 (MQ=255)

                                                           gTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACATTGTTGGt                              <  2:3406154/93‑1 (MQ=255)

                                                            tCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACATTGTTGGTa                             >  1:1052597/1‑93 (MQ=255)

                                                               tGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACATTGTTGGTAAcc                          <  2:2930059/93‑1 (MQ=255)

                                                                 aCCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACATTGTTGGTAACCAt                        <  1:1915544/93‑1 (MQ=255)

                                                                 aCCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACATTGTTGGTAACCAt                        <  1:3812701/93‑1 (MQ=255)

                                                                    ccGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACATTGTTGGTAACCATGCt                     <  2:2294838/93‑1 (MQ=255)

                                                                     cGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACATTGTTGGTAACCATGCtt                    <  1:1477330/93‑1 (MQ=255)

                                                                     cGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACATTGTTGGTAACCATGCtt                    >  1:1754338/1‑93 (MQ=255)

                                                                      gCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACATTGTTGGTAACCATGCttt                   >  1:590183/1‑93 (MQ=255)

                                                                      gCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACATTGTTGGTAACCATGCttt                   >  2:1589937/1‑93 (MQ=255)

                                                                      gCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACATTGTTGCGAACAATGCttt                   >  1:590182/1‑93 (MQ=255)

                                                                        tGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACATTGTTGGTAACCATGCTTTca                 >  1:2680149/1‑93 (MQ=255)

                                                                           cgcgCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACATTGTTGGTAACCATGCTTTcacca              <  1:2546280/93‑1 (MQ=255)

                                                                            gcgcAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACATTGTTGGTAACCATGCTTTcaccac             >  2:1350507/1‑93 (MQ=255)

                                                                            gcgcAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACATTGTTGGTAACCATGCTTTcaccac             >  2:3030785/1‑93 (MQ=255)

                                                                             cgcAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACATTGTTGGTAACCATGCTTTCaccacc            >  2:2229672/1‑93 (MQ=255)

                                                                              gcAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACATTGTTGGTAACCATGCTTTCACCACcg           <  2:4066008/93‑1 (MQ=255)

                                                                                 aCCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACATTGTTGGTAACCATGCTTTCACCACcgacg        >  1:748804/1‑93 (MQ=255)

                                                                                 aCCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACATTGTTGGTAACCATGCTTTCACCACcgacg        <  1:3416257/93‑1 (MQ=255)

                                                                                    gtgtTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACATTGTTGGTAACCATGCTTTCACCACCGACGAAc     >  1:2453646/1‑93 (MQ=255)

                                                                                    gtgtTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACATTGTTGGTAACCATGCTTTCACCACCGACGAAc     >  1:194326/1‑93 (MQ=255)

                                                                                       ttGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACATTGTTGGTAACCATGCTTTCACCACCGACGAACTGa  <  2:4003897/93‑1 (MQ=255)

                                                                                       ttGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACATTGTTGGTAACCATGCTTTCACCACCGACGAACTGa  <  1:1861791/93‑1 (MQ=255)

AGAAAGGTCTGGAAGACTTCTACTACAGCGTCGGTAAATATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGACCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACATTGTTGGTAACCATGCTTTCACCACCGACGAACTGA  >  NC_000913/198328‑198507

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