BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	4,572,227	G→T	17.9%	I213I (ATC→ATA)	hsdM ←	DNA methyltransferase M

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	4,572,227	0	G	T	17.9%	672.0 / 89.8	220	I213I (ATC→ATA)	hsdM	DNA methyltransferase M
Reads supporting (aligned to +/- strand): ref base G (78/102); new base T (19/21); total (97/123)
Fisher's exact test for biased strand distribution p-value = 7.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 222 total aligned reads displayed.

ACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGT  >  CP009273/4572078‑4572374

 cAATATGCGCCTTCGGCAGGTTTTCACCGTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGa                                                                                                                                                    >  2:1868840/1‑150 (MQ=255)

  aaTATGCGCCTTCGGCAGGTTTTCACCGTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGAt                                                                                                                                                   >  1:529998/1‑150 (MQ=255)

    tatGCGCCTTCGGCAGGTTTTCACCGTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATaa                                                                                                                                                 <  1:2578659/150‑1 (MQ=255)

          ccTTCGGCAGGTTTTCACCGTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAAcgcgc                                                                                                                                           >  1:1120249/1‑150 (MQ=255)

          ccTTCGGCAGGTTTTCACCGTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAAcgcgc                                                                                                                                           >  1:4079431/1‑150 (MQ=255)

          ccTTCGGCAGGTTTTCACCGTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAAcgcgc                                                                                                                                           <  1:2392200/150‑1 (MQ=255)

          ccTTCGGCAGGTTTTCACCGTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAAcgcgc                                                                                                                                           <  1:2392065/150‑1 (MQ=255)

            ttCGGCAGGTTTTCACCGTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCgg                                                                                                                                         <  2:684180/150‑1 (MQ=255)

                 cAGGTTTTCACCGTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGAt                                                                                                                                    >  1:4442418/1‑150 (MQ=255)

                   ggTTTTCACCGTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCt                                                                                                                                  <  2:833875/150‑1 (MQ=255)

                     ttttCACCGTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTgg                                                                                                                                <  2:1817381/150‑1 (MQ=255)

                      tttCACCGTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCACTATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGa                                                                                                                               <  2:238435/150‑1 (MQ=255)

                        tCACCGTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGaaa                                                                                                                             <  1:2011193/150‑1 (MQ=255)

                        tCACCGTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGaaa                                                                                                                             >  1:1510820/1‑150 (MQ=255)

                            cGTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCt                                                                                                                         >  1:3269711/1‑150 (MQ=255)

                             gTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTg                                                                                                                        <  1:1199103/150‑1 (MQ=255)

                             gTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTg                                                                                                                        <  1:3362915/150‑1 (MQ=255)

                             gTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTg                                                                                                                        <  1:3513250/150‑1 (MQ=255)

                             gTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTg                                                                                                                        >  2:3984597/1‑150 (MQ=255)

                                 cTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCgtg                                                                                                                    >  1:2216767/1‑150 (MQ=255)

                                 cTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCgtg                                                                                                                    <  1:922010/150‑1 (MQ=255)

                                      cAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCTATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGcc                                                                                                               >  2:1234213/1‑150 (MQ=255)

                                       agagTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCa                                                                                                              >  1:730471/1‑150 (MQ=255)

                                         agTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCTATAAACGCGCGGTGGATCTGGAAATCCTGCGTGGCGCCATc                                                                                                            >  1:1453564/1‑150 (MQ=255)

                                         agTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATc                                                                                                            >  1:3930296/1‑150 (MQ=255)

                                         agTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATc                                                                                                            <  2:1078211/150‑1 (MQ=255)

                                          gtgtTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCa                                                                                                           <  2:335131/150‑1 (MQ=255)

                                          gtgtTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCa                                                                                                           >  1:1275122/1‑150 (MQ=255)

                                          gtgtTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCa                                                                                                           <  2:4610197/150‑1 (MQ=255)

                                            gtTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAg                                                                                                         <  1:424451/150‑1 (MQ=255)

                                            gtTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAg                                                                                                         <  1:4501621/150‑1 (MQ=255)

                                               gCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGgtc                                                                                                      <  2:1603399/150‑1 (MQ=255)

                                                   aGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTcc                                                                                                  <  2:3898836/150‑1 (MQ=255)

                                                   aGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTcc                                                                                                  <  2:240730/150‑1 (MQ=255)

                                                    gACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCa                                                                                                 <  1:3682008/150‑1 (MQ=255)

                                                    gACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCa                                                                                                 <  2:4443665/150‑1 (MQ=255)

                                                    gACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCATCCa                                                                                                 <  2:11508/150‑1 (MQ=255)

                                                    gACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCATCCa                                                                                                 <  2:2216482/150‑1 (MQ=255)

                                                     aCGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAg                                                                                                <  2:554472/150‑1 (MQ=255)

                                                     aCGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAg                                                                                                >  2:3428026/1‑150 (MQ=255)

                                                      cGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGa                                                                                               >  1:1439733/1‑150 (MQ=255)

                                                           tGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCTATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCAtt                                                                                          >  1:490019/1‑150 (MQ=255)

                                                                cgccgTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCTATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGttt                                                                                     >  2:38235/1‑150 (MQ=255)

                                                                cgccgTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCTATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGttt                                                                                     >  2:386725/1‑150 (MQ=255)

                                                                cgccgTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCTATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGttt                                                                                     >  2:3546617/1‑150 (MQ=255)

                                                                     tGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGAc                                                                                <  1:2466543/150‑1 (MQ=255)

                                                                      ggTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACt                                                                               <  1:2876221/150‑1 (MQ=255)

                                                                      ggTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACt                                                                               >  1:1665431/1‑150 (MQ=255)

                                                                      ggTCGAGGTTGCCGTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATACACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACt                                                                               <  1:4188886/150‑1 (MQ=255)

                                                                           aGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCTATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACa                                                                          <  1:1163250/150‑1 (MQ=255)

                                                                            ggTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACAt                                                                         <  2:4555326/150‑1 (MQ=255)

                                                                              ttGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAg                                                                       <  2:4113446/150‑1 (MQ=255)

                                                                               tGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGc                                                                      <  1:2269054/150‑1 (MQ=255)

                                                                                gCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCg                                                                     >  2:2905193/1‑150 (MQ=255)

                                                                                 ccTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCgg                                                                    >  1:309997/1‑150 (MQ=255)

                                                                                   ttCAATATCGTGCATCATGCATTTCATCAGTGCCATTCGTCGGGTGCCGGGCACCAGTTCGAGGCCTATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTc                                                                  <  2:1650346/150‑1 (MQ=255)

                                                                                   ttCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCTATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTc                                                                  <  2:3159805/150‑1 (MQ=255)

                                                                                   ttCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCTATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTc                                                                  <  2:842043/150‑1 (MQ=255)

                                                                                   ttCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTc                                                                  >  1:3116508/1‑150 (MQ=255)

                                                                                   gtCAATAGCGTGCCGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAACTCCTGCGGGTCGCCATCCAGGTCGTCCAGATCAGGGGGTTGCGACGTAACATAGCGGTc                                                                  >  1:3397699/2‑150 (MQ=255)

                                                                                    tCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCg                                                                 >  2:2840718/1‑150 (MQ=255)

                                                                                      aaTATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCTATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGc                                                               >  1:2763895/1‑150 (MQ=255)

                                                                                        tatCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCTATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCtt                                                             <  2:122316/150‑1 (MQ=255)

                                                                                        tatCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCtt                                                             <  2:2029118/150‑1 (MQ=255)

                                                                                          tagTGCAGCTGGCAGTTCATCAGTGCCAGACGTCGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGGCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGGTtca                                                           <  1:765986/148‑1 (MQ=255)

                                                                                          tCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTtca                                                           <  1:2323828/150‑1 (MQ=255)

                                                                                            gTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCTATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTtcaat                                                         >  2:1621562/1‑150 (MQ=255)

                                                                                            gTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTtcaat                                                         <  2:4349966/150‑1 (MQ=255)

                                                                                                    ggcagTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTTCAATCAAAAAGc                                                 <  1:1302435/150‑1 (MQ=255)

                                                                                                    ggcagTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTTCAATCAAAAAGc                                                 >  1:688264/1‑150 (MQ=255)

                                                                                                     gcagTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTTCAATCAAAAAGcc                                                >  2:1937336/1‑150 (MQ=255)

                                                                                                          tcatcaGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCccgccg                                           >  2:3901788/1‑150 (MQ=255)

                                                                                                              caGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTAcc                                       >  2:2783757/1‑150 (MQ=255)

                                                                                                               aGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCt                                      <  2:1183551/150‑1 (MQ=255)

                                                                                                               aGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCt                                      >  2:2536929/1‑150 (MQ=255)

                                                                                                               aGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCt                                      <  2:247412/150‑1 (MQ=255)

                                                                                                                  gCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCTATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTgcc                                   >  2:4290831/1‑150 (MQ=255)

                                                                                                                  gCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTgcc                                   <  2:3624811/150‑1 (MQ=255)

                                                                                                                  gCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTgcc                                   >  2:3054342/1‑150 (MQ=255)

                                                                                                                   ccAGACGACGGGTGCCGGGCACCAGTTCGAGGCCTATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTgccg                                  >  2:1131138/1‑150 (MQ=255)

                                                                                                                   ccAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTgccg                                  >  2:4606428/1‑150 (MQ=255)

                                                                                                                    cAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTgccgc                                 >  2:2043836/1‑150 (MQ=255)

                                                                                                                     aGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTgccgcc                                >  1:3562434/1‑150 (MQ=255)

                                                                                                                       acgacgGGTGCCGGGCACCAGTTCGAGGCCTATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCgg                              <  1:10616/150‑1 (MQ=255)

                                                                                                                       acgacgGGTGCCGGGCACCAGTTCGAGGCCTATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCgg                              <  2:1372856/150‑1 (MQ=255)

                                                                                                                       acgacgGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCgg                              >  2:2019983/1‑150 (MQ=255)

                                                                                                                          acgGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTc                           <  2:2877787/150‑1 (MQ=255)

                                                                                                                           cgGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTcc                          <  1:842980/150‑1 (MQ=255)

                                                                                                                           cgGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTcc                          <  1:1860777/150‑1 (MQ=255)

                                                                                                                                     gCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACtt                >  1:2093445/1‑150 (MQ=255)

                                                                                                                                     gCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACtt                <  2:548737/150‑1 (MQ=255)

                                                                                                                                     gCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACtt                <  2:563818/150‑1 (MQ=255)

                                                                                                                                     gCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACtt                >  1:1530576/1‑150 (MQ=255)

                                                                                                                                     gCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACtt                <  2:716006/150‑1 (MQ=255)

                                                                                                                                        ccAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCAc             >  2:2756525/1‑150 (MQ=255)

                                                                                                                                        ccAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCAc             <  1:1687449/150‑1 (MQ=255)

                                                                                                                                          aGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCAcgc           >  2:2458604/1‑150 (MQ=255)

                                                                                                                                               gAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCAGCCGGGTCCTGCACCACTTCACGCGGCTg      <  1:2593835/150‑1 (MQ=255)

                                                                                                                                                aGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGc     <  2:4428774/150‑1 (MQ=255)

                                                                                                                                                   ccGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGt  >  2:1413822/1‑150 (MQ=255)

ACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCTACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAACGCGCGGTGGATCTGGAAATCCTGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTGCGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGT  >  CP009273/4572078‑4572374

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