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	94.9%	332.7 / 12.3	98	coding (739/1128 nt)	zapE	divisome ATPase
Reads supporting (aligned to +/- strand): ref base . (4/1); new base C (56/37); total (60/38)
Fisher's exact test for biased strand distribution p-value = 6.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.26e-01

Only 100 of 104 total aligned reads displayed.

GCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑C‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCG‑CCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGC  >  NC_000913/3379166‑3379353

gCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑c‑ccccc                                                                                           <  2:2047806/100‑1 (MQ=255)

    agtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑c‑cccccccgc                                                                                       >  2:814587/1‑97 (MQ=255)

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

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

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

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

       aaaaGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCCCCGCCAg                                                                                     >  2:1035313/1‑100 (MQ=255)

        aaaGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑C‑CCCCCCGCCAgcg                                                                                   >  2:1209719/1‑100 (MQ=255)

          agagACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCCCCGCCAgcgc                                                                                  >  1:821348/1‑100 (MQ=255)

             gaCCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCCCCGCCAGCGCCaa                                                                               >  1:1847039/1‑100 (MQ=255)

             gaCCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCCCCGCCAGCGCCaa                                                                               >  1:1587271/1‑100 (MQ=255)

             gaCCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCCCCGCCAGCGCCaa                                                                               >  2:1518227/1‑100 (MQ=255)

               ccGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCCCCGCCAGCGCCAAcc                                                                             >  2:1327491/1‑100 (MQ=255)

               ccGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCCCCGCCAGCGCCAAcc                                                                             <  2:868372/100‑1 (MQ=255)

                cGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCCCCGCCAGCGCCAACCa                                                                            >  2:660196/1‑100 (MQ=255)

                cGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCCCCGCCAGCGCCAACCa                                                                            >  1:1035618/1‑100 (MQ=255)

                 gCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCCCCGCCCGCGCCAACCAt                                                                           >  2:1558893/1‑100 (MQ=255)

                 gCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCCCCGCCAGCGCCAACCAt                                                                           >  1:770373/1‑100 (MQ=255)

                 gCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCCCCGCCAGCGCCAACCAt                                                                           >  1:2163412/1‑100 (MQ=255)

                  ccAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCCCCGCCAGCGCCAACCATa                                                                          >  1:516090/1‑100 (MQ=255)

                   cAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCCCCGCCAGCGCCAACCATAg                                                                         <  2:1299722/100‑1 (MQ=255)

                    aGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCCCCGCCAGCGCCAACCATAGt                                                                        <  1:1136805/100‑1 (MQ=255)

                     gCGTCTGGTTCTCGACGCCCATTGCTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCCCCGCCAGCGCCAACCATAGtt                                                                       >  2:950129/1‑100 (MQ=255)

                       gTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCCCCGCCAGCGCCAACCATAGTTTa                                                                     <  1:423065/100‑1 (MQ=255)

                         cTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCCCCGCCAGCGCCAACCATAGTTTATc                                                                   <  2:1748317/100‑1 (MQ=255)

                         cTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCCCCGCCAGCGCCAACCATAGTTTATc                                                                   <  1:418452/100‑1 (MQ=255)

                         cTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑C‑CCCCCCGCCAGCGCCAACCATAGTTTATcc                                                                  >  2:168070/1‑100 (MQ=255)

                          tGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCCCCGCCAGCGCCAACCATAGTTTATcc                                                                  >  2:601251/1‑100 (MQ=255)

                           ggTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCCCCGCCAGCTCCAACCATAGttn                                                                      <  1:1379287/95‑2 (MQ=255)

                           ggTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCCCCGCCAGCGCCAACCATAGTTTATCCa                                                                 >  2:1748686/1‑100 (MQ=255)

                            gTTTTCGACGTCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑CTCCCCTCTCCAGCGCCAACCATAGTTTATCAAt                                                                >  1:1326291/1‑100 (MQ=255)

                            gTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCCCCGCCAGCGCCAACCATAGTTTATCCAt                                                                >  2:1363550/1‑100 (MQ=255)

                              tctcGCCGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCCCCGCCAGCGCCAACCATAGTTTATCCAttt                                                              >  1:2098792/1‑100 (MQ=255)

                              nctcGACCCCCATTGTTGCTAATGGCCGATGGTTGATATCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCCCCGCCAGCGCCAACCATAGTTTCTCCAttt                                                              >  1:1195635/2‑100 (MQ=255)

                               ttccgatctcATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTg                                                             <  1:554890/91‑1 (MQ=255)

                                tcGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTgc                                                            <  1:162166/100‑1 (MQ=255)

                                 cGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTgcg                                                           >  1:356697/1‑100 (MQ=255)

                                 cGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTgcg                                                           >  1:992563/1‑100 (MQ=255)

                                 cGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTgcg                                                           >  2:1719550/1‑100 (MQ=255)

                                 cGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTgc                                                            >  1:1164237/1‑99 (MQ=255)

                                     gCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCG‑CCCg                                                      >  1:2021608/1‑100 (MQ=255)

                                      cccATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCG‑CCCgg                                                     >  2:628471/1‑100 (MQ=255)

                                      cccATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑C‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCG‑CCCggg                                                    >  2:612777/1‑100 (MQ=255)

                                       ccATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCG‑CCCggg                                                    >  2:398022/1‑100 (MQ=255)

                                        cATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTgagatcggaa                                                    >  2:554891/1‑91 (MQ=255)

                                        cATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCG‑CCCGGGt                                                   >  2:250847/1‑100 (MQ=255)

                                        cATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCCCCGCAAGCGCCAACCCTAGTTTATCCATTTGCG‑CCCGGGt                                                   >  1:768432/1‑100 (MQ=255)

                                         aTTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCG‑CCCGGGtt                                                  <  1:142764/100‑1 (MQ=255)

                                          ttgttgCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCG‑CCCGGGttt                                                 <  2:2200443/100‑1 (MQ=255)

                                          ttgttgCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCG‑CCCGGGttt                                                 >  1:1691572/1‑100 (MQ=255)

                                          ttgttgCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGAG‑CCCGGGttt                                                 >  2:2116772/1‑100 (MQ=255)

                                            gttgCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCG‑CCCGGGTTTCa                                               >  1:136328/1‑100 (MQ=255)

                                             ttgCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCCGGGTTTCa                                               <  2:712/100‑1 (MQ=255)

                                             ttgCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCG‑CCCGGGTTTCAt                                              <  1:1046018/100‑1 (MQ=255)

                                             ttgCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCG‑CCCGGGTTTCAt                                              <  1:1510454/100‑1 (MQ=255)

                                               gCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCG‑CCCGGGTTTCATCg                                            <  2:58328/100‑1 (MQ=255)

                                               gCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCG‑CCCGGGTTTCATCg                                            >  2:197561/1‑100 (MQ=255)

                                                cTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCG‑CCCGGGTTTCATCGt                                           >  1:342388/1‑100 (MQ=255)

                                                 tAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCG‑CCCGGGTTTCATCGTg                                          <  2:1803338/100‑1 (MQ=255)

                                                 tAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCG‑CCCGGGTTTCATCGTg                                          >  1:1248850/1‑100 (MQ=255)

                                                  aaTGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCG‑CCCGGGTTTCATCGTGa                                         >  2:1065612/1‑100 (MQ=255)

                                                  aaTGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCG‑CCCGGGTTTCATCGTGa                                         <  1:474372/100‑1 (MQ=255)

                                                    tGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCG‑CCCGGGTTTCATCGTGAAg                                       <  1:466411/100‑1 (MQ=255)

                                                    tGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCG‑CCCGGGTTTCATCGTGAAg                                       >  2:877506/1‑100 (MQ=255)

                                                    tGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCG‑CCCGGGTTTCATCGTGAAg                                       >  2:1783531/1‑100 (MQ=255)

                                                     ggCCGGTGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTCCC‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCG‑CCCGGGTTTCATCGTGa                                         <  1:150107/98‑1 (MQ=255)

                                                     ggCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCG‑CCCGGGTTTCATCGTGAAGt                                      >  2:625382/1‑100 (MQ=255)

                                                     ggCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCG‑CCCGGGTTTCATCGTGAAGt                                      <  1:1035312/100‑1 (MQ=255)

                                                        cGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCG‑CCCGGGTTTCATCGTGAAGTGGc                                   >  1:1604817/1‑100 (MQ=255)

                                                         gATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCG‑CCCGGGTTTCATCGTGAAGTGGCg                                  >  2:2095638/1‑100 (MQ=255)

                                                         gATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCG‑CCCGGGTTTCATCGTGAAGTGGCg                                  <  1:628470/100‑1 (MQ=255)

                                                          aTGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCG‑CCCGGGTTTCATCGTGAAGTGGCGa                                 <  1:950128/100‑1 (MQ=255)

                                                             gTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCG‑CCCGGGTTTCATCGTGAAGTGGCGAAAg                              >  2:513706/1‑100 (MQ=255)

                                                             gTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCG‑CCCGGGTTTCATCGTGAAGTGGCGAAAg                              <  2:659397/100‑1 (MQ=255)

                                                             gTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCG‑CCCGGGTTTCATCGTGAAGTGGCGAAAg                              >  1:1719786/1‑100 (MQ=255)

                                                             gTTGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCG‑CCCGCGTTTCATCGTGAAGTGGCGAAAg                              >  1:1768565/1‑100 (MQ=255)

                                                               tGATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCG‑CCCGGGTTTCATCGTGAAGTGGCGAAAGcc                            <  2:1587280/100‑1 (MQ=255)

                                                                gATTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCG‑CCCGGGTTTCATCGTGAAGTGGCGAAAGCca                           >  1:628845/1‑100 (MQ=255)

                                                                 aTTTCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCG‑CCCGGGTTTCATCGTGAAGTGGCGAAAGCcac                          <  2:821349/100‑1 (MQ=255)

                                                                    tCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCG‑CCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGa                       <  2:248832/100‑1 (MQ=255)

                                                                    tCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCG‑CCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGa                       >  1:2133406/1‑100 (MQ=255)

                                                                    tCTAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCG‑CCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGa                       >  2:2104147/1‑100 (MQ=255)

                                                                      tAACGTCGGTGAATTTTCTCGTTTC‑C‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCG‑CCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATg                     >  2:1889734/1‑100 (MQ=255)

                                                                      tAACGTCGGTGAATTCTCTCGTTTC‑C‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCG‑CCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATg                     <  2:2191005/100‑1 (MQ=255)

                                                                       aaCGTCGGTGAATTTTCTCGTTTC‑C‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCG‑CCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATgc                    >  1:1439370/1‑100 (MQ=255)

                                                                          gTCGGTGAATTTTCTCGTTTC‑C‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCG‑CCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGcc                 <  1:2228293/100‑1 (MQ=255)

                                                                          gTCGGTGAATTTTCTCGTTT‑‑C‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCG‑CCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCt                <  1:612776/100‑1 (MQ=255)

                                                                            cGGTGAATTTTCTCGTTTC‑C‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCG‑CCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtg               <  2:2149926/100‑1 (MQ=255)

                                                                             ggTGAATTTTCTCGTTTC‑C‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCG‑CCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtga              >  1:1523579/1‑100 (MQ=255)

                                                                              gTGAATTTTCTCGTTTC‑C‑CCCCCCGCCAGCGCCAACCATAGTTTTTCCATTTGCG‑CCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgag             >  1:1814937/1‑100 (MQ=255)

                                                                               tGAATTTTCTCGTTTC‑C‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCG‑CCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagt            >  1:2145867/1‑100 (MQ=255)

                                                                                  aTTTTCTCGTTTC‑C‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCG‑CCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtgag         >  2:1904625/1‑100 (MQ=255)

                                                                                  aTTTTCTCGTTTC‑C‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCG‑CCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtgag         <  1:510319/100‑1 (MQ=255)

                                                                                  aTTTTCTCGTTTC‑C‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCG‑CCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtgag         <  2:33493/100‑1 (MQ=255)

                                                                                    tttCTCGTTTC‑C‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCG‑CCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGagtgagtg       <  2:1115405/100‑1 (MQ=255)

                                                                                    tttCTCGTTTC‑C‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCG‑CCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGagtgagtg       <  2:1762330/100‑1 (MQ=255)

                                                                                      tctcGTTTC‑C‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCG‑CCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTa     >  2:2059826/1‑100 (MQ=255)

                                                                                       ctcGTTTC‑C‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCG‑CCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTAc    >  2:1180301/1‑100 (MQ=255)

                                                                                         cGTTTC‑C‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCG‑CCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGc  <  1:399558/100‑1 (MQ=255)

                                                                                         cGTTTC‑C‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCG‑CCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGc  >  2:16591/1‑100 (MQ=255)

GCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTT‑‑C‑CCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCG‑CCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGC  >  NC_000913/3379166‑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: `‑`