BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,042,457	C→A	91.2%	G30G (GGC→GGA)	yodB →	cytochrome b561 homolog

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,042,457	0	C	A	91.2%	363.5 / 23.3	123	G30G (GGC→GGA)	yodB	cytochrome b561 homolog
Reads supporting (aligned to +/- strand): ref base C (5/6); new base A (48/63); total (53/70)
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.99e-01

Only 100 of 123 total aligned reads displayed.

GGAATTTTTAAGGTTCGTGGGATAGCTTGACTGTGAAAATCACAGGAGCTACAAAAATGAACCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGCTGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATAAAACACCGTTATAGTGATCCTTCTATTGTGCCACCGCCAC  >  NC_000913/2042312‑2042596

 gAATTTTTAAGGTTCGTGGGATAGCTTGACTGTGAAAATCACAGGAGCTACAAAAATGAACCTATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGt                                                                                                                                         >  6:307200/1‑149 (MQ=255)

  aaTTTTTAAGGTTCGTGGGATAGCTTGACTGTGAAAATCACAGGAGCTACAAAAATGAACCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGtt                                                                                                                                        <  2:428604/149‑1 (MQ=255)

  aaTTTTTAAGGTTCGTGGGATAGCTTGACTGTGAAAATCACAGGAGCTACAAAAATGAACCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGtt                                                                                                                                        >  8:312983/1‑149 (MQ=255)

    tttttAATGTTCGTTGGATAGCTTGACTTAGAAAATCACAGCAGCTACAAAAATGAACCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGCTGGTTTc                                                                                                                                      <  6:199084/149‑1 (MQ=255)

    tttttAAGGTTCGTGGTATAGCTTGACTGGGACAATCACAGGAGCTACAAAACTTAACCGATTCTAAAAAAATCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTTGATGGTTTc                                                                                                                                      <  6:199081/149‑1 (MQ=255)

    tttttAAGGTTCGTGGGATAGCTTGACTGTGAAAATCACAGGAGCTACAAAAATGAACCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTc                                                                                                                                      <  5:188158/149‑1 (MQ=255)

    tttttAAGGTTCGTGGGATAGCTTGACTGTGAAAATCACAGGAGCTACAAAAATGAACCGATTCTCAAAAACTCAAATTTATTTACATTGGAAAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTc                                                                                                                                      <  5:150527/149‑1 (MQ=255)

    tttttAAGGTTCGTGGGATAGCGTGACTGTGAAAATCACAGCAGCTAAAAAAATGAACCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTGTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTc                                                                                                                                      >  1:388330/1‑149 (MQ=255)

      tttAAGGTTCGTGGGATAGCTTGACTGTGAAAATCACAGGAGCTACAAAAATGAACCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTcc                                                                                                                                     >  5:389934/1‑148 (MQ=255)

            gTTCGTGGGATAGCTTGACTGTGAAAATCACAGGAGCTACAAAAATGAACCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGgt                                                                                                                              <  1:296338/149‑1 (MQ=255)

             ttCGTGGGATAGCTTGACTGTGAAAATCACAGGAGCTACAAAAATGAACCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGCTGGTTTCCTAAAGgta                                                                                                                             <  4:254370/149‑1 (MQ=255)

               cGTGGGATAGCTTGACTGTGAAAATCACAGGAGCTACAAAAATGAACCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGgtag                                                                                                                            <  8:244018/148‑1 (MQ=255)

               cGTGGGATAGCTTGACTGTGAAAATCACAGGAGCTACAAAAATGAACCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGgtag                                                                                                                            <  1:372432/148‑1 (MQ=255)

               cGTGGGATAGCTTGACTGTGAAAATCACAGGAGCTACAAAAATGAACCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGgtag                                                                                                                            <  7:274947/148‑1 (MQ=255)

                   ggATAGCTTGACTGTGAAAATCACAGGAGCTACAAAATTGAACCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGgtag                                                                                                                            <  2:28415/144‑1 (MQ=255)

                   ggATAGCTTGACTGTGAAAATCACAGGAGCTACAAAATTGAACCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGgtag                                                                                                                            >  1:28415/1‑144 (MQ=255)

                   ggATAGCTTGACTGTGAAAATCACAGGAGCTACAAAAATGAACCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGgtagtagt                                                                                                                        <  6:257426/148‑1 (MQ=255)

                   ggATAGCTTGACTGTGAAAATCACAGGAGCTACAAAAATGAACCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCCTTTCGTTGCAATAACCTATGCCGCTATGGAACTCCGTGGATGGTTTCCTAAAGgtagtagt                                                                                                                        <  2:212181/148‑1 (MQ=255)

                    gATAGCTTGACTGTGAAAATCACAGGAGCTACAAAAATGAACCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATg                                                                                                                                           <  8:341240/128‑1 (MQ=255)

                    gATAGCTTGACTGTGAAAATCACAGGAGCTACAAAAATGAACCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATg                                                                                                                                           >  7:341240/1‑128 (MQ=255)

                    gATAGCTTGACTGTGAAAATCACAGGAGCTACAAAAATGAACCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGgtag                                                                                                                            >  4:147542/1‑143 (MQ=255)

                    gATAGCTTGACTGTGAAAATCACAGGAGCTACAAAAATGAACCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGgtag                                                                                                                            <  3:147542/143‑1 (MQ=255)

                    gATAGCTTGACTGTGAAAATCACAGGAGCTACAAAAATGAACCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTAc                                                                                                                      <  5:307200/149‑1 (MQ=255)

                    gATAGCTTGACTGTGAAAATCACAGGAGCTACAAAAATGAACCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTAc                                                                                                                      >  7:342259/1‑149 (MQ=255)

                        gCTTGACTGTGAAAATCACAGGAGCTACAAAAATGAACCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCTATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTa                                                                                                                   <  5:128910/148‑1 (MQ=255)

                        gCTTGACTGTGAAAATCACAGGAGCTACAAAAATGAACCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTa                                                                                                                   <  6:389934/148‑1 (MQ=255)

                        gCTTGACTGGGAAACTCACATGAGCTACAAAAATGACCCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTCTGCCGCGATGGAACTCCGTGGATGGTTGCCTAAAGGTAGTAGTACTTAt                                                                                                                  <  2:388330/149‑1 (MQ=255)

                          ttGACTGTGAAAATCACAGGAGCTACAAAAATGAACCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATct                                                                                                                >  1:378859/1‑149 (MQ=255)

                             aCTGTGAAAATCACAGGAGCTACAAAAATGAACCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATctgc                                                                                                              <  5:204405/148‑1 (MQ=255)

                             aCTGTGAAAATCACAGGAGCTACAAAAATGAACCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATctgc                                                                                                              <  7:140636/148‑1 (MQ=255)

                              cTGTGAAAATCACAGGAGCTACAAAAATGAACCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATctgct                                                                                                             >  4:335156/1‑148 (MQ=255)

                                 agataaTCACATTAGATAAAAAAATGAACCTATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCGTTTCGTTTCAATAACCTATGCCGCGATGGAACTCCGTGTATGGTTTCCTAAAGGTAGTAGTACTTATCTTCTGAt                                                                                                          <  6:273869/144‑1 (MQ=255)

                                   aaaaTCACAGGAGCTACAAAAATGAACCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGc                                                                                                        <  8:222247/148‑1 (MQ=255)

                                    aaaTCACAGGAGCTACAAAAATGAACCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCga                                                                                                      <  7:312983/149‑1 (MQ=255)

                                     aaTCACAGGAGCTACAAAAATGAACCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTGTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCgag                                                                                                     <  2:285158/149‑1 (MQ=255)

                                         acaGGAGCTACAAAAATGAACCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGaaa                                                                                                  >  2:304442/1‑148 (MQ=255)

                                         acaGGAGCTACAAAAATGAACCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGaaa                                                                                                  >  4:176757/1‑148 (MQ=255)

                                           aGGAGCTACAAAAATGAACCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGCTGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAacac                                                                                               <  8:359844/149‑1 (MQ=255)

                                             gAGCTACAAAAATGAACCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACAt                                                                                             <  7:2316/149‑1 (MQ=255)

                                              aGCTACAAAAATGAACCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACAtt                                                                                            <  6:394873/149‑1 (MQ=255)

                                                  aCAAAAATGAACCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACaa                                                                                        >  7:357795/1‑149 (MQ=255)

                                                    aaaaaTGAACCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACAtt                                                                                            <  8:173362/143‑1 (MQ=255)

                                                      aaaTGAACCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCg                                                                                    <  3:335156/149‑1 (MQ=255)

                                                                   ctcAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGtt                                                                        <  1:304442/148‑1 (MQ=255)

                                                                    tcAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACAt                                                                                             >  2:456559/1‑126 (MQ=255)

                                                                    tcAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACAt                                                                                             <  8:17660/126‑1 (MQ=255)

                                                                    tcAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACAt                                                                                             <  1:456559/126‑1 (MQ=255)

                                                                    tcAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACAt                                                                                             >  7:17660/1‑126 (MQ=255)

                                                                    tcAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTg                                                                      >  1:145159/1‑149 (MQ=255)

                                                                      aaaaaCTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGCTGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTggg                                                                    >  6:384311/1‑149 (MQ=255)

                                                                      aaaaaCTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTggg                                                                    <  3:123886/149‑1 (MQ=255)

                                                                      aaaaaCTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTggg                                                                    <  2:428981/149‑1 (MQ=255)

                                                                       aaaaCTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGTGTATATTCGTTTGGgt                                                                   >  4:47964/1‑149 (MQ=255)

                                                                        aaaCTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCTTTTGGgtg                                                                  <  7:106256/149‑1 (MQ=255)

                                                                            tCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAgg                                                                                                                               >  1:451511/1‑84 (MQ=255)

                                                                              aaaTTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGCTGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATg                                                            >  1:13254/1‑149 (MQ=255)

                                                                              aaaTTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATg                                                            >  6:456913/1‑149 (MQ=255)

                                                                               aaTTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGCTGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGt                                                           >  8:138221/1‑149 (MQ=255)

                                                                               aaTTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGt                                                           <  5:191091/149‑1 (MQ=255)

                                                                                atttatttACATTGGATAACGCTGCTTTTCGTTGCAATAACATATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTAATTATCTGCTGATGCGAGCAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGtt                                                          >  3:53577/1‑149 (MQ=255)

                                                                                 tttatttaCATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTa                                                                                           >  7:263376/1‑115 (MQ=255)

                                                                                 tttatttaCATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTc                                                                           >  6:234729/1‑131 (MQ=255)

                                                                                 tttatttaCATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTc                                                                           <  5:234729/131‑1 (MQ=255)

                                                                                  ttatttaCATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGttt                                                         >  7:120030/1‑148 (MQ=255)

                                                                                         cATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCt                                                 <  3:199914/149‑1 (MQ=255)

                                                                                           ttGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGa                                               >  4:322508/1‑149 (MQ=255)

                                                                                             ggATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCgag                                                                                                     <  8:298827/93‑1 (MQ=255)

                                                                                              gATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATta                                            <  2:378859/149‑1 (MQ=255)

                                                                                                tAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATtat                                           >  5:373562/1‑148 (MQ=255)

                                                                                                 aaCGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAAt                                                             <  3:97080/129‑1 (MQ=255)

                                                                                                     ctgctTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGCTGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATAAAaca                                     <  5:384311/149‑1 (MQ=255)

                                                                                                              gTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATAAAACACCGTtata                             <  3:322508/148‑1 (MQ=255)

                                                                                                               ttGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCgag                                                                                                     <  4:149208/75‑1 (MQ=255)

                                                                                                               ttGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCgag                                                                                                     >  3:149208/1‑75 (MQ=255)

                                                                                                                  cAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACATTGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGcc                                                  <  3:135611/123‑1 (MQ=255)

                                                                                                                  cAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAGTGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGcc                                                  >  4:135611/1‑123 (MQ=255)

                                                                                                                  cAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGTTATATTCGTTTGTGTGTTAATGTTTTCACGCCTGATTATAAAACACCGTTATAGTGa                         <  8:120030/148‑1 (MQ=255)

                                                                                                                  cAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATAAAACACCGTTATAGTGAt                        >  4:240002/1‑149 (MQ=255)

                                                                                                                   aaTAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATAAAACACCGTTATAGTGATc                       >  7:9511/1‑149 (MQ=255)

                                                                                                                    aTAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGTGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATAAAACACCGTTATAGTGATcc                      <  3:47964/149‑1 (MQ=255)

                                                                                                                     tAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATAAAACACCg                                  >  7:92340/1‑136 (MQ=255)

                                                                                                                     tAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATAAAACACCg                                  <  8:92340/136‑1 (MQ=255)

                                                                                                                       aCCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGACGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATAAAACACCGTTATAGTGATCCTTc                   <  7:121766/149‑1 (MQ=255)

                                                                                                                        ccTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATAAAACACCGTTATAGTGATCCTTc                   >  8:263911/1‑148 (MQ=255)

                                                                                                                         cTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATAAAACACCGTTATAGTGATCCTTCTa                 >  1:344858/1‑149 (MQ=255)

                                                                                                                         cTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATAAAACACCGTTATAGTGATCCTTCTa                 <  7:263911/149‑1 (MQ=255)

                                                                                                                           aTGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATAAAACACCGTTATAGTGATCCTTCTAtt               >  4:54153/1‑149 (MQ=255)

                                                                                                                            tGCCGCGATGGAACTCCGTGGCTGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATAAAACACCGTTATAGTGATCCTTCTAt                <  2:13254/147‑1 (MQ=255)

                                                                                                                            tGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATAAAACACCGTTATAGTGATCCTTCTAtt               <  6:373562/148‑1 (MQ=255)

                                                                                                                            tGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATAAAACACCGTTATAGTGATCCTTCTAtt               <  2:145159/148‑1 (MQ=255)

                                                                                                                             gCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGcc                                                  <  4:58874/112‑1 (MQ=255)

                                                                                                                             gCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGcc                                                  >  3:58874/1‑112 (MQ=255)

                                                                                                                              ccGCGATGGAACTCCTTGGATGGTTTCCTAAATGTAGTAGTACTTATATGCTGATGCTAGAAACACATTACAATGCTGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATAAAACACCGTTATAGTGATCCTTCTATtgt             <  4:53577/148‑1 (MQ=255)

                                                                                                                              ccGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATAAAACACCGTTATAGTGATCCTTCTATtgt             <  2:119231/148‑1 (MQ=255)

                                                                                                                                gcgATGGAACTCCGTGGCTGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATAAAACACCGTTATAGTGATCCTTCTATTGTgc           >  5:253498/1‑148 (MQ=255)

                                                                                                                                gcgATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATAAAACACCGTTATAGTGATCCTTCTATTGTgc           >  4:18812/1‑148 (MQ=255)

                                                                                                                                  gATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTgg                                                                     >  8:65581/1‑88 (MQ=255)

                                                                                                                                  gATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTgg                                                                     <  7:65581/88‑1 (MQ=255)

                                                                                                                                      gAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATAAAACACCGTTATAGTGATCCTTCTATTGTgccaccgc     <  8:357795/148‑1 (MQ=255)

                                                                                                                                        aCTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATAAAACACCGTTATAGTGATCCTTCTATTGTgccaccgccac  <  3:240002/149‑1 (MQ=255)

GGAATTTTTAAGGTTCGTGGGATAGCTTGACTGTGAAAATCACAGGAGCTACAAAAATGAACCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGCTGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATAAAACACCGTTATAGTGATCCTTCTATTGTGCCACCGCCAC  >  NC_000913/2042312‑2042596

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 14 ≤` `ATCG/ATCG` `< 27 ≤` `ATCG/ATCG` `< 36 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`