BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	4,572,800	A→T	45.2%	V22V (GTT→GTA)	hsdM ←	DNA methyltransferase M

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	4,572,800	0	A	T	45.2%	58.5 / inf	215	V22V (GTT→GTA)	hsdM	DNA methyltransferase M
Reads supporting (aligned to +/- strand): ref base A (56/62); new base T (51/46); total (107/108)
Fisher's exact test for biased strand distribution p-value = 4.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.26e-01

Only 100 of 217 total aligned reads displayed.

GTAGAACTGCAACTGCTCCTGGCCGATGCGGGATTTCAGGTCATCCCAGCGGTAACCTTCCGGCAGGTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCAT  >  CP009273/4572653‑4572949

gTAGAACTGCAACTGCTCCTGGCCGATGCGGGATTTCAGGTCATCCCAGCGGTAACCTTCCGGCAGGTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGATAc                                                                                                                                                     <  1:3380546/150‑1 (MQ=255)

gTAGAACTGCAACTGCTCCTGGCCGATGCGGGATTTCAGGTCATCCCAGCGGTAACCTTCCGGCAGGTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGAAAc                                                                                                                                                     >  2:4093397/1‑150 (MQ=255)

    aactgcaactgcTCCTGGCCGATGCGGGATTTCAGGTCATCCCAGCGGTAACCTTCCGGCAGGTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGAAAcgccg                                                                                                                                                 >  2:1661454/1‑150 (MQ=255)

      ctgcaactgcTCCTGGCCGATGCGGGATTTCAGGTCATCCCAGCGGTAACCTTCCGGCAGGTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGAAACgccgcc                                                                                                                                               <  1:773302/150‑1 (MQ=255)

       tgcaactgcTCCTGGCCGATGCGGGATTTCAGGTCATCCCAGCGGTAACCTTCCGGCAGGTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGATACGCCGCCa                                                                                                                                              >  1:1737856/1‑150 (MQ=255)

       tgcaactgcTCCTGGCCGATGCGGGATTTCAGGTCATCCCAGCGGTAACCTTCCGGCAGGTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGAAACGCCGCCa                                                                                                                                              >  2:882522/1‑150 (MQ=255)

         caactgcTCCTGGCCGATGCGGGATTTCAGGTCATCCCAGCGGTAACCTTCCGGCAGGTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGAAACGCCGCCATc                                                                                                                                            <  2:1425455/150‑1 (MQ=255)

          aactgcTCCTGGCCGATGCGGGATTTCAGGTCATCCCAGCGGTAACCTTCCGGCAGGTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGAAACGCCGCCATcg                                                                                                                                           >  2:3144976/1‑150 (MQ=255)

              gcTCCTGGCCGATGCGGGATTTCAGGTCATCCCAGCGGTAACCTTCCGGCAGGTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGATACGCCGCCATCGCGCa                                                                                                                                       <  2:711767/150‑1 (MQ=255)

              gcTCCTGGCCGATGCGGGATTTCAGGTCATCCCAGCGGTAACCTTCCGGCAGGTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGAAACGCCGCCATCGCGCa                                                                                                                                       >  2:4281928/1‑150 (MQ=255)

                tCCTGGCCGATGCGGGATTTCAGGTCATCCCAGCGGTAACCTTCCGGCAGGTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAgg                                                                                                                                     >  1:3329498/1‑150 (MQ=255)

                tCCTGGCCGATGCGGGATTTCAGGTCATCCCAGCGGTAACCTTCCGGCAGGTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAgg                                                                                                                                     >  1:3012506/1‑150 (MQ=255)

                 ccTGGCCGATGCGGGATTTCAGGTCATCCCAGCGGTAACCTTCCGGCAGGTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGATACGCCGCCATCGCGCAGGt                                                                                                                                    <  1:318129/150‑1 (MQ=255)

                 ccTGGCCGATGCGGGATTTCAGGTCATCCCAGCGGTAACCTTCCGGCAGGTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGt                                                                                                                                    >  2:1292688/1‑150 (MQ=255)

                   tGGCCGATGCGGGATTTCAGGTCATCCCAGCGGTAACCTTCCGGCAGGTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTg                                                                                                                                  <  1:4100636/150‑1 (MQ=255)

                     gCCGATGCGGGATTTCAGGTCATCCCAGCGGTAACCTTCCGGCAGGTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTc                                                                                                                                <  1:2487855/150‑1 (MQ=255)

                        gATGCGGGATTTCAGGTCATCCCAGCGGTAACCTTCCGGCAGGTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGca                                                                                                                             >  1:1353121/1‑150 (MQ=255)

                                    cAGGTCATCCCAGCGGTAACCTTCCGGCAGGTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGATACGCCGCCATCGCGCAGGTTGTCGcacagcttccacag                                                                                                                 >  2:1142114/1‑150 (MQ=255)

                                      ggTCATCCCAGCGGTAACCTTCCGGCAGGTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGATACGCCGCCATCGCGCAGGTTGTCGcacagcttccacagct                                                                                                               >  1:626122/1‑150 (MQ=255)

                                       gTCATCCCAGCGGTAACCTTCCGGCAGGTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGATACGCCGCCATCGCGCAGGTTGTCGcacagcttccacagctt                                                                                                              <  1:3756547/150‑1 (MQ=255)

                                        tCATCCCAGCGGTAACCTTCCGGCAGGTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGcacagcttccacagcttc                                                                                                             <  1:364898/150‑1 (MQ=255)

                                         cATCCCAGCGGTAACCTTCCGGCAGGTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGATACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTcg                                                                                                            >  1:3854280/1‑150 (MQ=255)

                                             ccAGCGGTAACCTTCCGTCAGGTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGATACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGAc                                                                                                        <  1:3954084/150‑1 (MQ=255)

                                                  ggTAACCTTCCGGCAGGTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGATACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGAt                                                                                                   >  2:3422169/1‑150 (MQ=255)

                                                   gTAACCTTCCGGCAGGTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATc                                                                                                  >  2:767031/1‑150 (MQ=255)

                                                    tAACCTTCCGGCAGGTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGATACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCg                                                                                                 >  2:908529/1‑150 (MQ=255)

                                                       ccTTCCGGCAGGTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGATACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTa                                                                                              >  2:1670230/1‑150 (MQ=255)

                                                          tCCGGCAGGTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTg                                                                                           >  2:4245067/1‑150 (MQ=255)

                                                           ccGGCAGGTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGATACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGt                                                                                          >  2:3942985/1‑150 (MQ=255)

                                                             ggCAGGTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGATACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTc                                                                                        >  1:4469963/1‑150 (MQ=255)

                                                               cAGGTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGATACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCAt                                                                                      >  2:4400984/1‑150 (MQ=255)

                                                               cAGGTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCAt                                                                                      <  2:4074631/150‑1 (MQ=255)

                                                                agttATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACTGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGATACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCAtt                                                                                     <  2:1814589/147‑1 (MQ=255)

                                                                aGGTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCAtt                                                                                     <  2:2373049/150‑1 (MQ=255)

                                                                 ggTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATtg                                                                                    <  2:2211175/150‑1 (MQ=255)

                                                                     ttCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGATACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAg                                                                                >  2:4628287/1‑150 (MQ=255)

                                                                      tCCGCTTCCTGACCGGTCTCTTTACTCATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGTATACGCCGCCATCGCTCAGTTTGTCGCACAGCTTCAACAGCTTCGCGACCATATCGTTATTGTTCATTGTGAGt                                                                               <  1:1868253/150‑1 (MQ=255)

                                                                           ttCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGt                                                                          >  2:2582441/1‑150 (MQ=255)

                                                                               tGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAAt                                                                      >  1:3031784/1‑150 (MQ=255)

                                                                               tGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAAt                                                                      >  1:2892148/1‑150 (MQ=255)

                                                                                 aCCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTa                                                                    <  2:2372103/150‑1 (MQ=255)

                                                                                   cGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGATACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAag                                                                  >  1:2298646/1‑150 (MQ=255)

                                                                                   cGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAag                                                                  <  2:1039744/150‑1 (MQ=255)

                                                                                   cGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAag                                                                  >  2:3881721/1‑150 (MQ=255)

                                                                                      tctcTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGATACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAagcag                                                               <  1:4354116/150‑1 (MQ=255)

                                                                                         ctTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCgg                                                            >  2:3837214/1‑150 (MQ=255)

                                                                                            tACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGccc                                                         >  1:3589312/1‑150 (MQ=255)

                                                                                             acacATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGATACGCCGCCATAGCGCAGGTTGTCGCACAGCTTCCTCAGATTCGAGACCAGATCGTTATTGTTCATTGTGTGTTCCGTAAATTAAGCAGCGGCCCa                                                        <  2:2375907/150‑1 (MQ=255)

                                                                                             acacATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCa                                                        <  2:2896929/150‑1 (MQ=255)

                                                                                              cacaTTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGATACACCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCaa                                                       <  1:2962002/150‑1 (MQ=255)

                                                                                              cacaTTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCaa                                                       >  2:2179428/1‑150 (MQ=255)

                                                                                                caTTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGATACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAAt                                                     >  2:633138/1‑150 (MQ=255)

                                                                                                 aTTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGATACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAAtt                                                    <  1:3407458/150‑1 (MQ=255)

                                                                                                     tCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGATACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATc                                                <  1:4429743/150‑1 (MQ=255)

                                                                                                     tCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGATACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATc                                                >  2:3671633/1‑150 (MQ=255)

                                                                                                      cAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCg                                               <  2:438838/150‑1 (MQ=255)

                                                                                                        aaaaCAGCAGCGAGGCGAGTTCATTGACGTGGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAg                                             >  1:2379875/1‑150 (MQ=255)

                                                                                                          aaCAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGATACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGcc                                           >  2:2792752/1‑150 (MQ=255)

                                                                                                             agcagcGAGGCGAGTTCATTGACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCa                                        >  1:11115/1‑150 (MQ=255)

                                                                                                              gcagcGAGGCGAGTTCATTGACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCag                                       <  1:2170227/150‑1 (MQ=255)

                                                                                                                agcGAGGCGAGTTCATTGACGTAGTTTTGATAGGATACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCagaa                                     >  2:3505578/1‑150 (MQ=255)

                                                                                                                agcGAGGCGAGTTCATTGACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCagaa                                     >  2:2243223/1‑150 (MQ=255)

                                                                                                                  cgaggcgagTTCATTGACGTAGTTTTGATAGGATACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCagaaga                                   >  1:1081106/1‑150 (MQ=255)

                                                                                                                      gcgagTTCATTGACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAgaagaagaa                               <  1:1095482/150‑1 (MQ=255)

                                                                                                                       cgagTTCATTGACGTAGTTTTGATAGGATACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGaaa                              <  1:2818850/150‑1 (MQ=255)

                                                                                                                        gagTTCATTGACGTAGTTTTGATAGGATACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAAt                             >  2:2367370/1‑150 (MQ=255)

                                                                                                                        gagTTCATTGACGTAGTTTTGATAGGATACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAAt                             <  1:2123927/150‑1 (MQ=255)

                                                                                                                        gagTTCATTGACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAAt                             >  1:2340487/1‑150 (MQ=255)

                                                                                                                         agTTCATTGACGTAGTTTTGATAGGATACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAAtt                            >  2:3320906/1‑150 (MQ=255)

                                                                                                                         agTTCATTGACGTAGTTTTGATAGGATACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAAtt                            <  1:4031766/150‑1 (MQ=255)

                                                                                                                              aTTGACGTAGTTTTGATAGGATACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGa                       >  1:3532486/1‑150 (MQ=255)

                                                                                                                              aTTGACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGa                       <  2:2663538/150‑1 (MQ=255)

                                                                                                                               ttGACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAg                      <  1:2710251/150‑1 (MQ=255)

                                                                                                                                tGACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCATATCGTTATTGTTCATTGTGAGTTCCGTAAATTTAGCAGCGGCCCAATTTCATCGAGCCGCAGAAGAAGAAATTGCCGAgg                     <  2:3820364/150‑1 (MQ=255)

                                                                                                                                 gACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAggg                    >  1:1901353/1‑150 (MQ=255)

                                                                                                                                 gACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAggg                    >  1:1331716/1‑150 (MQ=255)

                                                                                                                                   cGTAGTTTTGATAGGATACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTa                  <  1:2288843/150‑1 (MQ=255)

                                                                                                                                   cGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTa                  <  2:2315418/150‑1 (MQ=255)

                                                                                                                                   cGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTa                  <  2:4395278/150‑1 (MQ=255)

                                                                                                                                    gTAGTTTTGATAGGATACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTaa                 >  1:104736/1‑150 (MQ=255)

                                                                                                                                     tAGTTTTGATAGGATACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAat                <  2:3106772/150‑1 (MQ=255)

                                                                                                                                     tAGTTTTGATAGGATACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAat                <  1:4628439/150‑1 (MQ=255)

                                                                                                                                     tAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAat                >  1:4216415/1‑150 (MQ=255)

                                                                                                                                     tAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAat                >  1:4164403/1‑150 (MQ=255)

                                                                                                                                      aGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAata               <  2:623198/150‑1 (MQ=255)

                                                                                                                                       gTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAatat              >  1:1037706/1‑150 (MQ=255)

                                                                                                                                       gTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAatat              >  1:1683008/1‑150 (MQ=255)

                                                                                                                                       gTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAatat              >  2:3756025/1‑150 (MQ=255)

                                                                                                                                        ttttGATAGGATACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAAtata             >  2:4147143/1‑150 (MQ=255)

                                                                                                                                           tGATAGGATACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATacac          >  1:3369346/1‑150 (MQ=255)

                                                                                                                                           tGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATacac          >  1:4489695/1‑150 (MQ=255)

                                                                                                                                           tGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATacac          >  1:4075744/1‑150 (MQ=255)

                                                                                                                                            gATAGGATACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATAcaca         >  1:2893989/1‑150 (MQ=255)

                                                                                                                                              tAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACaaa       >  2:2380595/1‑150 (MQ=255)

                                                                                                                                              tAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACaaa       >  1:887463/1‑150 (MQ=255)

                                                                                                                                               aGGATACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACaaaa      <  1:95410/150‑1 (MQ=255)

                                                                                                                                               aGGATACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACaaaa      <  1:2275640/150‑1 (MQ=255)

                                                                                                                                                 gAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATc    <  2:1715245/150‑1 (MQ=255)

                                                                                                                                                 gAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGGTCCGTAAATTAAGCAGCGGCCCAAATTCCTCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATc    <  2:4040511/150‑1 (MQ=255)

                                                                                                                                                   aaCGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCAt  >  2:976402/1‑150 (MQ=255)

GTAGAACTGCAACTGCTCCTGGCCGATGCGGGATTTCAGGTCATCCCAGCGGTAACCTTCCGGCAGGTATTCCGCTTCCTGACCGGTCTCTTTACACATTTTCAAAAACAGCAGCGAGGCGAGTTCATTGACGTAGTTTTGATAGGAAACGCCGCCATCGCGCAGGTTGTCGCACAGCTTCCACAGCTTCGCGACCAGATCGTTATTGTTCATTGTGAGTTCCGTAAATTAAGCAGCGGCCCAAATTCATCGAGCCGCAGAAGAAGAAATTGCCGAGGGTAATATACACAAAATCAT  >  CP009273/4572653‑4572949

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: `‑`