BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,173,363	Δ2 bp	pseudogene (915‑916/1358 nt)	gatC ←	galactitol‑specific PTS enzyme IIC component

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,173,361	0	C	.	97.9%	588.2 / 5.1	144	pseudogene (918/1358 nt)	gatC	galactitol‑specific PTS enzyme IIC component
Reads supporting (aligned to +/- strand): ref base C (0/3); new base . (69/72); total (69/75)
Fisher's exact test for biased strand distribution p-value = 2.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*	NC_000913	2,173,362	0	C	.	98.6%	584.8 / 2.8	144	pseudogene (917/1358 nt)	gatC	galactitol‑specific PTS enzyme IIC component
Reads supporting (aligned to +/- strand): ref base C (0/2); new base . (69/72); total (69/75)
Fisher's exact test for biased strand distribution p-value = 4.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 170 total aligned reads displayed.

TTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATACCACC‑CC‑G‑CCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTG  >  NC_000913/2173222‑2173497

||

tttCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATaccac                                                                                                                                             <  1:1091788/140‑1 (MQ=255)

tttCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATaccac                                                                                                                                             <  2:2354740/140‑1 (MQ=255)

tttCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATaccac                                                                                                                                             <  1:151434/140‑1 (MQ=255)

tttCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATaccac                                                                                                                                             <  2:3291623/140‑1 (MQ=255)

tttCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATaccac                                                                                                                                             <  1:2551762/140‑1 (MQ=255)

 ttCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATaccacc                                                                                                                                            <  1:3059887/140‑1 (MQ=255)

  tCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATaccaccg                                                                                                                                           <  1:2908925/140‑2 (MQ=255)

  tCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATaccaccg                                                                                                                                           <  2:579220/140‑2 (MQ=255)

     aCGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGTATAAAAATCAGGC‑TTGCCGATaccaccgccg                                                                                                                                        <  1:609630/140‑5 (MQ=255)

        aTGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATaccaccgccgtat                                                                                                                                     <  1:2030885/140‑8 (MQ=255)

          gCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑‑CC‑G‑CCGTATGt                                                                                                                              <  2:1772602/140‑1 (MQ=255)

             cGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑‑CC‑G‑CCGTATGTccc                                                                                                                           >  1:2107702/1‑140 (MQ=255)

                  cccGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑‑CC‑G‑CCGTATGTCCcagcag                                                                                                                      <  2:3085703/139‑1 (MQ=255)

                  aCCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑‑CC‑G‑CCGTATGTCCcagcag                                                                                                                      <  1:3222822/140‑1 (MQ=255)

                   ccGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATTCCA‑‑‑CC‑G‑CCGTATGTCCCagcagc                                                                                                                     >  1:497877/1‑140 (MQ=255)

                   ccGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑‑CC‑G‑CCGTATGTCCCagcagc                                                                                                                     >  2:2137149/1‑140 (MQ=255)

                   ccGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑‑CC‑G‑CCGTATGTCCCagcagc                                                                                                                     >  1:991061/1‑140 (MQ=255)

                     gCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑‑CC‑G‑CCGTATGTCCCAGCAGCaa                                                                                                                   >  2:63541/1‑140 (MQ=255)

                       cATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑‑CC‑G‑CCGTATGTCCCAGCAGCAACg                                                                                                                 >  2:1900279/1‑140 (MQ=255)

                        aTCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑‑CC‑G‑CCGTATGTCCCAGCAGCAACGc                                                                                                                <  1:274650/140‑1 (MQ=255)

                         tCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑‑CC‑G‑CCGTATGTCCCAGCAGCAACGcc                                                                                                               <  1:3480864/140‑1 (MQ=255)

                         tCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑‑CC‑G‑CCGTATGTCCCAGCAACAACGcc                                                                                                               <  2:111784/140‑1 (MQ=255)

                                  agaagCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑‑CC‑G‑CCGTATGTCCCAGCAGCAACGCCGGATCAAgg                                                                                                      <  1:280131/140‑1 (MQ=255)

                                  agaagCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAACCAGGC‑TTGCCGATACCA‑‑‑CC‑G‑CCGTATGTCCCAGCAGCAACGCCGGATCAAgg                                                                                                      <  2:1331363/140‑1 (MQ=255)

                                   gaagCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑‑CC‑G‑CCGTATGTCCCAGCAGCAACGCCGGATCAAGGc                                                                                                     <  1:914910/140‑1 (MQ=255)

                                      gCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑‑CC‑G‑CCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCaa                                                                                                  <  1:2206653/140‑1 (MQ=255)

                                         gATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑‑CC‑G‑CCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCa                                                                                               >  2:1786784/1‑140 (MQ=255)

                                          aTGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑‑CC‑G‑CCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAg                                                                                              >  2:2510987/1‑140 (MQ=255)

                                            ggtggCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑‑CC‑G‑CCGTATGTCCCAGCAGCAACGCCGGATCAAg                                                                                                       <  1:32981/129‑1 (MQ=255)

                                             gtggCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATACCACCGCCTGTCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATc                                                                                                <  2:1893666/140‑1 (MQ=255)

                                              tggCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑‑CC‑G‑CCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAAc                                                                                          >  1:3095301/1‑140 (MQ=255)

                                              tggCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGAGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑‑CC‑G‑CCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAAc                                                                                          >  1:3227195/1‑140 (MQ=255)

                                                    gATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACACCAATTAAAATGGTGAGTGGAATAAAAATCAGGC‑TTGCCGATACCA‑‑‑CC‑G‑CCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGa                                                                                    <  2:1083474/140‑1 (MQ=255)

                                                       cGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑‑CC‑G‑CCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccg                                                                                 >  1:561682/1‑140 (MQ=255)

                                                        gCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑‑CC‑G‑CCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccgc                                                                                >  1:2169998/1‑140 (MQ=255)

                                                        gCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑‑CC‑G‑CCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccgc                                                                                >  1:1737114/1‑140 (MQ=255)

                                                        gCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑‑CC‑G‑CCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccgc                                                                                >  2:2743509/1‑140 (MQ=255)

                                                        gCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑‑CC‑G‑CCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccgc                                                                                <  2:22074/140‑1 (MQ=255)

                                                        gCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑‑CC‑G‑CCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccgc                                                                                >  1:2352927/1‑140 (MQ=255)

                                                               ggCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑‑CC‑G‑CCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACtt                                                                         <  2:2169998/140‑1 (MQ=255)

                                                               ggCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑‑CC‑G‑CCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACtt                                                                         <  2:2352927/140‑1 (MQ=255)

                                                                gcagccCCTGATTCCCCGGCACAAAAACAGCAATAAAAATGGTGAGTGGGAAAAAAATCAGGC‑TTGCCGATACCA‑‑‑CC‑G‑CCGTATGTACCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTc                                                                        <  2:833062/140‑1 (MQ=255)

                                                                 cagcaCCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑‑CC‑G‑CCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACt                                                                          >  2:3490330/1‑137 (MQ=255)

                                                                 cagcaCCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑‑CC‑G‑CCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCg                                                                       >  1:2253282/1‑140 (MQ=255)

                                                                  agcaCCTGATTCCCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑‑CC‑G‑CCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGc                                                                      <  2:2595269/140‑1 (MQ=255)

                                                                  agcaCCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑‑CC‑G‑CCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGc                                                                      >  2:3257621/1‑140 (MQ=255)

                                                                  agcaCCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑‑CC‑G‑CCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGc                                                                      >  2:1220918/1‑140 (MQ=255)

                                                                          aTTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑‑CC‑G‑CCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAAc                                                              <  2:41032/140‑1 (MQ=255)

                                                                               ccGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑‑CC‑G‑CCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTa                                                         >  2:2830989/1‑140 (MQ=255)

                                                                                cGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑‑CC‑G‑CCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTAc                                                        >  2:201327/1‑140 (MQ=255)

                                                                                 ggCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑‑CC‑G‑CCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACg                                                       >  1:725829/1‑140 (MQ=255)

                                                                                 ggCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑‑CC‑G‑CCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGCAACACCTGACCGCCGAACTTCGCCTGTAAACGACTACg                                                       >  1:193961/1‑140 (MQ=255)

                                                                                   cacacaAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑‑CC‑G‑CCGTATGT‑CCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGc                                                                      >  2:886593/1‑122 (MQ=255)

                                                                                   cacacaAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑‑CC‑G‑CCGTATGT‑CCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGc                                                                      <  1:886593/122‑1 (MQ=255)

                                                                                    acacaAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑‑CC‑G‑CCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGc                                                    >  2:1207075/1‑140 (MQ=255)

                                                                                    acacaAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑‑CC‑G‑CCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGc                                                    <  1:1786784/140‑1 (MQ=255)

                                                                                    acacaAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑‑CC‑G‑CCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGc                                                    >  2:1503714/1‑140 (MQ=255)

                                                                                    acacaAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑‑CC‑G‑CCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGc                                                    <  1:1423631/140‑1 (MQ=255)

                                                                                    acacaAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑‑CC‑G‑CCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGc                                                    <  1:2595354/140‑1 (MQ=255)

                                                                                    acacaAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑‑CC‑G‑CCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGc                                                    >  1:3486314/1‑140 (MQ=255)

                                                                                    acacaAACAGCAATTAAAAGGGTGAGTGGAATAAAAATCGGGC‑TTGCCGATCCCA‑‑‑CC‑G‑CCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGc                                                    <  2:2058445/140‑1 (MQ=255)

                                                                                      acaAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑‑CC‑G‑CCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCt                                                  <  1:2206104/140‑1 (MQ=255)

                                                                                          aCAGAAATAAAAAGGGTGAGTGGAATAAAAATCAGGC‑TTGCCGATCCCA‑‑‑CC‑G‑CCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCtt                                              <  2:1687821/140‑1 (MQ=255)

                                                                                                aTTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑‑CC‑G‑CCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGc                                                    >  2:384051/1‑128 (MQ=255)

                                                                                                aTTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑‑CC‑G‑CCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGc                                                    <  1:384051/128‑1 (MQ=255)

                                                                                                aTTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑‑CC‑G‑CCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGAt                                        >  1:1910898/1‑140 (MQ=255)

                                                                                                aTTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑‑CC‑G‑CCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGAt                                        <  2:2364100/140‑1 (MQ=255)

                                                                                                aTTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑‑CC‑G‑CCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGAt                                        <  2:2107702/140‑1 (MQ=255)

                                                                                                aTTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑‑CC‑G‑CCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGAt                                        >  1:2498881/1‑140 (MQ=255)

                                                                                                     aaTGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑‑CC‑G‑CCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCg                                   <  1:2120020/140‑1 (MQ=255)

                                                                                                       gggTGAGGGGAAAAAAAATCAGGC‑TTGCCGATAC‑‑CC‑CC‑G‑CCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGtt                                 <  2:561682/139‑1 (MQ=255)

                                                                                                        ggTGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑‑CC‑G‑CCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTa                                >  2:331785/1‑140 (MQ=255)

                                                                                                        ggTGAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑‑CC‑G‑CCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTa                                <  1:3046580/140‑1 (MQ=255)

                                                                                                          ggAGTGGGATAAAAACCAGGC‑TTGCCGAT‑CC‑CC‑CC‑G‑CCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTaaa                              <  2:2498881/139‑1 (MQ=255)

                                                                                                          ggAGTGGAAAAAAAACCAGGC‑TTGCCGAT‑CC‑CC‑CC‑G‑CCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTaaa                              <  2:1910898/139‑1 (MQ=255)

                                                                                                           gAGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑‑CC‑G‑CCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAc                             >  1:1461479/1‑140 (MQ=255)

                                                                                                            aGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑‑CC‑G‑CCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAcc                            >  2:2303122/1‑140 (MQ=255)

                                                                                                            aGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑‑CC‑G‑CCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAcc                            >  1:1816065/1‑140 (MQ=255)

                                                                                                            aGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑‑CC‑G‑CCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAcc                            >  1:3510838/1‑140 (MQ=255)

                                                                                                            aGTGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑‑CC‑G‑CCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAcc                            >  1:567323/1‑140 (MQ=255)

                                                                                                            aGTGGAATAAAAATCAGGC‑TTGCCGATACCA‑‑‑CC‑G‑CCGTATGTCCCAGCAGCAACACCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAcc                            >  2:3019554/1‑140 (MQ=255)

                                                                                                             gggggAAAAAAAT‑AGGCTTTGCCGAT‑CC‑CC‑CC‑G‑CCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAcca                           <  2:98261/138‑1 (MQ=255)

                                                                                                              tGGGATAAAAATCAGGC‑TTGCCGATACCA‑‑‑CC‑G‑CCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccat                          >  1:741345/1‑140 (MQ=255)

                                                                                                               ggaaTAAAAATCAGGC‑TTCCGGAT‑CC‑CC‑CC‑G‑CCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccatc                         <  2:567323/137‑1 (MQ=255)

                                                                                                                   tAAAAATCAGGC‑TTGCCGATACCA‑‑‑CC‑G‑CCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatg                     <  1:1394175/140‑1 (MQ=255)

                                                                                                                   tAAAAATCAGGC‑TTGCCGATACCA‑‑‑CC‑G‑CCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatg                     >  2:1394175/1‑140 (MQ=255)

                                                                                                                   tAAAAATCAGGC‑TTGCCGATACCA‑‑‑CC‑G‑CCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatg                     <  1:2510987/140‑1 (MQ=255)

                                                                                                                   aaaaaaTCAGGC‑TTGCCGATCCCA‑‑‑CC‑G‑CCGTATGT‑CCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatga                    <  1:613252/139‑1 (MQ=255)

                                                                                                                    aaaaaTCAGGC‑TTGCCGATACCA‑‑‑CC‑G‑CCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatga                    <  1:473972/140‑1 (MQ=255)

                                                                                                                    aaaaaTCAGGC‑TTGCCGATAC‑‑CT‑CC‑G‑CCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatga                    <  2:1414123/140‑1 (MQ=255)

                                                                                                                     aaaaTCAGGC‑TTGCCGATACCA‑‑‑CC‑G‑CCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatgat                   <  2:741345/140‑1 (MQ=255)

                                                                                                                       aaccAGGT‑TTGCGGATCCCA‑‑‑CC‑G‑CCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATgg                 <  2:2717988/137‑1 (MQ=255)

                                                                                                                        aTCAGGC‑TTGCCGATACCA‑‑‑CC‑G‑CCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATggg                >  2:3004962/1‑140 (MQ=255)

                                                                                                                        aTCAGGC‑TTGCCGATACCA‑‑‑CC‑G‑CCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATggg                >  1:262329/1‑140 (MQ=255)

                                                                                                                           aGGC‑TTGCCGATACCA‑‑‑CC‑G‑CCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGttt             >  2:212752/1‑140 (MQ=255)

                                                                                                                            ggC‑TTGCCGATACCA‑‑‑CC‑G‑CCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTa            >  1:2577909/1‑140 (MQ=255)

                                                                                                                                 tGCCGATACCA‑‑‑CC‑G‑CCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatg                     <  1:1896171/127‑1 (MQ=255)

                                                                                                                                  gCCGATACCA‑‑‑CC‑G‑CCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAAt          <  1:265341/137‑1 (MQ=255)

                                                                                                                                  gCCGATACCA‑‑‑CC‑G‑CCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAAt          >  2:265341/1‑137 (MQ=255)

                                                                                                                                   ccGATACCA‑‑‑CC‑G‑CCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCAcc      >  1:1569517/1‑140 (MQ=255)

||

TTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGC‑TTGCCGATACCACC‑CC‑G‑CCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTG  >  NC_000913/2173222‑2173497

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 15 ≤` `ATCG/ATCG` `< 27 ≤` `ATCG/ATCG` `< 33 ≤` `ATCG/ATCG` `< 40 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`