BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NZ_CP022686	3,935,341	T→G	7.9%	V699V (GTA→GTC)	lacZ ←	beta‑galactosidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP022686	3,935,341	0	T	G	7.9%	564.6 / 4.7	253	V699V (GTA→GTC)	lacZ	beta‑galactosidase
Reads supporting (aligned to +/- strand): ref base T (68/156); new base G (26/2); total (95/158)
Fisher's exact test for biased strand distribution p-value = 1.16e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 255 total aligned reads displayed.

CCGCTGGTGGTCAGTTGAGGGATGGCGTGGGACGCGGAGGGGAGTGTCACGCTGAGGTTTTCCGCCAGACGCCATTGCTGCCAGGCGCTGATGTGTCCGGCTTCTGACCATGCGGTCGCGTTTGGTTGCACTACGCGTACCGTTAGCCAGAGTTGTCCGGCGCTCTCCGGCTGCGGCAGTTCAGGCAGTTCAATCAACTGCTTACCTTGCGGGCCAACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTGCCATCCAGTGCCACCATCCAGTGC  >  NZ_CP022686/3935204‑3935479

ccGCTGGTGGTCAGTTGAGGGATGGCGTGGGACGCGGAGGGGAGTGTCACGCTGAGGTTTTCCGCCAGACGCCATTGCTGCCAGGCGCTGATGTGTCCGGCTTCTGACCATGCGGTCGCGGTTGGTTTCACAACGCGt                                                                                                                                            >  5:466127/1‑138 (MQ=255)

  gcgggtggtCAGCTGAGGGAGGGCGTGGGACGCGGAGGGGAGTGTCACGCTGAGGTGTTCCGCCAGACGCCATTGCTGCCAGGCGCTGATGTGTCCGGCTTCTGACCATGCGGTCGCGTTTGGTTGCACTACGCGTACCg                                                                                                                                        <  3:110303/137‑1 (MQ=255)

  gCTGGTGGTCAGTTGAGGGATGGCGTGGGACGCGGAGGGGAGTGTCACGCTGAGGTTTTCCGCCAGACGCCATTGCTGCCAGGCGCTGATGTGTCCGGCTTCTGACCATGCGGTCGCGTTTGGTTGCACTACGCGTACCg                                                                                                                                        <  8:517923/140‑1 (MQ=255)

  gCTGGTGGTCAGTTGAGGGATGGCGTGGGACGCGGAGGGGAGTGTCACGCTGAGGTTTTCCGCCAGACGCCATTGCTGCCAGGCGCTGATGTGTCCGGCTTCTGACCATGCGGGCGCGTTTGGTTGCACTACGCGTACCg                                                                                                                                        <  4:559816/140‑1 (MQ=255)

  gCTGGTGGCCAGTGGAGGGAGTGCGGGGGACGCGGAGGGGAGGGTCACGCTGAGGTTTTCCGCCAGACGCCATTGCTGCCAGGCGCTGATGTGTCCGGCTTCTGACCATGCGGTCGCGTTTGGTTGCACTACGCGTACCg                                                                                                                                        <  4:1338750/140‑1 (MQ=255)

   cTGGTGGGCAGTTGAGGGATGGCGTGGGACGCGGAGGGGAGTGTCACGCTGAGGTTTTCCGCCAGACGCCATTGCTGCCAGGCGCTGATGTGTCCGGCTTCTGACCATGCGGTCGCGTTTGGTTGCACTACGCGTACCGt                                                                                                                                       <  3:583650/140‑1 (MQ=255)

       ttgtaaGTTGAGGGATGGCGTGGGACGCGGAGGGTATTGGCACTCTGAGGTTTTCAGCCAGACGCCATTGCTGCCAGGCGATGATGTGTCCGGCTTCTGACCATGCGGTCGCGTTTGGTTGCACTACGCGTACCGTTa                                                                                                                                     <  2:380700/133‑1 (MQ=255)

       gggtCAGTTGAGGGATGGCGTGGGACGCGGAGGGGAGTGTCACGCTGAGGTTTTCCGCCAGACGCCATTGCTGCCAGGCGCTGATGTGTCCGGCTTCTGACCATGCGGTCGCGTTTGGTTGCACTACGCGTACCGTTa                                                                                                                                     <  6:1199990/137‑1 (MQ=255)

             gTTGAGGGATGGCGTGGGACGCGGAGGGGAGTGTCACGCTGAGGTTTTCCGCCAGACGCCATTGCTGCCAGGCGCTGATGTGTCCGGCTTCTGACCATGCGGTCGCGTTTTGTTGCACTACGCGTACCGTTAGCCaga                                                                                                                               >  2:879697/1‑138 (MQ=255)

                 aGGGATGGCGGGGGACGCGGAGGGGAGTGTCTCGCTGAGGTTTTCCGCCAGACGCCATTGCTGCCAGGCGCTGATGTGTCCGGCTTCTGACCATGCGGTCGCGTTTGGTTGCACTACGCGTACCGTTAGCCAGAGTTg                                                                                                                           <  1:645583/138‑1 (MQ=255)

                  gggATGGCGTGGGACGCGGAGGGGAGTGTCACGCTGAGGTTTTCCGCCAGACGCCATTGCTGCCAGGCGCTGATGTGTCCGGCTTCTGACCATGCGGTCGCGTTTGGTTGCACTACGCGTACCGTTAGCCAGAGTTGt                                                                                                                          >  6:1141108/1‑138 (MQ=255)

                    gctGGCGTGGGACGCGGAGGGGAGTGTCACGCTGAGGGTTTCCGCCAGACGCCATTGCTGCCAGGCGCTGATGTGTCCGGCTTCTGACCATGCGGTCGCGTTTGGTTGCACTACGCGTACCGTTAGCCAGAGTTGTCCgg                                                                                                                      <  4:140776/138‑1 (MQ=255)

                    gATGGCGTGGGACGCGGAGGGGAGTGTCACGCTGAGGTTTTCCGCCAGACGCCATTGCTGCCAGGCGCTGATGTGTCCGGCTTCTGACCATGCGGTCGCGTTTGGTTGCACTACGCGTACCGTTAGCCAGAGTTGTcc                                                                                                                        <  5:883079/138‑1 (MQ=255)

                        gCGTGGGACGCGGAGGGGAGTGTCACGCTGAGGTTTTCCGCCAGACGCCATTGCTGCCAGGCGCTGATGTGTCCGGCTTCTGACCATGCGGTCGCGTTTGGTTGCACTACGCGTACCGTTAGCCAGAGTTGTCCGgcgcg                                                                                                                  >  4:1142171/1‑139 (MQ=255)

                        gCGTGGGACGCGGAGGGGAGTGTCACGCTGAGGTTTTCCGCCAGACGCCATTGCTGCCAGGCGCTGATGTGTCCGGCTTCTGACCATGCGGTCGCGTTTGGTTGCACGACGCGTAACGTTAGCCAGAGTTGTCAGgcgcg                                                                                                                  >  4:586179/1‑139 (MQ=255)

                        gCGTGGGACGCGGAGGGGAGTGTCACGCTGAGGTTTTCCGCCAGACGCCATTGCTGCCAGGCGCTGATGTGTCCGGCTTCTGACCATGCGGTCGCGTGTGGTTGGACCAAGCGTACCGATATCCAGAGTTGTACGgcgcg                                                                                                                  >  7:545063/1‑139 (MQ=255)

                        gCGTGGGACGCGGAGGGGAGTGTCACGCTGAGGTTTTCCGCCAGACGCCATTGCTGCCAGGCGCTGATGTGTCCGGCTTCTGACCATGCGGTCGCGGTTGGTTGGACAACGCGGACCGTTAGGCAGAGTTGTACGgcc                                                                                                                    >  2:108530/1‑137 (MQ=255)

                        gCGTGGGACGCGGAGGGGAGTGTCACGCTGAGGTTTTCCGCCAGACGCCATTGCTGCCAGGCGCTGATGTGTCCGGCTTCTGACCATGCGGTCGCGATTGGTTGCACTACGCGTACCGGTAGCCAGAGTTGTCCGgcgcg                                                                                                                  >  8:21169/1‑139 (MQ=255)

                          gTGGGACGCGGAGGGGAGTGTCACGCTGAGGTTTTCCGCCAGACGCCATTGCTGCCAGGCGCTGATGTGTCCGGCTTCTGACCATGCGGTCGCGTTTGGTTGGACTAACCGTACCGTTATCCAGATTTGTCCGgcgca                                                                                                                  >  1:700153/1‑137 (MQ=255)

                          gTGGGACGCGGAGGGGAGTGTCACGCTGAGGTTTTCCGCCAGACGCCATTGCTGCCAGGCGCTGATGTGTCCGGCTTCTGACCATGCGGTCGCGGTTGGTTTCCCTACGCGGACCCTTATGCAGAATTTTGCGgcgca                                                                                                                  >  2:343769/1‑137 (MQ=255)

                          gTGGGACGCGGAGGGGAGTGTCACGCTGAGGTTTTCCGCCAGACGCCATTGCTGCCAGGCGCTGATGTGTCCGGCTTCTGACCATGCGGTCGCGGTTGGTTTCACTAACCGGACCGATTCCCAGAATTTTGCGgagcg                                                                                                                  >  1:1085351/1‑137 (MQ=255)

                               aCGCGGAGGGGAGTGTCACGCTGAGGTTTTCCGCCAGACGCCATTGCTGCCAGGCGCTGATGTGTCCGGCTTCTGAACATGCGGTCGCGGTTTGGTGGACAACGCGTACCGTTAGCCAGCGTTGTACGGCGCTCTccc                                                                                                             >  1:1324572/1‑137 (MQ=255)

                                  cgGAGGGGAGTGTCACGCTGAGGTTTTCCGCCAGACGCCATTGCTGCCAGGCGCTGATGTGTCCGGCTTCTGCCCATGCGGTCGCGTTTGGTTGCACATCGCGGACCGTTAGCCAGAGTTGTCCGGCGGTCTCCGGCTgg                                                                                                        >  7:770218/1‑139 (MQ=255)

                                   ggAGGGGAGTGTCACGCTGAGGTTTTCCGCCAGACGCCATTGCTGCCAGGCGCTGATGTGTCCGGCTTCTGACCATGCGGTCGCGTTTGGTTGCACTACGCGTACCGTTAGCCAGAGTTGTCCGGCGCTCTCCGGCTg                                                                                                         <  1:905826/138‑1 (MQ=255)

                                   ggAGGGGAGTGTCACGCTGAGGTTTTCCGCCAGACGCCATTGCTGCCAGGCGCTGATGTGTCCGGCTTCTGACCATGCGGTCGCGTTTGGTTGCACTACGCGTACCGTTAGCCAGAGTTGTCCGGCGCTCTCCGGCTGCg                                                                                                       <  4:1110103/140‑1 (MQ=255)

                                   ggAGGGGAGTGTCACGCTGAGGTTGTCCGCCAGACGCCATTGCTGCCAGGCGCTGATGTGTCCGGCTTCTGACCATGCGGTCGCGTTTGGTTGCACTACGCGTACCGTTAGCCAGAGTTGTCCGGCGCTCTCCGGCTGCg                                                                                                       <  4:748262/140‑1 (MQ=255)

                                   ggAGGGGAGTGCCACGCAGAGGGTTTCCGCCAGACGCCATTGCTGCCAGGCGCTGATGTGTCAGGCTTCTGACCATGCGGTCGCGTTTGGTTGCACTACGCGTACCGTTAGCCAGAGTTGTCCGGATCTCTCCGGCTg                                                                                                         <  6:1006270/138‑1 (MQ=255)

                                    gAGGGGAGTGTCACGCTGAGGTTTTCCGCCAGACGCCATTGCTGCCAGGCGCTGATGTGTCCGGCTTCTGACCATGCGGTCGCGTTTGGTTGCCCTACGCGTACCGTTAGCCAGAGTTGTCCGGCGCGCTCCGGCCGc                                                                                                        >  2:75266/1‑138 (MQ=255)

                                     aGGGGAGTGTCACGCTGAGGTTTTCCGCCAGACGCCATTGCTGCCAGGCGCTGATGTGTCCGGCTTCTGACCATGCGGTCGCGTTTGTTTGCACTACGCGTACCGTTTGCCAGACTTTTCCGGCGCTCACCCgcggcg                                                                                                       >  6:1045575/1‑138 (MQ=255)

                                        ggAGTGTCACTCTGAGGTTTTCCGCCAGACGCCATAGCTGCCAGGCGCTGATGTGTCCGGCTTCTGACCATGCTGTCGCGTTTGGTTGGACTACCCCTAACCATAGCCAGAGATGTCCGGCCCGAACCCGCTGCggcacg                                                                                                  >  8:1260117/1‑138 (MQ=255)

                                        ggAGTGTCACGCTGAGGGTTTCCGCCAGACGCCATTGCTGCCAGGCGCTGATGTGTCCGGCTTCTGACCATGCGGTCGCGTTTGGTTGCACTACGCGTACCGTTAGCCAGAGTTGTCCGGCGCTCTCCGGCTGCggca                                                                                                    <  5:813228/138‑1 (MQ=255)

                                            tgtCACGCTGAGGTTTTCCGCCAGACGCCATTGCTGCCAGGCGCTGATGTGTCCGGCTTCTGCCCATGCGGTCGCGTTTGGTTGCACTACGCGTACCGTTAGCCAGAGTTGTCCGGCGCTCTCCGGCTGCggcagttcag                                                                                              <  4:126863/140‑1 (MQ=255)

                                            tgtCACGCTGAGGTTTTCCGCCAGACGCCATTGCTGCCAGGCGCTGATGTGTCCGGCTTCTGACCATGCGGTCGCGTTTGGTTGCACAACTCGTACCGTTATCCAGAGTTTGCCGGCGCTCTCCGGAGGCggcagttcag                                                                                              >  7:1193371/1‑140 (MQ=255)

                                             gtCACGCTGAGGTTTTCCGCCAGACGCCATTGCTGCCAGGCGCTGATGTGTCCGGCTTCTGACCATGCGGTCGCGTTTGGTTGCACTACGCGTACCGTTAGCCAGAGTTGTCCGGCGCTCTCCGGCTGCggcagttca                                                                                               <  6:664758/138‑1 (MQ=255)

                                             gtCACGCTGAGGTTTTCCGCCAGACGCCATTGCTGCCAGGCGCTGATGTGTCCGGCTTCTGACCATGCGGTCGCGTTTGGTTGCACTACGCGTACCGTTAGCCAGAGTTGTCCGGCGCTCTCCGGCTGCggcagttca                                                                                               <  1:434951/138‑1 (MQ=255)

                                             gtCACGCTGAGGTATTCCGCCAGACGCCATTGCTGCCAGGCGGTGATGTGTCCGGCTTCTGACCATTCGGTCGCGTTTGGGTGCACTACGCGGACAGTTAGCCAGAGTTGTCCGGCGCTCTCCGGCTGCggcagttcagg                                                                                             <  7:1063098/140‑1 (MQ=255)

                                              tCACGCTGAGGTTTTCCGCCAGACGCCATTGCTGCCAGGCGCTGATGTGTCCGGCTTCTGACCATGCGGTCGCGGTTGGTTGGACTACGCGGACCGGTTGCCAGAGTTGTCCGGGGCGAGGCGGCGGCggccgttgca                                                                                              >  1:936442/1‑135 (MQ=255)

                                               cACGCTGAGGTTTTCCGCCAGACGCCATTGCTGCCAGGCGCTGATGTGTCCGGCTTCTGACCATGCGGTCGCGTTTGGTTGCACTACGCGTACCGTTAGCCAGAGTTGTCCGGCGCTCTCCGGCTTCggcagttcaggca                                                                                           <  3:85585/140‑1 (MQ=255)

                                                 cGCTGAGGTTTTCCGCCAGACGCCATTGCTGCCAGGCGCTGATGTGTCCGGCTTCTGACCATGCGGTCGCGTTTTGTTGCACTACGCGTACCGTTAGCCAGAGTTGTACGGCGCTCGCCCGCTGCgggcgttgaggccgt                                                                                         >  3:1121942/1‑140 (MQ=255)

                                                   cTGAGGTTTTCCGCCAGACGCCATTGCTGCCAGGCGCTGATGTGTCCGGCTTCTGACCATGCGGTCGCGGTTGGTTTCACTACCCCGACCCTTATCCAGAGTTGTGCGGCGCACTCCGGCTGCgggagttaaggcagt                                                                                         >  5:1222032/1‑138 (MQ=255)

                                                       ggTTTTCCGCCAGACGCCATTGCTGCCAGGCGCTGATGTGTCCGGCTTCTGACCATGCGGTCGCGTTTGGTTGCACTACGCGTACCGTTAGCCAGAGTTGTCCGGCGCTCTCCGGCTGCGGCAGTTCAGGCAGTtcaa                                                                                     <  6:1356916/138‑1 (MQ=255)

                                                       ggTTTTCCGCCAGACGCCATTGCTGCCAGGCGCTGATGTGTCCGGCTTCTGACCATGCGGTCGCGTTTGGTTGCACTACGCGTACCGTTAGCCAGAGTTGTCCGGCGCTCTCCGGCTGCGGCAGTTCAGGCAGTtcaa                                                                                     <  1:569431/138‑1 (MQ=255)

                                                        gtgttCCGCCAGACGCCATTGCTGCCAGGCGCTGATGTGTCCGGCTTCTGACCATGCGGTCGCGTTTGGTTGCACTACGCGTACCGTTAGCCAGAGTTGTCCGGCGCTCTCCGGCTGCGGCAGTTCAGGCAGTtcaat                                                                                    <  1:1372117/135‑1 (MQ=255)

                                                        gTTTTCCGCCCGACGCCATTGCTGCCAGGCGCTGATGTGTCCGGCTTCTGACCATGCGGTCGCGGTTGGTTGGCCTACGCGGGCCGGTTGCCCGAGTTTGGCGGGGGTCTACGGCGGGggcaggtgagggcgttcaca                                                                                    >  6:632219/1‑136 (MQ=255)

                                                        gTTTTCCGCCAGACGCCATTGCTGCCAGGCGCTGATGTGTCCTGCTTCTGACCATGCGGTCGCGTTTGGTTGCACTACGCGTACCGTTAGCCAGAGTTGTCCGGCGCTCTCCGGCTGCGGCAGTTCAGGCAGTtcaat                                                                                    <  1:246369/138‑1 (MQ=255)

                                                        gTTTTCCGCCAGACGCCATTGCTGCCAGGCGCTGATGTGTCCGGCTTCTGACCATGCGGTCGCGTTTGGTTGCACTACGCGTACCGTTAGCCAGAGTTGTCCGGCGCTCTCCGGCTGCGGCAGTTCAGGCAGTtcaatca                                                                                  <  4:599971/140‑1 (MQ=255)

                                                        gTTTTCCGCCAGACGCCATTGCTGCCAGGCGCTGATGTGTCCGGCTTCTGACCATGCGGTCGCGGTTGGTTGCACTACGCGGACCGTTATCCAGAGTTTTACGGCGGTCTCCGGCATCGGGAGTTACGGGAGGtcaat                                                                                    >  6:667109/1‑138 (MQ=255)

                                                         ttttCCGCCAGACGCCATTGCTGCCAGGCGCTGATGTGTCCGGCTTCTGACCATGCGGTCGCGTTTGGCTGCACTACGCGTACCGTTAGCCAGAGTTGTCCGGCGCTCTCCGGCTGCGGCAGTTCAGGCAGTtcaatc                                                                                   <  5:384044/138‑1 (MQ=255)

                                                          tttCCGCCAGACGCCATTGCTGCCAGGCGCTGATGTGTCCGGCTTCTGACCATGCGGTCGCGTTGGGTTGCACTACGCGTACCGTGAGCCAGAGTTGTCCGGCGCTCTCCGGCTGCGGCAGTTCAGGCAGTtcaatca                                                                                  <  5:214687/138‑1 (MQ=255)

                                                                    aCGCCATTGCTGCCAGGCGCTGATGTGTCCGGCTTCTGACCATGCGGTCGCGGTTGGTTGCACTACGCGTACCGTTAGCCAGAGTTGTCCGGCGCTCTCCGGCTGCGGCAGGTAAGGCAGTTAAATCAACTGCTTAcc                                                                        >  6:827744/1‑138 (MQ=255)

                                                                     cGCCATTGCTGCCAGGCGCTGATGTGGCCGGCTTCTGACCATGCGGTCGCGTTTGGTTGCACTACGCGTACCGTTAGCCAGAGTTGTCCGGCGCTCTCCGGCTGCGGCAGTTCAGGCAGTTCAATCAACTGCTTACCTTg                                                                     <  4:855252/140‑1 (MQ=255)

                                                                       ccATTGCTGCCAGGCGCTGATGTGTCCGGCTTCTGACCATGCGGTCGCGTTTGGTTGCACTACGCGTACCGTTAGCCAGAGTTGTCCGGCGCTCTCCGGCTGCGGCAGTTCAGGCAGTTCAATCAACTGCTTACCTTGCg                                                                   <  8:172271/140‑1 (MQ=255)

                                                                       ccATTGCTGCCAGGCGCTGATGTGTCCGGCTTCTGACCATGCGGTCGCGTTTGGTTGCACTACGCGTACCGTTAGCCAGAGTTGTCCGGCGCTCTCCGGCTGCGGCAGTTCAGGCAGTTCAATCAACTGCTTACCTTGCg                                                                   <  4:724817/140‑1 (MQ=255)

                                                                           tgctgcCAGGCGCTGATGTGTCCGGCTTCTGACCATGCGGTCGCGTTTGGTTGCACTACGCGTACCGTTAGCCAGAGTTGTCCGGCGCTCTCCGGCTGCGGCAGTTCAGGCAGTTCAATCAACTGCTTACCTTGCggg                                                                 <  2:772786/138‑1 (MQ=255)

                                                                           tgctgcCAGGCGCTGATGTGTCCGGCTTCTGACCATGCGGTCGCGTTTGGTTGCACTACGCGTACCGTTAGCCAGAGTTGTCCGGCGCTCTCCGGCTGCGGCAGTTCAGGCAGTTCAATCAACTGCTTACCTTGCGGGcc                                                               <  7:813563/140‑1 (MQ=255)

                                                                           tgctgcCAGGCGCTGATGTGTCCGGCTTCTGACCATGCGGTCGCGTTTGGTTGCACTACGCGTACCGTTAGCCAGAGTTGTCCGGCGCTCTCCGGCTGCGGCAGTTCAGGCAGTTCAATCAACTGCTTACCTTGCGGGcc                                                               <  7:288352/140‑1 (MQ=255)

                                                                            gctgcCAGGCGCTGATGTGTCCGGCTTCTGACCATGCGGTCGCGTTTGGTTGCACTACGCGTACCGTTAGCCAGAGTTGTCCGGCGCTCTCCGGCTGCGGCAGTTCAGGCAGTTCAATCAACTGCTTACCTTGCGGGCCa                                                              <  8:65671/140‑1 (MQ=255)

                                                                            gctgcCAGGCGCTGATGTGTCCGGCTTCTGACCATGCGGTCGCGTTTGGTTGCACTACGCGTACCGTTAGCCAGAGTTGTCCGGCGCTCTCCGGCTGCGGCAGTTCAGGCAGTTCAATCAACTGCTTACCTTGCGGGCCa                                                              <  3:1247219/140‑1 (MQ=255)

                                                                             ctgcCAGGCGCTGATGTGTCCGGCGTCTGACCATGCGGTCGCGTTTGGTTGCACTACGCGTACCGTTAGCCAGAGTTGTCCGGCGCTCTCCGGCTGCGGCAGTTCAGGCAGTTCAATCAACTGCTTACCTTGCGGGcc                                                               <  2:41384/138‑1 (MQ=255)

                                                                             ctgaCAGGCGCTGATGTGTCCGGCTTCTGACGATGCGGTCGCGTTTGGTTGAACTACGCGTACCGTTAGACAGAGTTGTCCGGCTCGCTCCGGCTGCGGCAGTTCAGGTAGTTCACTCAACTGCTTAACTTGCGGGcc                                                               <  6:1243000/138‑1 (MQ=255)

                                                                                  aGGCGCTGATGTGTCCGGCTTCTGACCATGCGGTCGCGTTTTGTTGGCCTACGCGTACCGTTAGCCAGCGTTTTACGGCGCGGTACCGCGGCGGCAGGTAAGGCAGTGTAATGAAAAGAGTACCATGCGGGCCGAcca                                                          >  6:967346/1‑136 (MQ=255)

                                                                                  aGGCGCTGATGTGTCCGGCTTCTGACCATGCGGTCGCGTTTGGTTGCACTACGCGTACCGTTAGCCAGAGTTGTCCGGCGCTCTCCGGCTGCGGCAGGTCAGGGACTTCAATAAACTGCTTACCTTGCGGGCCcacaa                                                          >  2:202864/1‑137 (MQ=255)

                                                                                    gcgcTGATGTGTCCTGCTTCTGACCATGCGGTAGCGTTTGTTTGCACTACGCGTACCGTTAGCCAGAGTTGTCCGGCGCTCTCAGGCTGCGGCAGTTCAGGGAGTTCAATCAACTGCTTACCTTGCGGGCCAACATCCag                                                      <  7:689868/140‑1 (MQ=255)

                                                                                    gcgcTGATGTGTCCGGCTTCTGACCATGCGGTCGCGTTTGGTTTCACTAAGCGTACCGTTAGCCAGAGTTGTACGGCGCTCTCCGGCGGCGGCAGTTGCCGCCGTTGAATAAACTGCGTACCATGGGGGCCAACATcc                                                        >  6:464617/1‑138 (MQ=255)

                                                                                     cgcTGATGTGTCCGGCTTCTGACCATGCGGTCGCGGTTGGTTGCACAACGCGGACCGTTAGCCAGAGTTGTCCGGCGCTCTCCGGCTGCGGCAGTTAAGGCAGTTCAATACACTGCATACCGTGCGGGCCCCCATccc                                                       >  2:605767/1‑137 (MQ=255)

                                                                                       cTGATGTGTCCGGCTTCTGACCATGCGGTCGCGTTTGGTTGCACTACGCGTACCGTTAGCCAGAGTTGTCCGGCGCTCTCCGGCTGCGGCAGTTCAGGCAGTTCAATCAACTGCTTACCTTGCGGGCCAACATCCaga                                                     <  5:1141108/138‑1 (MQ=255)

                                                                                          aTGTGTCCGGCTTCTGACCATGCGATCGCGTTTGGTTGCACTACGCGTACCGTTAGCCAGAGTTGTCCGGCGCTCTCCGGCTGCGGCAGTTCAGGCAGTTCAATCAACTGCTTACCTTGCGGGCCAACATCCAGAGGCAc                                                >  4:1239794/1‑140 (MQ=255)

                                                                                            gtgtCCGGCTTCTGACCATGCGGTCGCGGTTGGTTGCACTACGCGTACCGTTAGCCAGAGTTGTCCGGCGCTCTCCGGCTGCGGCAGTTCAGGCAGTTCAATCAACTGCTTACCTTGCGGGCCAACATCCAGAGGCACtt                                              <  4:394876/140‑1 (MQ=255)

                                                                                               tCCGGCTTCTGACCATGCGGTCGCGTTTGGTTGCCCTACGCCGACCGTTAGCCAGAGTTTTGCGGCGCTCTCCGGCGGCGGGACTTAACGCAGTTCAATCAACAGCTTACAGTGCGGGGCCACCTACAGCGGcccgtaac                                           >  4:231627/1‑133 (MQ=255)

                                                                                               tCCGGCTTCTGACCATGCGGTCGCGTTTGGTTGCACTACGCGTACCGTTAGCCAGAGTTGTCCGGCGCTCTCCGGCTGCGGCAGTTCAGGCAGTTCAATCAACTGCTTACCTTGCGGGCCAACATACAGAGGCACTTCAc                                           >  4:1038251/1‑140 (MQ=255)

                                                                                                 cGGCTTCTGACCATGCGGTCGCGTTTGGGTGCTCTACGCGTACCGTTAGCCAGAGTTGTCCGGCGCTCTCCGGCTGCGGCAGTTCAGGCAGTTCAATCAACTGCATACCTTGCGGGCCAACATCCAGAGGCACTTCACCg                                         <  3:589181/140‑1 (MQ=255)

                                                                                                      tCTGACCATGCGGTCGCGTTTGGTTGCACAACGCGTACCGTTAGCCCGAGGTGTCCGGCGCTCTCCGGCTGCGGCAGTTCAGGCAGTTGAAAACACAGCTTACCAAGCGGGCCCACCACCCGACGCACATAACCCCtt                                      >  2:848148/1‑138 (MQ=255)

                                                                                                       cTGACCATTCGGTCGCGGTTGGTTGCACTACGCGTACCGTTAGCCAGAGTTGTCCGGCGCTCTCCGGCTGCGGCAGTTCAGGCAGTTCAATCAACTGCTTACCTTGCGGGCCAACATCCAGAGGCACTTCACCGCTTGcc                                   <  7:113511/140‑1 (MQ=255)

                                                                                                                gggTCGCGTTTGGTTGCACTACGCGTACCGTTAGCCAGAGTTGTCCGGCGCTCTCCGGCTGCGGCAGTTCAGGCAGTTCAATCAACTGCTTACCTTGCGGGCCAACATCCAGAGGCACTTCACCGCTTGCCAGCGGCt                            <  6:240894/137‑1 (MQ=255)

                                                                                                                cGGTCGCGTTTGGTTGCACTACGCGTACCGTTAGCCAGAGTTGTCCGGCGCTCTCCGGCTGCGGCAGTTCAGGCAGTTCAATCAACTGCTTACCTTGCGGGCCAACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTg                          <  7:434887/140‑1 (MQ=255)

                                                                                                                        tttGGTTGCACTACGCGTACCGTTAGCCAGAGTTGTCCGGCGCTCTCCGGCTGCGGCAGTTCAGGCAGTTCAATCAACTGCTTACCTTGCGGGCCAACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTGCCATcc                    <  5:176827/138‑1 (MQ=255)

                                                                                                                        tttGGTTGCACTACGCGTACCGTTAGCCAGAGTTGTCCGGCGCTCTCCGGCTGCGGCAGTTCAGGCAGTTCAATCAACTGCTTACCTTGCGGGCCAACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTGCCATcc                    <  1:16654/138‑1 (MQ=255)

                                                                                                                        tttGGTTGCACTACGCGTACCGTTAGCCAGAGTTGTCCGGCGCTCTCCGGCTGCGGCAGTTCAGGCAGTTCAATCAACTGCTTACCTTGCGGGCCAACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTGCCATCCAg                  <  3:1305487/140‑1 (MQ=255)

                                                                                                                        tttGGTTGCACTACGCGTACCGTTAGCCAGAGTTGTCCGGCGCTCTCCGGCTGCGGCAGTTCAGGCAGTTCAATCAACTGCTTACCTTGCGGGCCAACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTGCCATCCAg                  <  4:106254/140‑1 (MQ=255)

                                                                                                                        tttGGTTGCACTACGAGTACCGTTAGCCAGAGTTGTCCGGCGCTCTCCGGCTGCGGCAGTTCAGGCAGTTCAATCAACTGCTTACCTTGCGGGCCAACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTGCCATcc                    <  2:574446/138‑1 (MQ=255)

                                                                                                                               gCACTACGCGTACCGTTAGCCAGAGGTGTCCGTCGCTCTCCGGCTGCGGCAGGTCAGGCAGTTCAATCAACTGCTTACCTTGCGGGCCAACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTGACATCGAGTGcca             <  6:599525/138‑1 (MQ=255)

                                                                                                                                  cTACGCGTACCGTTAGGCAGAGTTGTCCGGCGCTCTCCGGCTGCGGCAGTTCAGGCAGTTCAATCAACTGCTTACCTTGCGGGCCAACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTGCCATCCAGTGccacca          <  1:1039087/138‑1 (MQ=255)

                                                                                                                                  cTACGCGTACCGTTAGCCAGAGTTGTCCGGCGCTCTCCGGCTGCGGCAGTTCAGGCAGTTCAATCAACTGCTTACCTTGCGGGCCAACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTGCCATCCAGTGccccca          >  5:632966/1‑138 (MQ=255)

                                                                                                                                  cTACGCGTACCGTTAGCCAGAGTTGTCCGGCGCTCTCCGGCTGCGGCAGTTCAGGCAGTTCAATCAACTGCTTACCTTGCGGGCCAACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTGCCATCCAGTGccacca          <  6:718238/138‑1 (MQ=255)

                                                                                                                                  cTACGCGTACCGTTAGCCAGAGTTGTCCGGCGCTCTCCGGCTGCGGCAGTTCAGGCAGTTCAATCAACTGCTTACCTTGCGGGCCAACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTGCCATCCAGTGc               >  6:210219/1‑133 (MQ=255)

                                                                                                                                  cTACGCGTACCGTTAGCCAGAGTTGTCCGGCGCTCTCCGGCTGCGGCAGTTCAGGCAGTTCAATCAACTGCTTACCTTGCGGGCCAACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTGCCATCCAGTGCCACCATc        <  8:294862/140‑1 (MQ=255)

                                                                                                                                  cTACGCGTACCGTTAGCCAGAGTTGTCCGGCGCTCTCCGGCTGCGGCAGTTCAGGCAGTTCAATCAACTGCTTACCTTGCGGGCCAACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTGCCATCCAGTGCCACCATc        <  7:1018562/140‑1 (MQ=255)

                                                                                                                                  cTACGCGTACCGTTAGCCAGAGTTGTCCGGCGCTCTCCGGCTGCGGCAGTTCAGGCAGTTCAATCAACTGCTTACCTTGCGGGCCAACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTGCCATCCAGTGCCACCATc        <  3:914264/140‑1 (MQ=255)

                                                                                                                                  cTACGCGTACCGTTAGCCAGAGTTGTCCGGCGCTCTCCGGCTGCGGCAGTTCAGGCAGTTCAATCAACTGCTTACCTTGCGGGCCAACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTGCCATCCAGTGCCACCATc        <  8:1076698/140‑1 (MQ=255)

                                                                                                                                  cTACGCGTACCGTTAGCCAGAGTTGTCCGGCGCTCTCCGGCTGCGGCAGTTCAGGCAGTTAAATCAACTGCTTACCTTGCGGGCCAACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTGCCATCCAGTGccacca          >  2:1271648/1‑138 (MQ=255)

                                                                                                                                   tACGCGTACCGGTAGCCAGAGTTGTCCGGCGCTCTCCGGCTGCGGCAGTTCAGGCAGTTCAATCAACTGCTTACCTTGCGGGCCAACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTGCCATCCAGTGCCACCATcc       <  8:832568/140‑1 (MQ=255)

                                                                                                                                    actcgTTCCGTTAGCCATATGTGTCCGGCTCGCTCCGGCTGAGGAAGTTCAGGCATTTCAATCAACTGCTTACCTTGCGGGCCAACATCCAGAGGCACTTAACCGCTTGCCAGCGGCTTGCCATCCAGTGCCACCATCCa      <  4:311428/137‑1 (MQ=255)

                                                                                                                                     cgcgTACCGTTAGCCAGAGTTGTCCGGCGCTCTCCGGCTGCGGCAGTTCAGGCAGTTCAATCACCTGCTTACCTTGCGGGCCAACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTGCCATCCAGTGCCACCATCCAg     >  8:1174488/1‑140 (MQ=255)

                                                                                                                                     cgcgTACCGTTAGCCAGAGTTGTCCGGCGCTCTCCGGCTGCGGCAGTTCAGGCAGTTCAATCAACTGCTTACCTTGCGGGCCAACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTGCCATCCAGTGCCACCATCCAg     <  3:462874/140‑1 (MQ=255)

                                                                                                                                       cgTACCGTTAGCCAGAGTTGTCCGGCGCTCTCCGGCTGCggcagttca                                                                                               <  2:1443713/48‑1 (MQ=255)

                                                                                                                                       cgTACCGTTAGCCAGAGTTGTCCGGCGCTCTCCGGCTGCggcagttca                                                                                               >  1:1443713/1‑48 (MQ=255)

                                                                                                                                       cgTACCGTTAGCCAGAGTTGTCCGGCGCTCTCCGGCTGCGGCAGTTCAGGCAGTTCAATCAACTGCTTACCTTGCGGGCCAACATCCAGCGGCACTTCACCGCTTGCCAGCGGCTTGCCATCCAGTGCCACCCTCCAg     >  5:101377/1‑138 (MQ=255)

                                                                                                                                       cgTACCGTTAGCCAGAGTTGTCCGGCGCTCTCCGGCTGCGGCAGTTCAGGCAGTTCAATCAACTGCTTACCTTGCGGGCCAACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTGCCATCCAGTGCCACCATCCAg     >  6:378940/1‑138 (MQ=255)

                                                                                                                                        gTACCGTTAGCCAGAGTTGTCCGGCGCTCTCCGGCTGCGGCAGTTCAGGCAGTTCAATCAACTGCTTACCTTGCGGGCCAACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTGCCATCCAGTGCCACCATCCAGTGc  <  7:145874/140‑1 (MQ=255)

                                                                                                                                         tACCGTTAGCCAGAGTTGTCCGGCGCTCTCCGGCTGCCGCAGTTCAGGAAGTTCAATCAACTGCTTACATTGCGCGCCAAAATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTTCCCTCCAGTGCCACCATCCAGTg   >  6:765591/1‑138 (MQ=255)

CCGCTGGTGGTCAGTTGAGGGATGGCGTGGGACGCGGAGGGGAGTGTCACGCTGAGGTTTTCCGCCAGACGCCATTGCTGCCAGGCGCTGATGTGTCCGGCTTCTGACCATGCGGTCGCGTTTGGTTGCACTACGCGTACCGTTAGCCAGAGTTGTCCGGCGCTCTCCGGCTGCGGCAGTTCAGGCAGTTCAATCAACTGCTTACCTTGCGGGCCAACATCCAGAGGCACTTCACCGCTTGCCAGCGGCTTGCCATCCAGTGCCACCATCCAGTGC  >  NZ_CP022686/3935204‑3935479

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