BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,101,826	T→C	I198M (ATA→ATG)	csgD ←	CsgD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,101,826	0	T	C	81.5%	547.3 / 129.4	265	I198M (ATA→ATG)	csgD	CsgD
Reads supporting (aligned to +/- strand): ref base T (23/26); new base C (100/116); total (123/142)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.56e-01

Only 100 of 265 total aligned reads displayed.

CGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTATCTTCTTGAAAAGATTATAAAGATGCGTTTTAACCGTATTTTCGCTGATGAACAACGAACGAGCGATCTCGTTATTAGACGCGCCGATAC  >  NC_000913/1101735‑1101917

cgccACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTa                                                                                            <  2:80348/93‑1 (MQ=255)

cgccACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTa                                                                                            <  2:1997743/93‑1 (MQ=255)

 gccACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCCAt                                                                                           >  1:5740543/1‑93 (MQ=255)

   cacaATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCCAtct                                                                                         <  2:6968641/93‑1 (MQ=255)

   cacaATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCCAtct                                                                                         >  2:6739212/1‑93 (MQ=255)

    acaATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCCAtctt                                                                                        >  2:3306865/1‑93 (MQ=255)

     caATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCCAtcttc                                                                                       >  1:3006050/1‑93 (MQ=255)

     caATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCCAtcttc                                                                                       <  1:859101/93‑1 (MQ=255)

     caATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCCAtcttc                                                                                       <  2:2139035/93‑1 (MQ=255)

      aaTCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCCAtcttct                                                                                      >  1:4090151/1‑93 (MQ=255)

       aTCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCCATcttctt                                                                                     >  2:4363940/1‑93 (MQ=255)

       aTCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCCATcttctt                                                                                     >  2:332310/1‑93 (MQ=255)

        tCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTATCTTCTTg                                                                                    >  2:670378/1‑93 (MQ=255)

        tCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCCATCTTCTTg                                                                                    >  2:6166028/1‑93 (MQ=255)

         ccAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCCATCTTCTTGa                                                                                   <  1:6109535/93‑1 (MQ=255)

          cAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCCATCTTCTTGaa                                                                                  <  1:248010/93‑1 (MQ=255)

           aGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCCATCTTCTTGaaa                                                                                 <  1:110065/93‑1 (MQ=255)

           aGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCCATCTTCTTGaaa                                                                                 >  2:3940078/1‑93 (MQ=255)

            gCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCCATCTTCTTGaaaa                                                                                <  1:5496579/93‑1 (MQ=255)

            gCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCCATCTTCTTGaaaa                                                                                >  2:2947397/1‑93 (MQ=255)

              gTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCCATCTTCTTGAAAAGa                                                                              >  2:2976862/1‑93 (MQ=255)

               tAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCCATCTTCTTGAAAAGAt                                                                             <  2:5011019/93‑1 (MQ=255)

               tAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCCATCTTCTTGAAAAGAt                                                                             >  1:4655539/1‑93 (MQ=255)

                aaaTAACGTTTCATGGCTTTATCGCGTGAGGTTATCGTTTGCGCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTATCTTCTTGAAAAGAtt                                                                            <  2:4109986/93‑1 (MQ=255)

                aaaTAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCCATCTTCTTGAAAAGAtt                                                                            >  1:1627897/1‑93 (MQ=255)

                aaaTAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCCATCTTCTTGAAAAGAtt                                                                            >  2:7077521/1‑93 (MQ=255)

                 aaTAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCCATCTTCTTGAAAAGATta                                                                           <  2:3302876/93‑1 (MQ=255)

                 aaTAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCCATCTTCTTGAAAAGATta                                                                           <  1:6563584/93‑1 (MQ=255)

                  aTAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCCATCTTCTTGAAAAGATtat                                                                          <  1:2871346/93‑1 (MQ=255)

                    aaCGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTATCTTCTTGAAAAGATTATaa                                                                        <  2:2548933/93‑1 (MQ=255)

                    aaCGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCCATCTTCTTGAAAAGATTATaa                                                                        >  1:3972952/1‑93 (MQ=255)

                     aCGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCCATCTTCTTGAAAAGATTATaaa                                                                       <  1:77446/93‑1 (MQ=255)

                       gTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCCATCTTCTTGAAAAGATTATAAAGa                                                                     >  1:3437448/1‑93 (MQ=255)

                       gTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCCATCTTCTTGAAAAGATTATAAAGa                                                                     <  1:2876860/93‑1 (MQ=255)

                       gTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCCATCTTCTTGAAAAGATTATAAAGa                                                                     <  2:4219937/93‑1 (MQ=255)

                       gTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCCATCTTCTTGAAAAGATTATAAAGa                                                                     >  1:2265045/1‑93 (MQ=255)

                        tttCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTATCTTCTTGAAAAGATTATAAAGAt                                                                    >  1:1777192/1‑93 (MQ=255)

                          tCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCCATCTTCTTGAAAAGATTATAAAGATGc                                                                  >  2:891698/1‑93 (MQ=255)

                            aTGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCCATCTTCTTGAAAAGATTATAAAGATGCGt                                                                >  1:3881554/1‑93 (MQ=255)

                             tGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCCATCTTCTTGAAAAGATTATAAAGATGCGtt                                                               <  1:6166146/93‑1 (MQ=255)

                              ggCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCCATCTTCTTGAAAAGATTATAAAGATGCGttt                                                              <  2:1331710/93‑1 (MQ=255)

                                 tttATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTATCTTCTTGAAAAGATTATAAAGATGCGTTTTaa                                                           <  1:1346994/93‑1 (MQ=255)

                                 tttATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCCATCTTCTTGAAAAGATTATAAAGATGCGTTTTaa                                                           <  2:2360973/93‑1 (MQ=255)

                                   tATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCCATCTTCTTGAAAAGATTATAAAGATGCGTTTTAAcc                                                         >  2:4621885/1‑93 (MQ=255)

                                   tATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCCATCTTCTTGAAAAGATTATAAAGATGCGTTTTAAcc                                                         >  2:714244/1‑93 (MQ=255)

                                   tATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGACCATCTTCTTGAAAAGATTATAAAGATGCGTTTTAAcc                                                         >  1:4967409/1‑93 (MQ=255)

                                    aTCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCCATCTTCTTGAAAAGATTATAAAGATGCGTTTTAACCg                                                        <  1:6860583/93‑1 (MQ=255)

                                     tCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCCATCTTCTTGAAAAGATTATAAAGATGCGTTTTAACCGt                                                       >  1:5220698/1‑93 (MQ=255)

                                      cGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCCATCTTCTTGAAAAGATTATAAAGATGCGTTTTAACCGTa                                                      <  1:5785435/93‑1 (MQ=255)

                                         cTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCCATCTTCTTGAAAAGATTATAAAGATGCGTTTTAACCGTAttt                                                   >  2:2280506/1‑93 (MQ=255)

                                          tGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTATCTTCTTGAAAAGATTATAAAGATGCGTTTTAACCGTAtttt                                                  >  1:419464/1‑93 (MQ=255)

                                          tGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCCATCTTCTTGAAAAGATTATAAAGATGCGTTTTAACCGTAtttt                                                  <  1:2231771/93‑1 (MQ=255)

                                          tGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCCATCTTCTTGAAAAGATTATAAAGATGCGTTTTAACCGTAtttt                                                  <  1:1231847/93‑1 (MQ=255)

                                           gAGGTTATCGTTTGCCCAGGAAACCGCTTATGTCCGGTTTTTTACGGCCATCTTCTTGAAAAGATTATAAAGATGCGTTTTAACCGTATTTTc                                                 <  1:1473575/93‑1 (MQ=255)

                                            aGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCCATCTTCTTGAAAAGATTATAAAGATGCGTTTTAACCGTATTTTCg                                                >  1:6953396/1‑93 (MQ=255)

                                            aGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCCATCTTCTTGAAAAGATTATAAAGATGCGTTTTAACCGTATTTTCg                                                >  1:7253758/1‑93 (MQ=255)

                                             ggTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCCATCTTCTTGAAAAGATTATAAAGATGCGTTTTAACCGTATTTTCGc                                               >  2:385745/1‑93 (MQ=255)

                                                   cGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTATCTTCTTGAAAAGATTATAAAGATGCGTTTTAACCGTATTTTCGCtgatga                                         <  2:1498707/93‑1 (MQ=255)

                                                     tttGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTATCTTCTTGAAAAGATTATAAAGATGCGTTTTAACCGTATTTTCGCTGATGaac                                       >  1:6686228/1‑93 (MQ=255)

                                                          ccAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTATCTTCTTGAAAAGATTATAAAGATGCGTTTTAACCGTATTTTCGCTGATGAACaacga                                  >  1:1520068/1‑93 (MQ=255)

                                                          ccAGGAAACCGCTTGTGTCCGGTTTTTTACGGCCATCTTCTTGAAAAGATTATAAAGATGCGTTTTAACCGTATTTTCGCTGATGAACaacga                                  <  2:6489696/93‑1 (MQ=255)

                                                          ccAGGAAACCGCTTGTGTCCGGTTTTTTACGGCCATCTTCTTGAAAAGATTATAAAGATGCGTTTTAACCGTATTTTCGCTGATGAACaacga                                  <  1:3650534/93‑1 (MQ=255)

                                                          ccAGGAAACCGCTTGTGTCCGGTTTTTTACGGCCATCTTCTTGAAAAGATTATAAAGATGCGTTTTAACCGTATTTTCGCTGATGAACaacga                                  >  1:3132587/1‑93 (MQ=255)

                                                           cAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTATCTTCTTGAAAAGATTATAAAGATGCGTTTTAACCGTATTTTCGCTGATGAACaacgaa                                 <  1:7149470/93‑1 (MQ=255)

                                                           cAGGAAACCGCTTGTGTCCGGTTTTTTACGGCCATCTTCTTGAAAAGATTATAAAGATGCGTTTTAACCGTATTTTCGCTGATGAACaacgaa                                 <  2:7427372/93‑1 (MQ=255)

                                                            aGGAAACCGCTTGTGTCCGGTTTTTTACGGCTATCTTCTTGAAAAGATTATAAAGATGCGTTTTAACCGTATTTTCGCTGATGAACaacgaac                                >  1:3182562/1‑93 (MQ=255)

                                                             ggAAACCGCTTGTGTCCGGTTTTTTACGGCCATCTTCTTGAAAAGATTATAAAGATGCGTTTTAACCGTATTTTCGCTGATGAACaacgaacg                               >  2:558446/1‑93 (MQ=255)

                                                              gAAACCGCTTGTGTCCGGTTTTTTACGGCCATCTTCTTGAAAAGATTATAAAGATGCGTTTTAACCGTATTTTCGCTGATGAACAacgaacga                              <  1:6286185/93‑1 (MQ=255)

                                                               aaaCCGCTTGTGTCCGGTTTTTTACGGCCATCTTCTTGAAAAGATTATAAAGATGCGTTTTAACCGTATTTTCGCTGATGAACAACGAACGAg                             <  2:4429119/93‑1 (MQ=255)

                                                                 aCCGCTTGTGTCCGGTTTTTTACGGCTATCTTCTTGAAAAGATTATAAAGATGCGTTTTAACCGTATTTTCGCTGATGAACAACGAACGAGCg                           <  1:100025/93‑1 (MQ=255)

                                                                  ccGCTTGTGTCCGGTTTTTTACGGCCATCTTCTTGAAAAGATTATAAAGATGCGTTTTAACCGTATTTTCGCTGATGAACAACGAACGAGCGa                          >  1:4594173/1‑93 (MQ=255)

                                                                  ccGCTTGTGTCCGGTTTTTTACGGCCATCTTCTTGAAAAGATTATAAAGATGCGTTTTAACCGTATTTTCGCTGATGAACAACGAACGAGCGa                          >  2:7194678/1‑93 (MQ=255)

                                                                    gCTTGTGTCCGGTTTTTTACGGCTATCTTCTTGAAAAGATTATAAAGATGCGTTTTAACCGTATTTTCGCTGATGAACAACGAACGAGCGAtc                        <  2:7484045/93‑1 (MQ=255)

                                                                    gCTTGTGTCCGGTTTTTTACGGCTATCTTCTTGAAAAGATTATAAAGATGCGTTTTAACCGTATTTTCGCTGATGAACAACGAACGAGCGAtc                        >  1:3851475/1‑93 (MQ=255)

                                                                    gCTTGTGTCCGGTTTTTTACGGCCATCTTCTTGAAAAGATTATAAAGATGCGTTTTAACCGTATTTTCGCTGATGAACAACGAACGAGCGAtc                        <  1:6024047/93‑1 (MQ=255)

                                                                    gCTTGTGTCCGGTTTTTTACGGCCATCTTCTTGAAAAGATTATAAAGATGCGTTTTAACCGTATTTTCGCTGATGAACAACGAACGAGCGAtc                        >  2:1796639/1‑93 (MQ=255)

                                                                     cTTGTGTCCGGTTTTTTACGGCCATCTTCTTGAAAAGATTATAAAGATGCGTTTTAACCGTATTTTCGCTGATGAACAACGAACGAGCGAtct                       >  1:2623811/1‑93 (MQ=255)

                                                                     cTTGTGTCCGGTTTTTTACGGCCATCTTCTTGAAAAGATTATAAAGATGCGTTTTAACCGTATTTTCGCTGATGAACAACGAACGAGCGAtct                       >  2:2715652/1‑93 (MQ=255)

                                                                     cTTGTGTCCGGTTTTTTACGGCCATCTTCTTGAAAAGATTATAAAGATGCGTTTTAACCGTATTTTCGCTGATGAACAACGAACGAGCGAtct                       <  2:3071607/93‑1 (MQ=255)

                                                                      ttGTGTCCGGTTTTTTACGGCTATCTTCTTGAAAAGATTATAAAGATGCGTTTTAACCGTATTTTCGCTGATGAACAACGAACGAGCGAtctc                      >  2:2434240/1‑93 (MQ=255)

                                                                      ttGTGTCCGGTTTTATACGGCCATCTTCTTGAAAAGATTATAAAGATGCGTTTTAACCGTATTTTCGCTGATGAACAACGAACGAGCGAtctc                      >  2:6441891/1‑93 (MQ=255)

                                                                       tgtgTCCGGTTTTTTACGGCCATCTTCTTGAAAAGATTATAAAGATGCGTTTTAACCGTATTTTCGCTGATGAACAACGAACGAGCGATCTCg                     >  2:5801660/1‑93 (MQ=255)

                                                                        gtgtCCGGTTTTTTACGGCCATCTTCTTGAAAAGATTATAAAGATGCGTTTTAACCGTATTTTCGCTGATGAACAACGAACGAGCGATCTCGt                    <  2:6041792/93‑1 (MQ=255)

                                                                        gtgtCCGGTTTTTTACGGCCATCTTCTTGAAAAGATTATAAAGATGCGTTTTAACCGTATTTTCGCTGATGAACAACGAACGAGCGATCTCGt                    <  2:1244205/93‑1 (MQ=255)

                                                                         tgtCCGGTTTTTTACGGCCATCTTCTTGAAAAGATTATAAAGATGCGTTTTAACCGTATTTTCGCTGATGAACAACGAACGAGCGATCTCGtt                   <  1:2031263/93‑1 (MQ=255)

                                                                           tCCGGTTTTTTACGGCCATCTTCTTGAAAAGATTATAAAGATGCGTTTTAACCGTATTTTCGCTGATGAACAACGAACGAGCGATCTCGttat                 <  1:7451472/93‑1 (MQ=255)

                                                                             cGGTTTTTTACGGCTATCTTCTTGAAAAGATTATAAAGATGCGTTTTAACCGTATTTTCGCTGATGAACAACGAACGAGCGATCTCGttatta               <  1:663262/93‑1 (MQ=255)

                                                                             cGGTTTTTTACGGCCATCTTCTTGAAAAGATTATAAAGATGCGTTTTAACCGTATTTTCGCTGATGAACAACGAACGAGCGATCTCGttatta               <  1:2504556/93‑1 (MQ=255)

                                                                             cGGTTTTTTACGGCCATCTTCTTGAAAAGATTATAAAGATGCGTTTTAACCGTATTTTCGCTGATGAACAACGAACGAGCGATCTCGttatta               <  1:237072/93‑1 (MQ=255)

                                                                              ggTTTTTTACGGCCATCTTCTTGAAAAGATTATAAAGATGCGTTTTAACCGTATTTTCGCTGATGAACAACGAACGAGCGATCTCGTTATTAg              >  1:5330432/1‑93 (MQ=255)

                                                                               gTTTTTTACGGCCATCTTCTTGAAAAGATTATAAAGATGCGTTTTAACCGTATTTTCGCTGATGAACAACGAACGAGCGATCTCGTTATTAGa             <  2:3713113/93‑1 (MQ=255)

                                                                                 tttttACGGCCATCTTCTTGAAAAGATTATAAAGATGCGTTTTAACCGTATTTTCGCTGATGAACAACGAACGAGCGATCTCGTTATTAGAcg           <  2:3663122/93‑1 (MQ=255)

                                                                                   tttACGGCCATCTTCTTGAAAAGATTATAAAGATGCGTTTTAACCGTATTTTCGCTGATGAACAACGAACGAGCGATCTCGTTATTAGAcgcg         <  1:6354645/93‑1 (MQ=255)

                                                                                    ttACGGCCATCTTCTTGAAAAGATTATAAAGATGCGTTTTAACCGTATTTTCGCTGATGAACAACGAACGAGCGATCTCGTTATTAGACgcgc        <  1:4223188/93‑1 (MQ=255)

                                                                                    ttACGGCCATCTTCTTGAAAAGATTATAAAGATGCGTTTTAACCGTATTTTCGCTGATGAACAACGAACGAGCGATCTCGTTATTAGACgcgc        <  2:5281924/93‑1 (MQ=255)

                                                                                       cGGCCATCTTCTTGAAAAGATTATAAAGATGCGTTTTAACCGTATTTTCGCTGATGAACAACGAACGAGCGATCTCGTTATTAGACGCGCCGa     >  1:3445835/1‑93 (MQ=255)

                                                                                        ggCCATCTTCTTGAAAAGATTATAAAGATGCGTTTTAACCGTATTTTCGCTGATGAACAACGAACGAGCGATCTCGTTATTAGACGCGCCGAt    <  1:1267161/93‑1 (MQ=255)

                                                                                        ggCCATCTTCTTGAAAAGATTATAAAGATGCGTTTTAACCGTATTTTCGCTGATGAACAACGAACGAGCGATCTCGTTATTAGACGCGCCGAt    <  2:5961860/93‑1 (MQ=255)

                                                                                         gCTATCTTCTTGAAAAGATTATAAAGATGCGTTTTAACCGTATTTTCGCTGATGAACAACGAACGAGCGATCTCGTTATTAGACGCGCCGATa   >  2:375249/1‑93 (MQ=255)

                                                                                          cTATCTTCTTGAAAAGATTATAAAGATGCGTTTTAACCGTATTTTCGCTGATGAACAACGAACGAGCGATCTCGTTATTAGACGCGCCGATAc  <  1:1888170/93‑1 (MQ=255)

CGCCACAATCCAGCGTAAATAACGTTTCATGGCTTTATCGCCTGAGGTTATCGTTTGCCCAGGAAACCGCTTGTGTCCGGTTTTTTACGGCTATCTTCTTGAAAAGATTATAAAGATGCGTTTTAACCGTATTTTCGCTGATGAACAACGAACGAGCGATCTCGTTATTAGACGCGCCGATAC  >  NC_000913/1101735‑1101917

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