BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,949,344	G→A	Q174Q (CAG→CAA)	yecD →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,949,344	0	G	A	78.7%	1986.1 / inf	878	Q174Q (CAG→CAA)	yecD	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (82/105); new base A (385/305); total (467/411)
Fisher's exact test for biased strand distribution p-value = 3.86e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.19e-02

Only 100 of 945 total aligned reads displayed.

TCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTA  >  NC_000913/1949253‑1949436

 cGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAgcagca                                                                                            >  2:1674388/1‑93 (MQ=255)

  ggTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGACGACGCCTGTAGTGCCGCCAGCGCCGAGCAacac                                                                                           >  2:5939671/1‑93 (MQ=255)

  ggTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAacac                                                                                           <  2:5339748/93‑1 (MQ=255)

   gtgtTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAAcaca                                                                                          >  1:3626361/1‑93 (MQ=255)

    tgtTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACaa                                                                                         <  2:3631834/93‑1 (MQ=255)

     gtTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAAt                                                                                        >  1:549488/1‑93 (MQ=255)

     gtTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAAt                                                                                        >  2:2965824/1‑93 (MQ=255)

     gtTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAAt                                                                                        <  2:2251899/93‑1 (MQ=255)

      ttGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCAAGCAACACAATa                                                                                       >  1:5356470/1‑93 (MQ=255)

       tGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATaa                                                                                      >  1:2729844/1‑93 (MQ=255)

        gAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAAc                                                                                     >  2:4120953/1‑93 (MQ=255)

         aaTCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACa                                                                                    <  2:1950929/93‑1 (MQ=255)

          aTCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAg                                                                                   <  1:2454395/93‑1 (MQ=255)

           tCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAgc                                                                                  >  1:1438415/1‑92 (MQ=255)

            ccACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTa                                                                                 <  1:2926036/93‑1 (MQ=255)

            ccACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTa                                                                                 <  1:4540723/93‑1 (MQ=255)

            ccACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAgca                                                                                 <  1:1672501/93‑3 (MQ=255)

             cACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTAt                                                                                <  2:3204951/93‑1 (MQ=255)

                    cgcAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCat                                                                         >  2:2409699/1‑93 (MQ=255)

                     gcAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCata                                                                        <  1:384200/93‑1 (MQ=255)

                     gcAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCata                                                                        <  2:5011814/93‑1 (MQ=255)

                      cAATGCCTGGTAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCatat                                                                       <  1:2601891/93‑1 (MQ=255)

                      cAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCatat                                                                       >  2:5537091/1‑93 (MQ=255)

                       aaTGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATc                                                                      <  2:4209238/93‑1 (MQ=255)

                       aaTGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATc                                                                      <  1:4691338/93‑1 (MQ=255)

                        aTGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCt                                                                     >  1:4131166/1‑93 (MQ=255)

                         tGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTa                                                                    <  1:5282379/93‑1 (MQ=255)

                         tGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTa                                                                    <  2:5460320/93‑1 (MQ=255)

                         tGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTa                                                                    <  1:351660/93‑1 (MQ=255)

                         tGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTa                                                                    >  2:4626551/1‑93 (MQ=255)

                           ccTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTAcc                                                                  >  2:836837/1‑93 (MQ=255)

                            cggggACCTTGGCTTTAATCTGGTGATTGCCTAAGACGCCTGTAGTGCCGCCAGCCCCGAGCAACACAATAATAGCATTAATCATATCTTccc                                                                 <  1:491135/91‑1 (MQ=255)

                            cTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTAccc                                                                 >  1:4000982/1‑93 (MQ=255)

                            cTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTAccc                                                                 <  2:3727499/93‑1 (MQ=255)

                            cTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTAccc                                                                 >  1:4798124/1‑93 (MQ=255)

                             tGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCcg                                                                >  1:5243215/1‑93 (MQ=255)

                             tGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCcg                                                                >  2:1109606/1‑93 (MQ=255)

                             tGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCcg                                                                <  2:1571518/93‑1 (MQ=255)

                             tGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCACAGCCGAGAAACACAATAACAGCATTAATCATATCCACCcg                                                                >  1:581995/1‑93 (MQ=255)

                              gggAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCcgc                                                               >  2:5659860/1‑93 (MQ=255)

                              gggAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCcgc                                                               >  1:5583942/1‑93 (MQ=255)

                              gggAACTTGGCTTCAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCcgc                                                               >  2:4869870/1‑93 (MQ=255)

                                gAACTTGGCTTTAATCTGGGTGTTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCgcgc                                                             <  2:4765060/93‑1 (MQ=255)

                                 aaCTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCa                                                            <  2:4473375/93‑1 (MQ=255)

                                    ttGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCg                                                         >  1:5163445/1‑93 (MQ=255)

                                    ttGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCg                                                         >  2:2714766/1‑93 (MQ=255)

                                    ttGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCg                                                         >  1:5949609/1‑93 (MQ=255)

                                    ttGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCg                                                         >  1:604865/1‑93 (MQ=255)

                                     tGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGc                                                        >  1:4255869/1‑93 (MQ=255)

                                       gCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGccc                                                      <  1:5699021/93‑1 (MQ=255)

                                       gCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGccc                                                      >  1:2346347/1‑93 (MQ=255)

                                        cTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCg                                                     <  2:5140443/93‑1 (MQ=255)

                                        cTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCg                                                     <  2:3105195/93‑1 (MQ=255)

                                        cTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCg                                                     <  2:1202339/93‑1 (MQ=255)

                                        cTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCg                                                     <  1:2826519/93‑1 (MQ=255)

                                           tAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCgtgt                                                  >  2:6112922/1‑93 (MQ=255)

                                            aaTCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGtgtg                                                 <  1:5088634/93‑1 (MQ=255)

                                             aTCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGc                                                >  2:4561205/1‑93 (MQ=255)

                                             aTCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGc                                                >  2:401333/1‑93 (MQ=255)

                                              tCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCg                                               <  2:3840720/93‑1 (MQ=255)

                                               cTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGt                                              >  1:4400574/1‑93 (MQ=255)

                                               cTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGt                                              <  1:253624/93‑1 (MQ=255)

                                               cTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGt                                              <  1:681767/93‑1 (MQ=255)

                                                tGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTa                                             <  1:2058996/93‑1 (MQ=255)

                                                tGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTa                                             >  2:5224387/1‑93 (MQ=255)

                                                 ggTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAg                                            <  1:1968141/93‑1 (MQ=255)

                                                 ggTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAg                                            >  2:5789343/1‑93 (MQ=255)

                                                      ttGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTgg                                       >  2:2808102/1‑93 (MQ=255)

                                                       tGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGa                                      <  1:1202017/93‑1 (MQ=255)

                                                          cGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAaga                                   <  2:2051199/93‑1 (MQ=255)

                                                              gACGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATc                               <  2:2880829/93‑1 (MQ=255)

                                                               aCGCGTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATcc                              <  2:1317268/93‑1 (MQ=255)

                                                               aCGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATcc                              <  1:3026887/93‑1 (MQ=255)

                                                               aCGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATcc                              <  1:5027889/93‑1 (MQ=255)

                                                                cGCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCt                             <  1:5420928/93‑1 (MQ=255)

                                                                 gCCTGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTc                            >  2:1461688/1‑93 (MQ=255)

                                                                    tGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCttcac                         >  2:932212/1‑89 (MQ=255)

                                                                    tGTAGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGCGCGTAGCGTGGAAGAGATCCttcac                         >  2:258248/1‑89 (MQ=255)

                                                                       aGTGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCAcgcg                      >  2:1746531/1‑93 (MQ=255)

                                                                         tGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGtt                    >  1:3170717/1‑93 (MQ=255)

                                                                            cgccAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATg                 <  2:4364165/93‑1 (MQ=255)

                                                                             gccAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGa                <  1:3403970/93‑1 (MQ=255)

                                                                              ccAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGAt               <  2:690227/93‑1 (MQ=255)

                                                                                 gcgcCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTa            <  1:1007673/93‑1 (MQ=255)

                                                                                 gcgcCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTa            <  2:181909/93‑1 (MQ=255)

                                                                                   gcCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACa          <  1:21699/93‑1 (MQ=255)

                                                                                   gcCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACa          <  2:1402003/93‑1 (MQ=255)

                                                                                   gcCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACa          <  1:3333317/93‑1 (MQ=255)

                                                                                    ccGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACAt         <  1:1097367/93‑1 (MQ=255)

                                                                                    ccGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACAt         <  1:97899/93‑1 (MQ=255)

                                                                                    ccGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACAt         <  2:586890/93‑1 (MQ=255)

                                                                                    ccGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACAt         <  1:1354512/93‑1 (MQ=255)

                                                                                       agcagcACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCgg      >  2:3436520/1‑93 (MQ=255)

                                                                                       agcaacACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCgg      >  1:3712419/1‑93 (MQ=255)

                                                                                       agcaacACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCgg      <  1:155527/93‑1 (MQ=255)

                                                                                        gcagcaCAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGt     >  2:4252552/1‑93 (MQ=255)

                                                                                        gcaacaCAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGt     <  1:787003/93‑1 (MQ=255)

                                                                                        gcaacaCAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGt     <  1:1866064/93‑1 (MQ=255)

                                                                                        gcaacaCAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGCGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGt     <  2:78238/93‑1 (MQ=255)

                                                                                          agcaCAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCt   >  2:3302917/1‑93 (MQ=255)

TCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTA  >  NC_000913/1949253‑1949436

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