BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	561,800	A→G	T702A (ACC→GCC)	sfmD →	putative fimbrial usher protein SfmD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	561,800	0	A	G	99.8%	1680.7 / NA	403	T702A (ACC→GCC)	sfmD	putative fimbrial usher protein SfmD
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (173/229); total (173/230)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 403 total aligned reads displayed.

ATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTA  >  NC_000913/561702‑561897

 taACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATGCCa                                                                                                <  2:1093890/101‑1 (MQ=255)

  aaCTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATGCCaa                                                                                               <  1:1093450/101‑1 (MQ=255)

  aaCTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATGCCaa                                                                                               <  2:1094008/101‑1 (MQ=255)

   aCTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATGCCAAt                                                                                              >  2:1094580/1‑101 (MQ=255)

   aCTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATGCCAAt                                                                                              <  2:1093828/101‑1 (MQ=255)

    cTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATGCCAATg                                                                                             >  1:1094582/1‑101 (MQ=255)

     tACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATGCCAATGt                                                                                            <  2:1094002/101‑1 (MQ=255)

      aCGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATGCCAATGtt                                                                                           >  1:1094586/1‑101 (MQ=255)

       cGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATGCCAATGttt                                                                                          <  1:1094145/101‑1 (MQ=255)

        gACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATGCCAATGtttt                                                                                         <  1:1093978/101‑1 (MQ=255)

         aCCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATGCCAATGTTTTg                                                                                        >  1:1094592/1‑101 (MQ=255)

          ccGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATGCCAATGTTTTGa                                                                                       <  1:1093987/101‑1 (MQ=255)

          ccGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATGCCAATGTTTTGa                                                                                       >  2:1094596/1‑101 (MQ=255)

            gcgATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATGCCAATGtt                                                                                           <  2:1093926/95‑1 (MQ=255)

            gcgATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATGCCAATGTTTTGAtt                                                                                     >  1:1094599/1‑101 (MQ=255)

              gATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATGCCAATGTTTTGATTaa                                                                                   <  2:1094335/101‑1 (MQ=255)

               aTCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATGCCAATGTTTTGATTaaa                                                                                  >  1:1094605/1‑101 (MQ=255)

                tCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATGCCAATGTTTTGATTAAAg                                                                                 >  2:1094606/1‑101 (MQ=255)

                  aaCATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATGCCAATGTTTTGATTAAAgcg                                                                               >  2:1094613/1‑101 (MQ=255)

                  aaCATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATGCCAATGTTTTGATTAAAgcg                                                                               >  1:1094614/1‑101 (MQ=255)

                   aCATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATGCCAATGTTTTGATTAAAgcgc                                                                              <  2:1093901/101‑1 (MQ=255)

                    cATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATGCCAATGTTTTGATTAAAGCGcc                                                                             <  2:1093765/101‑1 (MQ=255)

                     aTGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATGCCAATGTTTTGATTAAAGCGCCt                                                                            <  2:1094023/101‑1 (MQ=255)

                      tGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATGCCAATGTTTTGATTAAAGCGCCTg                                                                           <  2:1094095/101‑1 (MQ=255)

                      tGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATGCCAATGTTTTGATTAAAGCGCCTg                                                                           <  1:1093959/101‑1 (MQ=255)

                      tGACGTTAACTGGCAGCTGTCTGGCGGGGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATGCCAATGTTTTGATTAAAGCGCCTg                                                                           >  1:1094622/1‑101 (MQ=255)

                       gACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATGCCAATGTTTTGATTAAAGCGCCTgg                                                                          <  1:1093788/101‑1 (MQ=255)

                        aCGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTATGGGATGCCAATGTTTTGATTAAAGCGCCTGgc                                                                         >  1:1094627/1‑101 (MQ=255)

                        aCGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATGCCAATGTTTTGATTAAAGCGCCTGgc                                                                         <  1:1093663/101‑1 (MQ=255)

                        aCGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATGCCAATGTTTTGATTAAAGCGCCTGgc                                                                         <  1:1094084/101‑1 (MQ=255)

                         cGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATGCCAATGTTTTGATTAAAGCGCCTGgcg                                                                        <  1:1093868/101‑1 (MQ=255)

                          gTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATGCCAATGTTTTGATTAAAGCGCCTGgcgc                                                                       >  2:1094629/1‑101 (MQ=255)

                            tAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATGCCAATGTTTTGATTAAAGCGCCTGGCGCAg                                                                     <  1:1093840/101‑1 (MQ=255)

                            tAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATGCCAATGTTTTGATTAAAGCGCCTGGCGCAg                                                                     <  2:1094241/101‑1 (MQ=255)

                             aaCTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATGCCAATGTTTTGATTAAAGCGCCTGGCGCAgg                                                                    >  1:1094633/1‑101 (MQ=255)

                                 ggCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATGCCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGgt                                                                >  2:1094637/1‑101 (MQ=255)

                                 ggCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATGCCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGgt                                                                <  1:1094190/101‑1 (MQ=255)

                                 ggCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATGCCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGgt                                                                <  2:1094058/101‑1 (MQ=255)

                                  gCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATGCCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGgtg                                                               <  2:1094033/101‑1 (MQ=255)

                                   cAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATGCCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGgtgt                                                              >  1:1094640/1‑101 (MQ=255)

                                   cAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATGCCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGgtgt                                                              <  2:1094268/101‑1 (MQ=255)

                                    aGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATGCCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTa                                                             <  2:1094306/101‑1 (MQ=255)

                                       tGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATGCCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGc                                                          >  1:1094645/1‑101 (MQ=255)

                                       tGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATGCCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGc                                                          <  2:1094074/101‑1 (MQ=255)

                                       tGTCTGGAGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATGCCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGc                                                          <  1:1094229/101‑1 (MQ=255)

                                        gTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATGCCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCa                                                         <  2:1093982/101‑1 (MQ=255)

                                         tCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATGCCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCAt                                                        <  2:1094442/101‑1 (MQ=255)

                                          cTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATGCCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCAtt                                                       <  1:1094119/101‑1 (MQ=255)

                                          cTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATGCCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCAtt                                                       <  2:1094206/101‑1 (MQ=255)

                                           tGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATGCCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTg                                                      >  2:1094651/1‑101 (MQ=255)

                                            ggCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATGCCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGa                                                     <  2:1094368/101‑1 (MQ=255)

                                            ggCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATGCCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGa                                                     >  1:1094652/1‑101 (MQ=255)

                                                  gtgGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATGCCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCa                                               >  2:1094659/1‑101 (MQ=255)

                                                   tgGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATGCCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCaa                                              >  1:1094661/1‑101 (MQ=255)

                                                      tCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATGCCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCaa                                              <  2:1094153/98‑1 (MQ=255)

                                                      tCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATGCCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACt                                           >  1:1094666/1‑101 (MQ=255)

                                                      tCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATGCCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACt                                           <  1:1094412/101‑1 (MQ=255)

                                                       cGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATGCCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTg                                          >  2:1094668/1‑101 (MQ=255)

                                                            aTGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATGCCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCAtt                                     <  2:1094051/101‑1 (MQ=255)

                                                              gAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATGCCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCAtttt                                   >  1:1094674/1‑101 (MQ=255)

                                                               aaaaTGGCATAACGCTGAGCCAGGCTTTAGGGGATGCCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGc                                        <  1:1094063/95‑1 (MQ=255)

                                                               aaaaTGGCATAACGCTGAGCCAGCGTTTAGGGGATGCCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCAt                                      <  2:1093859/97‑1 (MQ=255)

                                                                aaaTGGCATAACGCTGAGCCAGCCTTTAGGGGATGCCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTaa                                 <  1:1094168/101‑1 (MQ=255)

                                                                aaaTGGCATAACGCTGAGCCAGCCTTTAGGGGATGCCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTaa                                 <  2:1094187/101‑1 (MQ=255)

                                                                aaaTGGCATAACGCTGAGCCAGCCTTTAGGGGATGCCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTaa                                 >  2:1094679/1‑101 (MQ=255)

                                                                aaaTGGCATAACGCTGAGCCAGCCTTTAGGGGATGCCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTaa                                 <  2:1094256/101‑1 (MQ=255)

                                                                 aaTGGCATAACGCTGAGCCAGCCTTTAGGGGATGCCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAAc                                >  1:1094680/1‑101 (MQ=255)

                                                                 aaTGGCATAACGCTGAGCCAGCCTTTAGGGGATGCCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAAc                                <  1:1093924/101‑1 (MQ=255)

                                                                      cATAACGCTGAGCCAGCCTTTAGGGGATGCCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACt                           >  2:1094687/1‑101 (MQ=255)

                                                                      cATAACGCTGAGCCAGCCTTTAGGGGATGCCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACt                           <  1:1094338/101‑1 (MQ=255)

                                                                       aTAACGCTGAGCCAGCCTTTAGGGGATGCCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTg                          >  1:1094689/1‑101 (MQ=255)

                                                                       aTAACGCTGAGCCAGCCTTTAGGGGATGCCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTg                          <  1:1093945/101‑1 (MQ=255)

                                                                        tAACGCTGAGCCAGCCTTTAGGGGATGCCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTgg                         >  1:1094694/1‑101 (MQ=255)

                                                                           cGCTGAGCCAGCCTTTAGGGGATGCCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGgcgc                      <  1:1094130/101‑1 (MQ=255)

                                                                           cGCTGAGCCAGCCTTTAGGGGATGCCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGgcgc                      <  1:1094251/101‑1 (MQ=255)

                                                                            gCTGAGCCAGCCTTTAGGGGATGCCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGcgcg                     <  2:1093618/101‑1 (MQ=255)

                                                                            gCTGAGCCAGCCTTTAGGGGATGCCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGcgcg                     >  2:1094698/1‑101 (MQ=255)

                                                                              tGAGCCAGCCTTTAGGGGATGCCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGc                   >  1:1094705/1‑101 (MQ=255)

                                                                              tGAGCCAGCCTTTAGGGGATGCCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGc                   >  2:1094703/1‑101 (MQ=255)

                                                                               gAGCCAGCCTTTAGGGGATGCCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCt                  <  2:1094281/101‑1 (MQ=255)

                                                                               gAGCCAGCCTTTAGGGGATGCCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCt                  <  1:1093954/101‑1 (MQ=255)

                                                                                agccagccTTTAGGGGATGCCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTa                 >  2:1094709/1‑101 (MQ=255)

                                                                                agccagccTTTAGGGGATGCCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTa                 <  1:1094046/101‑1 (MQ=255)

                                                                                 gccagccTTTAGGGGATGCCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTAt                <  2:1094211/101‑1 (MQ=255)

                                                                                 gccagccTTTAGGGGATGCCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTAt                >  1:1094712/1‑101 (MQ=255)

                                                                                 gccagccTTTAGGGGATGCCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTAt                <  1:1094009/101‑1 (MQ=255)

                                                                                  ccagccTTTAGGGGATGCCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATg               <  2:1094112/101‑1 (MQ=255)

                                                                                   cagccTTTAGGGGATGCCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGc              <  2:1093961/101‑1 (MQ=255)

                                                                                   cagccTTTAGGGGATGCCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGc              <  1:1093904/101‑1 (MQ=255)

                                                                                    agccTTTAGGGGATGCCAATGTTTTGATTAAAGCGCCTgg                                                                          <  1:1094106/40‑1 (MQ=255)

                                                                                    agccTTTAGGGGATGCCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCg             >  1:1094718/1‑101 (MQ=255)

                                                                                     gccTTTAGGGGATGCCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCgg            >  2:1094723/1‑101 (MQ=255)

                                                                                       cTTTAGGGGATGCCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTg          >  1:1094727/1‑101 (MQ=255)

                                                                                        tttATGGGATGCCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGa         >  2:1094730/1‑101 (MQ=255)

                                                                                        tttAGGGGATGCCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCgg            <  1:1094287/98‑1 (MQ=255)

                                                                                        tttAGGGGATGCCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGa         >  1:1094733/1‑101 (MQ=255)

                                                                                          tAGGGGATGCCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATg       <  1:1094516/101‑1 (MQ=255)

                                                                                          tAGGGGATGCCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATg       >  2:1094736/1‑101 (MQ=255)

                                                                                           aGGGGATGCCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGc      >  2:1094743/1‑101 (MQ=255)

                                                                                            ggggATGCCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCt     >  1:1094744/1‑101 (MQ=255)

ATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTAACCGACTGGCGCGGCTATGCGGTGATGCTGTA  >  NC_000913/561702‑561897

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