BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,006,483	+C	coding (183/207 nt)	b4554 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,006,480	1	.	C	98.9%	602.6 / ‑0.7	176	coding (186/207 nt)	b4554	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (1/1); new base C (96/78); total (97/79)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.36e-01

Only 100 of 199 total aligned reads displayed.

GGCGCGCCACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGT‑CCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGG  >  NC_000913/4006388‑4006567

 gcgcgcCACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGT‑c                                                                                        <  1:2159984/93‑1 (MQ=255)

 gcgcgcCACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGT‑c                                                                                        >  1:2838416/1‑93 (MQ=255)

  cgcgcCACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGT‑cc                                                                                       >  1:1686913/1‑93 (MQ=255)

  cgcgcCACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGT‑cc                                                                                       >  1:4028546/1‑93 (MQ=255)

    cgcCACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGT‑cccc                                                                                     <  1:4949246/93‑2 (MQ=255)

    cgcCACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGT‑cccc                                                                                     >  1:97979/1‑92 (MQ=255)

     gcCACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGT‑cccca                                                                                    >  1:79009/1‑91 (MQ=255)

     gcCACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGT‑cccca                                                                                    >  1:3495517/1‑91 (MQ=255)

      ccACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGT‑ccccaa                                                                                   <  1:5509350/93‑4 (MQ=255)

      ccACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGT‑ccccaa                                                                                   <  1:6435612/93‑4 (MQ=255)

      ccACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGT‑ccccaa                                                                                   <  2:1390017/93‑4 (MQ=255)

       cACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGT‑ccccaac                                                                                  >  2:1345327/1‑89 (MQ=255)

         ctgctgGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGT‑CCCCAaccc                                                                                >  2:3030605/1‑93 (MQ=255)

          tgctgGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAAcccc                                                                                >  1:5928209/1‑93 (MQ=255)

          tgctgGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAAcccc                                                                                <  1:5121029/93‑1 (MQ=255)

          tgctgGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAAcccc                                                                                >  1:877006/1‑93 (MQ=255)

            ctgGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCAc                                                                              <  2:5319759/93‑1 (MQ=255)

              ggATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGc                                                                            >  1:817692/1‑93 (MQ=255)

                 ttCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCaa                                                                         >  1:222744/1‑93 (MQ=255)

                 ttCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTAAACACTTCAAACGCAATCACTGAGCGGGAAACTAGCGTCCCCAACCCCACGCCaa                                                                         >  1:2046955/1‑93 (MQ=255)

                  tCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAAc                                                                        >  1:1807985/1‑93 (MQ=255)

                    ccAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCa                                                                      >  1:56857/1‑93 (MQ=255)

                    acAGACACATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAACCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCa                                                                      <  1:138139/92‑1 (MQ=255)

                       gCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGca                                                                   >  1:744533/1‑93 (MQ=255)

                         cAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGcaca                                                                 >  2:4045314/1‑93 (MQ=255)

                         cAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGcaca                                                                 <  2:2611105/93‑1 (MQ=255)

                           aaTCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAat                                                               >  2:5513728/1‑93 (MQ=255)

                           aaTCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAat                                                               >  1:5209210/1‑93 (MQ=255)

                           aaTCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAat                                                               >  2:505136/1‑93 (MQ=255)

                           aaTCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAat                                                               >  2:1995968/1‑93 (MQ=255)

                             acAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAatat                                                             <  2:4928060/92‑1 (MQ=255)

                              cAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATg                                                            >  1:1657788/1‑93 (MQ=255)

                               aaGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGa                                                           >  1:62909/1‑93 (MQ=255)

                                 gTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAAt                                                         >  2:5957778/1‑93 (MQ=255)

                                 gTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAAt                                                         <  2:752402/93‑1 (MQ=255)

                                 gTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAAt                                                         <  1:554018/93‑1 (MQ=255)

                                   aaGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAACATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCg                                                       <  2:3430529/93‑1 (MQ=255)

                                    aGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGc                                                      <  1:6045384/93‑1 (MQ=255)

                                     gCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGcc                                                     >  2:4831328/1‑93 (MQ=255)

                                     gCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGcc                                                     <  2:2722508/93‑1 (MQ=255)

                                     gCCGCGCCGGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGcc                                                     <  1:4789480/93‑1 (MQ=255)

                                      ccGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCa                                                    <  2:146205/93‑1 (MQ=255)

                                       cgcgCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAg                                                   <  2:2564168/93‑1 (MQ=255)

                                       cgcgCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAg                                                   <  1:3698815/93‑1 (MQ=255)

                                        gcgcCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGt                                                  >  2:4649266/1‑93 (MQ=255)

                                         cgcCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTc                                                 >  1:6128761/1‑93 (MQ=255)

                                         cgcCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTc                                                 <  1:258031/93‑1 (MQ=255)

                                          gcCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTcc                                                >  2:4224681/1‑93 (MQ=255)

                                          gcCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTcc                                                <  2:1500749/93‑1 (MQ=255)

                                          gcCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTcc                                                <  1:1717470/93‑1 (MQ=255)

                                           ccTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCg                                               >  1:3016265/1‑93 (MQ=255)

                                           ccTGCCCACATCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCg                                               <  1:4543259/93‑1 (MQ=255)

                                            cTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCgg                                              <  2:5824035/93‑1 (MQ=255)

                                              gCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTc                                            >  2:5554375/1‑93 (MQ=255)

                                              gCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTc                                            <  1:1839906/93‑1 (MQ=255)

                                               cccACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCt                                           >  2:3906010/1‑93 (MQ=255)

                                                ccACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTg                                          >  1:5102442/1‑93 (MQ=255)

                                                ccACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTg                                          >  2:4162700/1‑93 (MQ=255)

                                                ccACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTg                                          >  2:5779254/1‑93 (MQ=255)

                                                ccACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTg                                          >  1:3001175/1‑93 (MQ=255)

                                                 cacttcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGa                                         <  2:44767/93‑1 (MQ=255)

                                                 cacttcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGa                                         <  1:263193/93‑1 (MQ=255)

                                                  acttcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGaa                                        >  2:1756140/1‑93 (MQ=255)

                                                  acttcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGaa                                        >  1:5478617/1‑93 (MQ=255)

                                                    ttcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGc                                      <  2:4965603/93‑1 (MQ=255)

                                                    ttcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGc                                      <  1:3187289/93‑1 (MQ=255)

                                                     tcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGcc                                     >  1:4125323/1‑93 (MQ=255)

                                                       aacacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCa                                   >  2:2178057/1‑93 (MQ=255)

                                                          acttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCa                                >  2:1039101/1‑93 (MQ=255)

                                                           cttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCaa                               <  1:3832896/93‑1 (MQ=255)

                                                            ttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAAt                              <  1:1457053/93‑1 (MQ=255)

                                                             tcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATa                             <  2:6099150/93‑1 (MQ=255)

                                                               aaaCGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCCACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGa                           >  1:2697074/1‑93 (MQ=255)

                                                               aaaCGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGa                           <  1:4980078/93‑1 (MQ=255)

                                                                aaCGCAATCACTGGGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGAc                          <  1:1140089/93‑1 (MQ=255)

                                                                aaCGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGAc                          >  2:5136330/1‑93 (MQ=255)

                                                                  cGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGc                        <  2:3377785/93‑1 (MQ=255)

                                                                   gCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGcc                       <  2:4412685/93‑1 (MQ=255)

                                                                   gCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGcc                       <  2:248600/93‑1 (MQ=255)

                                                                   gCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGcc                       >  2:1898973/1‑93 (MQ=255)

                                                                   gCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGcc                       <  2:5268269/93‑1 (MQ=255)

                                                                    cAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGAAATAGACGCCa                      >  1:2122214/1‑93 (MQ=255)

                                                                     aaTCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAc                     >  2:795386/1‑93 (MQ=255)

                                                                      aTCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAcc                    <  1:108688/93‑1 (MQ=255)

                                                                      aTCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAcc                    >  1:3570491/1‑93 (MQ=255)

                                                                       tCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAccg                   >  2:2297447/1‑93 (MQ=255)

                                                                       tCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAccg                   <  1:5339073/93‑1 (MQ=255)

                                                                          cTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAccgcgc                <  2:1254356/93‑1 (MQ=255)

                                                                           tGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAccgcgcc               >  2:2073059/1‑93 (MQ=255)

                                                                           tGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAccgcgcc               >  1:409790/1‑93 (MQ=255)

                                                                            gAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAccgcgccg              >  2:3180874/1‑93 (MQ=255)

                                                                            gAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAccgcgccg              <  2:3996415/93‑1 (MQ=255)

                                                                            gAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAccgcgccg              >  2:1881303/1‑93 (MQ=255)

                                                                            gAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAccgcgccg              >  2:4422080/1‑93 (MQ=255)

                                                                             aGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAccgcgccgc             >  2:639850/1‑93 (MQ=255)

                                                                              gCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAccgcgccgcg            >  1:530254/1‑93 (MQ=255)

                                                                                gggAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGAt          <  1:3392494/93‑1 (MQ=255)

                                                                                gggAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGAt          >  2:930210/1‑93 (MQ=255)

                                                                                 ggAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATa         >  2:2902163/1‑93 (MQ=255)

                                                                                    aaaTAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAAcc      >  2:3322019/1‑93 (MQ=255)

GGCGCGCCACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGT‑CCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGG  >  NC_000913/4006388‑4006567

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