BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,026,252	2 bp→GC	coding (110‑111/1164 nt)	fadA ←	thiolase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,026,252	0	C	G	100.0%	1055.2 / NA	331	T37S (ACG→AGC)	fadA	thiolase I
Reads supporting (aligned to +/- strand): ref base C (0/0); new base G (167/164); total (167/164)
*	NC_000913	4,026,253	0	G	C	99.7%	1061.5 / ‑4.9	333	T37S (ACG→AGC)	fadA	thiolase I
Reads supporting (aligned to +/- strand): ref base G (0/0); major base C (167/165); minor base A (0/1); total (167/166)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.69e-01

Only 100 of 347 total aligned reads displayed.

AAACCCTGCTCCAGCGTCTGCTGCACACAACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGCGTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCATCGGGGTGCGAATTGCA  >  NC_000913/4026166‑4026339

||

 aaCCCTGCTCCAGCGTCTGCTGCACACAACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGGTACGCGCCAGCAGGCTACGCa                                                                                  >  2:9583337/1‑93 (MQ=255)

     cTGCTCCAGCGTCTGCTGCACACAACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTaa                                                                              >  1:293773/1‑93 (MQ=255)

      tGCTCCAGCGTCTGCTGCACACAACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTaaa                                                                             <  1:6092805/93‑1 (MQ=255)

      tGCTCCAGCGTCTGCTGCACACAACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTaaa                                                                             <  1:2667086/93‑1 (MQ=255)

      tGCTCCAGCGTCTGCTCCACACACCCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTaaa                                                                             <  2:5374013/93‑1 (MQ=255)

       gCTCCAGCGTCTGCTGCACACAACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCCGCGGCGGGCTACGCGCCAGCAGGCTACGCATTAAAt                                                                            >  2:9765974/1‑93 (MQ=255)

         tCCAGCGTCTGCTGCACACACCCCCAGTAAATAGCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGa                                                                          <  2:5663000/93‑1 (MQ=255)

         tCCAGCGTCTGCTGCACACAACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGa                                                                          >  2:9524065/1‑93 (MQ=255)

         tCCAGCGTCTGCTGCACACAACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGa                                                                          <  1:5583053/93‑1 (MQ=255)

          ccAGCGTCTGCTGCACACAACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAg                                                                         <  2:7524755/93‑1 (MQ=255)

          ccAGCGTCTGCTGCACACAACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAg                                                                         <  2:2605856/93‑1 (MQ=255)

          ccAGCGTCTGCTGCACACAACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAg                                                                         <  2:6209041/93‑1 (MQ=255)

          ccAGCGTCTGCTGCACACAACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAg                                                                         <  1:4204941/93‑1 (MQ=255)

          ccAGCGTCTGCTGCACACAACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAg                                                                         >  1:3878798/1‑93 (MQ=255)

          ccAGCGTCTGCTGCACACAACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGGTACGCGCCAGCAGGCTACGCATTACATgcg                                                                         >  1:3824519/1‑91 (MQ=255)

            aGCGTCTGCTGCACACACCCCCAGTAAATATCGACGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCg                                                                       <  1:8966973/93‑1 (MQ=255)

            aGCGTCTGCTGCACACAACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCg                                                                       >  2:4059298/1‑93 (MQ=255)

            aGCGTCTGCTGCACACAACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCCGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCg                                                                       >  1:7596710/1‑93 (MQ=255)

             gCGTCTGCTGCACACAACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCgg                                                                      >  2:8399675/1‑93 (MQ=255)

              cGTCTGCTGCACACAACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGga                                                                     <  2:2129348/93‑1 (MQ=255)

              cGTCTGCTGCACACAACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGga                                                                     >  2:3608388/1‑93 (MQ=255)

              cGTCTGCCGCACACACCCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGga                                                                     <  1:308886/93‑1 (MQ=255)

               gTCTGCTGCACACAACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGgag                                                                    <  1:4350984/93‑1 (MQ=255)

               gTCTGCTGCACACAACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGgag                                                                    <  1:7429236/93‑1 (MQ=255)

                tCTGCTGCACACAACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGgaga                                                                   <  1:344767/93‑1 (MQ=255)

                tCTGCTGCACACAACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGCTACGCGCCAGCAGGCTACGACTTAAATGAGCggcac                                                                   >  1:5289928/1‑90 (MQ=255)

                  tgctgcACACAACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGgagaga                                                                 <  1:8418175/93‑1 (MQ=255)

                  cgctgcACACCACCCCAGTAAACACCGTCGCGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGgagaga                                                                 <  1:5334229/92‑1 (MQ=255)

                   gctgcACACAACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGAt                                                                <  2:3555747/93‑1 (MQ=255)

                   gctgcACACAACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGAt                                                                >  1:5774641/1‑93 (MQ=255)

                   gcagcACACAACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGAt                                                                <  1:4837588/90‑1 (MQ=255)

                     tgcACACAACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGAtct                                                              >  2:6873218/1‑93 (MQ=255)

                     tgcACACAACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGAtct                                                              >  2:5560055/1‑93 (MQ=255)

                       cacacaACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGAtcttc                                                            <  2:13717/93‑1 (MQ=255)

                       cacacaACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGAtcttc                                                            <  1:8027453/93‑1 (MQ=255)

                       cacacaACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGAtcttc                                                            <  1:581755/93‑1 (MQ=255)

                       cacacaACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGAtcttc                                                            <  2:7920127/93‑1 (MQ=255)

                        acacaACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGAtcttct                                                           <  1:1523412/93‑1 (MQ=255)

                         cacaACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTg                                                          >  1:4470253/1‑93 (MQ=255)

                         cacaACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGGTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTg                                                          >  2:7744005/1‑93 (MQ=255)

                          acaACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTgc                                                         <  1:7246585/93‑1 (MQ=255)

                           caACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTgca                                                        <  1:10022542/93‑1 (MQ=255)

                           caACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTgca                                                        <  1:3364411/93‑1 (MQ=255)

                             aCCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAACTGAGCGGAGAGATCTtgtgattg                                                      >  2:8478655/1‑89 (MQ=255)

                             aCCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTgcacg                                                      >  2:2320211/1‑93 (MQ=255)

                             aCCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTgcacg                                                      >  2:7282476/1‑93 (MQ=255)

                              ccccAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTgcacgc                                                     >  1:4945314/1‑93 (MQ=255)

                                ccAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTgcacgcac                                                   <  1:9376984/93‑1 (MQ=255)

                                 cAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGcacgcacg                                                  <  2:4574011/93‑1 (MQ=255)

                                  aGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGt                                                 >  1:4041689/1‑93 (MQ=255)

                                   gTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGtt                                                <  1:1407759/93‑1 (MQ=255)

                                   gTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGtt                                                <  1:8656830/93‑1 (MQ=255)

                                   gTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGtt                                                >  1:1982206/1‑93 (MQ=255)

                                   gTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGtt                                                >  2:4964876/1‑93 (MQ=255)

                                      aaTATCGTCGACGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACg                                             >  1:8691970/1‑93 (MQ=255)

                                        tatCGTCGAGGGCCGGCTCTTCCAGCGCCGGGGTACGCGCCAGCAGGCTTCGCATTAAATGAGCGAAGGGATCATTGTCATGCACGCTACgta                                           >  2:8964135/1‑91 (MQ=255)

                                        tatCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGaa                                           >  2:1783998/1‑93 (MQ=255)

                                         atCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGaaa                                          >  2:1490816/1‑93 (MQ=255)

                                          tcgtcgAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGaaaa                                         <  1:6907833/93‑1 (MQ=255)

                                          tcgtcgAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGaaaa                                         <  1:4147943/93‑1 (MQ=255)

                                           cgtcgAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAg                                        >  2:8129580/1‑93 (MQ=255)

                                            gtcgAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGc                                       <  1:2211013/93‑1 (MQ=255)

                                             tcgAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCa                                      >  2:910715/1‑93 (MQ=255)

                                              cgAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCAc                                     <  1:6254305/93‑1 (MQ=255)

                                              cgAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCAc                                     >  1:9756221/1‑93 (MQ=255)

                                              cgAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCAc                                     <  2:5172022/93‑1 (MQ=255)

                                               gAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCAcc                                    <  1:5177373/93‑1 (MQ=255)

                                               gAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCAcc                                    <  2:7444950/93‑1 (MQ=255)

                                                aGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCg                                   <  1:6789499/93‑1 (MQ=255)

                                                 gggCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGc                                  >  1:2543801/1‑93 (MQ=255)

                                                   gccgccGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGCGATCTTCTGCACGCACGTTACGAAAAGCACCGccc                                >  1:9820469/1‑93 (MQ=255)

                                                   gccgccGCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGccc                                <  2:3479162/93‑1 (MQ=255)

                                                    ccgccgCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCt                               >  1:7861438/1‑93 (MQ=255)

                                                    ccgccgCTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCt                               <  1:127124/93‑1 (MQ=255)

                                                      gccgcTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTc                             >  2:1220858/1‑93 (MQ=255)

                                                       ccgcTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCg                            >  1:9263401/1‑93 (MQ=255)

                                                       ccgcTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCg                            >  2:5990094/1‑93 (MQ=255)

                                                        cgcTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGa                           >  2:3310518/1‑93 (MQ=255)

                                                        cgcTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGa                           <  2:1894935/93‑1 (MQ=255)

                                                         gcTTCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGaa                          >  1:1841382/1‑93 (MQ=255)

                                                          cTTCCAGCGCCGGGTTACCCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGCAc                         >  1:7185911/1‑93 (MQ=255)

                                                           ttCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACg                        >  2:7385043/1‑93 (MQ=255)

                                                            tCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACgg                       >  2:6572731/1‑93 (MQ=255)

                                                            tCCAGCGCCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACgg                       <  2:2822448/93‑1 (MQ=255)

                                                                gcgcCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCa                   <  1:2364824/93‑1 (MQ=255)

                                                                gcgcCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCa                   >  1:9286509/1‑93 (MQ=255)

                                                                 cgcCGGGTTACGCACCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCAt                  <  1:8524664/93‑1 (MQ=255)

                                                                  gcCGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCATc                 <  2:9834365/93‑1 (MQ=255)

                                                                   ccGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCATCg                >  1:284356/1‑93 (MQ=255)

                                                                    cGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCATCgg               <  2:6762668/93‑1 (MQ=255)

                                                                    cGGGTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCATCgg               >  1:6564626/1‑93 (MQ=255)

                                                                      ggTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCATCgggg             >  2:5814444/1‑93 (MQ=255)

                                                                      ggTTACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCATCgggg             >  2:3106753/1‑93 (MQ=255)

                                                                        ttACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCATCGGGGTg           >  2:9337140/1‑93 (MQ=255)

                                                                         tACGCGCCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCATCGGGGTGc          >  2:6487061/1‑93 (MQ=255)

                                                                            gcgcCAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCATCGGGGTGCGaa       <  2:2953100/93‑1 (MQ=255)

                                                                               ccAGCAGGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCATCGGGGTGCGAATTg    >  1:3212022/1‑93 (MQ=255)

                                                                                cagcagGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCATCGGGGTGCGCATTGc   >  1:8858052/1‑93 (MQ=255)

                                                                                cagcagGCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCATCGGGGTGCGAATTGc   <  1:119758/93‑1 (MQ=255)

                                                                                 agcaggCTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCATCGGGGTGCGAATTGCa  <  2:9274759/93‑1 (MQ=255)

||

AAACCCTGCTCCAGCGTCTGCTGCACACAACCCCAGTAAATATCGTCGAGGGCCGCCGCTTCCAGCGCCGGGTTACGCGCCAGCAGCGTACGCATTAAATGAGCGGAGAGATCTTCTGCACGCACGTTACGAAAAGCACCGCCCTTCGAACGGCCCATCGGGGTGCGAATTGCA  >  NC_000913/4026166‑4026339

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