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	99.6%	753.3 / ‑5.1	235	D121G (GAC→GGC)	uxuA	mannonate hydrolase
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (112/122); total (113/122)
Fisher's exact test for biased strand distribution p-value = 4.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

Only 100 of 241 total aligned reads displayed.

ACCTGGCGCAGTGTGGCATTCGCACCGTGTGCTACAACTTCATGCCGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGACGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCGAAATGCATATCCTGAAACGCCCAGGCGCGGAAGCGGATTACA  >  NC_000913/4549476‑4549655

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

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

 ccTGGCGCAGTGTGGCATTCGCACCGTGTGCTACAACTTCATGCCGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCgg                                                                                        >  2:2128871/1‑93 (MQ=255)

  cTGGCGCAGTGTGGCATTCGCACCGTGTGCTACAACTTCATGCCGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGACGGc                                                                                       >  2:2700917/1‑93 (MQ=255)

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

    ggCGCAGTGTGGCATTCGCACCGTGTGCTACAACTTCATGCCGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTc                                                                                     <  2:6227366/93‑1 (MQ=255)

      cgcAGTGTGGCATTCGCACCGTGTGCTACAACTTCATGCCGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCa                                                                                   <  1:2983550/93‑1 (MQ=255)

       gcAGTGTGGCATTCGCACCGTGTGCTACAACTTCATGCCGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCaa                                                                                  <  2:3619040/93‑1 (MQ=255)

         aGTGTGGCATTCGCACCGTGTGCTACAACTTCATGCCGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAg                                                                                >  1:1074774/1‑93 (MQ=255)

          gtgtGGCATTCGCACCGTGTGCTACAACTTCATGCCGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGc                                                                               >  2:4678713/1‑93 (MQ=255)

            gtgGCATTCGCACCGTGTGCTACAACTTCATGCCGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGctc                                                                             >  2:811915/1‑93 (MQ=255)

            gtgGCATTCGCACCGTGTGCTACAACTTCATGCCGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGctc                                                                             >  1:5418391/1‑93 (MQ=255)

             tgGCATTCGCACCGTGTGCTACAACTTCATGCCGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGctct                                                                            <  1:569849/93‑1 (MQ=255)

              ggCATTCGCACCGTGTGCTACAACTTCATGCCGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTg                                                                           <  2:2870992/93‑1 (MQ=255)

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

              ggCATTCGCACCGTGTGCTACAACTTCATGCCGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTg                                                                           <  2:3984694/93‑1 (MQ=255)

               gCATTCGCACCGTGTGCTACAACTTCATGCCGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTgc                                                                          <  2:5801205/93‑1 (MQ=255)

               gCATTCGCACCGTGTGCTACAACTTCATGCCGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTgc                                                                          <  1:4793403/93‑1 (MQ=255)

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

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

                  ttCGCACCGTGTGCTACAACTTCATGCCGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCt                                                                       <  2:3899068/93‑1 (MQ=255)

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

                      cACCGTGTGCTACAACTTCATGCCGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGa                                                                   <  1:950722/93‑1 (MQ=255)

                      cACCGTGTGCTACAACTTCATGCCGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGa                                                                   <  1:4851054/93‑1 (MQ=255)

                       aCCGTGTGCTACAACTTCATGCCGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGAc                                                                  <  2:6079867/93‑1 (MQ=255)

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

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

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

                           tgtgCTACAACTTCATGCCGGTGCTCGACTGGACCCGTACTGTCCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGa                                                              >  2:3559806/1‑93 (MQ=255)

                            gtgCTACAACTTCATGCCGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGAt                                                             <  2:3804617/93‑1 (MQ=255)

                            gtgCTACAACTTCATGCCGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGAt                                                             <  2:5279149/93‑1 (MQ=255)

                             tgCTACAACTTCATGCCTGTGCTCGACTGGACCCGTACTGACCTCGAATACGCGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATc                                                            >  1:617639/1‑93 (MQ=255)

                             tgCTACAACTTCATGCCGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATc                                                            >  1:3070627/1‑93 (MQ=255)

                             tgCTACAACTTCATGCCGGTGCCCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATc                                                            <  1:4157573/93‑1 (MQ=255)

                             ggCTTCAACTTCATGCCGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATc                                                            <  2:3056398/92‑1 (MQ=255)

                              gCTACAACTTCATGCCGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCg                                                           >  1:2628985/1‑93 (MQ=255)

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

                                  cAACTTCATGCCGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAAtt                                                       >  1:6604150/1‑93 (MQ=255)

                                   aaCTTCATGCCGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTc                                                      <  2:1444013/93‑1 (MQ=255)

                                   aaCTTCATGCCGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTc                                                      <  2:2121013/93‑1 (MQ=255)

                                   aaCTTCATGCCGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTc                                                      <  2:4194388/93‑1 (MQ=255)

                                   aaCTTCATGCCGGTGCCCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTc                                                      >  1:470076/1‑93 (MQ=255)

                                    aCTTCATGCCGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCg                                                     >  2:6535730/1‑93 (MQ=255)

                                     cTTCATGCCGGTGCTCGACTGGGCCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGc                                                    <  2:3426980/93‑1 (MQ=255)

                                     cTTCATGCCGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCTAAGCTCTGCGCTTCGACCAGCTCGAATTCGc                                                    >  1:6553944/1‑93 (MQ=255)

                                      ttCATGCCGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCt                                                   >  2:371012/1‑93 (MQ=255)

                                      ttCATGCCGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCt                                                   <  1:4866835/93‑1 (MQ=255)

                                      ttCATGCCGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCt                                                   >  1:384550/1‑93 (MQ=255)

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

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

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

                                           gCCGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCAtt                                              >  2:5163940/1‑93 (MQ=255)

                                            ccGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTc                                             <  1:5010621/93‑1 (MQ=255)

                                             cGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCg                                            >  1:2827121/1‑93 (MQ=255)

                                                tGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCGaaa                                         <  1:5904232/93‑1 (MQ=255)

                                                 gCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCGAAAt                                        <  2:1542801/93‑1 (MQ=255)

                                                  cTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCGAAATg                                       <  1:826965/93‑1 (MQ=255)

                                                  cTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCGAAATg                                       <  2:1084971/93‑1 (MQ=255)

                                                  cTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCGAAATg                                       <  2:2441639/93‑1 (MQ=255)

                                                   tCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCGAAATGc                                      >  1:1251668/1‑93 (MQ=255)

                                                   tCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCGAAATGc                                      <  2:5728704/93‑1 (MQ=255)

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

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

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

                                                      aCTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCGAAATGCata                                   <  1:5015122/93‑1 (MQ=255)

                                                       cTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCGAAATGCatat                                  >  1:391714/1‑93 (MQ=255)

                                                       cTGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCGAAATGCatat                                  >  1:5487354/1‑93 (MQ=255)

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

                                                        tGGACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCGAAATGCATATc                                 <  2:5408045/93‑1 (MQ=255)

                                                         ggACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCGAAATGCATATcc                                >  1:5296178/1‑93 (MQ=255)

                                                         ggACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCGAAATGCATATcc                                >  1:6292965/1‑93 (MQ=255)

                                                         ggACCCGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCGAAATGCATATcc                                <  2:3408486/93‑1 (MQ=255)

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

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

                                                            cccGTACTGACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCGAAATGCATATCCTGa                             >  1:2147903/1‑93 (MQ=255)

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

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

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

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

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

                                                                    gACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCGAAATGCATATCCTGAAACGCCCa                     <  1:6394541/93‑1 (MQ=255)

                                                                    gACCTCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCGAAATGCATATCCTGAAACGCCCa                     <  1:46675/93‑1 (MQ=255)

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

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

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

                                                                        tCGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCGAAATGCATATCCTGAAACGCCCAGgcg                 >  1:5859911/1‑93 (MQ=255)

                                                                         cGAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCGAAATGCATATCCTGAAACGCCCAGgcgc                >  2:6385144/1‑93 (MQ=255)

                                                                          gAATACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCGAAATGCATATCCTGAAACGCCCAGGcgcg               >  2:3369943/1‑93 (MQ=255)

                                                                            aTACGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCGAAATGCATATCCTGAAACGCCCAGGCGCGGa             <  1:1386699/93‑1 (MQ=255)

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

                                                                              aCGTGCTGCCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCGAAATGCATATCCTGAAACGCCCAGGCGCGGGAg           <  1:1194172/93‑1 (MQ=255)

                                                                              aCGTGCTGCCAGGCGCCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCGAAATGCATATCCTGAAACGCCCAGGCGCGGAAg           <  2:2743566/93‑1 (MQ=255)

                                                                                 tgctgcCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCGAAATGCATATCCTGAAACGCCCAGGCGCGGAAGCgg        >  1:1763214/1‑93 (MQ=255)

                                                                                  gctgcCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCGAAATGCATATCCTGAAACGCCCAGGCGCGGAAGCGGa       <  2:6069836/93‑1 (MQ=255)

                                                                                  gctgcCAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCGAAATGCATATCCTGAAACGCCCAGGCGCGGAAGCGGa       <  2:1574788/93‑1 (MQ=255)

                                                                                    ggcCAGGCGGCTCCAAAGCTCTGCCCTTCGACCAGATCGAATTCGCTGCATTCGAAATGCATATCCTGAAACGCCCAGGCGCGGAAGCGGAtt     <  1:100571/92‑1 (MQ=255)

                                                                                      ccAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCGAAATGCATATCCTGAAACGCCCAGGCGCGGCAGCGGATTac   >  2:1123537/1‑93 (MQ=255)

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

                                                                                       cAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCGAAATGCATATCCTGAAACGCCCAGGCGCGGAAGCGGATTacc  >  1:5581919/1‑92 (MQ=255)

                                                                                       cAGGCGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCGAAATGCATATCCTGAAACGCCCAGGCGCGGAAGCGGATTaca  >  1:1905290/1‑93 (MQ=255)

ACCTGGCGCAGTGTGGCATTCGCACCGTGTGCTACAACTTCATGCCGGTGCTCGACTGGACCCGTACTGACCTCGAATACGTGCTGCCAGACGGCTCCAAAGCTCTGCGCTTCGACCAGATCGAATTCGCTGCATTCGAAATGCATATCCTGAAACGCCCAGGCGCGGAAGCGGATTACA  >  NC_000913/4549476‑4549655

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