BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,215,052	A→G	E122G (GAG→GGG)	aceA →	isocitrate lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,215,052	0	A	G	98.7%	459.4 / ‑3.4	148	E122G (GAG→GGG)	aceA	isocitrate lyase
Reads supporting (aligned to +/- strand): ref base A (2/0); new base G (62/84); total (64/84)
Fisher's exact test for biased strand distribution p-value = 1.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 154 total aligned reads displayed.

GATGGCAGGTAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGAGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCC  >  NC_000913/4214962‑4215141

gATGGCAGGTAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGGGc                                                                                         >  1:2299076/1‑93 (MQ=255)

  tGGCAGGTAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGGGCgg                                                                                       >  1:1307988/1‑93 (MQ=255)

   ggCAGGTAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGGGCGGa                                                                                      <  1:2403554/93‑1 (MQ=255)

    gCAGGTAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGGGCGGAt                                                                                     <  1:421200/93‑1 (MQ=255)

    gCAGGTAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGGGCGGAt                                                                                     <  1:421199/93‑1 (MQ=255)

    gCAGGTAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGGGCGGAt                                                                                     <  1:3914323/93‑1 (MQ=255)

       ggTAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGGGCGGATcaa                                                                                  >  2:426414/1‑93 (MQ=255)

         tAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGGGCGGATcaaca                                                                                <  2:2693713/93‑1 (MQ=255)

          tgcggcgGACGCTACCCTGGCGAACAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGGGCGGATcaacaa                                                                               <  2:3395150/92‑1 (MQ=255)

          aGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGGGCGGATcaacaa                                                                               >  2:2641265/1‑93 (MQ=255)

            cggcggACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGGGCGGATCAacaaca                                                                             <  1:2113513/93‑1 (MQ=255)

             ggcggACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGGGCGGATCAACAacac                                                                            >  2:3798452/1‑93 (MQ=255)

              gcggACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGGGCGGATCAACAACAcc                                                                           >  1:458656/1‑93 (MQ=255)

              gcggACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGGGCGGATCAACAACAcc                                                                           >  2:3888103/1‑93 (MQ=255)

               cggACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGGGCGGATCAACAACACCt                                                                          >  2:3260479/1‑93 (MQ=255)

               cggACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGGGCGGATCAACAACACCt                                                                          <  1:999675/93‑1 (MQ=255)

                 gACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGGGCGGATCAACAACACCTTc                                                                        <  2:948897/93‑1 (MQ=255)

                 gACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGGGCGGATCAACAACACCTTc                                                                        <  1:377840/93‑1 (MQ=255)

                 gACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGGGCGGATCAACAACACCTTc                                                                        >  2:1622694/1‑93 (MQ=255)

                   cGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGGGGTGGGGCGTATCAACAACACCTTCcg                                                                      <  2:4326602/93‑1 (MQ=255)

                    gCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGGGCGGATCAACAACACCTTCcgt                                                                     <  2:1865246/93‑1 (MQ=255)

                    gCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGGGCGGATCAACAACACCTTCcgt                                                                     <  1:4186071/93‑1 (MQ=255)

                     cTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGGGCGGATCAACAACACCTTCcgtc                                                                    <  1:1448375/93‑1 (MQ=255)

                      tAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGGGCGGATCAACAACACCTTCcgtcg                                                                   <  2:287911/93‑1 (MQ=255)

                      tAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGGGCGGATCAACAACACCTTCcgtcg                                                                   >  1:989368/1‑93 (MQ=255)

                      tAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGGGCGGATCAACAACACCTTCcgtcg                                                                   >  1:3469437/1‑93 (MQ=255)

                       aaCCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGGGCGGATCAACAACACCTTCcgtcgt                                                                  <  1:3554256/93‑1 (MQ=255)

                          cTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGGGCGGATCAACAACACCTTCCGTCGTGcc                                                               >  2:2546556/1‑93 (MQ=255)

                          cTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGGGCGGATCAACAACACCCTCCGTCGTGcc                                                               >  2:2777059/1‑93 (MQ=255)

                           tGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGGGCGGATCAACAACACCTTCCGTCGTGCCg                                                              >  2:4062449/1‑93 (MQ=255)

                           tGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGGGCGGATCAACAACACCTTCCGTCGTGCCg                                                              >  1:2625574/1‑93 (MQ=255)

                           tGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGGGCGGATCAACAACACCTTCCGTCGTGCCg                                                              >  2:3587553/1‑93 (MQ=255)

                           tGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGGGCGGATCAACAACACCTTCCGTCGTGCCg                                                              <  1:1593701/93‑1 (MQ=255)

                            ggcggcCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGGGCGGATCAACAACACCTTCCGTCGTGCCGa                                                             >  1:665852/1‑93 (MQ=255)

                               ggcCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGGGCGGATCAACAACACCTTCCGTCGTGCCGATCa                                                          >  1:2285120/1‑93 (MQ=255)

                                  cAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGAt                                                       <  1:1397136/93‑1 (MQ=255)

                                  cAGCATGTATCCAGATCATTCGCTCTATCCGGCAAACTCTGTTCCAGCTCTGGTAGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGAt                                                       <  1:747449/93‑1 (MQ=255)

                                   aGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATc                                                      >  2:3852130/1‑93 (MQ=255)

                                   aGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATc                                                      >  2:842636/1‑93 (MQ=255)

                                   aGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATc                                                      >  1:2041506/1‑93 (MQ=255)

                                   aGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATc                                                      >  1:1666653/1‑93 (MQ=255)

                                      aTGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCaa                                                   <  1:717016/93‑1 (MQ=255)

                                      aTGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCaa                                                   <  1:3501529/93‑1 (MQ=255)

                                       tGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAAt                                                  <  1:2594670/93‑1 (MQ=255)

                                       tGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAAt                                                  >  2:2178660/1‑93 (MQ=255)

                                       tGTATCCGGATCAGTCGCTATATCCGGCAAACTCGGTGCCAGCTGTGGTGGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAAt                                                  <  1:3880304/93‑1 (MQ=255)

                                        gTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATg                                                 <  2:2944728/93‑1 (MQ=255)

                                        gTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATg                                                 <  1:17574/93‑1 (MQ=255)

                                        gTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATg                                                 <  1:17573/93‑1 (MQ=255)

                                        gTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATg                                                 <  1:4096807/93‑1 (MQ=255)

                                         tATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGAGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATgg                                                >  1:3081747/1‑93 (MQ=255)

                                           tCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTc                                              <  2:659653/93‑1 (MQ=255)

                                           tCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTc                                              <  2:659652/93‑1 (MQ=255)

                                            ccGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTcc                                             <  2:3809397/93‑1 (MQ=255)

                                            ccGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTcc                                             <  2:1467416/93‑1 (MQ=255)

                                             cGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCcg                                            <  2:367542/93‑1 (MQ=255)

                                             cGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCcg                                            >  1:3604664/1‑93 (MQ=255)

                                              ggATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCcgc                                           <  2:1605775/93‑1 (MQ=255)

                                               gATCAGTCCATATATCCGGCAAACGCGGTGACAGCTGTGGTGGGGCGGACCACCAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCcgcg                                          <  2:757809/93‑1 (MQ=255)

                                               gATCAGTCACTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCcgcg                                          >  2:604782/1‑93 (MQ=255)

                                                   aGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCa                                      >  1:2093373/1‑93 (MQ=255)

                                                   aGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCa                                      >  1:321688/1‑93 (MQ=255)

                                                   aGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGCCCCCGGGCa                                      >  1:321686/1‑93 (MQ=255)

                                                      cGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTg                                   <  2:921307/93‑1 (MQ=255)

                                                       gCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGa                                  <  1:4682331/93‑1 (MQ=255)

                                                       gCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGa                                  <  1:528758/93‑1 (MQ=255)

                                                       gCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGa                                  <  2:4000365/93‑1 (MQ=255)

                                                       gCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGa                                  >  1:1460466/1‑93 (MQ=255)

                                                        ctctATCCGGCAAACTCGGTGCCAGCTGTGGTGGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAg                                 >  2:3103438/1‑93 (MQ=255)

                                                         tctATCCGGCAAACTCGGTGCCAGCTGTGGTGGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGc                                <  2:2383588/93‑1 (MQ=255)

                                                         tctATCCGGCAAACTCGGTGCCAGCTGTGGTGGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGGCCGCGGGCATTGAGc                                >  1:396994/1‑93 (MQ=255)

                                                          ctgtCCGGCAAACTCGGTGCCAGCTGTGGTGGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGcc                               <  2:492943/90‑1 (MQ=255)

                                                          ctATCCGGCAAACTCGGTGCCAGCTGTGGTGGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGcc                               <  1:2368069/93‑1 (MQ=255)

                                                          ctATCCGGCAAACTCGGTGCCAGCTGTGGTGGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGcc                               <  1:4242516/93‑1 (MQ=255)

                                                           tATCCGGCAAGCTCGGTGCCAGCTGTGGTGGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCg                              <  2:2608866/93‑1 (MQ=255)

                                                            aTCCGGCAAACTCGGTGCCAGCTGTGGTGGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCgg                             >  2:13951/1‑93 (MQ=255)

                                                            aTCCGGCAAACTCGGTGCCAGCTGTGGTGGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCgg                             >  1:1043774/1‑93 (MQ=255)

                                                             tCCGGCAAACTCGGTGCCAGCTGTGGTGGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCggg                            >  2:3279857/1‑93 (MQ=255)

                                                               cGGCAAACTCGGTGCCAGCTGTGGTGGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCg                          >  1:4011234/1‑93 (MQ=255)

                                                               cGGCAAACTCGGTGCCAGCTGTGGTGGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCg                          <  1:1062788/93‑1 (MQ=255)

                                                                ggCAAACTCGGTGCCAGCTGTGGTGGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAAAGGTCCGCGGGCATTGAGCCGGGCGa                         <  2:1274822/93‑1 (MQ=255)

                                                                 gCAAACTCGGTGCCAGCTGTGGTGGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGAt                        <  1:3186749/93‑1 (MQ=255)

                                                                 gCAAACTCGGTGCCAGCTGTGGTGGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGAt                        >  1:2252369/1‑93 (MQ=255)

                                                                   aaaCTCGGTGCCAGCTGTGGTTGGGCGGATCAACAACCCCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATcc                      <  2:1170919/93‑1 (MQ=255)

                                                                    aaCTCGGTGCCAGCTGTGGTGGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCcg                     >  1:4519479/1‑93 (MQ=255)

                                                                     aCTCGGTGCCAGCTGTGGTGGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCcgc                    >  1:429375/1‑93 (MQ=255)

                                                                     aCTCGGTGCCAGCTGTGGTGGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCcgc                    <  1:745297/93‑1 (MQ=255)

                                                                        cGGTGCCAGCTGTGGTGGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCt                 >  2:2186504/1‑93 (MQ=255)

                                                                          gTGCCAGCTGTGGTGGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTAt               <  1:2953891/93‑1 (MQ=255)

                                                                          gTGCCAGCTGTGGTGGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTAt               >  2:3996835/1‑93 (MQ=255)

                                                                           tGCCAGCTGTGGTGGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATg              <  1:1670133/93‑1 (MQ=255)

                                                                             ccAGCTGTGGTGGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTc            <  2:3490010/93‑1 (MQ=255)

                                                                               aGCTGTGGTGGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGa          >  1:4438247/1‑93 (MQ=255)

                                                                                   gtggtgGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTAc      <  2:2367084/93‑1 (MQ=255)

                                                                                   gtgctgGGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGTCGGGCAATCCGCGCTATGTCGATTAc      <  2:2538685/93‑1 (MQ=255)

                                                                                    tggtggGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACt     <  1:659569/93‑1 (MQ=255)

                                                                                     ggtggGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACtt    <  1:1490133/93‑1 (MQ=255)

                                                                                     ggtggGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACtt    >  2:923948/1‑93 (MQ=255)

                                                                                       tggGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTcc  <  1:726114/93‑1 (MQ=255)

                                                                                       tggGGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTcc  >  1:1766050/1‑93 (MQ=255)

GATGGCAGGTAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGAGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCC  >  NC_000913/4214962‑4215141

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