BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,379,266	(C)_7→8	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	98.0%	360.7 / 0.8	99	coding (739/1128 nt)	zapE	divisome ATPase
Reads supporting (aligned to +/- strand): ref base . (1/1); new base C (45/52); total (46/53)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.70e-01

Only 100 of 112 total aligned reads displayed.

CGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑C‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGC  >  NC_000913/3379160‑3379353

cGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGttt                                                                                                  >  2:926041/1‑100 (MQ=255)

  cacaGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑c‑c                                                                                              <  1:1060348/100‑1 (MQ=255)

  cacaGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑c‑c                                                                                              <  2:1153306/100‑1 (MQ=255)

   acaGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑c‑cc                                                                                             >  2:1764153/1‑100 (MQ=255)

   acaGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑c‑cc                                                                                             >  1:1353265/1‑100 (MQ=255)

    caGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑c‑ccc                                                                                            >  1:1230418/1‑100 (MQ=255)

    caGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑c‑ccc                                                                                            <  2:758022/100‑1 (MQ=255)

     aGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑c‑cccc                                                                                           >  1:2363955/1‑100 (MQ=255)

          agtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑c‑cccccccg                                                                                       >  2:2750731/1‑97 (MQ=255)

           gtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑c‑cccccccgcc                                                                                     >  1:1575453/1‑96 (MQ=255)

           gtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑c‑cccccccgcc                                                                                     <  1:1764164/100‑5 (MQ=255)

            taAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑c‑cccccccgcca                                                                                    >  1:584996/1‑95 (MQ=255)

             aaaaGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑CCCCCCGCCAg                                                                                    >  1:1331129/1‑100 (MQ=255)

              aaaGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑CCCCCCGCCAgc                                                                                   >  1:2356172/1‑100 (MQ=255)

              aaaGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑CCCCCCGCCAgc                                                                                   >  2:773205/1‑100 (MQ=255)

              aaaGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑CCCCCCGCCAgc                                                                                   >  1:2120321/1‑100 (MQ=255)

                agagACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑CCCCCCGCCAgcgc                                                                                 >  1:410418/1‑100 (MQ=255)

                 gagaCCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑CCCCCCGCCAGCGcc                                                                                <  1:1914865/100‑1 (MQ=255)

                 gagaCCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑CCCCCCGCCAGCGcc                                                                                <  2:1719571/100‑1 (MQ=255)

                  agaCCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑CCCCCCGCCAGCGCCa                                                                               <  2:1537219/100‑1 (MQ=255)

                   gaCCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑CCCCCCGCCAGCGCCaa                                                                              <  1:1883238/100‑1 (MQ=255)

                    aCCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑CCCCCCGCCAGCGCCAAc                                                                             <  2:2233199/100‑1 (MQ=255)

                    aCCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑CCCCCCGCCAGCGCCAAc                                                                             <  2:27897/100‑1 (MQ=255)

                     ccgcccgcGTCTGGGTCTCGACTCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑CCCCCCGCCAGCGCCAAc                                                                             <  2:2363941/99‑1 (MQ=255)

                         cAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑CCCCCCGCCAGCGCCAACCATAg                                                                        <  1:2434567/100‑1 (MQ=255)

                         cAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑CCCCCCGCCAGCGCCAACCATAg                                                                        <  1:1611307/100‑1 (MQ=255)

                         cAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑CCCCCCGCCAGCGCCAACCATAg                                                                        <  1:2360743/100‑1 (MQ=255)

                          aGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑CCCCCCGCCAGCGCCAACCATAGt                                                                       <  2:2562562/100‑1 (MQ=255)

                           gCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑CCCCCCGCCAGCGCCAACCATAGtt                                                                      <  1:311259/100‑1 (MQ=255)

                            cGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑CCCCCCGCCAGCGCCAAcagatcgga                                                                     >  1:2637512/1‑92 (MQ=255)

                             gTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑CCCCCCGCCAGCGCCAACCATAGTTTa                                                                    >  1:867678/1‑100 (MQ=255)

                             gTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑CCCCCCGCCAGCGCCAACCATAGTTTa                                                                    <  2:1421332/100‑1 (MQ=255)

                              tCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑CCCCCCGCCAGCGCCCACCATAGTTTAt                                                                   >  2:1959126/1‑100 (MQ=255)

                              attGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑CCCCCCGCCAGCGCCAACCATAGTTTAt                                                                   <  2:6706/98‑1 (MQ=255)

                                 ggTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑CCCCCCGCCAGCGCCAACCATAGTTTATCCa                                                                >  2:701324/1‑100 (MQ=255)

                                 ggTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑CCCCCCGCCAGCGCCAACCATACTTTATCCa                                                                >  1:449145/1‑100 (MQ=255)

                                  gTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑CCCCCCGCCAGCGCCAACCATAGTTTATCCAt                                                               <  2:1398030/100‑1 (MQ=255)

                                    tctcGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑CCCCCCGCCAGCGCCAACCATAGTTTATCCAttt                                                             <  2:1564559/100‑1 (MQ=255)

                                    tctcGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑CCCCCCGCCAGCGCCAACCATAGTTTATCCAttt                                                             <  2:932819/100‑1 (MQ=255)

                                    tctcGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑CCCCCCGCCAGCGCCAACCATAGTTTATCCAttt                                                             <  1:2248990/100‑1 (MQ=255)

                                      tcGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTgc                                                           >  2:2407861/1‑100 (MQ=255)

                                       cGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTgcg                                                          >  1:2526038/1‑100 (MQ=255)

                                          cGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGccc                                                       >  1:1685069/1‑100 (MQ=255)

                                            cccATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCgg                                                     <  1:2951936/100‑1 (MQ=255)

                                             ccATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCggg                                                    <  2:2883797/100‑1 (MQ=255)

                                               aTTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGtt                                                  <  2:1575444/100‑1 (MQ=255)

                                               aTTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑C‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGttt                                                 <  1:884816/100‑1 (MQ=255)

                                                ttgttgCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGttt                                                 <  2:1046005/100‑1 (MQ=255)

                                                ttgttgCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGttt                                                 <  1:2074221/100‑1 (MQ=255)

                                                 tgttgCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTc                                                <  2:1353258/100‑1 (MQ=255)

                                                  gttgCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCa                                               >  1:1094075/1‑100 (MQ=255)

                                                   ttgCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCAt                                              <  1:2307122/100‑1 (MQ=255)

                                                    tgCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATc                                             >  1:589527/1‑100 (MQ=255)

                                                    tgCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATc                                             <  1:2880874/100‑1 (MQ=255)

                                                    tgCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATc                                             <  1:1920585/100‑1 (MQ=255)

                                                     gCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCg                                            <  1:218223/100‑1 (MQ=255)

                                                     gCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGG‑TTCATCGt                                           >  1:1369410/1‑100 (MQ=255)

                                                      cTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGt                                           >  1:212804/1‑100 (MQ=255)

                                                      cTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGt                                           >  1:2156395/1‑100 (MQ=255)

                                                       tAATGGCCGATGGTTGATTTCTTACGTCGGTGAATTTTCTCGTTTCC‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTg                                          <  2:1331122/100‑1 (MQ=255)

                                                        aaTGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGa                                         <  1:1049621/100‑1 (MQ=255)

                                                        aaTGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGa                                         >  1:289126/1‑100 (MQ=255)

                                                        aaTGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGa                                         <  2:290531/100‑1 (MQ=255)

                                                        aaTGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGa                                         <  2:1230412/100‑1 (MQ=255)

                                                            gCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTg                                     >  1:1995771/1‑100 (MQ=255)

                                                            gCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTg                                     >  1:635535/1‑100 (MQ=255)

                                                             ccGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTgg                                    >  2:2020497/1‑100 (MQ=255)

                                                              cGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGc                                   <  1:1606823/100‑1 (MQ=255)

                                                                 tGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGaa                                <  1:2704137/100‑1 (MQ=255)

                                                                   gTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑CTCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAg                              >  2:2886940/1‑100 (MQ=255)

                                                                    ttGATTTCTAACGTCGGTGAATTTTCTCGTTTCC‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGc                             <  1:1785711/100‑1 (MQ=255)

                                                                        tttCTAACGTCGGTGAATTTTCTCGTTTCC‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCcaca                         >  2:1900380/1‑100 (MQ=255)

                                                                        tttCTAACGTCGGTGAATTTTCTCGTTTCC‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCcaca                         <  2:559900/100‑1 (MQ=255)

                                                                         ttCTAACGTCGGTGAATTTTCTCGTTTCC‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAg                        >  2:2673845/1‑100 (MQ=255)

                                                                          tCTAACGTCGGTGAATTTTCTCGTTTCC‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGa                       <  2:1995759/100‑1 (MQ=255)

                                                                           cTAACGTCGGTGAATTTTCTCGTTTCC‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGAt                      <  2:2120309/100‑1 (MQ=255)

                                                                             aaCGTCGGTGAATTTTCTCGTTTCC‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATgc                    >  2:2762205/1‑100 (MQ=255)

                                                                              aCGTCGGTGAATTTTCTCGTTTCC‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATgcg                   <  2:1569532/100‑1 (MQ=255)

                                                                               cGTCGGTGAATTTTCTCGTTTCC‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATgcgc                  >  2:2703915/1‑100 (MQ=255)

                                                                                gTCGGTGAATTTTCTCGTTTCC‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGcc                 <  1:701327/100‑1 (MQ=255)

                                                                                  cGGTGAATTTTCTCGTTTCC‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtg               >  2:1983424/1‑100 (MQ=255)

                                                                                     tGAATTTTCTCGTTTCC‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagt            <  1:1496606/100‑1 (MQ=255)

                                                                                      gAATTTTCTCGTTTCC‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtg           >  2:2381146/1‑100 (MQ=255)

                                                                                        aTTTTCTCGTTTCC‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtgag         >  1:1068523/1‑100 (MQ=255)

                                                                                        aTTTTCTCGTTTCC‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtgag         >  1:2802941/1‑100 (MQ=255)

                                                                                          tttCTCGTTTCC‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGagtgagtg       >  2:676837/1‑100 (MQ=255)

                                                                                           ttCTCGTTTCC‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAgtgt      >  2:1288915/1‑100 (MQ=255)

                                                                                            tctcGTTTCC‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTa     <  2:1702091/100‑1 (MQ=255)

                                                                                            tctcGTTTCC‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTa     >  2:2264781/1‑100 (MQ=255)

                                                                                            tctcGTTTCC‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTa     <  2:2131028/100‑1 (MQ=255)

                                                                                            tctcGTTTCC‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTa     <  2:1499694/100‑1 (MQ=255)

                                                                                             ctcGTTTCC‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTAc    >  1:1678942/1‑100 (MQ=255)

                                                                                              tcGTTTCC‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACg   <  2:2351351/100‑1 (MQ=255)

                                                                                               cGTTTCC‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGc  >  2:1012373/1‑100 (MQ=255)

                                                                                               cGTTTCC‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGc  >  2:2398361/1‑100 (MQ=255)

                                                                                               cGTTTCC‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGc  <  2:2356158/100‑1 (MQ=255)

                                                                                               cGTTTCC‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGc  >  1:1972351/1‑100 (MQ=255)

                                                                                               cGTTTCC‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGc  >  2:2044824/1‑100 (MQ=255)

                                                                                               cGTTTCC‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGc  >  1:1985202/1‑100 (MQ=255)

                                                                                               cGTTTCC‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGATTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGc  >  1:2219990/1‑100 (MQ=255)

CGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑C‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGC  >  NC_000913/3379160‑3379353

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