BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_pae_pgi	2,173,363	Δ2 bp	intergenic (‑1/+1)	gatC ← / ← gatC	pseudogene, galactitol‑specific enzyme IIC component of PTS;transport; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system galactitol‑specific enzyme IIC/pseudogene, galactitol‑specific enzyme IIC component of PTS;transport; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system galactitol‑specific enzyme IIC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_pae_pgi	2,173,361	0	C	.	99.0%	496.1 / ‑3.0	103	pseudogene (2/442 nt)	gatC	pseudogene, galactitol‑specific enzyme IIC component of PTS;transport; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system galactitol‑specific enzyme IIC
Reads supporting (aligned to +/- strand): ref base C (0/0); major base . (47/55); minor base T (0/1); total (47/56)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*	NC_000913_3_pae_pgi	2,173,362	0	C	.	99.0%	497.5 / ‑2.9	103	pseudogene (1/442 nt)	gatC	pseudogene, galactitol‑specific enzyme IIC component of PTS;transport; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system galactitol‑specific enzyme IIC
Reads supporting (aligned to +/- strand): ref base C (0/1); new base . (47/55); total (47/56)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 109 total aligned reads displayed.

GAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAG  >  NC_000913_3_pae_pgi/2173215‑2173504

||

gAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccgc                                                                                                                                              >  2:644100/1‑148 (MQ=255)

   cAGATTTCCACGATGCACGGCGACCGCCATCGCTACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccgccgt                                                                                                                                           <  4:605424/150‑6 (MQ=255)

     gATTTCCACGCTGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccgccgtat                                                                                                                                         >  1:152952/1‑143 (MQ=255)

      aTTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATg                                                                                                                                      >  2:544267/1‑150 (MQ=255)

         atcatgATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTcc                                                                                                                                   <  4:748215/145‑1 (MQ=255)

          ccACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTccc                                                                                                                                  <  2:359573/150‑1 (MQ=255)

           cACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCca                                                                                                                                 >  2:1254841/1‑150 (MQ=255)

           cACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCca                                                                                                                                 >  4:69497/1‑150 (MQ=255)

            aCGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGACACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCcag                                                                                                                                >  1:426918/1‑150 (MQ=255)

                tGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCagcagc                                                                                                                            <  3:154571/150‑1 (MQ=255)

                tGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCagcagc                                                                                                                            <  2:806320/150‑1 (MQ=255)

                   aCGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAAc                                                                                                                         <  2:481906/150‑1 (MQ=255)

                   aCGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAAc                                                                                                                         >  3:413306/1‑150 (MQ=255)

                   aCGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAAc                                                                                                                         >  1:481906/1‑150 (MQ=255)

                   aCGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAACAGCAAc                                                                                                                         >  1:1274308/1‑150 (MQ=255)

                     ggCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGc                                                                                                                       <  4:19744/150‑1 (MQ=255)

                       cGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCg                                                                                                                     >  4:371183/1‑150 (MQ=255)

                        gACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCgg                                                                                                                    <  3:48938/150‑1 (MQ=255)

                         aCCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCCAACGGCCGCACCTGACTCCCCGCCAAACAAAC‑GCTAATACAATCAGTGAGTCGGATAAAAACCCGGCTTGCCCAGACCATCTC‑CTG‑ATGTCCCAGCAGCAACGCCGGa                                                                                                                   <  2:1080095/150‑1 (MQ=255)

                         aCCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGa                                                                                                                   >  2:1197143/1‑150 (MQ=255)

                           cGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATc                                                                                                                 >  3:274868/1‑150 (MQ=255)

                           cGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATc                                                                                                                 <  2:779748/150‑1 (MQ=255)

                              cATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAg                                                                                                              >  1:549236/1‑150 (MQ=255)

                              cATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAg                                                                                                              >  2:676529/1‑150 (MQ=255)

                                   ccACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCaa                                                                                                         >  4:680370/1‑150 (MQ=255)

                                    cACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAAt                                                                                                        <  4:575488/150‑1 (MQ=255)

                                        aagaagCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAgg                                                                                                    <  4:737309/150‑1 (MQ=255)

                                          gaagCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGaa                                                                                                  <  4:554011/150‑1 (MQ=255)

                                           aagCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAAc                                                                                                 >  2:11587/1‑150 (MQ=255)

                                           aagCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAAc                                                                                                 >  4:587926/1‑150 (MQ=255)

                                             gCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTc                                                                                               <  2:1088925/150‑1 (MQ=255)

                                             gCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTc                                                                                               >  3:769825/1‑150 (MQ=255)

                                                     tggCAAGGTCGCCAAACGGCAGCAGCTGATTACCCGGGACACAAACAGCAATTAAAAT‑GGAGAGTGGGATAAAAATCGGGCTTGCCGATAGCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGTACTCCTGAccgc                                                                                       >  4:220237/1‑150 (MQ=255)

                                                     tggCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccgc                                                                                       <  3:207113/150‑1 (MQ=255)

                                                       gCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCGGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAcagccg                                                                                     >  2:602460/1‑150 (MQ=255)

                                                       gCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccgccg                                                                                     >  2:26193/1‑150 (MQ=255)

                                                       gCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccgccg                                                                                     >  2:472692/1‑150 (MQ=255)

                                                         aaGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGaa                                                                                   >  1:156936/1‑150 (MQ=255)

                                                           gATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACt                                                                                 >  2:1122648/1‑150 (MQ=255)

                                                             tCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGATCTCCTGACCGCCGAACTTc                                                                               <  2:1274308/150‑1 (MQ=255)

                                                             tCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTc                                                                               <  3:69497/150‑1 (MQ=255)

                                                               gCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGc                                                                             <  3:146296/150‑1 (MQ=255)

                                                               gCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGc                                                                             >  1:720753/1‑150 (MQ=255)

                                                               gCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGc                                                                             <  4:106654/150‑1 (MQ=255)

                                                               cccAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGc                                                                             <  4:106827/149‑1 (MQ=255)

                                                                ccAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGcc                                                                            <  2:508382/150‑1 (MQ=255)

                                                                 cAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCATATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCt                                                                           <  2:66897/150‑1 (MQ=255)

                                                                   aaCGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGt                                                                         <  1:26193/150‑1 (MQ=255)

                                                                      ggCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTaaa                                                                      <  1:134278/150‑1 (MQ=255)

                                                                       gcagcaCCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCTCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAAc                                                                     <  2:816812/150‑1 (MQ=255)

                                                                         agcaCCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGa                                                                   <  1:1197143/150‑1 (MQ=255)

                                                                         agcaCCTGATTACCCGACACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGa                                                                   <  2:426918/150‑1 (MQ=255)

                                                                                gATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAg                                                            >  3:712686/1‑150 (MQ=255)

                                                                                 aTTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGc                                                           >  2:387655/1‑150 (MQ=255)

                                                                                      ccGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCt                                                      >  4:695360/1‑150 (MQ=255)

                                                                                      ccGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCt                                                      >  4:740413/1‑150 (MQ=255)

                                                                                      ccGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCt                                                      >  3:729670/1‑150 (MQ=255)

                                                                                       cGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccgc                                                                                                                                              <  3:799588/63‑3 (MQ=255)

                                                                                       cGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccgc                                                                                                                                              >  4:799588/1‑61 (MQ=255)

                                                                                        ggCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTa                                                    >  3:721417/1‑150 (MQ=255)

                                                                                        ggCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTa                                                    >  1:407058/1‑150 (MQ=255)

                                                                                        ggCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTa                                                    >  1:1146042/1‑150 (MQ=255)

                                                                                         gCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCGAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAg                                                   >  1:782069/1‑150 (MQ=255)

                                                                                         gCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAg                                                   <  4:769825/150‑1 (MQ=255)

                                                                                           acacaAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCg                                                 >  4:725920/1‑150 (MQ=255)

                                                                                            cacaAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGa                                                >  1:547777/1‑150 (MQ=255)

                                                                                            cacaAACAGCAATTAAAAT‑AGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGa                                                <  2:1215739/150‑1 (MQ=255)

                                                                                               aaaCAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATgg                                             >  2:372796/1‑150 (MQ=255)

                                                                                                aaCAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccgccg                                                                                     >  3:502292/1‑109 (MQ=255)

                                                                                                aaCAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccgccg                                                                                     <  4:502292/109‑1 (MQ=255)

                                                                                                 aCAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCAGGATCAAGGCCAATGAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGc                                           <  1:1254841/150‑1 (MQ=255)

                                                                                                 aCAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGAGACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGc                                           <  1:11587/150‑1 (MQ=255)

                                                                                                      aaTTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGGCAAGCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTaa                                      <  3:371183/150‑1 (MQ=255)

                                                                                                       agtAAAAT‑GGTTGGTGGGGTAAGAATCAGGCTTGCAGATACCA‑‑CCGGCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTaaa                                     <  2:236705/148‑1 (MQ=255)

                                                                                                       aTTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTaaa                                     >  2:995494/1‑150 (MQ=255)

                                                                                                       aTTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTaaa                                     <  1:563274/150‑1 (MQ=255)

                                                                                                           aaaT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccat                                 >  4:430522/1‑150 (MQ=255)

                                                                                                            aaT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccatc                                <  1:372796/150‑1 (MQ=255)

                                                                                                                 gTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatg                            >  3:147995/1‑150 (MQ=255)

                                                                                                                   gAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatgat                          <  2:407058/150‑1 (MQ=255)

                                                                                                                    aGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatgatg                         <  4:546245/150‑1 (MQ=255)

                                                                                                                    aGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatgatg                         <  1:588947/150‑1 (MQ=255)

                                                                                                                    aGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatgatg                         >  2:780420/1‑150 (MQ=255)

                                                                                                                      tGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGAc                                                                  >  3:705087/1‑107 (MQ=255)

                                                                                                                      tGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGAc                                                                  <  4:705087/107‑1 (MQ=255)

                                                                                                                       gggATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGt                      <  2:782069/150‑1 (MQ=255)

                                                                                                                        ggATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGtt                     <  1:1129309/150‑1 (MQ=255)

                                                                                                                         gATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGttt                    <  2:696487/150‑1 (MQ=255)

                                                                                                                          aTAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCAAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATgggacag                   >  4:674688/1‑146 (MQ=255)

                                                                                                                           tAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTaa                  >  3:164183/1‑150 (MQ=255)

                                                                                                                            aaaaaTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACt                                                                 <  2:254867/102‑1 (MQ=255)

                                                                                                                            aaaaaTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACt                                                                 >  1:254867/1‑102 (MQ=255)

                                                                                                                            aaaaaTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAAt                 <  2:547777/150‑1 (MQ=255)

                                                                                                                              aaaTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCa               <  4:147995/150‑1 (MQ=255)

                                                                                                                               aaTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCAc              <  1:1097606/150‑1 (MQ=255)

                                                                                                                                 tCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCAccc            >  2:351912/1‑150 (MQ=255)

                                                                                                                                   aGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGCTGAACGTGAAACCATCCATGATGGGTTTAATCACCCGt          <  3:587926/150‑1 (MQ=255)

                                                                                                                                     gCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCCATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTgg        <  2:210591/150‑1 (MQ=255)

                                                                                                                                           cGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTag  <  3:680370/150‑1 (MQ=255)

                                                                                                                                             cgacacccccccCGTATGTCCACTCAGCGACGCCGCATCAAGGCCAATCAGGATCTACTGTCCGCCGAACGTCGCCTGTAAACGACTACGAGCCTGCTGAGCGATGGGCGTTAAACCCTCCATGATGGTTTTAATCACCCGTGGCATTag  <  4:538076/144‑1 (MQ=255)

||

GAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAAT‑GGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAG  >  NC_000913_3_pae_pgi/2173215‑2173504

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