BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	2,136,421	A→G	19.9%	D3D (GAT→GAC)	udk ←	uridine/cytidine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	2,136,421	0	A	G	19.9%	824.1 / 204.3	352	D3D (GAT→GAC)	udk	uridine/cytidine kinase
Reads supporting (aligned to +/- strand): ref base A (160/122); new base G (44/26); total (204/148)
Fisher's exact test for biased strand distribution p-value = 4.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

Only 100 of 355 total aligned reads displayed.

TGATCTTTGTAATAGCAGTCTTCGGGAATTACGCCGATGTGTTCATCACCGACTTGCTCACGCAATTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATCCGTCTGATACCAGTCGAATAGCGTTGCCGCGCGCTCAGAGTTAATTGTTGACAAAGAATTCC  >  CP009273/2136272‑2136570

  aTCTTTGTAATAGCAGTCTTCGGGAATTACGCCGATGTGTTCATCACCGACTTGCTCACGCACTTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATc                                                                                                                                                     >  1:625033/1‑150 (MQ=255)

  aTCTTTGTAATAGCAGTCTTCGGGAATTACGCCGATGTGTTCATCACCGACTTGCTCACGCAATTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGGTc                                                                                                                                                     <  1:5391266/150‑1 (MQ=255)

   tCTTTGTAATAGCAGTCTTCGGGAATTACGCCGATGTGTTCATCACCGACTTGCTCACGCAATTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGGTCa                                                                                                                                                    <  1:1201788/150‑1 (MQ=255)

   tCTTTGTAATAGCAGTCTTCGGGAATTACGCCGATGTGTTCATCACCGACTTGCTCACGCAATTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCa                                                                                                                                                    <  1:400893/150‑1 (MQ=255)

      ttGTAATAGCAGTCTTCGGGAATTACGCCGATGTGTTCATCACCGACTTGCTCACGCAATTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTc                                                                                                                                                 <  2:3780082/150‑1 (MQ=255)

       tGTAATAGCAGTCTTCGGGAATTACGCCGATGTGTTCATCACCGACTTGCTCACGCAATTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCa                                                                                                                                                >  2:1308017/1‑150 (MQ=255)

       tGTAATAGCAGTCTTCGGGAATTACGCCGATGTGTTCATCACCGACTTGCTCACGCAATTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCa                                                                                                                                                >  1:2048375/1‑150 (MQ=255)

         tAATAGCAGTCTTCGGGAATTACGCCGATGTGTTCATCACCGACTTGCTCACGCAATTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCata                                                                                                                                              >  1:2130262/1‑150 (MQ=255)

          aaTAGCAGTCTTCGGGAATTACGCCGATGTGTTCATCACCGACTTGCTCACGCAATTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCatat                                                                                                                                             >  2:1583589/1‑150 (MQ=255)

          aaTAGCAGTCTTCGGGAATTACGCCGATGTGTTCATCACCGACTTGCTCACGCAATTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCatat                                                                                                                                             >  1:1949469/1‑150 (MQ=255)

          aaTAGCAGTCTTCGGGAATTACGCCGATGTGTTCATCACCGACTTGCTCACGCAATTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACGGATCAGTCatat                                                                                                                                             >  2:2580111/1‑150 (MQ=255)

             aGCAGTCTTCGGGAATTACGCCGATGTGTTCATCACCGACTTGCTCACGCAATTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATAtt                                                                                                                                          >  2:3572916/1‑150 (MQ=255)

              gCAGTCTTCGGGAATTACGCCGATGTGTTCATCACCGACTTGCTCACGCAATTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATAttt                                                                                                                                         <  2:432609/150‑1 (MQ=255)

               cAGTCTTCGGGAATTACGCCGATGTGTTCATCACCGACTTGCTCACGCAATTCACGATACAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGCGACTGGTCAGTCCTTTAtttg                                                                                                                                        >  1:3530047/1‑149 (MQ=255)

                    ttCGGGAATTACGCCGATGTGTTCATCACCGACTTGCTCACGCAATTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGGTCAGTCATATATTTAGCGAc                                                                                                                                   >  1:4814091/1‑150 (MQ=255)

                    ttCGGGAATTACGCCGATGTGTTCATCACCGACTTGCTCACGCAATTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGAc                                                                                                                                   <  1:4150320/150‑1 (MQ=255)

                     tCGGGAATTACGCCGATGTGTTCATCACCGACTTGCTCACGCAATTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGGTCAGTCATATATTTAGCGAcc                                                                                                                                  <  2:278174/150‑1 (MQ=255)

                       gtgAATTACGCATATGTGTTCATCACCGAGTTGCTCACTCAATTCACGATAAAGGGTACTGGAAATAAGACTCTTGACGGAAGCCGATGCTCCAGCGATACAGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTg                                                                                                                                <  1:3392067/148‑1 (MQ=255)

                       gggAATTACGCCGATGTGTTCATCACCGACTTGCTCACGCAATTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGGTCAGTCATATATTTAGCGACCTg                                                                                                                                >  1:4227749/1‑150 (MQ=255)

                        ggAATTACGCCGATGTGTTCATCACCGACTTGCTCACGCAATTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGGTCAGTCATATATTTAGCGACCTGa                                                                                                                               >  1:1684071/1‑150 (MQ=255)

                        ggAATTACGCCGATGTGTTCATCACCGACTTGCTCACGCAATTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGa                                                                                                                               <  2:1175615/150‑1 (MQ=255)

                           aTTACGCCGATGTGTTCATCACCGACTTGCTCACGCAATTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTa                                                                                                                            <  2:3670530/150‑1 (MQ=255)

                             tACGCCGATGTGTTCATCACCGACTTGCTCACGCAATTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGGTCAGTCATATATTTAGCGACCTGATTAAc                                                                                                                          <  1:34678/150‑1 (MQ=255)

                             tACGCCGATGTGTTCATCACCGACTTGCTCACGCAATTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAAc                                                                                                                          <  1:2441060/150‑1 (MQ=255)

                                 ccGATGTGTTCATCACCGACTTGCTCACGCAATTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTgg                                                                                                                      >  1:2767916/1‑150 (MQ=255)

                                 ccGATGTGTTCATCACCGACTTGCTCACGCAATTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTgg                                                                                                                      <  2:942458/150‑1 (MQ=255)

                                  cGATGTGTTCATCACCGACTTGCTCACGCAATTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGa                                                                                                                     >  1:4921309/1‑150 (MQ=255)

                                   gATGTGTTCATCACCGACTTGCTCACGCAATTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGAt                                                                                                                    >  1:2459133/1‑150 (MQ=255)

                                     tgtgTTCATCACCGACTTGCTCACGCACTTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGt                                                                                                                  >  1:813473/1‑150 (MQ=255)

                                     tgtgTTCATCACCGACTTGCTCACGCAATTCACGATAAAGGGTACTGGCAATAATACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCAATGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGt                                                                                                                  >  1:1554938/1‑150 (MQ=255)

                                     tgtgTTCATCACCGACTTGCTCACGCAATTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGt                                                                                                                  >  2:1054426/1‑150 (MQ=255)

                                       tgtTCATCACCGACTTGCTCACGCAATTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTa                                                                                                                <  1:5522344/150‑1 (MQ=255)

                                       tgtTCATCACCGACTTGCTCACGCAATTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTa                                                                                                                <  1:2713541/150‑1 (MQ=255)

                                        gtTCATCACCGACTTGCTCACGCAATTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAg                                                                                                               >  1:128746/1‑150 (MQ=255)

                                         ttCATCACCGACTTGCTCACGCAATTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGGTCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAgg                                                                                                              >  1:1003888/1‑150 (MQ=255)

                                          tcatcaCCGACTTGCTCACGCAATTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGa                                                                                                             <  1:5078336/150‑1 (MQ=255)

                                            atcaCCGACTTGCTCACGCAATTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAg                                                                                                           >  2:4345523/1‑150 (MQ=255)

                                                    cTTGCTCACGCAATTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGGTCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGcgac                                                                                                   >  2:837768/1‑150 (MQ=255)

                                                     ttGCTCACGCAATTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGGTCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGcgacg                                                                                                  <  2:5172093/150‑1 (MQ=255)

                                                     ttGCTCACGCAATTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGcgacg                                                                                                  >  1:1511252/1‑150 (MQ=255)

                                                       gCTCACGCAATTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGaa                                                                                                >  1:4428855/1‑150 (MQ=255)

                                                        cTCACGCAATTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGGTCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAg                                                                                               <  2:4252527/150‑1 (MQ=255)

                                                         tCACGCAATTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAgc                                                                                              <  1:926018/150‑1 (MQ=255)

                                                          cACGCAATTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAgcg                                                                                             >  2:2088764/1‑150 (MQ=255)

                                                              cAATTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCaa                                                                                         >  1:2556210/1‑150 (MQ=255)

                                                                aTTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGGTCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAAc                                                                                       <  2:478808/150‑1 (MQ=255)

                                                                aTTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAAc                                                                                       <  2:5615067/150‑1 (MQ=255)

                                                                 ttCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAAcg                                                                                      >  1:645996/1‑150 (MQ=255)

                                                                   cACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAAcgcg                                                                                    >  2:77249/1‑150 (MQ=255)

                                                                      gATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCa                                                                                 >  2:2489426/1‑150 (MQ=255)

                                                                      gATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCa                                                                                 <  1:2218809/150‑1 (MQ=255)

                                                                       aTAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCaa                                                                                <  2:1225099/150‑1 (MQ=255)

                                                                       aTAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCaa                                                                                <  1:1386231/150‑1 (MQ=255)

                                                                        tAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAAt                                                                               >  1:4924203/1‑150 (MQ=255)

                                                                         aaaGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAAtt                                                                              <  2:3920636/150‑1 (MQ=255)

                                                                          aaGGGTACTGGCTTTAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTTGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATta                                                                             <  1:4995039/150‑1 (MQ=255)

                                                                          aaGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATta                                                                             <  1:3027220/150‑1 (MQ=255)

                                                                          aaGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATta                                                                             >  2:3462533/1‑150 (MQ=255)

                                                                            gggTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATtata                                                                           >  1:2641916/1‑150 (MQ=255)

                                                                                aCTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGa                                                                       <  2:310407/150‑1 (MQ=255)

                                                                                     cAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCa                                                                  >  1:1977216/1‑150 (MQ=255)

                                                                                            aCTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATCCGTCt                                                           <  1:1166322/150‑1 (MQ=255)

                                                                                            aCTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATCCGTCt                                                           >  1:5443933/1‑150 (MQ=255)

                                                                                             ctctTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGGTCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATCCGTCTg                                                          >  2:170184/1‑150 (MQ=255)

                                                                                             ctctTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATCCGTCTg                                                          <  2:1398059/150‑1 (MQ=255)

                                                                                               ctTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATCCGTCTGAt                                                        >  2:5105391/1‑150 (MQ=255)

                                                                                                 agCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGACGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATCCTTCTGATAc                                                      >  2:5259020/2‑150 (MQ=255)

                                                                                                   ccGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGGTCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATCCGTCTGATACCa                                                    >  1:991296/1‑150 (MQ=255)

                                                                                                     ggAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGGTCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATCCGTCTGATACCAGt                                                  >  1:1048802/1‑150 (MQ=255)

                                                                                                        aGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATCCGTCTGATACCAGTCGa                                               >  2:1857208/1‑150 (MQ=255)

                                                                                                         gCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATCCGTCTGATACCAGTCGaa                                              <  2:2391427/150‑1 (MQ=255)

                                                                                                           cGATGCGCCAGCGATACCGATAATGACTCACTGATGAGACTGATCAGTCATATATTTAGCGACCTTATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATCCGTCTGATACCAGTCGAATa                                            <  2:4392337/150‑1 (MQ=255)

                                                                                                            gATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGGTCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATCCGTCTGATACCAGTCGAATAg                                           >  2:2862875/1‑150 (MQ=255)

                                                                                                            gATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATCCGTCTGATACCAGTCGAATAg                                           <  2:5010652/150‑1 (MQ=255)

                                                                                                               gcgcCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATCCGTCTGATACCAGTCGAATAGCGt                                        <  1:5112127/150‑1 (MQ=255)

                                                                                                                 gcCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATCCGTCTGATACCAGTCGAATAGCGTTg                                      <  1:4751673/150‑1 (MQ=255)

                                                                                                                 gcCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATCCGTCTGATACCAGTCGAATAGCGTTg                                      <  2:3700320/150‑1 (MQ=255)

                                                                                                                  ccAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATCCGTCTGATACCAGTCGAATAGCGTTGc                                     >  1:3113675/1‑150 (MQ=255)

                                                                                                                       gATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATCCGTCTGATACCAGTCGAATAGCGTTGCcgcgc                                >  1:3033926/1‑150 (MQ=255)

                                                                                                                         tACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATCCGTCTGATACCAGTCGAATAGCGTTGCCgcgcgc                              >  1:2330627/1‑150 (MQ=255)

                                                                                                                            cGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATCCGTCTGATACCAGTCGAATAGCGTTGCCGCGCGCTCa                           <  1:3287775/150‑1 (MQ=255)

                                                                                                                            cGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATCCGTCTGATACCAGTCGAATAGCGTTGCCGCGCGCTCa                           <  2:4828511/150‑1 (MQ=255)

                                                                                                                             gATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATCCGTCTGATACCAGTCGAATAGCGTTGCCGCGCGCTCag                          <  2:4010613/150‑1 (MQ=255)

                                                                                                                             gATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATCCGTCTGATACCAGTCGAATAGCGTTGCCGCGCGCTCag                          >  1:5226272/1‑150 (MQ=255)

                                                                                                                              aTAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATCCGTCTGATACCATTCGAATAGCGTTGCCGCGCGCTCaga                         >  2:460719/1‑150 (MQ=255)

                                                                                                                              aTAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATCCGTCTGATACCAGTCGAATAGCGTTGCCGCGCGCTCaga                         >  1:3786759/1‑150 (MQ=255)

                                                                                                                                aaTGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATCCGTCTGATACCAGTCGAATAGCGTTGCCGCGCGCTCAGAGt                       >  1:2838239/1‑150 (MQ=255)

                                                                                                                                aaTGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATCCGTCTGATACCAGTCGAATAGCGTTGCCGCGCGCTCAGAGt                       >  1:4293579/1‑150 (MQ=255)

                                                                                                                                 aTGACGCACTGATGAGACTGGTCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATCCGTCTGATACCAGTCGAATAGCGTTGCCGCGCGCTCAGAGtt                      <  1:1430722/150‑1 (MQ=255)

                                                                                                                                  tGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATCCGTCTGATACCAGTCGAATAGCGTTGCCGCGCGCTCAGAGTTa                     >  1:3634851/1‑150 (MQ=255)

                                                                                                                                   gACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATCCGTCTGATACCAGTCGAATAGCGTTGCCGCGCGCTCAGAGTTaa                    <  1:3172495/150‑1 (MQ=255)

                                                                                                                                     cGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATCCGTCTGATACCAGTCGAATAGCGTTGCCGCGCGCTCAGAGTTAAtt                  <  2:4195732/150‑1 (MQ=255)

                                                                                                                                          tgatgaGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATCCGTCTGATACCAGTCGAATAGCGTTGCCGCGCGCTCAGAGTTAATTGTTGa             >  1:1816553/1‑150 (MQ=255)

                                                                                                                                           gatgaGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATCCGTCTGATACCAGTCGAATAGCGTTGCCGCGCGCTCAGAGTTAATTGTTGAc            >  1:4173744/1‑150 (MQ=255)

                                                                                                                                             tgaGACTGGTCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATCCGTCTGATACCAGTCGAATAGCGTTGCCGCGCGCTCAGAGTTAATTGTTGACaa          >  2:5529756/1‑150 (MQ=255)

                                                                                                                                             tgaGACTGATCAGTCATATATTTAGCGACCTGATTAACCTTGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATCCGTCTGATACCAGTCGAATAGCGTTGCCGCGCGCTCAGAGTTAATTGTTGACaa          <  2:2205303/150‑1 (MQ=255)

                                                                                                                                             tgaGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATCCGTCTGATACCAGTCGAATAGCGTTGCCGCGCGCTCAGAGTTAATTGTTGACaa          >  2:4849604/1‑150 (MQ=255)

                                                                                                                                                gaCTGGTCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACTCGGCAATTATAGGGATTTCATACGTCTGATACCAGTCGAATAGCGTTGCCTCGCGCTCAGAGTTAATTGTTGACAAAGa       >  2:195060/1‑150 (MQ=255)

                                                                                                                                                 aCTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATCCGTCTGATACCAGTCGAATAGCGTTGCCGCGCGCTCAGAGTTAATTGTTGACAAAGaa      >  2:3219977/1‑150 (MQ=255)

                                                                                                                                                  cTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATCCGTCTGATACCAGTCGAATAGCGTTGCCGCGCGCTCAGAGTTAATTGTTGACAAAGAAt     <  2:2048375/150‑1 (MQ=255)

TGATCTTTGTAATAGCAGTCTTCGGGAATTACGCCGATGTGTTCATCACCGACTTGCTCACGCAATTCACGATAAAGGGTACTGGCAATAAGACTCTTGCCGGAAGCCGATGCGCCAGCGATACCGATAATGACGCACTGATGAGACTGATCAGTCATATATTTAGCGACCTGATTAACCTGGATGTTAGGAAGGGGGCGACGAAGCGCCAAACGCGGCAATTATAGGGATTTCATCCGTCTGATACCAGTCGAATAGCGTTGCCGCGCGCTCAGAGTTAATTGTTGACAAAGAATTCC  >  CP009273/2136272‑2136570

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