BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,549,566	A→G	D121G (GAC→GGC)	uxuA →	mannonate hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,549,566	0	A	G	97.5%	616.9 / 3.3	200	D121G (GAC→GGC)	uxuA	mannonate hydrolase
Reads supporting (aligned to +/- strand): ref base A (1/4); new base G (89/106); total (90/110)
Fisher's exact test for biased strand distribution p-value = 3.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.56e-01

Only 100 of 207 total aligned reads displayed.

ACCTGGCGCAGTGTGGCATTCGCACCGTGTGCTACAACTTCATGCCGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGACGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCGAAATGCATATCCTGAAACGCCCAGGCGCGGAAGCGGATTACA  >  NC_000913/4549476‑4549655

aCCTGGCGCAGTGTGGCATTCGCACCGTGTGCTACAACTTCATGCCGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCg                                                                                         >  1:3512355/1‑93 (MQ=255)

 ccTGGCGCAGTGTGGCATTCGCACCGTGTGCTACAACTTCATGCCGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCgg                                                                                        <  2:3171815/93‑1 (MQ=255)

  cTGGCGCAGTGTGTCATTCGCACCGTGTGCTACAACTTCATGCCGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAggcggc                                                                                       <  1:2996409/93‑1 (MQ=255)

  cTGGCGCAGTGTGGCATTCGCACCGTGTGCTACAACTTCATGCCGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAggcggc                                                                                       <  2:825101/93‑1 (MQ=255)

   tGGCGCAGTGTGGCATTCGCACCGTGTGCTACAACTTCATGCCGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCt                                                                                      >  1:516771/1‑93 (MQ=255)

   tGGCGCAGTGTGGCATTCGCACCGTGTGCTACAACTTCATGCCGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCt                                                                                      <  2:507362/93‑1 (MQ=255)

   tGGCGCAGTGTGGCATTCGCACCGTGTGCTACAACTTCATGCCGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCt                                                                                      >  1:381868/1‑93 (MQ=255)

   tGGCGCAGTGTGGCATTCGCACCGTGTGCTACAACTTCATGCCGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCt                                                                                      <  1:2267279/93‑1 (MQ=255)

    ggCGCAGTGTGGCATTCGCACCGTGTGCTACAACTTCATGCTGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTc                                                                                     <  2:1370202/93‑1 (MQ=255)

    ggCGCAGTGTGGCATTCGCACCGTGTGCTACAACTTCATGCCGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTc                                                                                     >  2:1222719/1‑93 (MQ=255)

        cAGTGTGGCATTCGCACCGTGTGCTACAACTTCATGCCGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCaaa                                                                                 <  2:4543510/93‑1 (MQ=255)

            gtgGCATTCGCACCGTGTGCTACAACTTCATGCCGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGctc                                                                             <  1:863915/93‑1 (MQ=255)

              ggCATTCGCACCGTGTGCTACAACTTCATGCCGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTg                                                                           >  2:1084060/1‑93 (MQ=255)

               gCATTCGCACCGTGTGCTACAACTTCATGCCGGTGCTCAACTGGACCCGTACTGGCCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTgc                                                                          <  1:1038657/93‑1 (MQ=255)

                 aTTCGCACCGTGTGCTACAACTTCATGCCGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTgcgc                                                                        <  2:3364650/93‑1 (MQ=255)

                 aTTCGCACCGTGTGCTACAACTTCATGCCGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTgcgc                                                                        <  2:833300/93‑1 (MQ=255)

                 aTTCGCACCGTGTGCTACAACTTCATGCCGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTgcgc                                                                        <  1:4579528/93‑1 (MQ=255)

                 aTTCGCACCGTGTGCTACAACTTCATGCCGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTgcgc                                                                        >  1:3458758/1‑93 (MQ=255)

                 aTTCGCACCGTGTGCTACAACTTCATGCCGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTgcgc                                                                        >  2:1804218/1‑93 (MQ=255)

                  ttCGCACCGTGTGCTACAACTTCATGCCGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCt                                                                       >  2:3063806/1‑93 (MQ=255)

                   tCGCACCGTGTGCTACAACTTCATGCCGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCtt                                                                      <  2:3634868/93‑1 (MQ=255)

                   tCGCACCGTGTGCTACAACTTCATGCCGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGACGGCTCCAAAGCTCTGCGCtt                                                                      <  2:4010468/93‑1 (MQ=255)

                    cGCACCGTGTGCTACAACTTCATGCCGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTc                                                                     >  1:2703098/1‑93 (MQ=255)

                     gCACCGTGTGCTACAACTTCATGCCGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCg                                                                    >  2:2171742/1‑93 (MQ=255)

                       cccGTGTGCTACAACTTCATGCCGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGAc                                                                  <  1:583002/92‑1 (MQ=255)

                       aCCGTGTGCTACAACTTCATGCCGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGAc                                                                  >  2:1198281/1‑93 (MQ=255)

                        ccGTGTGCTACAACTTCATGCCGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGAcc                                                                 >  2:4358139/1‑93 (MQ=255)

                        ccGTGTGCTACAACTTCATGCCAGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGACGGCTCCAAAGCTCTGCGCTTCGAcc                                                                 <  2:3995371/93‑1 (MQ=255)

                         cGTGTGCTACAACTTCATGCCGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCa                                                                <  2:1737312/93‑1 (MQ=255)

                         cGTGTGCTACAACTTCATGCCGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCa                                                                >  2:127812/1‑93 (MQ=255)

                         cGTGTGCTACAACTTCATGCCGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCa                                                                <  2:3842634/93‑1 (MQ=255)

                          gtgtGCTACAACTTCATGCCGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAg                                                               >  1:2328513/1‑93 (MQ=255)

                           tgtgCTACAACTTCATGCCGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGa                                                              <  1:2617116/93‑1 (MQ=255)

                           tgtgCTACAACTTCATGCCGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGa                                                              <  1:49211/93‑1 (MQ=255)

                            gtgCTACAACTTCATGCCGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAgct                                                             <  2:3292077/93‑3 (MQ=255)

                                tACAACTTCATGCCGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGat                                                         >  1:1394056/1‑92 (MQ=255)

                                tACAACTTCATGCCGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGaa                                                         <  2:2069511/93‑1 (MQ=255)

                                tACAACTTCATGCCGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGaa                                                         >  2:2189515/1‑93 (MQ=255)

                                tACAACTTCATGCCGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGaa                                                         >  1:118702/1‑93 (MQ=255)

                                  cAACTTCATGCCGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAAtt                                                       <  1:2780287/93‑1 (MQ=255)

                                     cTTCATGCCGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGc                                                    >  1:4081852/1‑93 (MQ=255)

                                        cATGCCGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGc                                                 >  2:2611021/1‑93 (MQ=255)

                                        cATGCCGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGc                                                 <  2:2419628/93‑1 (MQ=255)

                                          tGCCGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCAt                                               <  1:4608400/93‑1 (MQ=255)

                                          tGCCGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCAt                                               >  1:88486/1‑93 (MQ=255)

                                          tGCCGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCAt                                               >  2:34412/1‑93 (MQ=255)

                                             cGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCg                                            <  1:1347357/93‑1 (MQ=255)

                                                tGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGCCCAGATCGAATTCGCTCCATTctaca                                         >  2:3877632/1‑89 (MQ=255)

                                                tGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCGaaa                                         >  1:1324378/1‑93 (MQ=255)

                                                 gCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCAGCGCTTCGCCCATCTCGAATTCGCTGCATTCGAAAt                                        >  1:3236580/1‑93 (MQ=255)

                                                    cGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCGAAATGCa                                     >  2:4243696/1‑93 (MQ=255)

                                                    cGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCGAAATGCa                                     <  1:4341165/93‑1 (MQ=255)

                                                     gACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCGAAATGCat                                    <  2:3364094/93‑1 (MQ=255)

                                                     gACTGGACCCGCACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCGAAATGCat                                    >  2:2679635/1‑93 (MQ=255)

                                                      aCTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCGAAATGCata                                   >  1:1284742/1‑93 (MQ=255)

                                                      aCTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCGAAATGCata                                   >  1:89484/1‑93 (MQ=255)

                                                       cTGGACCCGTACTGCCCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCGAAATGCatat                                  <  2:4142556/93‑1 (MQ=255)

                                                       cTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCGAAATGCatat                                  <  1:450250/93‑1 (MQ=255)

                                                        tGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCGAAATGCATATc                                 >  2:2496509/1‑93 (MQ=255)

                                                         ggACCCGTACTGACCTCGAATACGTGCTGCCAGACGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCGAAATGCATATcc                                <  1:1579972/93‑1 (MQ=255)

                                                          gACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCGAAATGCATATCCt                               <  2:1318600/93‑1 (MQ=255)

                                                          gACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCGAAATGCATATCCt                               >  2:224240/1‑93 (MQ=255)

                                                           aCCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCGAAATGCATATCCTg                              <  2:1645807/93‑1 (MQ=255)

                                                            cccGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCGAAATGCATATCCTGa                             <  2:580825/93‑1 (MQ=255)

                                                             ccGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCGAAATGCATATCCTGaa                            >  1:909188/1‑93 (MQ=255)

                                                              cGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCGAAATGCATATCCTGaaa                           >  2:2888644/1‑93 (MQ=255)

                                                               gTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCGAAATGCATATCCTGAAAc                          >  2:3250790/1‑93 (MQ=255)

                                                               gTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCGAAATGCATATCCTGAAAc                          <  1:729667/93‑1 (MQ=255)

                                                                tACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCGAAATGCATATCCTGAAACg                         >  2:3535856/1‑93 (MQ=255)

                                                                tACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCGAAATGCATATCCTGAAACg                         >  1:4706259/1‑93 (MQ=255)

                                                                tACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCGAAATGCATATCCTGAAACg                         <  1:3387563/93‑1 (MQ=255)

                                                                   tGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCGAAATGCATATCCTGAAACGccc                      >  2:570405/1‑93 (MQ=255)

                                                                   tGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCGAAATGCATATCCTGAAACGccc                      >  2:1508409/1‑93 (MQ=255)

                                                                     aCCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGTTTCGACCAGATCGAATTCGCTGCATTCGAAATGCATATCCTGAAACGCCCAg                    <  2:4096930/93‑1 (MQ=255)

                                                                     aCCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCGAAATGCATATCCTGAAACGCCCAg                    <  1:4641391/93‑1 (MQ=255)

                                                                     aCCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCGAAATGCATATCCTGAAACGCCCAg                    <  2:633133/93‑1 (MQ=255)

                                                                      ccTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCGAAATGCATATCCTGAAACGCCCAgg                   <  2:1965699/93‑1 (MQ=255)

                                                                      ccTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCGAAATGCATATCCTGAAACGCCCAgg                   <  1:295506/93‑1 (MQ=255)

                                                                      ccTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCGAAATGCATATCCTGAAACGCCCAgg                   <  1:4005959/93‑1 (MQ=255)

                                                                       cTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCGAAATGCATATCCTGAAACGCCCAGgc                  <  2:2854440/93‑1 (MQ=255)

                                                                       cTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCGAAATGCATATCCTGAAACGCCCAGgc                  <  1:3227718/93‑1 (MQ=255)

                                                                           aaTACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCGAAATGCATATCCTGAAACGCCCAGGCGCgg              >  2:2691028/1‑93 (MQ=255)

                                                                           aaTACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCGAAATGCATATCCTGAAACGCCCAGGCGCgg              >  2:3721731/1‑93 (MQ=255)

                                                                           aaTACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCGAAATGCATATCCTGAAACGCCCAGGCGCgg              >  1:3866570/1‑93 (MQ=255)

                                                                            aTACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCGAAATGCATATCCTGAAACGCCCAGGCGCGGa             >  1:3732448/1‑93 (MQ=255)

                                                                             tACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCGAAATGCATATCCTGAAACGCCCAGGCGCGGaa            >  1:2485967/1‑93 (MQ=255)

                                                                             tACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCGAAATGCATATCCTGAAACGCCCAGGCGCGGaa            <  2:1715988/93‑1 (MQ=255)

                                                                             tACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCGAAATGCATATCCTGAAACGCCCAGGCGCGGaa            <  1:4673571/93‑1 (MQ=255)

                                                                              aCGTGCTGCCCGGCGTCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCGAAATGCATATCCTGAAACGCCCAGGCGCGGAAg           <  1:1627025/93‑1 (MQ=255)

                                                                               cGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCGAAATGCATATCCTGAAACGCCCAGGCGCGGAAGc          <  1:2063603/93‑1 (MQ=255)

                                                                                gTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCGAAATGCATATCCTGAAACGCCCAGGCGCGGAAGCg         <  1:4291246/93‑1 (MQ=255)

                                                                                gTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCGAAATGCATATCCTGAAACGCCCAGGCGCGGAAGCg         <  1:4749934/93‑1 (MQ=255)

                                                                                gTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCGAAATGCATATCCTGAAACGCCCAGGCGCGGAAGCg         >  2:1041530/1‑93 (MQ=255)

                                                                                   ctgcCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCGAAATGCATATCCTGAAACGCCCAGGCGCGGAAGCGGAt      <  1:2180915/93‑1 (MQ=255)

                                                                                    tgcCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCGAAATGCATATCCTGAAACGCCCAGGCGCGGAAGCGGAtt     <  2:1459920/93‑1 (MQ=255)

                                                                                    tgcCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCGAAATGCATATCCTGAAACGCCCAGGCGCGGAAGCGGAtt     <  1:3005381/93‑1 (MQ=255)

                                                                                    tgcCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCGAAATGCATATCCTGAAACGCCCAGGCGCGGAAGCGGAtt     >  2:2778823/1‑93 (MQ=255)

                                                                                    tgcCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCGAAATGCATATCCTGAAACGCCCAGGCGCGGAAGCGGAtt     >  1:1244174/1‑93 (MQ=255)

                                                                                     gcCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCGAAATGCATATCCTGAAACGCCCAGGCGCGGAAGCGGATTa    <  1:4470310/93‑1 (MQ=255)

                                                                                      ccAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCGAAATGCATATCCTGAAACGCCCAGGCGCGGAAGCGGATTac   <  2:109938/93‑1 (MQ=255)

ACCTGGCGCAGTGTGGCATTCGCACCGTGTGCTACAACTTCATGCCGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGACGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCGAAATGCATATCCTGAAACGCCCAGGCGCGGAAGCGGATTACA  >  NC_000913/4549476‑4549655

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