BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,554,088	+C	coding (2/624 nt)	yhgKJ ←	RNA phosphage cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,554,088	1	.	C	99.6%	828.2 / ‑5.1	255	coding (2/624 nt)	yhgKJ	RNA phosphage cyclase
Reads supporting (aligned to +/- strand): ref base . (0/1); new base C (150/104); total (150/105)
Fisher's exact test for biased strand distribution p-value = 4.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 274 total aligned reads displayed.

GGCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCG‑AGGCGGCGCCGACGGGTTATCGGTCGCCG‑CGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTT  >  NC_000913/3553997‑3554174

 gCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCgca                                                                                        <  2:4732313/93‑3 (MQ=255)

 gCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCgca                                                                                        >  1:7526065/1‑91 (MQ=255)

 gCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCgca                                                                                        <  2:5348677/93‑3 (MQ=255)

 gCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCgca                                                                                        <  2:1235870/93‑3 (MQ=255)

  cAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCgcag                                                                                       >  2:5526090/1‑90 (MQ=255)

  cAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCgcag                                                                                       >  2:1504537/1‑90 (MQ=255)

   aGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCgcagg                                                                                      <  2:2144392/93‑5 (MQ=255)

   aGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCgcagg                                                                                      <  1:1800601/93‑5 (MQ=255)

   aGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCCGCGGATAAAATCCgcagg                                                                                      >  1:1093146/1‑89 (MQ=255)

     ggTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAggcg                                                                                    <  1:5484334/93‑1 (MQ=255)

        aaaaaCCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGgcggcg                                                                                 >  2:2829453/1‑93 (MQ=255)

        aaaaaCCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGgcggcg                                                                                 <  1:1395359/93‑1 (MQ=255)

        aaaaaCCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATACGCAGgcggcg                                                                                 >  1:2704112/1‑93 (MQ=255)

         aaaaCCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGgcgc                                                                                <  2:1040601/93‑1 (MQ=255)

         aaaaCCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGgcgc                                                                                <  1:1344017/93‑1 (MQ=255)

           aaCCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCg                                                                              >  2:2988875/1‑93 (MQ=255)

            aCCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGa                                                                             <  2:2418023/93‑1 (MQ=255)

             ccGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGCTAAAATCCGCAGGCGGCGCCGAc                                                                            >  1:4702129/1‑93 (MQ=255)

             ccGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGAc                                                                            <  1:3856739/93‑1 (MQ=255)

              cGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACg                                                                           <  2:4557803/93‑1 (MQ=255)

              cGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACg                                                                           >  1:3700843/1‑93 (MQ=255)

               gTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACgg                                                                          <  2:5987507/93‑1 (MQ=255)

               gTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACgg                                                                          <  1:1518190/93‑1 (MQ=255)

                tGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACggg                                                                         >  2:1756687/1‑93 (MQ=255)

                 ggCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGt                                                                        <  2:2549957/93‑1 (MQ=255)

                  gCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGtt                                                                       <  2:794428/93‑1 (MQ=255)

                  gCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGtt                                                                       >  1:6203835/1‑93 (MQ=255)

                     tAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATc                                                                    >  1:4220513/1‑93 (MQ=255)

                     tAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATc                                                                    <  1:7001309/93‑1 (MQ=255)

                     tAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATc                                                                    <  2:7438011/93‑1 (MQ=255)

                     tAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATc                                                                    >  2:6358360/1‑93 (MQ=255)

                     tAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATc                                                                    >  2:4367233/1‑93 (MQ=255)

                      aaCAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCg                                                                   >  2:6305821/1‑93 (MQ=255)

                       aCAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCgg                                                                  <  1:5287541/93‑1 (MQ=255)

                       aCAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCgg                                                                  <  2:7357069/93‑1 (MQ=255)

                       aCAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCgg                                                                  >  2:2599820/1‑93 (MQ=255)

                        cAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGt                                                                 <  2:7041406/93‑1 (MQ=255)

                         aGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGtg                                                                >  1:6115654/1‑92 (MQ=255)

                             tGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGt‑gccg‑c                                                          <  2:6624290/93‑1 (MQ=255)

                             tGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGt‑gccg‑c                                                          <  2:1297113/93‑1 (MQ=255)

                             tGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGt‑gccg‑c                                                          >  2:955033/1‑93 (MQ=255)

                              ggTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGt‑gccg‑cc                                                         >  1:3912919/1‑92 (MQ=255)

                              ggTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGt‑gccg‑cc                                                         >  2:4858598/1‑92 (MQ=255)

                                 ttGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGt‑gccg‑ccgct                                                      >  1:4259140/1‑89 (MQ=255)

                                   gCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCa                                                     <  2:5835890/93‑1 (MQ=255)

                                    cTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCAc                                                    <  1:5623135/93‑1 (MQ=255)

                                    cTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCAc                                                    >  2:6091267/1‑93 (MQ=255)

                                    cTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCAc                                                    >  2:3344014/1‑93 (MQ=255)

                                    cTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCAc                                                    >  1:3432080/1‑93 (MQ=255)

                                     tgatgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACt                                                   >  1:1031745/1‑93 (MQ=255)

                                       atgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCCCCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTc                                                 <  1:478115/93‑1 (MQ=255)

                                       atgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTc                                                 <  1:6303648/93‑1 (MQ=255)

                                       atgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTc                                                 >  1:7202493/1‑93 (MQ=255)

                                       atgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGGTATCGGT‑GCCGCCGCTCACTTc                                                 >  1:447061/1‑93 (MQ=255)

                                        tgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCa                                                >  2:4683171/1‑93 (MQ=255)

                                        tgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCa                                                >  1:571397/1‑93 (MQ=255)

                                         gaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCaa                                               <  1:7401157/93‑1 (MQ=255)

                                         gaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCaa                                               >  2:2002110/1‑93 (MQ=255)

                                         gaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCCGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCaa                                               >  1:6522059/1‑93 (MQ=255)

                                          aaTTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAac                                              >  1:6796810/1‑93 (MQ=255)

                                          aaTTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAac                                              >  1:6177475/1‑93 (MQ=255)

                                           aTTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAaca                                             >  1:5109581/1‑93 (MQ=255)

                                            ttCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAacac                                            >  1:2926158/1‑93 (MQ=255)

                                            ttCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAacac                                            <  1:2513756/93‑1 (MQ=255)

                                            ttCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAacac                                            >  1:1624520/1‑93 (MQ=255)

                                               cTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACAcgcg                                         >  1:7306228/1‑93 (MQ=255)

                                                tATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGa                                        >  2:4296525/1‑93 (MQ=255)

                                                 aTTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGCCGGGTTATCGG‑GGCCGCCGCTCACTTCAACACGCGaa                                       >  1:1231960/1‑93 (MQ=255)

                                                 aTTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGaa                                       <  1:4814846/93‑1 (MQ=255)

                                                  ttttCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAg                                      >  1:4582926/1‑93 (MQ=255)

                                                  ttttCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAg                                      >  2:3591591/1‑93 (MQ=255)

                                                   tttCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAgg                                     <  1:856181/93‑1 (MQ=255)

                                                        ccAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCAt                                >  2:227856/1‑93 (MQ=255)

                                                          agcagcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCg                              <  2:5746508/93‑1 (MQ=255)

                                                           gcagcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCgg                             <  1:2058644/93‑1 (MQ=255)

                                                           gcagcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCgg                             >  2:6197328/1‑93 (MQ=255)

                                                           gcagcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCgg                             >  2:1174527/1‑93 (MQ=255)

                                                           gcagcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCgg                             >  1:7145384/1‑93 (MQ=255)

                                                              gcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCaa                          >  1:6851230/1‑93 (MQ=255)

                                                              gcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCaa                          >  1:4848221/1‑93 (MQ=255)

                                                              gcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCaa                          >  2:5225650/1‑93 (MQ=255)

                                                              gcGGCTCCAGCACCCCGCGGATAAAATCCGCAGGCGGAGCCGCCGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCaa                          >  1:359583/1‑93 (MQ=255)

                                                                ggCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAAc                        >  2:3147080/1‑93 (MQ=255)

                                                                ggCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAAc                        <  2:687609/93‑1 (MQ=255)

                                                                  cTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACca                      <  2:345756/93‑1 (MQ=255)

                                                                    ccAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACcaca                    <  2:6839901/93‑1 (MQ=255)

                                                                      aGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAgc                  <  1:3209224/93‑1 (MQ=255)

                                                                      aGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAgc                  <  1:6266200/93‑1 (MQ=255)

                                                                      aGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAgc                  <  1:3132455/93‑1 (MQ=255)

                                                                        cACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAgcgc                >  1:1890896/1‑93 (MQ=255)

                                                                              gcggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgcag          >  2:2700572/1‑93 (MQ=255)

                                                                              gcggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgcag          <  1:2208982/93‑1 (MQ=255)

                                                                                ggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgcagca        >  2:1596749/1‑93 (MQ=255)

                                                                                 gATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCAc       >  2:6690696/1‑93 (MQ=255)

                                                                                 gATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCAc       >  2:3976884/1‑93 (MQ=255)

                                                                                 gATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCAc       >  2:7294432/1‑93 (MQ=255)

                                                                                 gATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCAc       >  2:5264939/1‑93 (MQ=255)

                                                                                  aTAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCAcc      >  1:5916176/1‑93 (MQ=255)

                                                                                  aTAAAATCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCAcc      <  1:6370927/93‑1 (MQ=255)

                                                                                      aaTCCGCAGGCGGCGCCGACGGGTTATCGGT‑GCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGttt  <  2:2845061/93‑1 (MQ=255)

GGCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCG‑AGGCGGCGCCGACGGGTTATCGGTCGCCG‑CGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTT  >  NC_000913/3553997‑3554174

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