BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,943,196	T→C	K223K (AAA→AAG)	ygdL ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,943,196	0	T	C	99.7%	902.1 / ‑5.3	284	K223K (AAA→AAG)	ygdL	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (141/142); total (141/143)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.25e-01

Only 100 of 293 total aligned reads displayed.

AACCGCAACAAACCCAAAGGTGGCGGTCACCATCGTTGCCGCGCCAAATCCAGATGCACAATCCATCCGCTTCGGCCCTTCTGCCGTGGCTTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTT  >  NC_000913/2943106‑2943285

aaCCGCAACAAACCCAAAGGTGGCGGTCACCATCGTTGCCGCGCCAAATCCAGATGCACAATCCATCCGCTTCGGCCCTTCTGCCGTGGCctt                                                                                         <  2:3874617/93‑1 (MQ=255)

aaCCGCAACAAACCCAAAGGTGGCGGTCACCATCGTTGCCGCGCCAAATCCAGATGCACAATCCATCCGCTTCGGCCCTTCTGCCGTGGCctt                                                                                         <  1:6266913/93‑1 (MQ=255)

 aCCGCAACAAACCCAAAGGTGGCGGTCACCATCGTTGCCGCGCCAAATCCAGATGCACAATCCATCCGCTTCGGCCCTTCTGCCGTGGCCTTc                                                                                        >  2:2508665/1‑93 (MQ=255)

   cGCAACAAACCCAAAGGTGGCGGTCACCATCGTTGCCGCGCCAAATCCAGATGCACAATCCATCCGCTTCGGCCCTTCTGCCGTGGCCTTCAt                                                                                      >  1:4970514/1‑93 (MQ=255)

   cGCAACAAACCCAAAGGTGGCGGTCACCATCGTTGCCGCGCCAAATCCAGATGCACAATCCATCCGCTTCGGCCCTTCTGCCGTGGCCTTCAt                                                                                      <  1:7328433/93‑1 (MQ=255)

   cGCAACAAACCCAAAGGTGGCGGTCACCATCGTTGCCGCGCCAAATCCAGATGCACAATCCATCCGCTTCGGCCCTTCTGCCGTGGCCTTCAt                                                                                      <  1:2808404/93‑1 (MQ=255)

    gCAACAAACCCAAAGGTGGCGGTCACCATCGTTGCCGCGCCAAATCCAGATGCACAATCCATCCGCTTCGGCCCTTCTGCCGTGGCCTTCATc                                                                                     >  2:5741811/1‑93 (MQ=255)

    gCAACAAACCCAAAGGTGGCGGTCACCATCGTTGCCGCGCCAAATCCAGATGCACAATCCATCCGCTTCGGCCCTTCTGCCGTGGCCTTCATc                                                                                     >  1:5621900/1‑93 (MQ=255)

    gCAACAAACCCAAAGGTGGCGGTCACCATCGTTGCCGCGCCAAATCCAGATGCACAATCCATCCGCTTCGGCCCTTCTGCCGTGGCCTTCATc                                                                                     <  2:4556875/93‑1 (MQ=255)

       acaaACCCAAAGGTGGCGGTCACCATCGTTGCCGCGCCAAATCCAGATGCACAATCCATCCGCTTCGGCCCTTCTGCCGTGGCCTTCATCGca                                                                                  >  1:2013475/1‑93 (MQ=255)

        caaACCCAAAGGTGGCGGTCACCATCGTTGCCGCGCCAAATCCAGATGCACAATCCATCCGCTTCGGCCCTTCTGCCGTGGCCTTCATCGcac                                                                                 >  1:6048341/1‑93 (MQ=255)

        caaACCCAAAGGTGGCGGTCACCATCGTTGCCGCGCCAAATCCAGATGCACAATCCATCCGCTTCGGCCCTTCTGCCGTGGCCTTCATCGcac                                                                                 <  1:7714384/93‑1 (MQ=255)

          aaCCCAAAGGTGGCGGTCACCATCGTTGCCGCGCCAAATCCAGATGCACAATCCATCCGCTTCGGCCCTTCTGCCGTGGCCTTCATCGcacac                                                                               >  2:2372302/1‑93 (MQ=255)

          aaCCCAAAGGTGGCGGTCACCATCGTTGCCGCGCCAAATCCAGATGCACAATCCATCCGCTTCGGCCCTTCTGCCGTGGCCTTCATCGcacac                                                                               <  1:7068123/93‑1 (MQ=255)

           aCCCAAAGGTGGCGGTCACCATCGTTGCCGCGCCAAATCCAGATGCACAATCCATCCGCTTCGGCCCTTCTGCCGTGGCCTTCATCGcacaca                                                                              <  2:7870895/93‑1 (MQ=255)

            cccAAAGGTGGCGGTCACCATCGTTGCCGCGCCAAATCCAGATGCACAATCCATCCGCTTCGGCCCTTCTGCCGTGGCCTTCATCGCacacac                                                                             <  2:6642245/93‑1 (MQ=255)

            cccAAAGGTGGCGGTCACCATCGTTGCCGCGCCAAATCCAGATGCACAATCCATCCGCTTCGGCCCTTCTGCCGTGGCCTTCATCGCacacac                                                                             <  2:3135624/93‑1 (MQ=255)

               aaaGGTGGCGGTCACCATCGTTGCCGCGCCAAATCCAGATGCACAATCCATCCGCTTCGGCCCTTCTGCCGTGGCCTTCATCGCACACaccgt                                                                          >  1:1093200/1‑93 (MQ=255)

               aaaGGTGGCGGTCACCATCGTTGCCGCGCCAAATCCAGATGCACAATCCATCCGCTTCGGCCCTTCTGCCGTGGCCTTCATCGCACACaccgt                                                                          >  1:2247530/1‑93 (MQ=255)

                 aGGTGGCGGTCACCATCGTTGCCGCGCCAAATCCAGATGCACAATCCATCCGCTTCGGCCCTTCTGCCGTGGCCTTCATCGCACACaccgtac                                                                        <  2:4935983/93‑1 (MQ=255)

                   gTGGCGGTCACCATCGTTGCCGCGCCAAATCCAGATGCACAATCCATCCGCTTCGGCCCTTCTGCCGTGGCCTTCATCGCACACaccgtaccg                                                                      >  1:3872564/1‑93 (MQ=255)

                   gTGGCGGTCACCATCGTTGCCGCGCCAAATCCAGATGCACAATCCATCCGCTTCGGCCCTTCTGCCGTGGCCTTCATCGCACACaccgtaccg                                                                      <  1:5348775/93‑1 (MQ=255)

                    tGGCGGTCACCATCGTTGCCGCGCCAAATCCAGATGCACAATCCATCCGCTTCGGCCCTTCTGCCGTGGCCTTCATCGCACACaccgtaccgt                                                                     >  1:6557028/1‑93 (MQ=255)

                      gCGGTCACCATCGTTGCCGCGCCAAATCCAGATGCACAATCCATCCGCTTCGGCCCTTCTGCCGTGGCCTTCATCGCACACACCGTACCGTCt                                                                   <  2:593700/93‑1 (MQ=255)

                        ggTCACCATCGTTGCCGCGCCAAATCCAGATGCACAATCCATCCGCTTCGGCCCTTCTGCCGTGGCCTTCATCGCACACACCGTACCGTCTGa                                                                 >  1:6930523/1‑93 (MQ=255)

                        ggTCACCATCGTTGCCGCGCCAAATCCAGATGCACAATCCATCCGCTTCGGCCCTTCTGCCGTGGCCTTCATCGCACACACCGTACCGTCTGa                                                                 >  2:7719359/1‑93 (MQ=255)

                        ggTCACCATCGTTGCCGCGCCAAATCCAGATGCACAATCCATCCGCTTCGGCCCTTCTGCCGTGGCCTTCATCGCACACACCGTACCGTCTGa                                                                 <  2:4615929/93‑1 (MQ=255)

                        ggTCACCATCGTTGCCGCGCCAAATCCAGATGCACAATCCATCCGCTTCGGCCCTTCTGCCGTGGCCTTCATCGCACACACCGTACCGTCTGa                                                                 <  1:2874584/93‑1 (MQ=255)

                        ggTCACCATCGTTGCCGCGCCAAATCCAGATGCACAATCCATCCGCTTCGGCCCTTCTGCCGTGGCCTTCATCGCACACACCGTACCGTCTGa                                                                 <  2:6340792/93‑1 (MQ=255)

                          tCACCATCGTTGCCGCGCCAAATCCAGATGCACAATCCATCCGCTTCGGCCCTTCTGCCGTGGCCTTCATCGCACACACCGTACCGTCTGACt                                                               >  1:6401039/1‑93 (MQ=255)

                             ccATCGTTGCCGCGCCAAATCCAGATGCACAATCCATCCGCTTCGGCCCTTCTGCCGTGGCCTTCATCGCACACACCGTACCGTCTGACTGCg                                                            <  1:4733595/93‑1 (MQ=255)

                              cATCGTTGCCGCGCCAAATCCAGATGCACAATCCATCCGCTTCGGCCCTTCTGCCGTGGCCTTCATCGCACACACCGTACCGTCTGACTGCgg                                                           <  1:764865/93‑1 (MQ=255)

                               aTCGTTGCCGCGCCAAATCCAGATGCACAATCCATCCGCTTCGGCCCTTCTGCCGTGGCCTTCATCGCACACACCGTACCGTCTGACTGCggg                                                          >  2:1988206/1‑93 (MQ=255)

                                tCGTTGCCGCGCCAAATCCAGATGCACAATCCATCCGCTTCGGCCCTTCTGCCGTGGCCTTCATCGCACACACCGTACCGTCTGACTGCGGGt                                                         >  2:2164607/1‑93 (MQ=255)

                                 cGTTGCCGCGCCAAATCCAGATGCACAATCCATCCGCTTCGGCCCTTCTGCCGTGGCCTTCATCGCACACACCGTACCGTCTGACTGCGGGTa                                                        >  2:1542836/1‑93 (MQ=255)

                                 cGTTGCCGCGCCAAATCCAGATGCACAATCCATCCGCTTCGGCCCTTCTGCCGTGGCCTTCATCGCACACACCGTACCGTCTGACTGCGGGTa                                                        <  2:442196/93‑1 (MQ=255)

                                  gTTGCCGCGCCAAATCCAGATGCACAATCCATCCGCTTCGGCCCTTCTGCCGTGTCCTTCATCGCATACACCGTATCGTCTGACTGCGGGGac                                                       >  2:7036418/1‑93 (MQ=255)

                                  gTTGCCGCGCCAAATCCAGATGCACAATCCATCCGCTTCGGCCCTTCTGCCGTGGCCTTCATCGCACACACCGTACCGTCTGACTGCGGGTac                                                       >  1:2503509/1‑93 (MQ=255)

                                   ttGCCGCGCCAAATCCAGATGCACAATCCATCCGCTTCGGCCCTTCTGCCGTGGCCTTCATCGCACACACCGTACCGTCTGACTGCGGGTaca                                                      >  1:3237766/1‑93 (MQ=255)

                                   ttGCCGCGCCAAATCCAGATGCACAATCCATCCGCTTCGGCCCTTCTGCCGTGGCCTTCATCGCACACACCGTACCGTCTGACTGCGGGTaca                                                      >  1:4215378/1‑93 (MQ=255)

                                     gCCGCGCCAAATCCAGATGCACAATCCATCCGCTTCGGCCCTTCTGCCGTGGCCTTCATCGCACACACCGTACCGTCTGACTGCGGGTACAcc                                                    >  2:5221046/1‑93 (MQ=255)

                                      ccGCGCCAAATCCAGATGCACAATCCATCCGCTTCGGCCCTTCTGCCGTGGCCTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCa                                                   >  1:556941/1‑93 (MQ=255)

                                        gcgcCAAATCCAGATGCACAATCCATCCGCTTCGGCCCTTCTGCCGTGGCCTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAgc                                                 >  2:3418048/1‑93 (MQ=255)

                                           ccAAATCCAGATGCACAATCCATCCGCTTCGGCCCTTCTGCCGTGGCCTTCATCGCACACACCGTACCGTCTGACTGCGGGCACACCAGCGCt                                              >  1:3435669/1‑93 (MQ=255)

                                            cAAATCCAGATGCACAATCCATCCGCTTCGGCCCTTCTGCCGTGGCCTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCtt                                             <  1:6216582/93‑1 (MQ=255)

                                            cAAATCCAGATGCACAATCCATCCGCTTCGGCCCTTCTGCCGTGGCCTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCtt                                             >  2:2054188/1‑93 (MQ=255)

                                             aaaTCCAGATGCACAATCCATCCGCTTCGGCCCTTCTGCCGTGGCCTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTc                                            >  2:11005/1‑93 (MQ=255)

                                              aaTCCAGATGCACAATCCATCCGCTTCGGCCCTTCTGCCGTGGCCTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCa                                           <  2:7681224/93‑1 (MQ=255)

                                              aaTCCAGATGCACAATCCATCCGCTTCGGCCCTTCTGCCGTGGCCTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCa                                           <  2:206306/93‑1 (MQ=255)

                                                tCCAGATGCACAATCCATCCGCTTCGGCCCTTCTGCCGTGGCCTTCATCGCACCCACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGt                                         <  2:6988164/93‑1 (MQ=255)

                                                tCCAGATGCACAATCCATCCGCTTCGGCCCTTCTGCCGTGGCCTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGt                                         >  2:2248241/1‑93 (MQ=255)

                                                 ccAGATGCACAATCCATCCGCTTCGGCCCTTCTGCCGTGGCCTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTa                                        <  1:5059991/93‑1 (MQ=255)

                                                 ccAGATGCACAATCCATCCGCTTCGGCCCTTCTGCCGTGGCCTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCACTa                                        <  2:4153205/93‑1 (MQ=255)

                                                  cAGATGCACAATCCATCCGCTTCGGCCCTTCTGCCGTGGCCTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAg                                       <  2:7468880/93‑1 (MQ=255)

                                                    gATGCACAATCCATCCGCTTCGGCCCTTCTGCCGTGGCCTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGaa                                     >  2:1730180/1‑93 (MQ=255)

                                                     aTGCACAATCCATCCGCTTCGGCCCTTCTGCCGTGGCCTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGaaa                                    <  1:4000305/93‑1 (MQ=255)

                                                      tGCACAATCCATCCGCTTCGGCCCTTCTGCCGTGGCCTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGaaaa                                   <  2:1189964/93‑1 (MQ=255)

                                                      tGCACAATCCATCCGCTTCGGCCCTTCTGCCGTGGCCTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGaaaa                                   <  1:6655073/93‑1 (MQ=255)

                                                       gCACAATCCATCCGCTTCGGCCCTTCTGCCGTGGCCTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAac                                  <  1:4528814/93‑1 (MQ=255)

                                                       gCACAATCCATCCGCTTCGGCCCTTCTGCCGTGGCCTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAac                                  <  2:5410124/93‑1 (MQ=255)

                                                        cacaATCCATCCGCTTCGGCCCTTCTGCGGTGGCCTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAaca                                 <  1:4349601/93‑1 (MQ=255)

                                                        cacaATCCATCCGCTTCGGCCCTTCTGCCGTGTCCTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAaca                                 <  2:7255220/93‑1 (MQ=255)

                                                        cacaATCCATCCGCTTCGGCCCTTCTGCCGTGGCCTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAaca                                 <  1:3024359/93‑1 (MQ=255)

                                                        cacaATCCATCCGCTTCGGCCCTTCTGCCGTGGCCTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAaca                                 <  2:6169468/93‑1 (MQ=255)

                                                        cacaATCCATCCGCTTCGGCCCTTCTGCCGTGGCCTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAaca                                 <  1:6575620/93‑1 (MQ=255)

                                                        cacaATCCATCCGCTTCGGCCCTTCTGCCGTGGCCTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAaca                                 >  2:722991/1‑93 (MQ=255)

                                                         acaATCCATCCGCTTCGGCCCTTCTGCCGTGGCCTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAacac                                <  1:5114435/93‑1 (MQ=255)

                                                          caATCCATCCGCTTCGGCCCTTCTGCCGTGGCCTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACg                               <  2:2565204/93‑1 (MQ=255)

                                                             tccatccGCTTCGGCCCTTCTGCCGTGGCCTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGcca                            >  2:4651737/1‑91 (MQ=255)

                                                             tccatccGCTTCGGCCCTTCTGCCGTGGCCTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCaa                            <  2:3996477/93‑1 (MQ=255)

                                                               catccGCTTCGGCCCTTCTGCCGTGGCCTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATc                          <  2:2430068/93‑1 (MQ=255)

                                                                ctccGCTTCGGCCCTTCTGCCGTGGCCTTCAGCGCACACACCGTACCGTCTGACTCCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATcc                         <  2:2655521/92‑1 (MQ=255)

                                                                atccGCTTCGGCCCTTCTGCCGTGGCCTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATcc                         >  1:1675645/1‑93 (MQ=255)

                                                                 tccGCTTCGGCCCTTCTGCCGTGGCCTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATccc                        >  2:5692767/1‑92 (MQ=255)

                                                                 tccGCTTCGGCCCTTCTGCCGTGGCCTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCa                        <  1:7840088/93‑1 (MQ=255)

                                                                 tccGCTTCGGCCCTTCTGCCGTGGCCTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCa                        <  1:1293654/93‑1 (MQ=255)

                                                                 tccGCTTCGGCCCTTCTGCCGTGGCCTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCa                        <  1:6841257/93‑1 (MQ=255)

                                                                  cctcTTCGGCCCTTCTGCCGTGGCCTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCAc                       <  1:2345642/90‑1 (MQ=255)

                                                                   cGCTTCGGCCCTTCTGCCGTGGCCTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACg                      >  1:5798595/1‑93 (MQ=255)

                                                                     cTTCGGCCCTTCTGCCGTGGCCTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGcc                    >  1:1129006/1‑93 (MQ=255)

                                                                      ttCGGCCCTTCTGCCGTGGCCTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCg                   >  2:5844169/1‑93 (MQ=255)

                                                                         ggCCCTTCTGCCGTGGCCTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGc                >  1:4444388/1‑93 (MQ=255)

                                                                          gCCCTTCTGCCGTGGCCTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCt               <  2:2871657/93‑1 (MQ=255)

                                                                          gCCCTTCTGCCGTGGCCTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCt               <  2:2793597/93‑1 (MQ=255)

                                                                           cccTTCTGCCGTGGCCTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCtt              >  1:3698457/1‑93 (MQ=255)

                                                                            ctttCTGCCGTGGCCTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTa             <  2:430363/91‑1 (MQ=255)

                                                                             cTTCTGCCGTGGCCTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTAc            <  1:1875487/93‑1 (MQ=255)

                                                                             cTTCTGCCGTGGCCTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTAc            >  1:1626871/1‑93 (MQ=255)

                                                                              ttCTGCCGTGGCCTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCTACGCCGAGCTTAcc           >  1:7188785/1‑93 (MQ=255)

                                                                              ttCTGCCGTGGCCTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCa                        >  2:560658/1‑80 (MQ=255)

                                                                              ttCTGCCGTGGCCTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTAcc           >  1:4943692/1‑93 (MQ=255)

                                                                                cTGCCGTGGCCTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCtt         <  2:4718558/93‑1 (MQ=255)

                                                                                cTGCCGTGGCCTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCtt         <  2:5115798/93‑1 (MQ=255)

                                                                                 tGCCGTGGCCTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCttt        >  1:3478604/1‑93 (MQ=255)

                                                                                   ccGTGGCCTTCATCGCACACACCGTACCGTCTGACTGCTGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTAc      <  1:984572/93‑1 (MQ=255)

                                                                                   ccGTGGCCTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTAc      <  1:5706557/93‑1 (MQ=255)

                                                                                    cGTGGCCTTCATCGCCCACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTagt     <  1:2465510/93‑3 (MQ=255)

                                                                                     gTGGCCTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTg    >  2:2329601/1‑93 (MQ=255)

                                                                                      tGGCCTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGt   <  1:156012/93‑1 (MQ=255)

                                                                                      gggCCTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGt   <  1:2677601/92‑1 (MQ=255)

AACCGCAACAAACCCAAAGGTGGCGGTCACCATCGTTGCCGCGCCAAATCCAGATGCACAATCCATCCGCTTCGGCCCTTCTGCCGTGGCTTTCATCGCACACACCGTACCGTCTGACTGCGGGTACACCAGCGCTTCAGTAGAAAACACGCAATCCACGCCGAGCTTACCTTTACTGTT  >  NC_000913/2943106‑2943285

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