BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,056,815	T→C	M352V (ATG→GTG)	yeeO ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,056,815	0	T	C	86.8%	507.3 / 68.2	212	M352V (ATG→GTG)	yeeO	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (11/17); new base C (77/107); total (88/124)
Fisher's exact test for biased strand distribution p-value = 8.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.43e-01

Only 100 of 218 total aligned reads displayed.

GTAGAAGCAGAGCCGAGCGCACTTCCGGGTAAGTTGATAAGAGCCGCAATTGAAAACGCGATAAAATTTCCGGCAATAACACTGGTCCCCATCCCGGCAACGAACATTTGGGTTAATAACCGACCACTGGTAAATAACACTGATTCGACACTCGCGGGAATACCAATCCCCATGACTTCCCAGA  >  NC_000913/2056724‑2056907

 tAGAAGCAGAGCCGAGCGCACTTCCGGGTAAGTTGATAAGAGCCGCAATTGAAAACGCGATAAAATTTCCGGCAATAACACTGGTCCCCAccc                                                                                            <  2:6027629/93‑1 (MQ=255)

 tAGAAGCAGAGCCGAGCGCACTTCCGGGTAAGTTGATAAGAGCCGCAATTGAAAACGCGATAAAATTTCCGGCAATAACACTGGTCCCCAccc                                                                                            <  1:1894153/93‑1 (MQ=255)

  aGAAGCAGAGCCGAGCGCACTTCCGGGTAAGTTGATAAGAGCCGCAATTGAAAACGCGATAAAATTTCCGGCAATAACACTGGTCCCCAcccc                                                                                           <  1:1055307/93‑1 (MQ=255)

   gAAGCAGAGCCGAGCGCACTTCCGGGTAAGTTGATAAGAGCCGCAATTGAAAACGCGATAAAATTTCCGGCAATAACACTGGTCCCCATCCCg                                                                                          <  2:382603/93‑1 (MQ=255)

      gCAGAGCCGAGCGCACTTCCGGGTAAGTTGATAAGAGCCGCAATTGAAAACGCGATAAAATTTCCGGCAATAACACTGGTCCCCACCCCGGCa                                                                                       >  2:7457212/1‑93 (MQ=255)

      gCAGAGCCGAGCGCACTTCCGGGTAAGTTGATAAGAGCCGCAATTGAAAACGCGATAAAATTTCCGGCAATAACACTGGTCCCCACCCCGGCa                                                                                       >  1:7342981/1‑93 (MQ=255)

       cAGAGCCGAGCGCACTTCCGGGTAAGTTGATAAGAGCCGCAATTGAAAACGCGATAAAATTTCCGGCAATAACACTGGTCCCCACCCCGGCaa                                                                                      >  2:1735316/1‑93 (MQ=255)

         gagCCGAGCGCACTTCCGGGTAAGTTGATAAGAGCCGCAATTGAAAACGCGATAAAATTTCCGGCAATAACACTGGTCCCCACCCCGGCAACg                                                                                    >  1:586205/1‑93 (MQ=255)

         gagCCGAGCGCACTTCCGGGTAAGTTGATAAGAGCCGCAATTGAAAACGCGATAAAATTTCCGGCAATAACACTGGTCCCCACCCCGGCAACg                                                                                    >  2:673618/1‑93 (MQ=255)

           gCCGAGCGCACTTCCGGGTAAGTTGATAAGAGCCGCAATTGAAAACGCGATAAAATTTCCGGCAATAACACTGGTCCCCATCCCGGCAACGaa                                                                                  >  1:2169934/1‑93 (MQ=255)

           gCCGAGCGCACTTCCGGGTAAGTTGATAAGAGCCGCAATTGAAAACGCGATAAAATTTCCGGCAATAACACTGGTCCCCACCCCGGCAACGaa                                                                                  <  1:7424575/93‑1 (MQ=255)

             cGAGCGCACTTCCGGGTAAGTTGATAAGAGCCGCAATTGAAAACGCGATAAAATTTCCGGCAATAACACTGGTCCCCACCCCGGCAACGAACa                                                                                >  1:5155484/1‑93 (MQ=255)

               aGCGCACTTCCGGGTAAGTTGATAAGAGCCGCAATTGAAAACGCGATAAAATTTCCGGCAATAACACTGGTCCCCACCCCGGCAACGAACAtt                                                                              <  1:5905085/93‑1 (MQ=255)

               aGCGCACTTCCGGGTAAGTTGATAAGAGCCGCAATTGAAAACGCGATAAAATTTCCGGCAATAACACTGGTCCCCACCCCGGCAACGAACAtt                                                                              <  1:837620/93‑1 (MQ=255)

               aGCGCACTTCCGGGTAAGTTGATAAGAGCCGCAATTGAAAACGCGATAAAATTTCCGGCAATAACACTGGTCCCCACCCCGGCAACGAACAtt                                                                              <  2:1032445/93‑1 (MQ=255)

                gcgcACTTCCGGGTAAGTTGATAAGAGCCGCAATTGAAAACGCGATAAAATTTCCGGCAATAACACTGGTCCCCACCCCGGCAACGAACAttt                                                                             >  2:4388479/1‑93 (MQ=255)

                gcgcACTTCCGGGTAAGTTGATAAGAGCCGCAATTGAAAACGCGATAAAATTTCCGGCAATAACACTGGTCCCCACCCCGGCAACGAACAttt                                                                             <  2:3148040/93‑1 (MQ=255)

                 cgcACTTCCGTGTAAGTTGATAAGAGCCTCAATTGAAAACGCGATAAAATTTCCGGCAATAACACCGGTCCCCACCCCGGCAACGAACATTTg                                                                            <  1:2075989/93‑1 (MQ=255)

                 cgcACTTCCGGGTAAGTTGATAAGAGCCGCAATTGAAAACGCGATAAAATTTCCGGCAATAACACTGGTCCCCACCCCGGCAACGAACATTTg                                                                            <  1:3550763/93‑1 (MQ=255)

                  gcACTTCCGGGTAAGTTGATAAGAGCCGCAATTGAAAACGCGATAAAATTTCCGGCAATAACACTGGTCCCCATCCCGGCAACGAACATTTgg                                                                           <  1:475900/93‑1 (MQ=255)

                  gcACTTCCGGGTAAGTTGATAAGAGCCGCAATTGAAAACGCGATAAAATTTCCGGCAATAACACTGGTCCCCACCCCGGCAACGAACATTTgg                                                                           <  2:6380095/93‑1 (MQ=255)

                   cACTTCCGGGTAAGTTGATAAGAGCCGCAATTGAAAACGCGATAAAATTTCCGGCAATAACACTGGTCCCCATCCCGGCAACGAACATTTggg                                                                          <  1:2580316/93‑1 (MQ=255)

                   cACTTCCGGGTAAGTTGATAAGAGCCGCAATTGAAAACGCGATAAAATTTCCGGCAATAACACTGGTCCCCATCCCGGCAACGAACATTTggg                                                                          >  1:6688346/1‑93 (MQ=255)

                   cACTTCCGGGTAAGTTGATAAGAGCCGCAATTGAAAACGCGATAAAATTTCCGGCAATAACACTGGTCCCCACCCCGGCAACGAACATTTggg                                                                          >  2:3374302/1‑93 (MQ=255)

                    aCTTCCGGGTAAGTTGATAAGAGCCGCAATTGAAAACGCGATAAAATTTCCGGCAATAACACTGGTCCCCACCCCGGCAACGAACATTTGGGt                                                                         <  2:2977563/93‑1 (MQ=255)

                      ttCCGGGTAAGTTGATAAGAGCCGCAATTGAAAACGCGATAAAATTTCCGGCAATAACACTGGTCCCCATCCCGGCAACGAACATTTGGGTta                                                                       >  1:1722987/1‑93 (MQ=255)

                      ttCCGGGTAAGTTGATAAGAGCCGCAATTGAAAACGCGATAAAATTTCCGGCAATAACACTGGTCCCCACCCCGGCAACGAACATTTGGGTta                                                                       >  1:2003118/1‑93 (MQ=255)

                      ttCCGGGTAAGTTGATAAGAGCCGCAATTGAAAACGCGATAAAATTTCCGGCAATAACACTGGTCCCCACCCCGGCAACGAACATTTGGGTta                                                                       >  2:167381/1‑93 (MQ=255)

                        ccGGGTAAGTTGATAAGAGCCGCAATTGAAAACGCGATAAAATTTCCGGCAATAACACTGGTCCCCACCCCGGCAACGAACATTTGGGTtaat                                                                     <  2:5949732/93‑1 (MQ=255)

                          gggTAAGTTGATAAGAGCCGCAATTGAAAACGCGATAAAATTTCCGGCAATAACACTGGTCCCCACCCCGGCAACGAACATTTGGGTtaatga                                                                   <  1:3986702/93‑3 (MQ=255)

                          gggTAAGTTGATAAGAGCCGCAATTGAAAACGCGATAAAATTTCCGGCAATAACACTGGTCCCCACCCCGGCAACGAACATTTGGGTtaataa                                                                   >  2:6486544/1‑93 (MQ=255)

                          gggTAAGTTGATAAGAGCCGCAATTGAAAACGCGATAAAATTTCCGGCAATAACACTGGTCCCCACCCCGGCAACGAACATTTGGGTtaataa                                                                   >  2:4611482/1‑93 (MQ=255)

                          gggTAAGTTGATAAGAGCCGCAATTGAAAACGCGATAAAATTTCCGGCAATAACACTGGTCCCCACCCCGGCAACGAACATTTGGGTtaataa                                                                   <  1:886265/93‑1 (MQ=255)

                           ggTAAGTTGATAAGAGCCGCAATTGAAAACGCGATAAAATTTCCGGCAATAACACTGGTCCCCACCCCGGCAACGAACATTTGGGTTAATAAc                                                                  >  1:3488156/1‑93 (MQ=255)

                           ggTAAGTTGATAAGAGCCGCAATTGAAAACGCGATAAAATTTCCGGCAATAACACTGGTCCCCACCCCGGCAACGAACATTTGGGTTAATAAc                                                                  <  1:5046255/93‑1 (MQ=255)

                           ggTAAGTTGATAAGAGCCGCAATTGAAAACGCGATAAAATTTCCGGCAATAACACTGGTCCCCACCCCGGCAACGAACATTTGGGTTAATAAc                                                                  >  1:1084616/1‑93 (MQ=255)

                           ggTAAGTTGATAAGAGCCGCAATTGAAAACGCGATAAAATTTCCGGCAATAACACTGGTCCCCACCCCGGCAACGAACATTTGGGTTAATAAc                                                                  >  1:49908/1‑93 (MQ=255)

                            gTAAGTTGATAAGAGCCGCAATTGAAAACGCGATAAAATTTCCGGCAATAACACTGGTCCCCATCCCGGCAACGAACATTTGGGTTAATAAcc                                                                 <  2:3956932/93‑1 (MQ=255)

                             tAAGTTGATAAGAGCCGCAATTGAAAACGCGATAAAATTTCCGGCAATAACACTGGTCCCCACCCCGGCAACGAACATTTGGGTTAATAACCg                                                                >  1:854002/1‑93 (MQ=255)

                             tAAGTTGATAAGAGCCGCAATTGAAAACGCGATAAAATTTCCGGCAATAACACTGGTCCCCACCCCGGCAACGAACATTTGGGTTAATAACCg                                                                <  1:2438214/93‑1 (MQ=255)

                              aaGTTGATAAGAGCCGCAATTGAAAACGCGATAAAATTTCCGGCAATAACACTGGTCCCCACCCCGGCAACGAACATTTGGGTTAATAACCGa                                                               <  2:6288427/93‑1 (MQ=255)

                               aGTTGATAAGAGCCGCAATTGAAAGCGCGATAAAATTTCCGGCAATAACACTGGTCCCCACCCCGGCAACGAACATTTGGGTTAATAACCGAc                                                              <  2:4785531/93‑1 (MQ=255)

                               aGTTGATAAGAGCCGCAATTGAAAACGCGATAAAATTTCCGGCAATAACACTGGTCCCCACCCCGGCAACGAACATTTGGGTTAATAACCGAc                                                              <  2:2886412/93‑1 (MQ=255)

                               aGTTGATAAGAGCCGCAATTGAAAACGCGATAAAATTTCCGGCAATAACACTGGTCCCCACCCCGGCAACGAACATTTGGGTTAATAACCGAc                                                              <  1:5201065/93‑1 (MQ=255)

                                gTTGATAAGAGCCGCAATTGAAAACGCGATAAAATTTCCGGCAATAACACTGGTCCCCACCCCGGCAACGAACATTTGGGTTAATAACCGAcc                                                             >  1:1151497/1‑93 (MQ=255)

                                 ttGATAAGAGCCGCAATTGAAAACGCGATAAAATTTCCGGCAATAACACTGGTCCCCACCCCGGCAACGAACATTTGGGTTAATAACCGACCa                                                            <  2:6949521/93‑1 (MQ=255)

                                  tGATAAGAGCCGCAATTGAAAACGCGATAAAATTTCCGGCAATAACACTGGTCCCCACCCCGGCAACGAACATTTGGGTTAATAACCGACCAc                                                           <  2:4520823/93‑1 (MQ=255)

                                    aTAAGAGCCGCAATTGAAAACGCGATAAAATTTCCGGCAATAACACTGGTCCCCACCCCGGCAACGAACATTTGGGTTAATAACCGACCACTg                                                         >  1:3826042/1‑93 (MQ=255)

                                    aTAAGAGCCGCAATTGAAAACGCGATAAAATTTCCGGCAATAACACTGGTCCCCACCCCGGCAACGAACATTTGGGTTAATAACCGACCACTg                                                         <  1:519067/93‑1 (MQ=255)

                                     tAAGAGCCGCAATTGAAAACGCGATAAAATTTCCGGCAATAACACTGGTCCCCATCCCGGCAACGAACATTTGGGTTAATAACCGACCACTgg                                                        <  1:6482996/93‑1 (MQ=255)

                                     tAAGAGCCGCAATTGAAAACGCGATAAAATTTCCGGCAATAACACTGGTCCCCATCCCGGCAACGAACATTTGGGTTAATAACCGACCACTgg                                                        <  2:7179112/93‑1 (MQ=255)

                                     tAAGAGCCGCAATTGAAAACGCGATAAAATTTCCGGCAATAACACTGGTCCCCACCCCGGCAACGAACATTTGGGTTAATAACCGACCACTgg                                                        >  2:7054294/1‑93 (MQ=255)

                                       agagCCGCAATTGAAAACGCGATAAAATTTCCGGCAATAACACTGGTCCCCACCCCGGCAACGAACATTTGGGTTAATAACCGACCACTGGTa                                                      >  2:4041769/1‑93 (MQ=255)

                                         agCCGCAATTGAAAACGCGATAAAATTTCCGGCAATAACACTGGTCCCCACCCCGGCAACGAACATTTGGGTTAATAACCGACCACTGGTaaa                                                    <  1:5766306/93‑1 (MQ=255)

                                         agCCGCAATTGAAAACGCGATAAAATTTCCGGCAATAACACTGGTCCCCACCCCGGCAACGAACATTTGGGTTAATAACCGACCACTGGTaaa                                                    <  2:5862805/93‑1 (MQ=255)

                                          gCCGCAATTGAAAACGCGATAAAATTTCCGGCAATAACACTGGTCCCCACCCCGGCAACGAACATTTGGGTTAATAACCGACCACTGGTAAAt                                                   <  2:1976713/93‑1 (MQ=255)

                                           ccGCAATTGAAAACGCGATAAAATTTCCGGCAATAACACTGGTCCCCACCCCGGCAACGAACATTTGGGTTAATAACCGACCACTGGTAAATa                                                  >  2:5013400/1‑93 (MQ=255)

                                            cGCAATTGAAAACGCGATAAAATTTCCGGCAATAACACTGGTCCCCACCCCGGCAACGAACATTTGGGTTAATAACCGACCACTGGTAAATaa                                                 <  2:1278551/93‑1 (MQ=255)

                                               aaTTGAAAACGCGATAAAATTTCCGGCAATAACACTGGTCCCCACCCCGGCAACGAACATTTGGGTTAATAACCGACCACTGGTAAATAacac                                              >  1:5996386/1‑93 (MQ=255)

                                                aTTGAAAACGCGATAAAATTTCCGGCAATAACACTGGTCCCCACCCCGGCAACGAACATTTGGGTTAATAACCGACCACTGGTAAATAACACt                                             >  2:2167490/1‑93 (MQ=255)

                                                aTTGAAAACGCGATAAAATTTCCGGCAATAACACTGGTCCCCACCCCGGCAACGAACATTTGGGTTAATAACCGACCACTGGTAAATAACACt                                             <  1:1214017/93‑1 (MQ=255)

                                                 ttGAAAACGCGATAAAATTTCCGGCAATAACACTGGTCCCCATCCCGGCAACGAACATTTGGGTTAATAACCGACCACTGGTAAATAACACTg                                            >  2:4968074/1‑93 (MQ=255)

                                                 ttGAAAACGCGATAAAATTTCCGGCAATAACACTGGTCCCCACCCCGGCAACGAACATTTGGGTTAATAACCGACCACTGGTAAATAACACTg                                            >  1:6425596/1‑93 (MQ=255)

                                                 ttGAAAACGCGATAAAATTTCCGGCAATAACACTGGTCCCCACCCCGGCAACGAACATTTGGGTTAATAACCGACCACTGGTAAATAACACTg                                            >  2:6505750/1‑93 (MQ=255)

                                                 ttGAAAACGCGATAAAATTTCCGGCAATAACACTGGTCCCCACCCCGGCAACGAACATTTGGGTTAATAACCGACCACTGGTAAATAACACTg                                            >  2:1075569/1‑93 (MQ=255)

                                                     aaaCGCGATAAAATTTCCGGCAATAACACTGGTCCCCACCCCGGCAACGAACATTTGGGTTAATAACCGACCACTGGTAAATAACACTGATTc                                        <  1:2539369/93‑1 (MQ=255)

                                                      aaCGCGATAAAATTTCCGGCAATAACACTGGTCCCCACCCCGGCAACGAACATTTGGGTTAATAACCGACCACTGGTAAATAACACTGATTCg                                       <  1:2937621/93‑1 (MQ=255)

                                                        cgcgATAAAATTTCCGGCAATAACACTGGTCCCCACCCCGGCAACGAACATTTGGGTTAATAACCGACCACTGGTAAATAACACTGATTCGac                                     >  1:6998304/1‑93 (MQ=255)

                                                           gATAAAATTTCCGGCAATAACACTGGTCCCCACCCCGGCAACGAACATTTGGGTTAATAACCGACCCCTGGTAAATAACACTGATTCGACACt                                  >  1:2881117/1‑93 (MQ=255)

                                                           gATAAAATTTCCGGCAATAACACTGGTCCCCACCCCGGCAACGAACATTTGGGTTAATAACCGACCACTGGTAAATAACACTGATTCGACACt                                  <  2:5414669/93‑1 (MQ=255)

                                                           gATAAAATTTCCGGCAATAACACTGGTCCCCACCCCGGCAACGAACATTTGGGTTAATAACCGACCACTGGTAAATAACACTGATTCGACACt                                  >  1:4947027/1‑93 (MQ=255)

                                                           gATAAAATTTCCGGCAATAACACTGGTCCCCACCCCGGCAACGAACATTTGGGTTAATAACCGACCACTGGTAAATAACACTGATTCGACACt                                  >  1:4252732/1‑93 (MQ=255)

                                                           gATAAAATTTCCGGCAATAACACTGGTCCCCACCCCGGCAACGAACATTTGGGTTAATAACCGACCACTGGTAAATAACACTGATTCGACACt                                  >  2:2574681/1‑93 (MQ=255)

                                                            aTAAAATTTCCGGCAATAACACTGGTCCCCATCCCGGCAACGAACATTTGGGTTAATAACCGACCACTGGTAAATAACACTGATTCGACACTc                                 <  1:6004582/93‑1 (MQ=255)

                                                            aTAAAATTTCCGGCAATAACACTGGTCCCCACCCCGGCAACGAACATTTGGGTTAATAACCGACCACTGGTAAATAACACTGATTCGACACTc                                 >  1:148565/1‑93 (MQ=255)

                                                             tAAAATTTCCGGCAATAACACTGGTCCCCACCCCGGCAACGAACATTTGGGTTAATAACCGACCACTGGTAAATAACACTGATTCGACACTcg                                >  1:1672448/1‑93 (MQ=255)

                                                              aaaaTTTCCGGCAATAACACTGGTCCCCACCCCGGCAACGAACATTTGGGTTAATAACCGACCACTGGTAAATAACACTGATTCGACACTcgc                               <  2:309104/93‑1 (MQ=255)

                                                              aaaaTTTCCGGCAATAACACTGGTCCCCACCCCGGCAACGAACATTTGGGTTAATAACCGACCACTGGTAAATAACACTGATTCGACACTcgc                               <  2:6255243/93‑1 (MQ=255)

                                                                aaTTTCCGGCAATAACACTGGTCCCCACCCCGGCAACGAACATTTGGGTTAATAACCGACCACTGGTAAATAACACTGATTCGACACTCGCgg                             <  1:4737887/93‑1 (MQ=255)

                                                                aaTTTCCGGCAATAACACTGGTCCCCACCCCGGCAACGAACATTTGGGTTAATAACCGACCACTGGTAAATAACACTGATTCGACACTCGCgg                             <  2:2124018/93‑1 (MQ=255)

                                                                   ttCCGGCAATAACACTGGTCCCCACCCCGGCAACGAACATTTGGGTTAATAACCGACCACTGGTAAATAACACTGATTCGACACTCGCGGGaa                          <  2:4288877/93‑1 (MQ=255)

                                                                     ccGGCAATAACACTGGTCCCCACCCCGGCAACGAACATTTGGGTTAATAACCGACCACTGGTAAATAACACTGATTCGACACTCGCGGGAATa                        <  2:5517729/93‑1 (MQ=255)

                                                                       ggCAATAACACTGGTCCCCACCCCGGCAACGAACATTTGGGTTAATAACCGACCACTGGTAAATAACACTGATTCGACACTCGCGGGAATAcc                      <  2:458803/93‑1 (MQ=255)

                                                                       ggCAATAACACTGGTCCCCACCCCGGCAACGAACATTTGGGTTAATAACCGACCACTGGTAAATAACACTGATTCGACACTCGCGGGAATAcc                      <  2:7371650/93‑1 (MQ=255)

                                                                       ggCAATAACACTGGTCCCCACCCCGGCAACGAACATTTGGGTTAATAACCGACCACTGGTAAATAACACTGATTCGACACTCGCGGGAATAcc                      <  2:612398/93‑1 (MQ=255)

                                                                        gCAATAACACTGGTCCCCACCCCGGCAACGAACATTTGGGTTAATAACCGACCACTGGTAAATAACACTGATTCGACACTCGCGGGAATACCa                     >  2:4265204/1‑93 (MQ=255)

                                                                         cAATAACACTGGTCCCCACCCCGGCAACGAACATTTGGGTTAATAACCGACCACTGGTAAATAACACTGATTCGACACTCGCGGGAATACCaa                    <  1:6910342/93‑1 (MQ=255)

                                                                           aTAACACTGGTCCCCACCCCGGCAACGAACATTTGGGTTAATAACCGACCACTGGTAAATAACACTGATTCGACACTCGCGGGAATACCAATc                  <  2:2506097/93‑1 (MQ=255)

                                                                            tAACACTGGTCCCCATCCCGGCAACGAACATTTGGGTTAATAACCGACCACTGGTAAATAACACTGATTCGACACTCGCGGGAATACCAATcc                 >  2:1204217/1‑93 (MQ=255)

                                                                            tAACACTGGTCCCCACCCCGGCAACGAACATTTGGGTTAATAACCGACCACTGGTAAATAAGACTGATTCGACACTCGCGGGAATACCAATcc                 >  1:3905708/1‑93 (MQ=255)

                                                                             aaCACTGGTCCCCACCCCGGCAACGAACATTTGGGTTAATAACCGACCACTGGTAAATAACACTGATTCGACACTCGCGGGAATACCAATccc                >  1:3207764/1‑93 (MQ=255)

                                                                              acacTGGTCCCCACCCCGGCAACGAACATTTGGGTTAATAACCGACCACTGGTAAATAACACTGATTCGACACTCGCGGGAATACCAATcccc               >  1:6517284/1‑93 (MQ=255)

                                                                               cacTGGTCCCCATCCCGGCAACGAACATTTGGGTTAATAACCGACCACTGGTAAATAACACTGATTCGACACTCGCGGGAATACCAATCCCCa              <  1:345325/93‑1 (MQ=255)

                                                                               cacTGGTCCCCACCCCGGCAACGAACATTTGGGTTAATAACCGACCACTGGTAAATAACACTGATTCGACACTCGCGGGAATACCAATCCCCa              <  2:2027211/93‑1 (MQ=255)

                                                                                 cTGGTCCCCACCCCGGCAACGAACATTTGGGTTAATAACCGACCACTGGTAAATAACACTGATTCGACACTCGCGGGAATACCAATCCCCATg            >  2:3656094/1‑93 (MQ=255)

                                                                                    gTCCCCACCCCGGCAACGAACATTTGGGTTAATAACCGACCACTGGTAAATAACACTGATTCGACACTCGCGGGAATACCAATCCCCATGACt         <  1:7172257/93‑1 (MQ=255)

                                                                                    gTCCCCACCCCGGCAACGAACATTTGGGTTAATAACCGACCACTGGTAAATAACACTGATTCGACACTCGCGGGAATACCAATCCCCATGACt         <  1:339207/93‑1 (MQ=255)

                                                                                    gTCCCCACCCCGGCAACGAACATTTGGGTTAATAACCGACCACTGGTAAATAACACTGATTCGACACTCGCGGGAATACCAATCCCCATGACt         <  2:1390405/93‑1 (MQ=255)

                                                                                     tCCCCACCCCGGCAACGAACATTTGGGTTAATAACCGACCACTGGTAAATAACACTGATTCGACACTCGCGGGAATACCAATCCCCATGACtt        <  2:286184/93‑1 (MQ=255)

                                                                                           tCCCGGCAACGAACATTTGGGTTAATAACCGACCACTGGTAAATAACACTGATTCGACACTCGCGGGAATACCAATCCCCATGACTTCCCAGa  >  1:3759151/1‑93 (MQ=255)

GTAGAAGCAGAGCCGAGCGCACTTCCGGGTAAGTTGATAAGAGCCGCAATTGAAAACGCGATAAAATTTCCGGCAATAACACTGGTCCCCATCCCGGCAACGAACATTTGGGTTAATAACCGACCACTGGTAAATAACACTGATTCGACACTCGCGGGAATACCAATCCCCATGACTTCCCAGA  >  NC_000913/2056724‑2056907

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