BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	568,823	A→G	G233G (GGA→GGG)	ybcK →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	568,823	0	A	G	99.1%	690.4 / ‑4.9	213	G233G (GGA→GGG)	ybcK	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (2/0); new base G (120/91); total (122/91)
Fisher's exact test for biased strand distribution p-value = 5.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

Only 100 of 216 total aligned reads displayed.

CAGGATGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGACCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAAGAGGGAAAGGGATAAGTGA  >  NC_000913/568733‑568912

 aGGATGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCt                                                                                        >  1:1636402/1‑93 (MQ=255)

  ggATGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCtt                                                                                       <  1:2062105/93‑1 (MQ=255)

   gATGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTc                                                                                      >  1:669811/1‑93 (MQ=255)

   gATGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTc                                                                                      >  1:2208346/1‑93 (MQ=255)

    aTGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCt                                                                                     <  1:1710139/93‑1 (MQ=255)

     tGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTg                                                                                    <  1:3909696/93‑1 (MQ=255)

      ggAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGt                                                                                   >  1:803604/1‑93 (MQ=255)

      ggAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGt                                                                                   >  2:1171377/1‑93 (MQ=255)

        aaaGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTaa                                                                                 <  1:3400228/93‑1 (MQ=255)

         aaGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAAt                                                                                <  1:4255238/93‑1 (MQ=255)

         aaGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAAt                                                                                >  1:952238/1‑93 (MQ=255)

          aGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAAtt                                                                               >  2:4502506/1‑93 (MQ=255)

           ggCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTg                                                                              <  2:3896064/93‑1 (MQ=255)

           ggCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTg                                                                              >  1:428967/1‑93 (MQ=255)

             cgcTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGaa                                                                            <  1:27050/93‑1 (MQ=255)

              gcTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGaaa                                                                           >  2:1851341/1‑93 (MQ=255)

               cTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGaaaa                                                                          <  1:1467275/93‑1 (MQ=255)

                tCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGaaaaa                                                                         >  1:2613002/1‑93 (MQ=255)

                  aTTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAatt                                                                       <  1:2355668/93‑1 (MQ=255)

                  aTTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAatt                                                                       <  1:1941499/93‑1 (MQ=255)

                   ttGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAatta                                                                      >  1:2747255/1‑93 (MQ=255)

                   ttGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAatta                                                                      <  2:3663542/93‑1 (MQ=255)

                     gAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAattatt                                                                    >  1:3632969/1‑93 (MQ=255)

                     gAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAattatt                                                                    <  2:1446421/93‑1 (MQ=255)

                     gAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAattatt                                                                    >  1:3153111/1‑93 (MQ=255)

                      aaTGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAAttatta                                                                   <  1:108270/93‑1 (MQ=255)

                      aaTGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAAttatta                                                                   <  1:1199801/93‑1 (MQ=255)

                      aaTGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAAttatta                                                                   >  1:2596594/1‑93 (MQ=255)

                      aaTGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAAttatta                                                                   <  1:3120091/93‑1 (MQ=255)

                      aaTGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAAttatta                                                                   <  1:3944681/93‑1 (MQ=255)

                      aaTGAAATCGCCAAGTTTTTTAATGATCATACTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAAttatta                                                                   <  2:3970930/93‑1 (MQ=255)

                       aTGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAg                                                                  >  2:3220264/1‑93 (MQ=255)

                        tGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGc                                                                 <  2:1844437/93‑1 (MQ=255)

                          cAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGa                                                               >  1:484039/1‑93 (MQ=255)

                            aTAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAAt                                                             >  1:4419662/1‑93 (MQ=255)

                             tAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATa                                                            <  2:2937107/93‑1 (MQ=255)

                             tAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATa                                                            >  2:2511300/1‑93 (MQ=255)

                             tAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATa                                                            <  2:2695/93‑1 (MQ=255)

                                 cAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGACCTTCTGTAATTGAAAAATTATTAGCGAATAAAGc                                                        >  2:4031920/1‑93 (MQ=255)

                                  aaGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAACGct                                                       >  2:4359614/1‑93 (MQ=255)

                                  aaGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGct                                                       >  1:3103497/1‑93 (MQ=255)

                                  aaGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGct                                                       >  1:1977772/1‑93 (MQ=255)

                                   aGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGctc                                                      <  2:2002284/93‑1 (MQ=255)

                                   aGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGctc                                                      <  1:3251019/93‑1 (MQ=255)

                                   aGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGctc                                                      <  1:3032544/93‑1 (MQ=255)

                                    gTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGctct                                                     >  1:2253374/1‑93 (MQ=255)

                                    gTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGctct                                                     >  2:732099/1‑93 (MQ=255)

                                     tATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTCTTAGCGAATAAAGCTCTg                                                    >  1:2732849/1‑93 (MQ=255)

                                     tATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTg                                                    >  2:2264642/1‑93 (MQ=255)

                                      aTTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGa                                                   >  1:1014454/1‑93 (MQ=255)

                                       tttAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGAt                                                  >  2:1911701/1‑93 (MQ=255)

                                        ttAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATa                                                 <  2:4602904/93‑1 (MQ=255)

                                         tAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAg                                                >  2:397976/1‑93 (MQ=255)

                                            aTGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTa                                             <  2:43797/93‑1 (MQ=255)

                                            aTGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTa                                             <  2:3314993/93‑1 (MQ=255)

                                            aTGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGAGCATGGGGGCATTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTa                                             >  1:2348277/1‑93 (MQ=255)

                                             tGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTAt                                            >  2:1968178/1‑93 (MQ=255)

                                               aTCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTAttt                                          >  1:4221069/1‑93 (MQ=255)

                                                 cATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGc                                        >  2:638002/1‑93 (MQ=255)

                                                  aTGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCg                                       >  1:1254120/1‑93 (MQ=255)

                                                  aTGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCg                                       >  1:4261968/1‑93 (MQ=255)

                                                   tGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGt                                      >  2:696386/1‑93 (MQ=255)

                                                      tGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTAcc                                   >  2:2659005/1‑93 (MQ=255)

                                                       gTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCt                                  <  1:539869/93‑1 (MQ=255)

                                                       gTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCt                                  >  1:4481065/1‑93 (MQ=255)

                                                        tAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCtt                                 >  2:3004313/1‑93 (MQ=255)

                                                          aaGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTct                               >  1:4454850/1‑92 (MQ=255)

                                                          aaGAATTTCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCa                               >  2:4289983/1‑93 (MQ=255)

                                                                ttCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGt                         >  1:1414029/1‑93 (MQ=255)

                                                                ttCTCAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGt                         <  2:2021386/93‑1 (MQ=255)

                                                                  ctcAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGc                       >  2:1245842/1‑93 (MQ=255)

                                                                   tcAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCa                      >  2:3031475/1‑93 (MQ=255)

                                                                    cAGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCaa                     >  2:4334199/1‑93 (MQ=255)

                                                                     aGGAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAag                    >  1:1808739/1‑93 (MQ=255)

                                                                      ggAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAaga                   >  2:2387554/1‑93 (MQ=255)

                                                                      ggAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAaga                   >  1:148118/1‑93 (MQ=255)

                                                                      ggAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAaga                   <  1:2724464/93‑1 (MQ=255)

                                                                       gAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAagag                  <  1:2883380/93‑1 (MQ=255)

                                                                       gAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAagag                  <  1:3166613/93‑1 (MQ=255)

                                                                       gAAAAGAAAGTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAagag                  <  1:3444980/93‑1 (MQ=255)

                                                                         aaagaaagTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAAGAggg                <  2:1622888/93‑1 (MQ=255)

                                                                         aaagaaagTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAAGAggg                >  1:3564207/1‑93 (MQ=255)

                                                                          aagaaagTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAAGAGGGa               <  1:3193192/93‑1 (MQ=255)

                                                                           agaaagTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAAGAGGGaa              >  2:2843794/1‑93 (MQ=255)

                                                                           agaaagTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAAGAGGGaa              >  2:112242/1‑93 (MQ=255)

                                                                            gaaagTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAAGAGGGaaa             >  2:1672179/1‑93 (MQ=255)

                                                                            gaaagTGCAGGGGGGCCTTCTGTAATGAAAAAATTAATAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAAGAGGGaaa             <  2:486938/93‑1 (MQ=255)

                                                                              aagTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAAGAGGGAAAgg           >  1:4737651/1‑93 (MQ=255)

                                                                               agTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAAGAGGGAAAggg          <  1:3715723/93‑1 (MQ=255)

                                                                               agTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAAGAGGGAAAggg          <  2:1544069/93‑1 (MQ=255)

                                                                                gTGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAAGAGGGAAAGGGa         <  2:857117/93‑1 (MQ=255)

                                                                                 tGCATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAAGAGGGAAAGGGAt        <  2:3770446/93‑1 (MQ=255)

                                                                                   cATGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAAGAGGGAAAGGGATaa      <  2:1080923/93‑1 (MQ=255)

                                                                                     tGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAAGAGGGAAAGGGATAAGt    >  2:3274329/1‑93 (MQ=255)

                                                                                     tGGGGGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAAGAGGGAAAGGGATAAGt    >  2:476883/1‑93 (MQ=255)

                                                                                      ggggGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAAGAGGGAAAGGGATAAGTg   >  2:2571676/1‑93 (MQ=255)

                                                                                      ggggGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAAGAGGGAAAGGGATAAGTg   >  2:2138109/1‑93 (MQ=255)

                                                                                      ggggGCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAAGAGGGAAAGGGATAAGTg   >  1:2175493/1‑93 (MQ=255)

                                                                                       ggggCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAAGAGGGAAAGGGATAAGTGa  <  2:1819245/93‑1 (MQ=255)

                                                                                       ggggCCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAAGAGGGAAAGGGATAAGTGa  <  1:1073701/93‑1 (MQ=255)

CAGGATGGAAAGGCGCTCATTGAATGCAATAGCCAAGTATTTAAATGATCATGCTGTAAAGAATTTCTCAGGAAAAGAAAGTGCATGGGGACCTTCTGTAATTGAAAAATTATTAGCGAATAAAGCTCTGATAGGTATTTGCGTACCTTCATATCGTGCAAGAGGGAAAGGGATAAGTGA  >  NC_000913/568733‑568912

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