BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,935,895	A→G	T15A (ACC→GCC)	rbsR →	RbsR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,935,895	0	A	G	99.2%	828.4 / ‑3.3	238	T15A (ACC→GCC)	rbsR	RbsR
Reads supporting (aligned to +/- strand): ref base A (2/0); new base G (129/107); total (131/107)
Fisher's exact test for biased strand distribution p-value = 5.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

Only 100 of 245 total aligned reads displayed.

CCGTGGCGTGAAGAGATCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAAT  >  NC_000913/3935805‑3935984

ccGTGGCGTGAAGAGATCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGcc                                                                                         >  1:1890189/1‑93 (MQ=255)

  gTGGCGTGAAGAGATCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTc                                                                                       >  2:4587824/1‑93 (MQ=255)

  gTGGCGTGAAGAGATCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTc                                                                                       >  2:4868169/1‑93 (MQ=255)

  gTGGCGTGAAGAGATCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTc                                                                                       <  1:4680936/93‑1 (MQ=255)

   gggCGTGAAGAGATCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCa                                                                                      <  1:4393155/92‑1 (MQ=255)

    ggCGTGAAGAGATCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCaa                                                                                     >  1:4586648/1‑93 (MQ=255)

     gCGTGAAGAGATCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAAc                                                                                    >  2:4082132/1‑93 (MQ=255)

     gCGTGAAGAGATCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAAc                                                                                    <  1:6151103/93‑1 (MQ=255)

      cGTGAAGAGATCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACa                                                                                   <  1:515927/93‑1 (MQ=255)

       gTGAAGAGATCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGTCTCAACAg                                                                                  >  2:617562/1‑93 (MQ=255)

       gTGAAGAGATCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAg                                                                                  <  2:6080664/93‑1 (MQ=255)

       gTGAAGAGATCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGGTTCTGCCTCAAccg                                                                                  >  2:5418726/1‑91 (MQ=255)

        tGAAGAGATCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGt                                                                                 <  2:3520653/93‑1 (MQ=255)

          aaGAGATCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGttt                                                                               >  2:2041569/1‑93 (MQ=255)

           agagATCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTtc                                                                              >  2:3925276/1‑93 (MQ=255)

            gagaTCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTtct                                                                             >  1:983261/1‑93 (MQ=255)

             agaTCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTtctc                                                                            <  1:2575502/93‑1 (MQ=255)

              gaTCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCa                                                                           <  2:452497/93‑1 (MQ=255)

              gaTCGACGCATTTTGAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCa                                                                           <  2:3181500/93‑1 (MQ=255)

               aTCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCAc                                                                          >  1:3510046/1‑93 (MQ=255)

                tCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACg                                                                         >  1:5613303/1‑93 (MQ=255)

                 cGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGt                                                                        <  2:710244/93‑1 (MQ=255)

                 cGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGt                                                                        <  1:5098102/93‑1 (MQ=255)

                  gACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGtt                                                                       <  2:2309483/93‑1 (MQ=255)

                   aCGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTa                                                                      <  2:2761277/93‑1 (MQ=255)

                   aCGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTa                                                                      >  2:6253240/1‑93 (MQ=255)

                     gCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATc                                                                    >  1:3092611/1‑93 (MQ=255)

                      cATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCa                                                                   >  1:3120043/1‑93 (MQ=255)

                       aTTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCaa                                                                  >  2:1294627/1‑93 (MQ=255)

                        tttttAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAAt                                                                 >  2:2453339/1‑93 (MQ=255)

                         ttttAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATa                                                                <  1:2915409/93‑1 (MQ=255)

                           ttAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCCCGTTATCAATaaa                                                              <  1:6356850/93‑1 (MQ=255)

                           ttAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATaaa                                                              <  2:2752519/93‑1 (MQ=255)

                             aGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGa                                                            >  2:2926436/1‑93 (MQ=255)

                             aGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGa                                                            <  1:5740004/93‑1 (MQ=255)

                              gACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGAt                                                           <  1:3345009/93‑1 (MQ=255)

                              gACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGAt                                                           >  2:4965492/1‑93 (MQ=255)

                              gACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGAt                                                           >  2:5712542/1‑93 (MQ=255)

                              gACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGAt                                                           >  1:2412669/1‑93 (MQ=255)

                              gACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGAt                                                           <  1:5628053/93‑1 (MQ=255)

                              gACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGAt                                                           <  2:3143500/93‑1 (MQ=255)

                                 aGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGc                                                        <  1:1634002/93‑1 (MQ=255)

                                  ggCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCt                                                       <  2:5533423/93‑1 (MQ=255)

                                  ggCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCt                                                       <  2:4563716/93‑1 (MQ=255)

                                   gCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCtt                                                      >  2:3852784/1‑93 (MQ=255)

                                      gagGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGt                                                   >  1:4071246/1‑93 (MQ=255)

                                       agGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTc                                                  <  1:1705185/93‑1 (MQ=255)

                                       agGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTc                                                  >  1:1905574/1‑93 (MQ=255)

                                       agGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTc                                                  <  2:4268797/93‑1 (MQ=255)

                                       agGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTc                                                  <  2:3726827/93‑1 (MQ=255)

                                       agGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTc                                                  <  1:711597/93‑1 (MQ=255)

                                        ggTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCa                                                 <  1:2291199/93‑1 (MQ=255)

                                          tGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGt                                               <  2:395664/93‑1 (MQ=255)

                                           gACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTg                                              <  2:5007498/93‑1 (MQ=255)

                                            aCGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGa                                             >  1:2890428/1‑93 (MQ=255)

                                             cGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGaa                                            >  1:3851115/1‑93 (MQ=255)

                                             cGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGaa                                            >  1:3601365/1‑93 (MQ=255)

                                             cGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGaa                                            <  2:2483880/93‑1 (MQ=255)

                                              gCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAg                                           >  1:4286879/1‑93 (MQ=255)

                                              gCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAg                                           <  1:1434714/93‑1 (MQ=255)

                                                ttGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCg                                         >  2:5615454/1‑93 (MQ=255)

                                                  ggCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGAt                                       <  2:5170253/93‑1 (MQ=255)

                                                    cTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTa                                     <  1:4418618/93‑1 (MQ=255)

                                                     tACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTAc                                    >  1:3047211/1‑93 (MQ=255)

                                                      aCAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTAcc                                   <  1:5802389/93‑1 (MQ=255)

                                                      aCAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTAcc                                   >  1:1101221/1‑93 (MQ=255)

                                                        aaTGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGc                                 >  2:1078768/1‑93 (MQ=255)

                                                         aTGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGcc                                >  1:1280793/1‑93 (MQ=255)

                                                         aTGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGcc                                >  1:5245678/1‑93 (MQ=255)

                                                         aTGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGcc                                >  1:2242585/1‑93 (MQ=255)

                                                         aTGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGcc                                <  1:2072064/93‑1 (MQ=255)

                                                           gAAAGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCaa                              <  2:1959150/93‑1 (MQ=255)

                                                            aaaGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCaaa                             <  1:4968054/93‑1 (MQ=255)

                                                             aaGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAg                            >  2:1368627/1‑93 (MQ=255)

                                                             aaGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAg                            >  1:941548/1‑93 (MQ=255)

                                                             aaGATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAg                            <  1:696991/93‑1 (MQ=255)

                                                               gATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGtt                          >  2:1870539/1‑93 (MQ=255)

                                                               gATGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGtt                          >  1:4000706/1‑93 (MQ=255)

                                                                 tGTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGa                        <  1:1321422/93‑1 (MQ=255)

                                                                  gTTGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGaa                       <  2:1240106/93‑1 (MQ=255)

                                                                    tGCCCGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAgc                     >  2:533634/1‑93 (MQ=255)

                                                                       ccGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAgcggc                  <  1:6029894/93‑1 (MQ=255)

                                                                        cGCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAgcggcg                 >  2:1609246/1‑93 (MQ=255)

                                                                         gCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGa                >  2:3641025/1‑93 (MQ=255)

                                                                         gCCTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGa                >  2:6257368/1‑93 (MQ=255)

                                                                         gCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGa                >  2:792764/1‑93 (MQ=255)

                                                                          ccTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGAt               >  1:3694317/1‑93 (MQ=255)

                                                                           cTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGAtt              >  1:5456091/1‑93 (MQ=255)

                                                                           cTGGCGGGCGTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGAtt              <  2:5796436/93‑1 (MQ=255)

                                                                             ggcgggcgTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTaa            >  2:4827562/1‑93 (MQ=255)

                                                                               cgggcgTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAATGTTGAAGCGGCGATTAAAg          <  1:6319058/93‑1 (MQ=255)

                                                                               cgggcgTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGGTGAAGCGGCGATTAAAg          >  1:4172280/1‑93 (MQ=255)

                                                                                gggcgTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGa         <  1:532916/93‑1 (MQ=255)

                                                                                 ggcgTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGaa        >  2:2228054/1‑93 (MQ=255)

                                                                                    gTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACta     >  1:3764211/1‑92 (MQ=255)

                                                                                    gTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTc     >  1:742354/1‑93 (MQ=255)

                                                                                    gTTTCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGAGTAAAGAACTc     >  1:5873766/1‑93 (MQ=255)

                                                                                      ttCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCaa   >  2:1551886/1‑93 (MQ=255)

                                                                                       tCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAAt  <  1:4550988/93‑1 (MQ=255)

                                                                                       tCTGCCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAAt  >  1:1587379/1‑93 (MQ=255)

CCGTGGCGTGAAGAGATCGACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAAT  >  NC_000913/3935805‑3935984

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