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	93.8%	442.3 / 22.1	130	coding (739/1128 nt)	zapE	divisome ATPase
Reads supporting (aligned to +/- strand): ref base . (6/2); new base C (59/62); total (66/64)
Fisher's exact test for biased strand distribution p-value = 2.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 148 total aligned reads displayed.

CGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCG‑TTT‑‑CC‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAG  >  NC_000913/3379160‑3379355

cGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCG‑ttt                                                                                                     >  1:2111710/1‑100 (MQ=255)

 gCACAGCGTAGTAAAACAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCG‑TTT‑‑c                                                                                                  <  2:1919386/100‑1 (MQ=255)

  cacaGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCG‑TTT‑‑cc                                                                                                 <  2:838554/100‑1 (MQ=255)

  cacaGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCG‑TTT‑‑cc                                                                                                 >  1:2429129/1‑100 (MQ=255)

   acaGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCG‑TTT‑‑cc‑c                                                                                               <  1:404043/100‑1 (MQ=255)

    caGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTACTGGCCGATGGTTGCTTTCTCACGTCGGTGAATTTTCTCG‑TTT‑‑cc‑cc                                                                                              >  2:1883765/1‑100 (MQ=255)

    caGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCG‑TTT‑‑cc‑cc                                                                                              >  2:2200849/1‑100 (MQ=255)

     aGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCG‑TTT‑‑cc‑ccc                                                                                             <  1:119230/100‑1 (MQ=255)

     aGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCG‑TTT‑‑cc‑ccc                                                                                             <  1:1587863/100‑1 (MQ=255)

           gtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCG‑TTT‑‑cc‑ccccccgcc                                                                                       >  1:1009752/1‑96 (MQ=255)

           gtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCG‑TTT‑‑cc‑ccccccccc                                                                                       >  2:2826294/1‑96 (MQ=255)

            taAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCG‑TTT‑‑cc‑cccccc                                                                                          >  2:683161/1‑95 (MQ=255)

             aaaaGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCG‑TTT‑‑cc‑ccccccccc                                                                                       >  2:962924/1‑94 (MQ=255)

                agagACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCG‑TTTC‑CC‑CCCCCGCCAgcgc                                                                                   <  2:923097/100‑1 (MQ=255)

                  agaCCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCG‑TTTC‑CC‑CCCCCGCCAGCGCCa                                                                                 <  2:1988433/100‑1 (MQ=255)

                   gaCCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCG‑TTTC‑CC‑CCCCCGCCAGCGCCaa                                                                                >  2:2861872/1‑100 (MQ=255)

                    aCCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCG‑TTTC‑CC‑CCCCCGCCAGCGCCAAc                                                                               >  1:1389850/1‑100 (MQ=255)

                     ccGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTTT‑‑CC‑CCCCCGCCAGCGCCAAcc                                                                              >  2:882352/1‑100 (MQ=255)

                     ccGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCG‑TTTC‑CC‑CCCCCGCCAGCGCCAAcc                                                                              <  1:387778/100‑1 (MQ=255)

                     ccGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCG‑TTTC‑CC‑CCCCCGCCAGCGCCAAcc                                                                              >  2:2930089/1‑100 (MQ=255)

                      cGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCG‑TTTCCCC‑CCCCCGCCAGCGCCAAcc                                                                              >  2:2762966/1‑100 (MQ=255)

                      cGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCG‑TTTC‑CC‑CCCCCGCCAGCGCCAACCa                                                                             >  1:1007430/1‑100 (MQ=255)

                      cGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCG‑TTTC‑CC‑CCCCCGCCAGCGCCAACCa                                                                             <  1:2861886/100‑1 (MQ=255)

                      cGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCG‑TTTC‑CC‑CCCCCGCCAGCGCCAACCa                                                                             >  1:2215609/1‑100 (MQ=255)

                       gCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCG‑TTTC‑CC‑CCCCCGCCAGCGCCAACCAt                                                                            >  1:591117/1‑100 (MQ=255)

                        ncAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCG‑TTTC‑CC‑CCCCCGCCAGCGCCAACCATa                                                                           >  2:1848959/2‑100 (MQ=255)

                         cAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCG‑TTTC‑CC‑CCCCCGCCAGCGCCAACCATAg                                                                          <  1:2142274/100‑1 (MQ=255)

                              tCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCG‑TTTC‑CC‑CCCCCGCCAGCGCCAACCATAGTTTAt                                                                     <  1:520033/100‑1 (MQ=255)

                              tCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCG‑TTTC‑CC‑CCCCCGCCAGCGCCAACCATAGTTTAt                                                                     >  1:1878313/1‑100 (MQ=255)

                              tCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCG‑TTTC‑CC‑CCCCCGCCAGCGCCAACCATAGTTTAt                                                                     >  1:1878285/1‑100 (MQ=255)

                               cTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCG‑TTTC‑CC‑CCCCCGCCAGCGCCAACCATAGTTTATc                                                                    <  1:1378204/100‑1 (MQ=255)

                                tGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCG‑TTTC‑CC‑CCCCCGCCAGCGCCAACCATAGTTTATcc                                                                   <  2:1365665/100‑1 (MQ=255)

                                tGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCG‑TTTC‑CC‑CCCCCGCCAGCGCCAACCATAGTTTATcc                                                                   >  1:2037524/1‑100 (MQ=255)

                                tGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCACG‑TTTG‑CC‑CCCCCGCCAGCGCCAACCATAGTTTATcc                                                                   >  2:2202480/1‑100 (MQ=255)

                                    tctcGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCG‑TTTC‑CC‑CCCCCGCCAGCGCCAACCATAGTTTATCCAttt                                                               >  1:2672815/1‑100 (MQ=255)

                                    tctcGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCG‑TTTC‑CC‑CCCCCGCCAGCGCCAACCATAGTTTATCCAttt                                                               <  1:1523964/100‑1 (MQ=255)

                                     ctcGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCG‑TTTC‑CC‑CCCCCGCCCGCGCCAACCATAGTTTATCCAttt                                                               >  1:2447350/1‑99 (MQ=255)

                                     ctcGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCG‑TTTC‑CC‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTg                                                              <  1:2643724/100‑1 (MQ=255)

                                     ctcGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCG‑TTTC‑CC‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTg                                                              <  1:3241766/100‑1 (MQ=255)

                                       cGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCG‑TTTC‑CC‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTgcg                                                            >  1:2477303/1‑100 (MQ=255)

                                           gCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCG‑TTTC‑CC‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGGGCCCg                                                        >  1:1467726/1‑100 (MQ=255)

                                            cccATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCG‑TTTC‑CC‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCgg                                                       >  1:3256800/1‑100 (MQ=255)

                                            cccATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCG‑TTTC‑CC‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCgg                                                       <  1:2492265/100‑1 (MQ=255)

                                             ccATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCG‑TTTC‑CC‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCggg                                                      <  1:2691844/100‑1 (MQ=255)

                                              cATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCG‑TTTC‑CC‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGt                                                     >  1:2808897/1‑100 (MQ=255)

                                               aTTGTTGCTAATGGCCGATGGTTGATTTCTTACGTCGGTGAATTTTCTCG‑TTTC‑CC‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGAGtt                                                    <  2:73318/100‑1 (MQ=255)

                                                  gttgCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCG‑TTTC‑CC‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCa                                                 >  2:1250767/1‑100 (MQ=255)

                                                  gttgCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCG‑TTTC‑CC‑CCCCCGCCAGCGCCAACCAAAGTTTATCCATTTGCGCCCGGGTTTCa                                                 >  1:1047887/1‑100 (MQ=255)

                                                     gCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCG‑TTTC‑CC‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCg                                              <  1:1856458/100‑1 (MQ=255)

                                                      cTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCG‑TTTC‑CC‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGt                                             <  1:333287/100‑1 (MQ=255)

                                                      cTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCG‑TTTC‑CC‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGt                                             >  2:2768254/1‑100 (MQ=255)

                                                       tAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCG‑TTTC‑CC‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTg                                            <  2:155714/100‑1 (MQ=255)

                                                        aaTGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCG‑TTTC‑CC‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGa                                           <  2:2136205/100‑1 (MQ=255)

                                                          tGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCG‑TTTC‑CC‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAg                                         <  1:1141685/100‑1 (MQ=255)

                                                           ggCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCG‑TTTC‑CC‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGt                                        >  2:2101891/1‑100 (MQ=255)

                                                           ggCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCG‑TTTC‑CC‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGt                                        >  2:1349214/1‑100 (MQ=255)

                                                             ccGATGGTTGATTTCTAACGTCGGTGAATTTTCTCG‑TTTC‑CC‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTgg                                      >  2:2336021/1‑100 (MQ=255)

                                                             ccGATGGTTGATTTCTAACGTCGGTGAATTTTCTCG‑TTTC‑CC‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTgg                                      >  1:1710106/1‑100 (MQ=255)

                                                             ccGATGGTTGATTTCTAACGTCGGTGAATTTTCTCG‑TTTC‑CC‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTgg                                      >  2:53916/1‑100 (MQ=255)

                                                             ccGATGGTTGATTTCTAACGTCGGTGAATTTTCTCG‑TTTC‑CC‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTgg                                      >  1:2925045/1‑100 (MQ=255)

                                                              cGATGGTTGATTTCTAACGTCGGTGAATTTTCTCG‑TTTC‑CC‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTgg                                      >  2:126766/1‑99 (MQ=255)

                                                              cGATGGTTGATTTCTAACGTCGGTGAATTTTCTCG‑TTTC‑CC‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGc                                     <  1:986551/100‑1 (MQ=255)

                                                              cGATGGTTGATTTCTAACGTCGGTGAATTTTCTCG‑TTTC‑CC‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGc                                     <  1:2486013/100‑1 (MQ=255)

                                                                aTGGTTGATTTCTAACGTCGGTGAATTTTCTCG‑TTTCCCC‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCg                                    <  2:522580/100‑1 (MQ=255)

                                                                   gTTGATTTCTAACGTCGGTGAATTTTCTCG‑TTTC‑CC‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAg                                <  2:1499029/100‑1 (MQ=255)

                                                                    ttGATTTCTAACGTCGGTGAATTTTCTCG‑TTTC‑CC‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGc                               >  2:1383264/1‑100 (MQ=255)

                                                                    ttGATTTCTAACGTCGGTGAATTTTCTCG‑TTTC‑CC‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGc                               <  1:3069225/100‑1 (MQ=255)

                                                                    ttGATTTCTAACGTCGGTGAATTTTCTCG‑TTTC‑CC‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTCCGCCCGGGTTTCATCGTGAAGTGGCGAAAGc                               <  1:2200752/100‑1 (MQ=255)

                                                                    ttGATTTCTAACGTCGGTGAATTTTCTCG‑TTT‑‑CCACCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGaa                                  >  2:499107/1‑97 (MQ=255)

                                                                     tGATTTCTAACGTCGGTGAATTTTCTCG‑TTTC‑CC‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGcc                              >  2:1540827/1‑100 (MQ=255)

                                                                     tGATTTCTAACGTCGGTGAATTTTCTCG‑TTTC‑CC‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGcc                              <  1:509429/100‑1 (MQ=255)

                                                                      gATTTCTAACGTCGGTGAATTTTCTCG‑TTTC‑CC‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCca                             >  1:1675000/1‑100 (MQ=255)

                                                                       aTTTCTAACGTCGGTGAATTTTCTCG‑TTTC‑CC‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCcac                            <  2:1347149/100‑1 (MQ=255)

                                                                        tttCTAACGTCGGTGAATTTTCTCG‑TTTC‑CC‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCcaca                           <  2:1463048/100‑1 (MQ=255)

                                                                        tttCTAACGTCGGTGAATTTTCTCG‑TTTC‑CC‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCcaca                           <  1:2200858/100‑1 (MQ=255)

                                                                        tttCTAACGTCGGTGAATTTTCTCG‑TTTC‑CC‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCcaca                           <  1:242607/100‑1 (MQ=255)

                                                                            tAACGTCGGTGAATTTTCTCG‑TTTC‑CC‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATg                       <  1:74477/100‑1 (MQ=255)

                                                                             aaCGTCGGTGAATTTTCTCG‑TTTC‑CC‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATgc                      >  2:2908893/1‑100 (MQ=255)

                                                                              aCGTCGGTGAATTTTCTCG‑TTTC‑CC‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATgcg                     <  1:2854198/100‑1 (MQ=255)

                                                                               cGTCGGTGAATTTTCTCG‑TTTC‑CC‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATgcgc                    >  2:2719727/1‑100 (MQ=255)

                                                                                gTCGGTGAATTTTCTCG‑TTTC‑CC‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGcc                   <  2:1009750/100‑1 (MQ=255)

                                                                                  cGGTGAATTTTCTCG‑TTTC‑CC‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtg                 <  2:3232257/100‑1 (MQ=255)

                                                                                    gTGAATTTTCTCG‑TTTC‑CC‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgag               <  2:2309234/100‑1 (MQ=255)

                                                                                    gTGAATTTTCTCG‑TTTC‑CC‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgag               >  1:1584769/1‑100 (MQ=255)

                                                                                    gTGAATTTTCTCG‑TTTC‑CC‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgag               <  2:241718/100‑1 (MQ=255)

                                                                                    gTGAATTTTCTCG‑TTTC‑CC‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgag               <  2:241643/100‑1 (MQ=255)

                                                                                      gAATTTTCTCG‑TTTCCCC‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagt              <  2:24023/100‑1 (MQ=255)

                                                                                      gAATTTTCGCG‑TTTC‑CC‑CCCCCGCCAGCGCCAACCAGAGTTTATCCATTTGCGCCCGGGTGTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtg             <  1:3166542/100‑1 (MQ=255)

                                                                                       aaTTTTCTCG‑TTTC‑CC‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtg             <  1:1506521/99‑1 (MQ=255)

                                                                                        aTTTTCTCG‑TTTC‑CC‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtgag           >  2:340540/1‑100 (MQ=255)

                                                                                        aTTTTCTCG‑TTTC‑CC‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtgag           <  2:49887/100‑1 (MQ=255)

                                                                                        aTTTTCTCG‑TTTC‑CC‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtgag           <  1:820109/100‑1 (MQ=255)

                                                                                        aTTTTCTCG‑TTTC‑CC‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtgag           <  2:678622/100‑1 (MQ=255)

                                                                                        aTTTTCTCG‑TTTC‑CC‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtgag           >  1:1749302/1‑100 (MQ=255)

                                                                                        aTTTTCTCG‑TTT‑‑CC‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTgagtgagt          >  1:2994686/1‑100 (MQ=255)

                                                                                         ttttCTCG‑TTTC‑CC‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTgagtgagt          >  1:2144133/1‑100 (MQ=255)

                                                                                              tcG‑TTTC‑CC‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACg     <  1:1108068/100‑1 (MQ=255)

                                                                                              tcG‑TTTC‑CC‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACg     >  2:762268/1‑100 (MQ=255)

                                                                                               cG‑TTTC‑CC‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGc    <  1:4144/100‑1 (MQ=255)

                                                                                                g‑TTT‑‑CC‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGCAGTGGCGCAAGCCACCGATGCCCCTGCGTGAGTGTACGCAg  >  2:2652841/1‑100 (MQ=255)

CGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCG‑TTT‑‑CC‑CCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAG  >  NC_000913/3379160‑3379355

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: `‑`