BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP025268	4,331,820	(G)_7→6	pseudogene (473/817 nt)	rdcB →	YjcZ‑like family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP025268	4,331,814	0	G	.	99.2%	1705.1 / 6.4	361	pseudogene (467/817 nt)	rdcB	YjcZ‑like family protein
Reads supporting (aligned to +/- strand): ref base G (3/0); new base . (145/210); total (150/211)
Fisher's exact test for biased strand distribution p-value = 6.98e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 364 total aligned reads displayed.

GTTAACGCAACGAGGCCTACTGGTGCAGGCAGGTGATCAGAGCGAGCAAATTTTTACCCGCTGGCGTGCTGGGGCGTATAACGCCTGGTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCTGGGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGACACAACATCTGGCAGAGAGTATCAATGCGGCACCGACCACGCGTCACTATTACC  >  NZ_CP025268/4331683‑4331949

gTTAACGCAACGAGGCCTACTGGTGCAGGCAGGTGATCAGAGCGAGCAAATTTTTACCCGCTGGCGTGCTGGGGCGTATAACGCCTGGTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCTggggggcat                                                                                                                                 <  1:3919295/140‑4 (MQ=255)

gTTAACGCAACGAGGCCTACTGGTGCAGGCAGGTGATCAGAGCGAGCAAATTTTTACCCGCTGGCGTGCTGGGGCGTATAACGCCTGGTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCTggggggcat                                                                                                                                 >  2:4630168/1‑137 (MQ=255)

   aaCGCAACGAGGCCTACTGGTGCAGGCAGGTGATCAGAGCGAGCAAATTTTTACCCGCTGGCGTGCTGGGGCGTATAACGCCTGGTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCT‑GGGGGGCATTTg                                                                                                                             >  1:1495433/1‑140 (MQ=255)

     cGCAACGAGGCCTACTGGTGCAGGCAGGTGATCAGAGCGAGCAAATTTTTACCCGCTGGCGTGCTGGGGCGTATAACGCCTGGTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCT‑GGGGGGCATTTGGc                                                                                                                           <  2:2128561/140‑1 (MQ=255)

     cGCAACGAGGCCTACTGGTGCAGGCAGGTGATCAGAGCGAGCAAATTTTTACCCGCTGGCGTGCTGGGGCGTATAACGCCTGGTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCT‑GGGGGGCATTTGGc                                                                                                                           <  2:3938773/140‑1 (MQ=255)

         aCGAGGCCTACTGGTGCAGGCAGGTGATCAGAGCGAGCAAATTTTTACCCGCTGGCGTGCTGGGGCGTATAACGCCTGGTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCT‑GGGGGGCATTTGGCGAtg                                                                                                                       >  2:335592/1‑140 (MQ=255)

            aGGCCTACTGGTGCAGGCAGGTGATCAGAGCTAGCAAATTTTTACCCGCTGGCGTGCTGGGGCGTATAACGCCTGGTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCT‑GGGGGGCATTTGGCGAtgcct                                                                                                                    <  1:4131340/140‑1 (MQ=255)

             ggCCTACTGGTGCAGGCAGGTGATCAGAGCGAGCAAATTTTTACCCGCTGGCGTGCTGGGGCGTATAACGCCTGGTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCT‑GGGGGGCATTTGGCGAtgcctg                                                                                                                   <  1:2013368/140‑1 (MQ=255)

                 tACTGGTGCAGGCAGGTGATCAGAGCGAGCAAATTTTTACCCGCTGGCGTGCTGGGGCGTATAACGCCTGGTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCT‑GGGGGGCATTTGGCGATGCCTGCCGt                                                                                                               <  1:4293215/140‑1 (MQ=255)

                        gcaggcagGTGATCAGAGCGAGCAAATTTTTACCCGCTGGCGTGCTGGGGCGTATAACGCCTGGTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCT‑GGGGGGCATTTGGCGATGCCTGCCGTCTGGGaa                                                                                                        <  1:468301/140‑1 (MQ=255)

                         ccggcaggTGATCAGAGCGAGCAAATTTTTACCCGCTTTCGTGCTGTGGCGTATCACGCCTGGTCGTTGCCGGGGCGCCGTTTTATCGTTCGGGAGGAGTTGCGCT‑GGGTAGCATTGTGCGATTCCTGCCGTCTGGTAAg                                                                                                       >  2:2185566/3‑140 (MQ=255)

                         caggcaggTGATCAGAGCGAGCAAATTTTTACCCGCTGGCGTGCTGGGGCGTATAACGCCTGGTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCT‑GGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAg                                                                                                       >  1:4973048/1‑140 (MQ=255)

                           ggcaggTGATCAGAGCGAGCAAATTTTTACCCGCTGGCGTGCTGGGGCGTATAACGCCTGGTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCT‑GGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGcc                                                                                                     >  1:3970533/1‑140 (MQ=255)

                           ggcaggTGATCAGAGCGAGCAAATTTTTACCCGCTGGCGCGCTGGGGCGTATAACGCCTGGTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCT‑GTGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGcc                                                                                                     >  2:3825611/1‑140 (MQ=255)

                            gcaggTGATCAGAGCGAGCAAAGTTTTACCCGCTGGCGTACTGGGGCGTATAACGCCTGGTCGTTGCCTGGGCGCTGGTTTATCGTTCTGGAGGAGTTGCGCT‑GGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGccc                                                                                                    >  1:2347946/1‑140 (MQ=255)

                                gTGATCAGAGCGAGCAAATTTTTACCCGCTGGCGTGCTGGGGCGTATAACGCCTGGTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCT‑GGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCaa                                                                                                >  1:768235/1‑140 (MQ=255)

                                 tGATCAGAGCGAGCAAATTTTTACCCGCTGGCGTGCTGGGGCGTATAACGCCTGGTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCT‑GGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAg                                                                                               >  2:4259307/1‑140 (MQ=255)

                                    tCAGAGCGAGCAAATTTTTACCCGCTGGCGTGCTGGGGCGTATAACGCCTGGTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCT‑GGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCgg                                                                                            >  1:3483673/1‑140 (MQ=255)

                                     cAGAGCGAGCAAATTTTTACCCGCTGGCGTGCTGGGGCGTATAACGCCTGGTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCT‑GGGGGGCATTTGGCGATGCCTGCCGTCTGTTAAGCCCGCAAGCGGt                                                                                           >  2:40089/1‑140 (MQ=255)

                                     cAGAGCGAGCAAATTTTTACCCGCTGGCGTGCTGGGGCGTATAACGCCTGGTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCT‑GGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGt                                                                                           <  2:3915273/140‑1 (MQ=255)

                                      agagCGAGCAAATTTTTACCCGCTGGCGTGCTGGGGCGTATAACGCCTGGTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCT‑GGGGGGCATTTGGCGATGCCTGCCGTCCGGGAAGCCCGCagccgcgg                                                                                          >  2:3694455/1‑133 (MQ=255)

                                        agtgagcAAATTTTTACCCGCTTGCGGGCTGGCAGGTATAACGCCTGGTCGTGGCCTGGGCGCTGTGTTATCGTTCTGGAGGACTGTCGCT‑GGGAGGCATTTGGCGATGCCTGCCGTCTGGGAAGCGCGCAAGCGGGGGc                                                                                        >  2:1357032/4‑140 (MQ=255)

                                         gcgagcAAATTTTTACCCGCTGGCGTGCTGGGGCGTATAACGCCTGGTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCT‑GGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCg                                                                                       <  2:2477159/140‑1 (MQ=255)

                                         gcgagcAAATTTTTACCCGCTGGCGTGCTGGGGCGTATAACGCCTGGTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCT‑GGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCg                                                                                       <  1:2297690/140‑1 (MQ=255)

                                          cgagcAAATTTTTACCCGCTGGCGTGCTGGGGCGTATAACGCCTGGTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCT‑GGGGGGCAGTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCgg                                                                                      >  2:12759/1‑139 (MQ=255)

                                             gcAAATTTTTACCCGCTGGCGTGCTGGGGCGTATAACGCCTGGTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCT‑GGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCgttgt                                                                                   >  1:982725/1‑140 (MQ=255)

                                             gcAAATTTTTACCCGCTGGCGTGCTGGGGCGTATAACGCCTGGTCGTTGCCTGGGAGCTGTTTTATCGTTCTGGAGGAGTTGCGCT‑GGGGGGCATTTGGAGATGCCTGCCGTCTAGCAAGCCCGTAAGCGGTGGCgttgt                                                                                   >  2:1737569/1‑140 (MQ=255)

                                                      tACCCGCTGGCGTGCTGGGGCGTATAACGCCTGGTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCT‑GGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTg                                                                          <  2:292537/140‑1 (MQ=255)

                                                       aCCCGCTGGCGTGCTGGGGCGTATAACGCCTGGGCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCT‑GGTGGGCATGTGGCGATGCCTGCCGTCTGGTAAGCCCGCAAGCGGTGGCGGTGTTGCTGGGTGa                                                                         >  1:3242841/1‑140 (MQ=255)

                                                             tGGCGTGCTGGGGCGTATAACGCCTGGTCTTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCT‑GGGGGGCATTTGGCGATGCCTGACGTCTGTGAAGCCCGCAAGCGGTGGCGTTGTTGATGGGTGATTTGCt                                                                   >  2:764501/1‑140 (MQ=255)

                                                             tGGCGTGCTGGGGCGTATAACGCCTGGTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCT‑GGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCt                                                                   >  1:1002220/1‑140 (MQ=255)

                                                              ggCGTGCTGGGGCGTATAACGCCTGGTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCT‑GGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTc                                                                  <  1:2518927/140‑1 (MQ=255)

                                                                   gCTGGGGCGTATAACGCCTGGTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCT‑GGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTGGTTGCTGGGTGATTTGCTCGAGaa                                                             <  1:938332/140‑1 (MQ=255)

                                                                   gCTGGGGCGTATAACGCCTGGTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCT‑GGGGGGCATTTGGCGAGGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGaa                                                             <  1:4596155/140‑1 (MQ=255)

                                                                      ggggCGTATAACGCCTGGTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCT‑GGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGc                                                          <  2:2160477/140‑1 (MQ=255)

                                                                          cgcctagcGCCCGGTCGTTGCCTGGGCGATGTTTTATCGTTCTGGAGCAGTTGCGCT‑GGGGGGAATTTGGCGATGCTCGCCGTCTAGGCAGCCCGCAAGCGGTGGCGTTGTCCCCGGGTGATTTGCTCGAGAAAGCGaca                                                      <  2:2287570/133‑1 (MQ=255)

                                                                            tataACGCCTGGTCGTTGCCTGGGCGCTGTTTTATCGTTTTGGAGGAGTTGCGCT‑GGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGAgaca                                                    >  1:291533/1‑140 (MQ=255)

                                                                            tataACGCCTGGTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCT‑GGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGAcaca                                                    >  1:5099401/1‑140 (MQ=255)

                                                                            tataACGCCTGGTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCT‑GGGGAGCATTTGGCGATGCCAGCCGTGTGGGAAGCCCGGAAGCGGGGTCGTTGTTGCTGGCTGATTAGCTTGAGAGAGCGacatc                                                    >  2:1086863/1‑138 (MQ=255)

                                                                              taACGCCTGGTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCT‑GGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGACacaac                                                  <  2:4807158/140‑1 (MQ=255)

                                                                              taACGCCTGGTCGTTGCCTGGGCGCTGTTGTATCGTTGTGAAGGAGTTGCG‑GCGGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCTGCAAGCGGGGGCAATGTAGCTGGAAGAGGTGCTCGAGAAAACGACacaac                                                  >  2:1938407/1‑140 (MQ=255)

                                                                                aCGCCTGGTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCT‑GGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGACACAACAt                                                <  2:2413697/140‑1 (MQ=255)

                                                                                aCGCCTGGTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCT‑GGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGACACAACAt                                                <  1:819229/140‑1 (MQ=255)

                                                                                 cGCCTGGTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCT‑GGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGACACAACATc                                               >  1:3337368/1‑140 (MQ=255)

                                                                                 cGCCTGGTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCT‑GGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGACACAACATc                                               >  1:3037128/1‑140 (MQ=255)

                                                                                   ccTGGTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCT‑GGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGACACAACATCTg                                             >  2:294971/1‑140 (MQ=255)

                                                                                   ccTGGTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCT‑GGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGACACAACATCTg                                             >  2:4950226/1‑140 (MQ=255)

                                                                                   ccTGGTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCT‑GGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGACACAACATCTg                                             >  1:4933179/1‑140 (MQ=255)

                                                                                   ccTGGTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGC‑GGGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCACGAGAAAGCGACACAACATCTg                                             >  2:2830604/1‑140 (MQ=255)

                                                                                      ggTCGTTGCCTGGGCGCTGTTTTCTCGTTCTGGAGGGGTTGCGCT‑GGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGAGTTGTTGCTGGGTGATTTGCGCGAAACAGCGACACAACATCTGACa                                          >  2:2497348/1‑140 (MQ=255)

                                                                                      ggTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCT‑GGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGACACAACATCTGGCa                                          >  1:635739/1‑140 (MQ=255)

                                                                                      ggTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCT‑GGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGACACAACATCTGGCa                                          <  2:2356495/140‑1 (MQ=255)

                                                                                      ggTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCT‑GGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGACACAACATCTGGCa                                          <  2:4908237/140‑1 (MQ=255)

                                                                                      ggTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCT‑GGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGACACAACATCTGGCa                                          >  1:4795323/1‑140 (MQ=255)

                                                                                      ggTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCT‑GGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGACACAACATCTGGCa                                          <  1:1685365/140‑1 (MQ=255)

                                                                                      ggTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCT‑GGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGACACAACATCTGGCa                                          >  2:562706/1‑140 (MQ=255)

                                                                                      ggTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCT‑GGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGACACAACATCTGGCa                                          <  2:1486897/140‑1 (MQ=255)

                                                                                      ggTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCT‑GGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGACACAACATCTGGCa                                          <  1:1344524/140‑1 (MQ=255)

                                                                                      ggTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCT‑GGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGACACAACATCTGGCa                                          >  2:1985505/1‑140 (MQ=255)

                                                                                      ggTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCT‑GGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCCTTGTTGCTGGGTGATTTGCTCGAGAAAGCGACACAACATCTGGCa                                          >  2:87797/1‑140 (MQ=255)

                                                                                      ggTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGC‑GGGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGACACAACATCTGGCa                                          >  2:580503/1‑140 (MQ=255)

                                                                                        tCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCT‑GGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGACACAACATCTGGCaga                                        <  2:4526565/140‑1 (MQ=255)

                                                                                        tCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCT‑GGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGACACAACATCTGGCaga                                        <  1:562735/140‑1 (MQ=255)

                                                                                        tCGTTGCCTGGGCGCTGTTGTATCGTTCTGGAGGAGTTGCGCT‑GGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGACACAACATCTGGCaga                                        <  1:2105872/140‑1 (MQ=255)

                                                                                        tCGTTGCCTGGGCGCTGTGTTATCGTTCTGGAGGAGTTGCGCT‑GGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGACACAACATCTGGCaga                                        <  1:2223274/140‑1 (MQ=255)

                                                                                         cGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCT‑GGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGACACAACATCTGGCagag                                       <  1:3680493/140‑1 (MQ=255)

                                                                                              ccTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCT‑GGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGCGAAAGCGACACAACATCTGGCAGAGAGTAt                                  <  1:3618546/140‑1 (MQ=255)

                                                                                              ccTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCT‑GGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGACACAACATCTGGCAGAGAGTAt                                  <  1:2150935/140‑1 (MQ=255)

                                                                                              ccTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCT‑GGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGACACAACATCTGGCAGAGAGTAt                                  <  2:4676550/140‑1 (MQ=255)

                                                                                              ccTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCT‑GGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGACACAACATCTGGCAGAGAGTAt                                  <  2:3464427/140‑1 (MQ=255)

                                                                                              ccTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCT‑GGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGACACAACATCTGGCAGAGAGTAt                                  <  2:3183927/140‑1 (MQ=255)

                                                                                              ccTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCT‑GGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCACAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGACACAACATCTGGCAGAGAGTAt                                  <  1:2821102/140‑1 (MQ=255)

                                                                                                  ggCGCTGTTTTATCGTTCTGGAGGAGTTGCGCT‑GGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGACACAACATCTGGCAGAGAGTATCAAt                              <  1:3810253/140‑1 (MQ=255)

                                                                                                  ggCGCTGTTTTATCGTTCTGGAGGAGTTGCGCT‑GGGGGGCATTTGGCGATGCCTGCCGGCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGACACAACATCTGGCAGAGAGTATCAAt                              >  2:3652910/1‑140 (MQ=255)

                                                                                                   gcgcTGTTTTATCGTTCTGGAGGAGTTGCGCT‑GGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGACACAACATCTGGCAGAGAGTATCAATg                             <  1:1300009/140‑1 (MQ=255)

                                                                                                   gcgcTGTTTTATCGTTCTGGAGGAGTTGCGCT‑GGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGACACAACATCTGGCAGAGAGTATCAATg                             <  1:1629195/140‑1 (MQ=255)

                                                                                                   gcgcTGTTTTAGATTTCTGGAGGAGTGGTGCT‑GGGCAGCATCAGGTGTTTCCAGCCGTCTGGGAATGCAGCAAGCAGTGGGCTTGTTGGTGGGTGATTAGCTCGAGAAAGTGACACAACTTCTGGCAGAGAGTATCAATg                             <  2:670539/140‑1 (MQ=255)

                                                                                                     gcTGTTTTATCGTTCTGGAGGAGTTGCGCT‑GGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGACACAACATCTGGCAGAGAGTATCAATGCg                           <  1:3530508/140‑1 (MQ=255)

                                                                                                      cTGTTTTATCGTTCTGGAGGAGTTGCGCT‑GGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATCTGCTCGAGAAAGCGACACAACATCTGGCAGAGAGTATCAATGCgg                          >  1:4004297/1‑140 (MQ=255)

                                                                                                        gTTTTATCGTTCTGGAGGAGTTGCGCT‑GGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGACACAACATCTGGCAGAGAGTATCAATGCGGCa                        <  1:3123899/140‑1 (MQ=255)

                                                                                                             aTCGTTCTGGAGGAGTTGCGCT‑GGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGACACAACATCTGGCAGAGAGTATCAATGCGGCACCGAc                   <  1:465769/140‑1 (MQ=255)

                                                                                                               cGTTCTGGAGGAGTTGCGCT‑GGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGACACAACATCTTGCAGAGAGTATCAATGCGGCCCCGACCa                 <  2:1130857/140‑1 (MQ=255)

                                                                                                               cGTTCTGGAGGAGTTGCGCT‑GGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGACACAACATCTGGCAGAGAGTATCAATGCGGCACCGACCa                 <  1:74380/140‑1 (MQ=255)

                                                                                                               cGTTCTGGAGGAGTTGCGCT‑GGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGACACAACATCTGGCAGAGAGTATCAATGCGGCACCGACCa                 <  1:4497362/140‑1 (MQ=255)

                                                                                                                  tCTGGAGGAGTTGCGCT‑GGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGACACAACATCTGGCAGAGAGTATCAATGCGGCACCGACCAcgc              <  1:103311/140‑1 (MQ=255)

                                                                                                                  tCTGGAGGAGTTGCGCT‑GGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGACACAACATCTGGCAGAGAGTATCAATGCGGCACCGACCAcgc              <  1:2616670/140‑1 (MQ=255)

                                                                                                                   cTGGCGGGGCTGCGCT‑GGGGGGCCTTTGGCGATGCCTGCCGTCTGGGAAGGCAGCCAGAGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAGGCGACACAAAATCTCGCAGAGCGTATCAAGGCGGAACAGAGTAcgcg             >  2:5035854/1‑140 (MQ=255)

                                                                                                                   cTGGAGGAGTTGCGCT‑GGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAGGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGACACAACATCTGGCAGAGAGTATCAATGCGGCACCGACCAcgcg             >  2:4182826/1‑140 (MQ=255)

                                                                                                                   cTGGAGGAGTTGCGCT‑GGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGACACAACATCTGGCAGAGAGTATCAATGCGGCACCGACCAcgcg             >  1:1902276/1‑140 (MQ=255)

                                                                                                                    tGGAGGAGTTGCGCT‑GGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGACACAACATCTGGCAGAGAGTATCAATGCGGCACCGACCACGCGt            <  1:3369723/140‑1 (MQ=255)

                                                                                                                    tGGAGGAGTTGCGCT‑GGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGACACAACATCTGGCAGAGAGTATCAATGCGGCACCGACCACGCGt            <  2:3046505/140‑1 (MQ=255)

                                                                                                                    tGGAGGAGTTGCGCT‑GGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGACACAACATCTGGCAGAGAGTATCAATGCGGCACCGACCACGCGt            <  1:1776157/140‑1 (MQ=255)

                                                                                                                    tGGAGGAGTTGCGCT‑GGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGACACAACATCTGGCAGAGAGTATCAATGCGGCACCGACCACGCGt            <  2:2753552/140‑1 (MQ=255)

                                                                                                                    tGGAGGAGTTGCGCT‑GGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGACACAACATCTGGCAGAGAGTATCAATGCGGCACCGACCACGCGt            <  1:524140/140‑1 (MQ=255)

                                                                                                                    tGGAGGAGTTGCGCT‑GGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGACACAACATCTGGCAGAGAGTATCAATGCGGCACCGACCACGCGt            <  2:1284875/140‑1 (MQ=255)

                                                                                                                      gaggagTTGCGCT‑GGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGACACAACATCTGGCAGAGAGTATCAATGCGGCACCGACCACGCGTCa          <  2:4760319/140‑1 (MQ=255)

                                                                                                                      gaggagTTGCGCT‑GGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGACACAACATCTGGCAGAGAGTATCAATGCGGCACCGACCACGCGTCa          >  2:619990/1‑140 (MQ=255)

                                                                                                                        ggagTTGCGCT‑GGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGACACAACATCTGGCAGAGAGTATCAATGCGGCACCGACCACGCGTCACt        >  2:1554548/1‑140 (MQ=255)

                                                                                                                           gTTGCGCT‑GGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGACACAACATCTGGCAGAGAGTATCAATGCGGCACCGACCACGCGTCACTAtt     >  2:4988202/1‑140 (MQ=255)

                                                                                                                              gcgcT‑GGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGACACAACATCTGGCAGAGAGTATCAATGCGGCACCGACCACGCGTCACTATTAcc  >  2:2685869/1‑140 (MQ=255)

GTTAACGCAACGAGGCCTACTGGTGCAGGCAGGTGATCAGAGCGAGCAAATTTTTACCCGCTGGCGTGCTGGGGCGTATAACGCCTGGTCGTTGCCTGGGCGCTGTTTTATCGTTCTGGAGGAGTTGCGCTGGGGGGGCATTTGGCGATGCCTGCCGTCTGGGAAGCCCGCAAGCGGTGGCGTTGTTGCTGGGTGATTTGCTCGAGAAAGCGACACAACATCTGGCAGAGAGTATCAATGCGGCACCGACCACGCGTCACTATTACC  >  NZ_CP025268/4331683‑4331949

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