BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,369,183	A→G	Q194R (CAG→CGG)	mopA →	GroEL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,369,183	0	A	G	98.1%	482.0 / ‑0.7	155	Q194R (CAG→CGG)	mopA	GroEL
Reads supporting (aligned to +/- strand): ref base A (2/1); new base G (67/85); total (69/86)
Fisher's exact test for biased strand distribution p-value = 5.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.89e-01

Only 100 of 163 total aligned reads displayed.

AAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCAGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGAAACTGGCGCAGTAGAACTGGAAAGCCCGTTCATCCTGCTGGCTGA  >  NC_000913/4369093‑4369273

aaGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCGGt                                                                                          <  2:2694831/93‑1 (MQ=255)

  gCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCGGTTc                                                                                        <  1:710085/93‑1 (MQ=255)

  gCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCGGTTc                                                                                        <  1:4488941/93‑1 (MQ=255)

   cGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCGGTTCg                                                                                       <  1:2894177/93‑1 (MQ=255)

        gACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCGGTTCGACCGt                                                                                  >  1:2625679/1‑93 (MQ=255)

        gACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCGGTTCGACCGt                                                                                  >  1:2001473/1‑93 (MQ=255)

        gACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCGGTTCGACCGt                                                                                  <  1:1529291/93‑1 (MQ=255)

         aCAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCGGTTCGACCGTg                                                                                 <  1:4557933/93‑1 (MQ=255)

         aCAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCGGTTCGACCGTg                                                                                 >  1:4333348/1‑93 (MQ=255)

         aCAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCGGTTCGACCGTg                                                                                 >  1:3712709/1‑93 (MQ=255)

          cAAAGTCGTTAAAGAAGGCGTTATCACCGTTGTAGACCGTACCGGTCTGCAGGACGAACTGGACGGGGGTGAAGGTATGCAGTTCGACCGTgg                                                                                >  2:2546573/1‑93 (MQ=255)

           aaaGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCGGTTCGACCGTGGc                                                                               >  1:1223731/1‑93 (MQ=255)

           aaaGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCGGTTCGACCGTGGc                                                                               <  1:3553042/93‑1 (MQ=255)

             aGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCGGTTCGACCGTGGCTa                                                                             <  2:4120029/93‑1 (MQ=255)

             aGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCGGTTCGACCGTGGCTa                                                                             >  2:3866866/1‑93 (MQ=255)

             aGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCGGTTCGACCGTGGCTa                                                                             >  2:2915686/1‑93 (MQ=255)

              gTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCGGTTCGACCGTGGCTAc                                                                            >  1:724088/1‑93 (MQ=255)

              gTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCAGTTCGACCGTGGCTAc                                                                            >  1:4023274/1‑93 (MQ=255)

               tCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCGGTTCGACCGTGGCTAcc                                                                           <  1:2350464/93‑1 (MQ=255)

               gcGGTAGAGAAGGCGTTATCACCGTTGAAGCCGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCGGTTCGACCGTGGCTAcc                                                                           <  1:2403596/92‑1 (MQ=255)

                 ggTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCGGTTCGACCGTGGCTACCTg                                                                         <  1:342483/93‑1 (MQ=255)

                  gTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCGGTTCGACCGTGGCTACCTGt                                                                        >  2:172922/1‑93 (MQ=255)

                   tAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCGGTTCGACCGTGGCTACCTGtc                                                                       <  2:1226301/93‑1 (MQ=255)

                    aaaGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCGGTTCGACCGTGGCTACATGtct                                                                      >  1:1378834/1‑93 (MQ=255)

                     aagaagGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCGGTTCGACCGTGGCTACCTGtctc                                                                     <  1:1478569/93‑1 (MQ=255)

                      agaagGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCGGTTCGACCGTGGCTACCTGTCTcc                                                                    <  1:2532483/93‑1 (MQ=255)

                      agaagGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCGGTTCGACCGTGGCTACCTGTCTcc                                                                    >  2:277246/1‑93 (MQ=255)

                      agaagGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCGGTTCGACCGTGGCTACCTGTCTcc                                                                    >  1:3140832/1‑93 (MQ=255)

                      agaagGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCGGTTCGACCGTGGCTACCTGTCTcc                                                                    <  2:47938/93‑1 (MQ=255)

                        aagGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCGGTTCGACCGTGGCTACCTGTCTCCtt                                                                  <  2:1908474/93‑1 (MQ=255)

                        aagGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCGGTTCGACCGTGGCTACCTGTCTCCtt                                                                  <  2:3412264/93‑1 (MQ=255)

                          ggCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCGGTTCGACCGTGGCTACCTGTCTCCTTAc                                                                <  1:287791/93‑1 (MQ=255)

                          ggCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCGGTTCGACCGTGGCTACCTGTCTCCTTAc                                                                <  2:3744768/93‑1 (MQ=255)

                          ggCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCGGTTCGACCGTGGCTACCTGTCTCCTTAc                                                                <  2:3818985/93‑1 (MQ=255)

                           gCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCGGTTCGACCGTGGCTAACTGTCTCCTTACt                                                               >  2:1556373/1‑93 (MQ=255)

                             gTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCGGTTCGACCGTGGCTACCTTTCTCCTTACtta                                                             >  2:623295/1‑92 (MQ=255)

                              ttATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCGGTTCGACCGTGGCTACCTGTCTCCTTACTtca                                                            <  2:2354096/93‑1 (MQ=255)

                              ttATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCGGTTCGACCGTGGCTACCTGTCTCCTTACTtca                                                            >  2:2283505/1‑93 (MQ=255)

                                 tCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCGGTTCGACCGTGGCTACCTGTCTCCTTACTtcatca                                                         <  2:24095/93‑1 (MQ=255)

                                 tCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCGGTTCGACCGTGGCTACCTGTCTCCTTACTtcatca                                                         >  1:251723/1‑93 (MQ=255)

                                  cACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCGGTTCGACCGTGGCTACCTGTCTCCTTACTTCATcaa                                                        <  1:3775653/93‑1 (MQ=255)

                                   aCCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCGGTTCGACCCTGGCTCCCTGTCTCCTTACTtcttcacc                                                       >  1:914834/1‑91 (MQ=255)

                                    ccGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCGGTTCGACCGTGGCTACCTGTCTCCTTACTTCATcaaca                                                      >  1:463511/1‑93 (MQ=255)

                                    ccGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCGGTTCGACCGTGGCTACCTGTCTCCTTACTTCATcaaca                                                      >  2:1475099/1‑93 (MQ=255)

                                     cGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCGGTTCGACCGTGGCTACCTGTCTCCTTACTTCATcaacaa                                                     <  1:2295418/93‑1 (MQ=255)

                                     cGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCGGTTCGACCGTGGCTACCTGTCTCCTTACTTCATcaacaa                                                     <  1:2038623/93‑1 (MQ=255)

                                       ttGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCGGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGc                                                   >  1:785571/1‑93 (MQ=255)

                                           aGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCGGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGa                                               >  1:2746693/1‑93 (MQ=255)

                                           aGACGGGTCCGGTCTGCAGGAATAACTGGACGTGGTTGTAGGTATGCGGTTCGACCGTGGCTACCTGTCTCCTTACGGCATCAACAAGCCGGa                                               <  2:2858838/93‑1 (MQ=255)

                                             aCGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCGGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGaaa                                             <  2:1312093/93‑1 (MQ=255)

                                             aCGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCGGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGaaa                                             <  2:694580/93‑1 (MQ=255)

                                              cGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCGGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGAAAc                                            <  2:886209/93‑1 (MQ=255)

                                              cGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCGGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGAAAc                                            <  2:702986/93‑1 (MQ=255)

                                              cGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCGGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGAAAc                                            >  1:669087/1‑93 (MQ=255)

                                                 tACCGGTCTGCAGGACGAAATGGACGTGGTTGAAGGTATGCGGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGAAACTgg                                         <  1:902848/93‑1 (MQ=255)

                                                  aCCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCGGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGAAACTGgc                                        >  2:1882091/1‑93 (MQ=255)

                                                    cGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCGGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGAAACTGgcgc                                      >  2:3965528/1‑93 (MQ=255)

                                                    cGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCGGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGAAACTGgcgc                                      <  1:389296/93‑1 (MQ=255)

                                                     ggTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCGGTTCGACCGTGGCTACCTGTCTCCTTACTTCTTCAACACGCCGGAAACTGGCCCa                                     >  2:1948342/1‑93 (MQ=255)

                                                      gTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCGGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGAAACTGGCGCAg                                    >  1:3429017/1‑93 (MQ=255)

                                                       tCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCGGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGAAACTGGCGCAGt                                   >  1:1678884/1‑93 (MQ=255)

                                                          gCAGGACGAACTGGACGTGGTTGAAGGTATGCGGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGAAACTGGCGCAGTAGa                                <  1:2497839/93‑1 (MQ=255)

                                                          gCAGGACGAACTGGACGTGGTTGAAGGTATGCGGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGAAACTGGCGCAGTAGa                                >  2:327103/1‑93 (MQ=255)

                                                           cAGGACGAACCGGACGTGGTTGAAGGTATGCGGTACGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGAAACTGGCGCAGTAGaa                               <  2:4423459/93‑1 (MQ=255)

                                                            aGGACGAACTGGACGTGGTTGAAGGTATGCGGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGAAACTGGCGCAGTAGAAc                              <  1:367587/93‑1 (MQ=255)

                                                             ggACGAACTGGACGTGGTTGAAGGTATGCGGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGAAACTGGCGCAGTAGAACt                             <  2:3972378/93‑1 (MQ=255)

                                                             ggACGAACTGGACGTGGTTGAAGGTATGCGGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGAAACTGGCGCAGTAGAACt                             >  1:2611430/1‑93 (MQ=255)

                                                              gACGAACTGGACGTGGTTGAAGGTATGCGGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGAAACTGGCGCAGTAGAACTg                            >  1:4054142/1‑93 (MQ=255)

                                                              gACGAACTGGACGTGGTTGAAGGTATGCGGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGAAACTGGCGCAGTAGAACTg                            <  2:4240920/93‑1 (MQ=255)

                                                                cGAACTGGACGTGGTTGAAGGTATGCGGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGAAACTGGCGCAGTAGAACTGGa                          <  1:1042073/93‑1 (MQ=255)

                                                                cGAACTGGACGTGGTTGAAGGTATGCGGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGAAACTGGCGCAGTAGAACTGGa                          <  2:2856531/93‑1 (MQ=255)

                                                                  aaCTGGACGTGGTTGAAGGTATGCGGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGAAACTGGCGCAGTAGAACTGGaaa                        >  2:452623/1‑93 (MQ=255)

                                                                  aaCTGGACGTGGTTGAAGGTATGCGGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGAAACTGGCGCAGTAGAACTGGaaa                        <  1:1046615/93‑1 (MQ=255)

                                                                   aCTGGACGTGGTTGAAGGTATGCGGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGAAACTGGCGCAGTAGAACTGGAAAg                       <  2:2674141/93‑1 (MQ=255)

                                                                     tGGACGTGGTTGAAGGTATGCGGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGAAACTGGCGCAGTAGAACTGGAAAGcc                     >  2:2522953/1‑93 (MQ=255)

                                                                     tGGACGTGGTTGAAGGTATGCGGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGAAACTGGCGCAGTAGAACTGGAAAGcc                     >  2:2078926/1‑93 (MQ=255)

                                                                      ggACGTGGTTGAAGGTATGCGGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGAAACTGGCGCAGTAGAACTGGAAAGccc                    <  1:2174900/93‑1 (MQ=255)

                                                                       gACGTGGTTGAAGGTATGCGGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGAAACTGGCGCAGTAGAACTGGAAAGCCCg                   >  2:4442988/1‑93 (MQ=255)

                                                                       gACGTGGTTGAAGGTATGCGGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGAAACTGGCGCAGTAGAACTGGAAAGCCCg                   >  1:1572838/1‑93 (MQ=255)

                                                                         cGTGGTTGAAGGTATGCGTTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGAAACTGGCGCAGTAGAACTGGAAAGCCCGtt                 <  1:2216422/93‑1 (MQ=255)

                                                                         cGTGGTTGAAGGTATGCGGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGAAACTGGCGCAGTAGAACTGGAAAGCCCGtt                 <  1:1869809/93‑1 (MQ=255)

                                                                          gTGGTTGAAGGTATGCGGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGAAACTGGCGCAGTAGAACTGGAAAGCCCGTTc                >  2:4692714/1‑93 (MQ=255)

                                                                           tGGTTGAAGGTATGCGGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGAAACTGGCGCAGTAGAACTGGAAAGCCCGTTCa               <  1:4724772/93‑1 (MQ=255)

                                                                            ggTTGAAGGTATGCGGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGAAACTGGCGCAGTAGAACTGGAAAGCCCGTTCAt              <  1:3364251/93‑1 (MQ=255)

                                                                             gTTGAAGGTATGCGGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGAAACTGGCGCAGTAGAACTGGAAAGCCCGTTCATc             <  1:1721653/93‑1 (MQ=255)

                                                                             gTTGAAGGTATGCGGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGAAACTGGCGCAGTAGAACTGGAAAGCCCGTTCATc             >  1:1293613/1‑93 (MQ=255)

                                                                             gTTGAAGGTATGCGGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGAAACTGGCGCAGTAGAACTGGAAAGCCCGTTCATc             <  1:3256101/93‑1 (MQ=255)

                                                                              ttGAAGGTATGCGGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGAAACTGGCGCAGTAGAACTGGAAAGCCCGTTCATcc            <  2:1166786/93‑1 (MQ=255)

                                                                              ttGAAGGTATGCGGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGAAACTGGCGCAGTAGAACTGGAAAGCCCGTTCATcc            <  2:2565771/93‑1 (MQ=255)

                                                                               tGAAGGTATGCGGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGAAACTGGCGCAGTAGAACTGGAAAGCCCGTTCATCct           <  1:2649741/93‑1 (MQ=255)

                                                                                gAAGGTATGCGGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGAAACTGGCGCAGTAGAACTGGAAAGCCCGTTCATCctg          >  1:1618793/1‑93 (MQ=255)

                                                                                 aaGGTATGCGGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGAAACTGGCGCAGTAGAACTGGAAAGCCCGTTCATCctgc         <  2:3796669/93‑1 (MQ=255)

                                                                                  aGGTATGCGGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGAAACTGGCGCAGTAGAACTGGAAAGCCCGTTCATCctgct        <  2:1750724/93‑1 (MQ=255)

                                                                                   ggTATGCGGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGAAACTGGCGCAGTAGAACTGGAAAGCCCGTTCATCctgctg       >  2:1866649/1‑93 (MQ=255)

                                                                                    gTATGCGGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGAAACTGGCGCAGTAGAACTGGAAAGCCCGTTCATCCTgctgg      <  2:1403777/93‑1 (MQ=255)

                                                                                    gTATGCGGTTCGACCGTGGCTACCTGTATCCTTACTTCATCAACAAGCCGGAAACTGGCGCAGTAGAACTGGAAAGCCCGTTCATCCTgctgg      <  2:3661070/93‑1 (MQ=255)

                                                                                     tATGCGGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGAAACTGGCGCAGTAGAACTGGAAAGCCCGTTCATCCTgctggc     >  1:15022/1‑93 (MQ=255)

                                                                                      aTGCGGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGAAACTGGCGCAGTAGAACTGGAAAGCCCGTTCATCCTgctggct    <  1:3748601/93‑1 (MQ=255)

                                                                                       tGCGGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGAAACTGGCGCAGTAGAACTGGAAAGCCCGTTCATCCTgctggctg   >  2:3745546/1‑93 (MQ=255)

                                                                                        gCAGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGAAACTGGCGCAGTAGAACTGGAAAGCCCGTTCATCCTGCTGGCTGa  >  2:4105420/1‑93 (MQ=255)

AAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCAGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACAAGCCGGAAACTGGCGCAGTAGAACTGGAAAGCCCGTTCATCCTGCTGGCTGA  >  NC_000913/4369093‑4369273

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