BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,379,266:1	(C)_7→8	100%	coding (732/1128 nt)	zapE ←	divisome ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,379,259	1	.	C	94.6%	323.8 / 11.3	93	coding (739/1128 nt)	zapE	divisome ATPase
Reads supporting (aligned to +/- strand): ref base . (2/3); new base C (46/42); total (48/45)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 104 total aligned reads displayed.

GCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑T‑‑CC‑CC‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGC  >  NC_000913/3379161‑3379353

gCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTGACGTCGGTGAATTTTCTCGTT‑T‑‑c                                                                                                 >  2:1142644/1‑100 (MQ=255)

gCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑T‑‑c                                                                                                 <  2:567723/100‑1 (MQ=255)

 cacaGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑T‑‑cc                                                                                                >  2:1717084/1‑100 (MQ=255)

 cacaGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑T‑‑cc                                                                                                <  1:1361934/100‑1 (MQ=255)

  acaGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑T‑‑cc‑c                                                                                              >  2:1842909/1‑100 (MQ=255)

  acaGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑T‑‑cc‑c                                                                                              >  1:1975413/1‑100 (MQ=255)

  acaGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑T‑‑c                                                                                                 >  2:2177601/1‑98 (MQ=255)

   caGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑T‑‑cc‑cc                                                                                             <  2:2205543/100‑1 (MQ=255)

    aGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑T‑‑cc‑cgc                                                                                            >  1:1215173/1‑98 (MQ=255)

        tagtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑T‑‑cc‑cc‑ccccg                                                                                       >  1:2157155/1‑98 (MQ=255)

        tagtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑T‑‑cc‑cc‑ccccg                                                                                       <  2:1567361/100‑3 (MQ=255)

            aaaaGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑CC‑CC‑CCCGCCAg                                                                                    >  2:1553567/1‑100 (MQ=255)

            aaaaGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑CC‑CC‑CCCGCCAg                                                                                    >  2:1825471/1‑100 (MQ=255)

            aaaaGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑T‑‑CC‑CC‑CCCGCCAgc                                                                                   >  1:846616/1‑100 (MQ=255)

               agagACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑CC‑CC‑CCCGCCAgcgc                                                                                 <  1:867838/100‑1 (MQ=255)

                gagaCCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑CC‑CC‑CCCGCCAGCGcc                                                                                >  2:720387/1‑100 (MQ=255)

                 agaCCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑CC‑CC‑CCCGCCAGCGCCa                                                                               >  1:415548/1‑100 (MQ=255)

                 agaCCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑CC‑CC‑CCCGCCAGCGCCa                                                                               >  2:2258265/1‑100 (MQ=255)

                   cccGCCAGCGTCTGGTTCTCGACGCGCATTGTTGCTACTGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑CC‑CC‑CCCGCCAGCGCCa                                                                               <  2:813346/97‑1 (MQ=255)

                   aCCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑CC‑CC‑CCCGCCAGCGCCAAc                                                                             <  2:628733/100‑1 (MQ=255)

                    ccGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑CC‑CC‑CCCGCCCGCGCCAAcc                                                                            >  2:2216849/1‑100 (MQ=255)

                    ccGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑CC‑CC‑CCCGCCAGCGCCAAcc                                                                            >  1:1166306/1‑100 (MQ=255)

                      gCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑CC‑CC‑CCCGCCAGCGCCAACCa                                                                           >  1:1467687/1‑99 (MQ=255)

                      gCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑CC‑CC‑CCCGCCAGCGCCAACCAt                                                                          >  1:2021511/1‑100 (MQ=255)

                       ccAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑CC‑CC‑CCCGCCAGCGCCAACCATa                                                                         >  1:1709484/1‑100 (MQ=255)

                        cAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑CC‑CC‑CCCGCCAGCGCCAACCATAg                                                                        <  2:2260877/100‑1 (MQ=255)

                          gCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGCTGGTTGATTTCTAACGTCGGTGACTTTTCTCGTT‑TC‑CC‑CC‑CCCGCCAGCGCCAACCATAGtt                                                                      <  1:1282976/100‑1 (MQ=255)

                           cGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑CC‑CC‑CCCGCCAGCGCCAACCATAGttt                                                                     <  2:1336900/100‑1 (MQ=255)

                            gTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGTATTTTCTCGTT‑TCCCC‑CC‑CCCGCCAGCGCCAACCATAGttt                                                                     <  1:1876972/100‑1 (MQ=255)

                            gTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑CC‑CC‑CCCGCCAGCGCCAACCATAGTTTa                                                                    <  1:239388/100‑1 (MQ=255)

                             tCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑CC‑CC‑CCCGCCAGCGCCAACCATAGTTTAt                                                                   <  2:473535/100‑1 (MQ=255)

                              cTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑CC‑CC‑CCCGCCAGCGCCAACCATAGTTTATc                                                                  >  2:2141208/1‑100 (MQ=255)

                                ggTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑CC‑CC‑CCCGCCAGCGCCAACCATAGTTTATCCa                                                                >  1:1420593/1‑100 (MQ=255)

                                ggTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑CC‑CC‑CCCGCCAGCGCCAACCATAGTTTATCCa                                                                >  2:455137/1‑100 (MQ=255)

                                 gTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑CC‑CC‑CCCGCCAGCGCCAACCATAGTTTATCCAt                                                               >  2:1344748/1‑100 (MQ=255)

                                    ctcGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑CC‑CC‑CCCGCCAGCGCCAACCATAGTTTATTCAttt                                                             >  2:1621793/1‑99 (MQ=255)

                                     tcGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑CC‑CC‑CCCGCCAGCGCCCACAATAGTTTATCCATTTg                                                            >  2:208642/1‑99 (MQ=255)

                                     tcGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑CC‑CC‑CCCGCCAGCGCCAACCATAGTTTATCCATTTgc                                                           >  2:515980/1‑100 (MQ=255)

                                     tcGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑CC‑CC‑CCCGCCAGCACCAACCATAGTTTATCCAttt                                                             >  2:337967/1‑98 (MQ=255)

                                       ggcGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑CC‑CC‑CCCGCCAGCGCCAACCATAGTTTATCCATTTgcgc                                                         >  1:1586812/3‑100 (MQ=255)

                                       gACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑CC‑CC‑CCCGCCAGCGCCAACCATAGTTTATCCATTTgcgc                                                         >  1:450181/1‑100 (MQ=255)

                                        aCGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑CC‑CC‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGcc                                                        <  1:144663/100‑1 (MQ=255)

                                        aCGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑CC‑CC‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGcc                                                        <  1:579510/100‑1 (MQ=255)

                                         cGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGCATTTTCTCGTT‑TC‑CC‑CC‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGccc                                                       >  2:2293580/1‑100 (MQ=255)

                                          gCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑CC‑CC‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCg                                                      <  1:1806443/100‑1 (MQ=255)

                                           cccATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑CC‑CC‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCgg                                                     >  1:2140997/1‑100 (MQ=255)

                                            ccATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑CC‑CC‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCggg                                                    <  1:209439/100‑1 (MQ=255)

                                            ccATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑CC‑CC‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCggg                                                    >  2:1576051/1‑100 (MQ=255)

                                              aTTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑CC‑CC‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGtt                                                  <  1:713937/100‑1 (MQ=255)

                                              aTTGTTGCGAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCTTT‑TC‑CC‑CC‑CCCTCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGt                                                   <  1:1055486/99‑1 (MQ=255)

                                                tgttgCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑CC‑CC‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTc                                                >  1:1449587/1‑100 (MQ=255)

                                                tgttgCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑CC‑CC‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTc                                                <  1:2100624/100‑1 (MQ=255)

                                                tgttgCTAATGGCCGATGGTTGATTTCTAACGGCGGGGAATTTTCTCGTT‑TC‑CC‑CC‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTc                                                <  2:362247/100‑1 (MQ=255)

                                                 gttgCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑CC‑CC‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCa                                               <  1:455140/100‑1 (MQ=255)

                                                   tgCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑CC‑CC‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATc                                             <  2:1383037/100‑1 (MQ=255)

                                                    gCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑CC‑CC‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCg                                            <  1:1892829/100‑1 (MQ=255)

                                                     cTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑CC‑CC‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGt                                           <  2:1213304/100‑1 (MQ=255)

                                                      tAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑CC‑CC‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTg                                          >  2:1694180/1‑100 (MQ=255)

                                                      tAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑CC‑CC‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTg                                          >  1:729299/1‑100 (MQ=255)

                                                         tGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑CC‑CC‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAg                                       >  1:1033949/1‑100 (MQ=255)

                                                          ggCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑CC‑CC‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGCAGt                                      >  2:2273106/1‑100 (MQ=255)

                                                           gCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑CC‑CC‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTg                                     >  2:2188268/1‑100 (MQ=255)

                                                           gCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑CC‑CC‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTg                                     <  1:1408957/100‑1 (MQ=255)

                                                           gCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑CC‑CC‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTg                                     >  1:1383804/1‑100 (MQ=255)

                                                              gATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑CC‑CC‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCg                                  <  2:418152/100‑1 (MQ=255)

                                                              gATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑CC‑CC‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCg                                  <  1:1576061/100‑1 (MQ=255)

                                                               aTGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑CC‑CC‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCTTGAAGTGGCGa                                 <  2:1808184/100‑1 (MQ=255)

                                                                tGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑T‑‑CCTCC‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGaa                                <  1:658031/100‑1 (MQ=255)

                                                                   ttGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑CC‑CC‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGc                             <  1:276461/100‑1 (MQ=255)

                                                                    tGATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑CC‑CC‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGctttcgcc                            <  1:2177620/100‑8 (MQ=255)

                                                                     gATTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑CC‑CC‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCca                           <  1:498103/100‑1 (MQ=255)

                                                                      aTTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑CC‑CC‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGctttcgccac                          <  1:1516120/100‑10 (MQ=255)

                                                                      aTTTCTAACGTCGGTGAATTTTCTCGTT‑TC‑CC‑CC‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTAAAGTGGCGAAAGCcac                          <  1:1717096/100‑1 (MQ=255)

                                                                       tttCTAACGTCGGTGAATTTTCTCGTT‑TC‑CC‑CC‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCcaca                         <  1:1596787/100‑1 (MQ=255)

                                                                        ttCTAACGTCGGTGAATTTTCTCGTT‑TC‑CC‑CC‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAg                        >  1:60973/1‑100 (MQ=255)

                                                                         tCTAACGTCGGTGAATTTTCTCGTT‑TC‑CC‑CC‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGa                       <  1:1445570/100‑1 (MQ=255)

                                                                          cTAACGTCGGTGAATTTTCTCGTT‑TC‑CC‑CC‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGAt                      >  1:34433/1‑100 (MQ=255)

                                                                           tAACGTCGGTGAATTTTCTCGTT‑TC‑CC‑CC‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATg                     <  2:897968/100‑1 (MQ=255)

                                                                            aaCGTCGGTGAATTTTCTCGTT‑TC‑CC‑CC‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATgc                    <  1:2022734/100‑1 (MQ=255)

                                                                             aCGTCGGTGAATTTTCTCGTT‑TC‑CC‑CC‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATgcg                   <  1:2141226/100‑1 (MQ=255)

                                                                              cGTCGGTGAATTTTCTCGTT‑TC‑CC‑CC‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATgcgc                  >  1:119909/1‑100 (MQ=255)

                                                                              cGTCGGTGAATTTTCTCGTT‑TC‑CC‑CC‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATgcgc                  >  2:646969/1‑100 (MQ=255)

                                                                              cGTCGGTGAATTTTCTCGTT‑TC‑CC‑CC‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATgcg                   <  1:1321679/99‑1 (MQ=255)

                                                                               ggcgcgaaaTTTTCTCGTT‑T‑‑CC‑CCTCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGcc                 <  1:2173049/93‑1 (MQ=255)

                                                                               gTCGGTGAATTTTCTCGTT‑TC‑CC‑CC‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGcc                 <  2:1614484/100‑1 (MQ=255)

                                                                               gTCGGTGAATTTTCTCGTT‑TC‑CC‑CC‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGcc                 <  1:1653449/100‑1 (MQ=255)

                                                                               gTCGGTGAATTTTCTCGTT‑TC‑CC‑CC‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGcc                 >  1:1480058/1‑100 (MQ=255)

                                                                               gTCGGTGAATTTTCTCGTT‑TC‑CC‑CC‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGcc                 <  1:1217/100‑1 (MQ=255)

                                                                                tCGGTGAATTTTCTCGTT‑TC‑CC‑CC‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCt                >  2:1588371/1‑100 (MQ=255)

                                                                                  ggTGAATTTTCTCGTT‑TC‑CC‑CC‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGGAAGCCACAGATGCGCCtga              >  1:497090/1‑100 (MQ=255)

                                                                                  ggTGAATTTTCTCGTT‑TC‑CC‑CC‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtga              >  1:2092947/1‑100 (MQ=255)

                                                                                  ggTGAATTTTCTCGTT‑TC‑CC‑CC‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtg               >  2:448216/1‑99 (MQ=255)

                                                                                  ggTGAATTTTCTCGTT‑T‑‑CC‑CC‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgag             >  1:407845/1‑100 (MQ=255)

                                                                                     gAATTTTCTCGTT‑TC‑CC‑CC‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtg           <  1:22053/100‑1 (MQ=255)

                                                                                      aaTTTTCTCGTT‑TC‑CC‑CC‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtga          >  1:854918/1‑100 (MQ=255)

                                                                                       aTTTTCTCGTTCC‑‑CC‑CC‑CCCGCCAGCGCCAACCAGAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtgag         <  2:2027929/100‑1 (MQ=255)

                                                                                       aTTTTCTCGTT‑TCCCC‑CC‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCCTCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtga          <  2:1265130/100‑1 (MQ=255)

                                                                                             tcGTT‑TC‑CC‑CC‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACg   >  1:870674/1‑100 (MQ=255)

                                                                                              cGTT‑TC‑CC‑CC‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGc  >  2:2132149/1‑100 (MQ=255)

                                                                                              cGTT‑TC‑CC‑CC‑CCCGCCAGCGCCAAAAATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGAGAAAGCCACAGATGCGCCTGAATGAGTGTACGc  >  2:1004134/1‑100 (MQ=255)

GCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑T‑‑CC‑CC‑CCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGC  >  NC_000913/3379161‑3379353

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