BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,146,005	T→C	W104R (TGG→CGG)	b2069 →	YegD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,146,005	0	T	C	97.3%	557.4 / 2.5	185	W104R (TGG→CGG)	b2069	YegD
Reads supporting (aligned to +/- strand): ref base T (3/2); new base C (101/79); total (104/81)
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.58e-01

Only 100 of 192 total aligned reads displayed.

TAACCGCGAAGAAGATATCGATGTTACGGCGAAAAGCGTGCAGTTCGGTCTTTCCTCACTGGCACAGTACATTGATGATCCAGAAGAAGTGTGGTTTGTGAAATCACCAAAATCGTTCCTCGGTGCCAGCGGCTTAAAACCGCAGCAGGTAGCGCTGTTTGAGGATCTGGTCTGCGCAATGA  >  NC_000913/2145914‑2146095

     gcgAAGAAGATATCGATGTTACGGCGAAAAGCGTGCAGTTCGGTCTTTCCTCACTGGCACAGTACATTGATGATCCAGAAGAAGTGCGGTTtg                                                                                      <  1:4543303/93‑1 (MQ=255)

       gaagaaGATATCGATGTTACGGCGAAAAGCGTGCAGTTCGGTCTTTCCTCACTGGCACAGTACATTGATGATCCAGAAGAAGTGCGGTTtgtg                                                                                    >  1:1589752/1‑93 (MQ=255)

        aagaagATATCGATGTTACGGCGAAAAGCGTGCAGTTCGGTCTTTCCTCACTGGCACAGTACATTGATGATCCAGAAGAAGTGCGGTTTGTGa                                                                                   >  2:2215043/1‑93 (MQ=255)

        aagaagATATCGATGTTACGGCGAAAAGCGTGCAGTTCGGTCTTTCCTCACTGGCACAGTACATTGATGATCCAGAAGAAGTGCGGTTTGTGa                                                                                   <  2:3350714/93‑1 (MQ=255)

         agaagaTATCGATGTTACGGCGAAAAGCGTGCAGTTCGGTCTTTCCTCACTGGCACAGTACATTGATGATCCAGAAGAAGTGCGGTTTGTGaa                                                                                  >  1:4624416/1‑93 (MQ=255)

         agaagaTATCGATGTTACGGCGAAAAGCGTGCAGTTCGGTCTTTCCTCACTGGCACAGTACATTGATGATCCAGAAGAAGTGCGGTTTGTGaa                                                                                  >  1:51101/1‑93 (MQ=255)

          gaagaTATCGATGTTACGGCGAAAAGCGTGCAGTTCGGTCTTTCCTCACTGGCACAGTACATTGATGATCCAGAAGAAGTGCGGTTTGTGaaa                                                                                 >  1:2282080/1‑93 (MQ=255)

          gaagaTATCGATGTTACGGCGAAAAGCGTGCAGTTCGGTCTTTCCTCACTGGCACAGTACATTGATGATCCAGAAGAAGTGCGGTTTGTGaaa                                                                                 >  2:1668660/1‑93 (MQ=255)

           aagaTATCGATGTTACGGCGAAAAGCGTGCAGTTCGGTCTTTCCTCACTGGCACAGTACATTGATGATCCAGAAGAAGTGCGGTTTGTGAAAt                                                                                <  2:2473750/93‑1 (MQ=255)

           aagaTATCGATGTTACGGCGAAAAGCGTGCAGTTCGGTCTTTCCTCACTGGCACAGTACATTGATGATCCAGAAGAAGTGCGGTTTGTGAAAt                                                                                <  2:3320523/93‑1 (MQ=255)

            agaTATCGATGTTACGGCGAAAAGCGTGCAGTTCGGTCTTTCCTCACTGGCACAGTACATTGATGATCCAGAAGAAGTGCGGTTTGTGAAATc                                                                               >  2:2932147/1‑93 (MQ=255)

            agaTATCGATGTTACGGCGAAAAGCGTGCAGTTCGGTCTTTCCTCACTGGCACAGTACATTGATGATCCAGAAGAAGTGCGGTTTGTGAAATc                                                                               >  1:4337505/1‑93 (MQ=255)

             gaTATCGATGTTACGGCGAAAAGCGTGCAGTTCGGTCTTTCCTCACTGGCACAGTACATTGATGATCCAGAAGAAGTGCGGTTTGTGAAATCa                                                                              >  2:2079694/1‑93 (MQ=255)

             gaTATCGATGTTACGGCGAAAAGCGTGCAGTTCGGTCTTTCCTCACTGGCACAGTACATTGATGATCCAGAAGAAGTGCGGTTTGTGAAATCa                                                                              >  2:648631/1‑93 (MQ=255)

              atatCGATGTTACGGCGAAAAGCGTGCAGTTCGGTCTTTCCTCACTGGCACAGTACATTGATGATCCAGAAGAAGTGCGGTTTGTGAAATCAc                                                                             >  2:4607290/1‑93 (MQ=255)

               tatCGATGTTACGGCGAAAAGCGTGCAGTTCGGTCTTTCCTCACTGGCACAGTACATTGATGATCCAGAAGAAGTGCGGTTTGTGAAATCAcc                                                                            >  2:2770883/1‑93 (MQ=255)

               tatCGATGTTACGGCGAAAAGCGTGCAGTTCGGTCTTTCCTCACTGGCACAGTACATTGATGATCCAGAAGAAGTGCGGTTTGTGAAATCAcc                                                                            >  2:388709/1‑93 (MQ=255)

                 tCGATGTTACGGCGAAAAGCGTGCAGTTCGGTCTTTCCTCACTGGCACAGTACATTGATGATCCAGAAGAAGTGCGGTTTGTGAAATCACCaa                                                                          >  2:3456234/1‑93 (MQ=255)

                   gATGTTACGGCGAAAAGCGTGCAGTTCGGTCTTTCCTCACTGGCACAGTACATTGATGATCCAGAAGAAGTGCGGTTTGTGAAATCACCaaaa                                                                        >  2:503842/1‑93 (MQ=255)

                    aTGTTACGGCGAAAAGCGTGCAGTTCGGTCTTTCCTCACTGGCACAGTACATTGATGATCCAGAAGAAGTGCGGTTTGTGAAATCACCAAAAt                                                                       <  2:4393327/93‑1 (MQ=255)

                    aTGTTACGGCGAAAAGCGTGCAGTTCGGTCTTTCCTCACTGGCACAGTACATTGATGATCCAGAAGAAGTGCGGTTTGTGAAATCACCAAAAt                                                                       >  1:199342/1‑93 (MQ=255)

                    aTGTTACGGCGAAAAGCGTGCAGTTCGGTCTTTCCTCACTGGCACAGTACATTGATGATCCAGAAGAAGTGCGGTTTGTGAAATCACCAAAAt                                                                       >  2:2832056/1‑93 (MQ=255)

                    aTGTTACGGCGAAAAGCGTGCAGTTCGGTCTTTCCTCACTGGCACAGTACATTGATGATCCAGAAGAAGTGCGGTTTGTGAAATCACCAAAAt                                                                       <  2:4241359/93‑1 (MQ=255)

                    aTGTTACGGCGAAAAGCGTGCAGTTCGGTCTTTCCTCACTGGCACAGTACATTGATGATCCAGAAGAAGTGCGGTTTGTGAAATCACCAAAAt                                                                       <  1:3881900/93‑1 (MQ=255)

                     tGTTACGGCGAAAAGCGTGCAGTTCGGTCTTTCCTCACTGGCACAGTACATTGATGATCCAGAAGAAGTGTGGTTTGTGAAATCACCAAAATc                                                                      <  2:3718206/93‑1 (MQ=255)

                     tGTTACGGCGAAAAGCGTGCAGTTCGGTCTTTCCTCACTGGCACAGTACATTGATGATCCAGAAGAAGTGCGGTTTGTGAAATCACCAAAATc                                                                      >  2:2420552/1‑93 (MQ=255)

                     tGTTACGGCGAAAAGCGTGCAGTTCGGTCTTTCCTCACTGGCACAGTACATTGATGATCCAGAAGAAGTGCGGTTTGTGAAATCACCAAAATc                                                                      <  1:388427/93‑1 (MQ=255)

                      gTTACGGCGAAAAGCGTGCAGTTCGGTCTTTCCTCACTGGCACAGTACATTGATGATCCAGAAGAAGTGCGGTTTGTGAAATCACCAAAATCg                                                                     >  1:2614354/1‑93 (MQ=255)

                       ttACGGCGAAAAGCGTGCAGTTCGGTCTTTCCTCACTGGCACAGTACATTGATGATCCAGAAGAAGTGCGGTTTGTGAAATCACCAAAATCGt                                                                    <  2:1119076/93‑1 (MQ=255)

                          cGGCGAAAAGCGTGCAGTTCGGTCTTTCCTCACTGGCACAGTACATTGATGATCCAGAAGAAGTGCGGTTTGTGAAATCACCAAAATCGTTcc                                                                 >  2:4355718/1‑93 (MQ=255)

                           ggCGAAAAGCGTGCAGTTCGGTCTTTCCTCACTGGCACAGTACATTGATGATCCAGAAGAAGTGCGGTTTGTGAAATCACCAAAATCGTTCCt                                                                <  1:4218973/93‑1 (MQ=255)

                           cgCGAAAAGCGTGCAGTTCGGTCTTTCCTCACTAGCACAGTACATTGATGATCCAGAAGAAGTGCGGTTTGTGAAATCACCAAAATCGTTCCt                                                                <  1:1165734/92‑1 (MQ=255)

                             cGAAAAGCGTGCAGTTCGGTCTTTCCTCACTGGCACAGTACATTGATGATCCAGAAGAAGTGCGGTTTGTGAAATCACCAAAATCGTTCCTCg                                                              <  1:477830/93‑1 (MQ=255)

                              gAAAAGCGTGCAGTTCGGTCTTTCCTCACTGGCACAGTACATTGATGATCCAGAAGAAGTGCGGTTTGTGAAATCACCAAAATCGTTCCTCgg                                                             >  1:2910289/1‑93 (MQ=255)

                              gAAAAGCGTGCAGTTCGGTCTTTCCTCACTGGCACAGTACATTGATGATCCAGAAGAAGTGCGGTTTGTGAAATCACCAAAATCGTTCCTCgg                                                             >  1:2904885/1‑93 (MQ=255)

                              gAAAAGCGTGCAGTTCGGTCTTTCCTCACTGGCACAGTACATTGATGATCCAGAAGAAGTGCGGTTTGTGAAATCACCAAAATCGTTCCTCgg                                                             >  2:441773/1‑93 (MQ=255)

                               aaaaGCGTGCAGTTCGGTCTTTCCTCACTGGCACAGTACATTGATGATCCAGAAGAAGTGCGGTTTGTGAAATCACCAAAATCGTTCCTCGGt                                                            >  2:4011681/1‑93 (MQ=255)

                               aaaaGCGTGCAGTTCGGTCTTTCCTCACTGGCACAGTACATTGATGATCCAGAAGAAGTGCGGTTTGTGAAATCACCAAAATCGTTCCTCGGt                                                            >  2:3634085/1‑93 (MQ=255)

                               aaaaGCGTGCAGTTCGGTCTTTCCTCACTGGCACAGTACATTGATGATCCAGAAGAAGTGCGGTTTGTGAAATCACCAAAATCGTTCCTCGGt                                                            <  1:4303198/93‑1 (MQ=255)

                                     gTGCAGTTCGGTCTTTCCTCACTGGCACAGTACATTGATGATCCAGAAGAAGTGCGGTTTGTGAAATCACCAAAATCGTTCCTCGGTGCCAGc                                                      >  1:3100713/1‑93 (MQ=255)

                                      tGCAGTTCGGTCTTTCCTCACTGGCACAGTACATTGATGATCCAGAAGAAGTGCGGTTTGTGAAATCACCAAAATCGTTCCTCGGTGCCAGCg                                                     >  1:1404002/1‑93 (MQ=255)

                                      tGCAGTTCGGTCTTTCCTCACTGGCACAGTACATTGATGATCCAGAAGAAGTGCGGTTTGTGAAATCACCAAAATCGTTCCTCGGTGCCAGCg                                                     >  1:1242171/1‑93 (MQ=255)

                                         aGTTCGGTCTTTCCTCACTGGCACAGTACATTGATGATCCAGAAGAAGTGCGGTTTGTGAAATCACCAAAATCGTTCCTCGGTGCCAGCGGCt                                                  <  2:2968609/93‑1 (MQ=255)

                                         aGTTCGGTCTTTCCTCACTGGCACAGTACATTGATGATCCAGAAGAAGTGCGGTTTGTGAAATCACCAAAATCGTTCCTCGGTGCCAGCGGCt                                                  <  2:642259/93‑1 (MQ=255)

                                         aGTTCGGTCTTTCCTCACTGGCACAGTACATTGATGATCCAGAAGAAGTGCGGTTTGTGAAATCACCAAAATCGTTCCTCGGTGCCAGCGGCt                                                  <  1:2800007/93‑1 (MQ=255)

                                           ttCGGTCTTTCCTCACTGGCACAGTACATTGATGATCCAGAAGAAGTGCGGTTTGTGAAATCACCAAAATCGTTCCTCGGTGCCAGCGGCTTa                                                >  2:2397201/1‑93 (MQ=255)

                                             cGGTCTTTCCTCACTGGCACAGTACATTGATGATCCAGAAGAAGTGCGGTTTGTGAAATCACCAAAATCTTTCCTCGGTGCCAGCGGCTTaaa                                              <  2:3096316/93‑1 (MQ=255)

                                             cGGTCTTTCCTCACTGGCACAGTACATTGATGATCCAGAAGAAGTGCGGTTTGTGAAATCACCAAAATCGTTCCTCGGTGCCAGCGGCTTaaa                                              <  1:4182877/93‑1 (MQ=255)

                                              ggTCTTTCCTCACTGGCACAGTACATTGATGATCCAGAAGAAGTGCGGTTTGTGAAATCACCAAAATCGTTCCTCGGTGCCAGCGGCTTaaaa                                             >  1:3643020/1‑93 (MQ=255)

                                               gTCTTTCCTCACTGGCACAGTACATTGATGATCCAGAAGAAGTGCGGTTTGTGAAATCACCAAAATCGTTCCTCGGTGCCAGCGGCTTAAAAc                                            >  1:2016753/1‑93 (MQ=255)

                                               gTCTTTCCTCACTGGCACAGTACATTGATGATCCAGAAGAAGTGCGGTTTGTGAAATCACCAAAATCGTTCCTCGGTGCCAGCGGCTTAAAAc                                            <  2:4151002/93‑1 (MQ=255)

                                               gTCTTTCCTCACTGGCACAGTACATTGATGATCCAGAAGAAGTGCGGTTTGTGAAATCACCAAAATCGTTCCTCGGTGCCAGCGGCTTAAAAc                                            <  1:4113951/93‑1 (MQ=255)

                                                tCTTTCCTCACTGGCACAGTACATTGATGATCCAGAAGAAGTGCGGTTTGTGAAATCACCAAAATCGTTCCTCGGTGCCAGCGGCTTAAAAcc                                           <  2:1550778/93‑1 (MQ=255)

                                                  tttCCTCACTGGCACAGTACATTGATGATCCAGAAGAAGTGCGGTTTGTGAAATCACCAAAATCGTTCCTCGGTGCCAGCGGCTTAAAACCgc                                         >  2:3602875/1‑93 (MQ=255)

                                                   ttCCTCACTGGCACAGTACATTGATGATCCAGAAGAAGTGCGGTTTGTGAAATCACCAAAATCGTTCCTCGGTGCCAGCGGCTTACAACCgca                                        >  2:513554/1‑93 (MQ=255)

                                                   ttCCTCACTGGCACAGTACATTGATGATCCAGAAGAAGTGCGGTTTGTGAAATCACCAAAATCGTTCCTCGGTGCCAGCGGCTTAAAACCgca                                        >  1:2876694/1‑93 (MQ=255)

                                                   ttCCTCACTGGCACAGTACATTGATGATCCAGAAGAAGTGCGGTTTGTGAAATCACCAAAATCGTTCCTCGGTGCCAGCGGCTTAAAACCgca                                        <  2:3274586/93‑1 (MQ=255)

                                                      cTCACTGGCACAGTACATTGATGATCCAGAAGAAGTGCGGTTTGTGAAATCACCAAAATCGTTCCTCGGTGCCAGCGGCTTAAAACCgcagca                                     <  2:562199/93‑1 (MQ=255)

                                                      cTCACTGGCACAGTACATTGATGATCCAGAAGAAGTGCGGTTTGTGAAATCACCAAAATCGTTCCTCGGTGCCAGCGGCTTAAAACCgcagca                                     >  1:926199/1‑93 (MQ=255)

                                                        cACTGGCACAGTACATTGATGATCCAGAAGAAGTGCGGTTTGTGAAATCACCAAAATCGTTCCTCGGTGCCAGCGGCTTAAAACCGCAGCAgg                                   >  1:3695754/1‑93 (MQ=255)

                                                         aCTGGCACAGTACATTGATGATCCAGAAGAAGTGCGGTTTGTGAAATCACCAAAATCGTTCCTCGGTGCCAGCGGCTTAAAACCGCAGCAGGt                                  >  1:3743411/1‑93 (MQ=255)

                                                         aCTGGCACAGTACATTGATGATCCAGAAGAAGTGCGGTTTGTGAAATCACCAAAATCGTTCCTCGGTGCCAGCGGCTTAAAACCGCAGCAGGt                                  >  1:3235004/1‑93 (MQ=255)

                                                         aCTGGCACAGTACATTGATGATCCAGAAGAAGTGCGGTTTGTGAAATCACCAAAATCGTTCCTCGGTGCCAGCGGCTTAAAACCGCAGCAGGt                                  <  1:3180424/93‑1 (MQ=255)

                                                          cTGGCACAGTACATTGATGATCCAGAAGAAGTGCGGTTTGTGAAATCACCAAAATCGTTCCTCGGTGCCAGCGGCTTAAAACCGCAGCAGGTa                                 <  2:231329/93‑1 (MQ=255)

                                                              cacaGTACATTGATGATCCAGAAGAAGTGCGGTTTGTGAAATCACCAAAATCGTTCCTCGGTGCCAGCGGCTTAAAACCGCAGCAGGTAgcgc                             <  1:73638/93‑1 (MQ=255)

                                                                caGTACATTGATGATCCAGAAGAAGTGCGGTTTGTGAAATCACCAAAATCGTTCCTCGGTGCCAGCGGCTTAAAACCGCAGCAGGTAGCGCTg                           >  2:889835/1‑93 (MQ=255)

                                                                caGTACATTGATGATCCAGAAGAAGTGCGGTTTGTGAAATCACCAAAATCGTTCCTCGGTGCCAGCGGCTTAAAACCGCAGCAGGTAGCGCTg                           >  1:337267/1‑93 (MQ=255)

                                                                caGTACATTGATGATCCAGAAGAAGTGCGGTTTGTGAAATCACCAAAATCGTTCCTCGGTGCCAGCGGCTTAAAACCGCAGCAGGTAGCGCTg                           >  2:1613339/1‑93 (MQ=255)

                                                                 aGTACATTGATGATCCAGAAGAAGTGCGGTTTGTGAAATCACCAAAATCGTTCCTCGGTGCCAGCGGCTTAAAACCGCAGCAGGTAGCGCTGt                          <  2:2506934/93‑1 (MQ=255)

                                                                  gTACATTGATGATCCAGAAGAAGTGCGGTTTGTGAAATCACCAAAATCGTTCCTCGGTGCCAGCGGCTTAAAACCGCAGCAGGTAGCGCTGtt                         >  2:1243245/1‑93 (MQ=255)

                                                                  gTACATTGATGATCCAGAAGAAGTGCGGTTTGTGAAATCACCAAAATCGTTCCTCGGTGCCAGCGGCTTAAAACCGCAGCAGGTAGCGCTGtt                         <  2:4460003/93‑1 (MQ=255)

                                                                   tACATTGATGATCCAGAAGAAGTGCGGTTTGTGAAATCACCAAAATCGTTCCTCGGTGCCAGCGGCTTAAAACCGCAGCAGGTAGCGCTGttt                        >  2:1023605/1‑93 (MQ=255)

                                                                   tACATTGATGATCCAGAAGAAGTGCGGTTTGTGAAATCACCAAAATCGTTCCTCGGTGCCAGCGGCTTAAAACCGCAGCAGGTAGCGCTGttt                        >  1:3805565/1‑93 (MQ=255)

                                                                    aCATTGATGATCCAGAAGAAGTGCGGTTTGTGAAATCACCAAAATCGTTCCTCGGTGCCAGCGGCTTAAAACCGCAGCAGGTAGCGCTGTTTg                       >  1:4002704/1‑93 (MQ=255)

                                                                     cATTGATGATCCAGAAGAAGTGCGGTTTGTGAAATCACCAAAATCGTTCCTCGGTGCCAGCGGCTTAAAACCGCAGCAGGTAGCGCTGTTTGa                      >  2:2426301/1‑93 (MQ=255)

                                                                     cATTGATGATCCAGAAGAAGTGCGGTTTGTGAAATCACCAAAATCGTTCCTCGGTGCCAGCGGCTTAAAACCGCAGCAGGTAGCGCTGTTTGa                      >  2:2131856/1‑93 (MQ=255)

                                                                      aTTGATGATCCAGAAGAAGTGTGGTTTGTGAAATCACCAAAATCGTTCCTCGGTGCCAGCGGCTTAAAACCGCAGCAGGTAGCGCTGTTTGAg                     <  1:338828/93‑1 (MQ=255)

                                                                      aTTGATGATCCAGAAGAAGTGCGGTTTGTGAAATCACCAAAATCGTTCCTCGGTGCCAGCGGCTTAAAACCGCAGCAGGTAGCGCTGTTTGAg                     >  1:1524366/1‑93 (MQ=255)

                                                                       ttGATGATCCAGAAGAAGTGCGGTTTGTGAAATCACCAAAATCGTTCCTCGGTGCCAGCGGCTTAAAACCGCAGCAGGTAGCGCTGTTTGAgg                    >  1:4465994/1‑93 (MQ=255)

                                                                       ttGATGATCCAGAAGAAGTGCGGTTTGTGAAATCACCAAAATCGTTCCTCGGTGCCAGCGGCTTAAAACCGCAGCAGGTAGCGCTGTTTGAgg                    >  1:4140101/1‑93 (MQ=255)

                                                                       ttGATGATCCAGAAGAAGTGCGGTTTGTGAAATCACCAAAATCGTTCCTCGGTGCCAGCGGCTTAAAACCGCAGCAGGTAGCGCTGTTTGAgg                    >  1:2164460/1‑93 (MQ=255)

                                                                        tgatgaTCCAGAAGAAGTGCGGTTTGTGAAATCACCAAAATCGTTCCTCGGTGCCAGCGGCTTAAAACCGCAGCAGGTAGCGCTGTTTGAGGa                   <  2:744868/93‑1 (MQ=255)

                                                                         gatgatCCAGAAGAAGTGCGGTTTGTGAAATCACCAAAATCGTTCCTCGGTGCCAGCGGCTTAAAACCGCAGCAGGTAGCGCTGTTTGAGGAt                  >  1:2346189/1‑93 (MQ=255)

                                                                         gatgatCCAGAAGAAGTGCGGTTTGTGAAATCACCAAAATCGTTCCTCGGTGCCAGCGGCTTAAAACCGCAGCAGGTAGCGCTGTTTGAGGAt                  <  2:3781673/93‑1 (MQ=255)

                                                                               ccAGAAGAAGTGCGGTTTGTGAAATCACCAAAATCGTTCCTCGGTGCCAGCGGCTTAAAACCGCAGCAGGTAGCGCTGTTTGAGGATCTGGTc            <  1:223443/93‑1 (MQ=255)

                                                                               ccAGAAGAAGTGCGGTTTGTGAAATCACCAAAATCGTTCCTCGGTGCCAGCGGCTTAAAACCGCAGCAGGTAGCGCTGTTTGAGGATCTGGTc            <  1:114547/93‑1 (MQ=255)

                                                                                cAGAAGAAGTGCGGTTTGTGAAATCACCAAAATCGTTCCTCGGTGCCAGCGGCTTAAAACCGCAGCAGGTAGCGCTGTTTGAGGATCTGGTCt           >  1:987628/1‑93 (MQ=255)

                                                                                 agaagaAGTGCGGTTTGTGAAATCACCAAAATCGTTCCTCGGTGCCAGCGGCTTAAAACCGCAGCAGGTAGCGCTGTTTGAGGATCTGGTCTg          >  1:1935659/1‑93 (MQ=255)

                                                                                  gaagaaGTGCGGTTTGTGAAATCACCAAAATCGTTCCTCGGTGCCAGCGGCTTAAAACCGCAGCAGGTAGCGCTGTTTGAGGATCTGGTCTgc         >  2:4200975/1‑93 (MQ=255)

                                                                                  gaagaaGTGCGGTTTGTGAAATCACCAAAATCGTTCCTCGGTGCCAGCGGCTTAAAACCGCAGCAGGTAGCGCTGTTTGAGGATCTGGTCTgc         >  1:1837564/1‑93 (MQ=255)

                                                                                  gaagaaGTGCGGTTTGTGAAATCACCAAAATCGTTCCTCGGTGCCAGCGGCTTAAAACCGCAGCAGGTAGCGCTGTTTGAGGATCTGGTCTgc         >  2:2899787/1‑93 (MQ=255)

                                                                                    agaagTGCGGTTTGTGAAATCACCAAAATCGTTCCTCGGTGCCAGCGGCTTAAAACCGCAGCAGGTAGCGCTGTTTGAGGATCTGGTCTgcgc       >  1:249206/1‑93 (MQ=255)

                                                                                      aagTGCGGTTTGTGAAATCACCAAAATCGTTCCTCGGTGCCAGCGGCTTAAAACCGCAGCAGGTAGCGCTGTTTGAGGATCTGGTCTGCGCaa     <  1:2666653/93‑1 (MQ=255)

                                                                                       agTGCGGTTTGTGAAATCACCAAAATCGTTCCTCGGTGCCAGCGGCTTAAAACCGCAGCAGGTAGCGCTGTTTGAGGATCTGGTCTGCGCAat    >  1:8222/1‑93 (MQ=255)

                                                                                       agTGCGGTTTGTGAAATCACCAAAATCGTTCCTCGGTGCCAGCGGCTTAAAACCGCAGCAGGTAGCGCTGTTTGAGGATCTGGTCTGCGCAat    <  2:541534/93‑1 (MQ=255)

                                                                                        gtgcGGTTTGTGAAATCACCAAAATCGTTCCTCGGTGCCAGCGGCTTAAAACCGCAGCAGGTAGCGCTGTTTGAGGATCTGGTCTGCGCAatg   >  1:4019593/1‑93 (MQ=255)

                                                                                        gtgcGGTTTGTGAAATCACCAAAATCGTTCCTCGGTGCCAGCGGCTTAAAACCGCAGCAGGTAGCGCTGTTTGAGGATCTGGTCTGCGCAatg   >  1:630254/1‑93 (MQ=255)

                                                                                        gtgcGGTTTGTGAAATCACCAAAATCGTTCCTCGGTGCCAGCGGCTTAAAACCGCAGCAGGTAGCGCTGTTTGAGGATCTGGTCTGCGCAatg   >  1:1767253/1‑93 (MQ=255)

                                                                                        gtgcGGTTTGTGAAATCACCAAAATCGTTCCTCGGTGCCAGCGGCTTAAAACCGCAGCAGGTAGCGCTGTTTGAGGATCTGGTCTGCGCAatg   <  2:2426875/93‑1 (MQ=255)

                                                                                         tgtgGTTTGTGAAATCACCAAAATCGTTCCTCGGTGCCAGCGGCTTAAAACCGCAGCAGGTAGCGCTGTTTGAGGATCTGGTCTGCGCAatga  <  2:2750154/93‑1 (MQ=255)

TAACCGCGAAGAAGATATCGATGTTACGGCGAAAAGCGTGCAGTTCGGTCTTTCCTCACTGGCACAGTACATTGATGATCCAGAAGAAGTGTGGTTTGTGAAATCACCAAAATCGTTCCTCGGTGCCAGCGGCTTAAAACCGCAGCAGGTAGCGCTGTTTGAGGATCTGGTCTGCGCAATGA  >  NC_000913/2145914‑2146095

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: `‑`