BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,866,592:1	+C	100%	coding (960/993 nt)	rpoS ←	RNA polymerase, sigma S (sigma 38) factor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,866,588	1	.	C	93.0%	493.8 / 34.2	158	coding (964/993 nt)	rpoS	RNA polymerase, sigma S (sigma 38) factor
Reads supporting (aligned to +/- strand): ref base . (6/5); new base C (77/66); total (86/72)
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 170 total aligned reads displayed.

CATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑G‑C‑CC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACA  >  NC_000913/2866489‑2866682

cATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑g                                                                                                   <  1:1333006/100‑1 (MQ=255)

 aTCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGAGCTTCGATATTC‑A‑G‑c                                                                                                 >  1:2559994/1‑100 (MQ=255)

    ggACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑G‑c‑cc‑c                                                                                            <  2:514280/100‑1 (MQ=255)

    ggACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑G‑c‑cc‑c                                                                                            >  1:750393/1‑100 (MQ=255)

     gACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑G‑c‑cc‑cc                                                                                           >  1:924305/1‑99 (MQ=255)

     gACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGATTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑G‑c‑cc‑cc                                                                                           >  1:3915701/1‑99 (MQ=255)

         ttttATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑G‑c‑cc‑cctgcg                                                                                       <  1:1130673/100‑6 (MQ=255)

         ttttATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑G‑c‑cc‑cctgcg                                                                                       >  2:1187878/1‑95 (MQ=255)

         ttttATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑G‑c‑cc‑cctgcg                                                                                       >  1:825094/1‑95 (MQ=255)

            tATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GCC‑CC‑CTGCGttt                                                                                     >  1:1808749/1‑100 (MQ=255)

            tATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GCC‑CC‑CTGCGttt                                                                                     >  2:562646/1‑100 (MQ=255)

               tgtgCACATAAAAGGCCAGCCTCGCTTGAGACTTGTCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GCC‑CC‑CTGCGTTTGCa                                                                                  <  1:499816/100‑1 (MQ=255)

               tgtgCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GCC‑CC‑CTGCGTTTGCa                                                                                  <  1:1778342/100‑1 (MQ=255)

                gtgCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GCC‑CC‑CTGCGTTTGCAg                                                                                 <  2:2968271/100‑1 (MQ=255)

                gtgCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GCC‑CC‑CTGCGTTTGCAg                                                                                 <  2:3896967/100‑1 (MQ=255)

                    acaGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GCC‑CC‑CTGCGTTTGCAGGAtt                                                                             >  2:2186390/1‑100 (MQ=255)

                     caGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGAAACAGCGATTAAATATTC‑A‑GCC‑CC‑CTACGTTTGCAGGAttt                                                                            >  1:3236250/1‑100 (MQ=255)

                      aGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GCC‑CC‑CTGCGTTTGCacgaaatc                                                                           >  1:3416647/1‑93 (MQ=255)

                      aGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GCC‑CC‑CTGCGTTTGCAGGATTTc                                                                           >  1:1551249/1‑100 (MQ=255)

                       gAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GCC‑CC‑CTGCGTTTGCAGGATTTcg                                                                          <  1:2421670/100‑1 (MQ=255)

                        aaaaGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GCC‑CC‑CTGCGTTTGCAGGATTTcg                                                                          <  2:3211113/99‑1 (MQ=255)

                         aaaGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GCC‑CC‑CTGCGTTTGCAGGATTTcgcg                                                                        <  1:2374246/100‑1 (MQ=255)

                          aaGGCCATCCTCGCTTGAGACTGTCTTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GCC‑CC‑CTGCGTTTGCAGGATTTCgcgc                                                                       <  1:313809/100‑1 (MQ=255)

                           aGGTCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GCC‑CC‑CTGCGTTTGCAGGATTTcgcg                                                                        <  1:1868488/98‑1 (MQ=255)

                           aGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGCAACATCACTTAGATATTC‑A‑A‑CACA‑CTGCATTAGCAGGATTACGCGCa                                                                      >  1:1846819/1‑100 (MQ=255)

                            ggCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑A‑C‑CCACTGCATTTGCAGGATTTCGCGCaa                                                                     >  1:2507894/1‑100 (MQ=255)

                            ggCCAGCCTCGCTGGAGACTGGCCTTTCTCACAGATCCTTACTTACTAGCGGAACAACGATTCGAAATTC‑A‑GCC‑CC‑CTTCGTTTGCAGGATTTCGAGCaa                                                                     >  1:1327592/1‑100 (MQ=255)

                             gCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GCC‑CC‑CTGCGTTTGCAGGATTTCGCGCaaa                                                                    >  2:3030385/1‑100 (MQ=255)

                              ccAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GCC‑CC‑CTGCGTTTGCAGGATTTCGCGCAAAc                                                                   >  1:607043/1‑100 (MQ=255)

                              ccAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GCC‑CC‑CTGCGTTTGCAGGATTTCGCGCAAAc                                                                   >  1:974230/1‑100 (MQ=255)

                              ccAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GCC‑CC‑CTGCGTTTGCAGGATTTCGCGCAAAc                                                                   <  2:1285657/100‑1 (MQ=255)

                              ccAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GCC‑CC‑CTGCGTTTGCAGGATTTCGCGCAAAc                                                                   >  1:1970862/1‑100 (MQ=255)

                                aGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GCC‑CC‑CTGCGTTTGCAGGATTTCGCGCAAACgg                                                                 <  2:2773390/100‑1 (MQ=255)

                                  ccTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GCC‑CC‑CTGCGTTTGCAGGATTTCGCGCAAACGgcg                                                               >  1:125103/1‑100 (MQ=255)

                                     cGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GCC‑CC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAg                                                            >  1:1672000/1‑100 (MQ=255)

                                     cGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GCC‑CC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAg                                                            <  2:1786125/100‑1 (MQ=255)

                                      gCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GCC‑CC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAgg                                                           <  2:1808734/100‑1 (MQ=255)

                                       cTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GCC‑CC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGc                                                          >  2:3620191/1‑100 (MQ=255)

                                        ttGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GCC‑CC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGcc                                                         >  2:2829366/1‑100 (MQ=255)

                                        ttGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GCC‑CC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGcc                                                         <  1:3758058/100‑1 (MQ=255)

                                         tGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GCC‑CC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCt                                                        >  2:1933578/1‑100 (MQ=255)

                                         tGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GCC‑CC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCt                                                        >  2:3408357/1‑100 (MQ=255)

                                         tGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GCC‑CC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCt                                                        >  1:3305069/1‑100 (MQ=255)

                                         tGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GCC‑CC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCt                                                        >  2:3209599/1‑100 (MQ=255)

                                          gagaCTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GCC‑CC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCtt                                                       >  2:252459/1‑100 (MQ=255)

                                           agaCTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GCC‑CC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTc                                                      <  1:3253460/100‑1 (MQ=255)

                                           agaCTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GAC‑CC‑CTGCGTTTGCAAGATTTCGCGAAAAAGGCGAAg                                                            >  1:3523407/1‑94 (MQ=255)

                                             aCTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GCC‑CC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCaa                                                    >  2:1404225/1‑100 (MQ=255)

                                              cTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GCC‑CC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAAc                                                   <  1:2754891/100‑1 (MQ=255)

                                              cTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GCC‑CC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAAc                                                   <  2:3093838/100‑1 (MQ=255)

                                               tGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GCC‑CC‑CTGCGTTTGCAGTATTTCGCGCAAACGGCGCAGGCCTTCAAc                                                   <  2:938359/99‑1 (MQ=255)

                                               tGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑G‑C‑CC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCt                                                 <  1:1121019/100‑1 (MQ=255)

                                               tGGCCTTTCTGACAGATGCTTACTTACCCGCGGAACAGCGCTTCGATATTC‑C‑GCC‑CC‑CTGCGTTTGCAGGATTTCGCTCAAACGGCGCAGGCCTTCAAcc                                                  <  1:869950/100‑1 (MQ=255)

                                                ggCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GCC‑CC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCt                                                 >  2:61717/1‑100 (MQ=255)

                                                ggCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GCC‑CC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCt                                                 <  2:825084/100‑1 (MQ=255)

                                                 gCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GCC‑CC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTg                                                >  1:2397161/1‑100 (MQ=255)

                                                  ccTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GCC‑CC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGa                                               >  1:3175905/1‑100 (MQ=255)

                                                    tttCTGACAGATGCTTATTTACTTGGGGAAGAGCGCTTCGTTTTTC‑A‑GCC‑CC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGaa                                              <  1:542028/99‑1 (MQ=255)

                                                         gACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GCC‑CC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGc                                        <  1:1645365/100‑1 (MQ=255)

                                                          aCAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GCC‑CC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCg                                       <  1:2071562/100‑1 (MQ=255)

                                                           cAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GCC‑CC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGa                                      >  2:1466011/1‑100 (MQ=255)

                                                           cAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GCC‑CC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGa                                      <  1:2255231/100‑1 (MQ=255)

                                                           cAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GCC‑CC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGa                                      <  2:2937346/100‑1 (MQ=255)

                                                           cAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GCC‑CC‑CTGCGTTTGCAGGATTTAGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGa                                      >  1:2334987/1‑100 (MQ=255)

                                                             gATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GCC‑CC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAac                                    <  1:1320674/100‑1 (MQ=255)

                                                              aTGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GCC‑CC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAaca                                   <  1:3364833/100‑1 (MQ=255)

                                                               tGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GCC‑CC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAacac                                  >  1:2928379/1‑100 (MQ=255)

                                                               tGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑GCC‑CC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAacac                                  >  1:2928844/1‑100 (MQ=255)

                                                               tGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑G‑CGCC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCg                                       >  2:1955728/1‑95 (MQ=255)

                                                                 cttacttaCTCGCGGAACAGCGCTTCGATATTC‑A‑GCC‑CC‑CTGCGTTTTCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGt                                >  1:1413586/1‑100 (MQ=255)

                                                                     cttactCGCGGAACAGCGCTTCGATATTC‑A‑GCC‑CC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCAc                            >  2:3167674/1‑100 (MQ=255)

                                                                        actCGCGGAACAGCGCTTCGATATTC‑A‑GCC‑CC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACggg                         >  2:109001/1‑100 (MQ=255)

                                                                         ctCGCGGAACAGCGCTTCGATATTC‑A‑GCC‑CC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCTACCTGAATCTGGCGAACACGTTCACGGGt                        >  2:833869/1‑100 (MQ=255)

                                                                         ctCGCGGAACAGCGCTTCGATATTC‑A‑GCC‑CC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACggg                         <  2:2057146/99‑1 (MQ=255)

                                                                         ctCGCGGAACAGCGCTTCGATATTC‑A‑GCC‑CC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGt                        >  1:170868/1‑100 (MQ=255)

                                                                         ctCGCGGAACAGCGCTTCGATATTC‑A‑GCC‑CC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGt                        >  1:1772188/1‑100 (MQ=255)

                                                                           cgcgGAACAGCGCTTCGATATTC‑A‑GCC‑CC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGa                      >  2:980073/1‑100 (MQ=255)

                                                                           cgcgAAACAGCGCTTCGATATTC‑A‑GCC‑CC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTg                       <  2:3297157/99‑1 (MQ=255)

                                                                            gtgttactgcgcTTCGATATTC‑ATT‑T‑CC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAg                     <  1:873867/92‑1 (MQ=255)

                                                                            gcgGAACAGCGCTTCGATATTCAA‑C‑C‑CC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCCCGGGTGAg                     >  2:1174121/1‑100 (MQ=255)

                                                                            gcgGAACAGCGCTTCGATATTC‑A‑G‑C‑CCTCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCGGAATCTGGCGAACACGTTCACGGGTGAg                     <  1:3221937/100‑1 (MQ=255)

                                                                             cgGAACAGCGCTTCGATATTC‑A‑GCC‑CC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAgg                    <  1:2625056/100‑1 (MQ=255)

                                                                               gAACAGCGCTTCGATATTC‑A‑GCC‑CC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGcc                  >  1:3764679/1‑100 (MQ=255)

                                                                                  cAGCGCTTCGATATTC‑A‑GCC‑CC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAAt               >  2:1983784/1‑100 (MQ=255)

                                                                                   aGCGCTTCGATATTC‑A‑GCC‑CC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAAtt              >  2:1680779/1‑100 (MQ=255)

                                                                                    gcgcTTCGATATTC‑A‑GCC‑CC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGTAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAAttt             <  1:1282709/100‑1 (MQ=255)

                                                                                      gcTTCGATATTC‑A‑GCC‑CC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCa           <  1:3312194/100‑1 (MQ=255)

                                                                                       cctCGATATTC‑A‑GCC‑CC‑CTGCGTTTGCAGGATCTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCAc          <  2:147296/98‑1 (MQ=255)

                                                                                       cTTCGATATTC‑A‑GCC‑CC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCAc          >  2:3296754/1‑100 (MQ=255)

                                                                                        ttCGATATTC‑A‑GCC‑CC‑CTGCGTTTGCAGGATTTCTCTCAAACGGCGCAGGCCTTCAACCTTAATCTGGCGAACACGTTCACGTGTGAGGCCAATTTCACg         <  2:1834346/100‑1 (MQ=255)

                                                                                         tCGATATTC‑A‑GCC‑CC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGa        >  2:3781472/1‑100 (MQ=255)

                                                                                         tCGATATTC‑A‑GCC‑CC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGa        >  2:447789/1‑100 (MQ=255)

                                                                                         tCGATATTC‑A‑GCC‑CC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGa        <  2:750384/100‑1 (MQ=255)

                                                                                         tCGATATTC‑A‑G‑C‑CC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGAc       >  2:3798701/1‑100 (MQ=255)

                                                                                         ncGATATTC‑A‑GCC‑CC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGa        >  1:3028385/2‑100 (MQ=255)

                                                                                          cGATATTC‑A‑GTC‑CC‑CTTCGTTTGCAGGATTTGGCGCAAACGGCGCATGCCTTCAATCTGGATCTGGCGAACACGTTCACGGTTGTGGCCAAGTGCATGAc       >  2:3526684/1‑100 (MQ=255)

                                                                                           gATATTC‑A‑GCC‑CC‑CTGCGATTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCAGGAGTGAGGCCAATTTCACGa        >  2:459047/1‑98 (MQ=255)

                                                                                             tatTC‑A‑GCC‑CC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCTTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTa    <  1:1260881/100‑1 (MQ=255)

                                                                                              atTC‑A‑GCC‑CC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTAc   >  2:2323933/1‑100 (MQ=255)

                                                                                               ttC‑A‑GCC‑CC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACa  >  1:2299953/1‑100 (MQ=255)

CATCGGACCTTTTATTGTGCACAGAAAAGGCCAGCCTCGCTTGAGACTGGCCTTTCTGACAGATGCTTACTTACTCGCGGAACAGCGCTTCGATATTC‑A‑G‑C‑CC‑CTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACA  >  NC_000913/2866489‑2866682

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