BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	3,116,216	G→T	5.2%	R317S (CGT→AGT)	glcB ←	malate synthase G

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	3,116,216	0	G	T	5.2%	1129.8 / 28.4	297	R317S (CGT→AGT)	glcB	malate synthase G
Reads supporting (aligned to +/- strand): ref base G (158/123); new base T (9/7); total (167/130)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 302 total aligned reads displayed.

CATGACGCCATCAAGAATGCCTTCCGGGATTTCATTGCCTTCGCTGTCCCAAATCACAGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGTTGCGGTACAGCAGGATTTTATCTTCCGCATCAACCGCCGCGACCGAATCTTCGCAGTCGAGAATGG  >  CP009273/3116067‑3116365

cATGACGCCATCAAGAATGCCTTCCGGGATTTCATTGCCTTCGCTGTCCCAAATCACAGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACg                                                                                                                                                       >  2:225769/1‑150 (MQ=255)

  tGACGCCATCAAGAATGCCTTCCGGGATTTCATTGCCTTCGCTGTCCCAAATCACAGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGAt                                                                                                                                                     <  2:2096249/150‑1 (MQ=255)

  tGACGCCATCAAGAATGCCTTCCGGGATTTCATTGCCTTCGCTGTCCCAAATCACAGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGAt                                                                                                                                                     >  1:4835770/1‑150 (MQ=255)

  tGACGCCATCAAGAATGCCTTCCGGGATTTCATTGCCTTCGCTGTCCCAAATCACAGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGAt                                                                                                                                                     <  2:5119057/150‑1 (MQ=255)

   gACGCCATCAAGAATGCCTTCCGGGATTTCATTGCCTTCGCTGTCCCAAATCACAGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATc                                                                                                                                                    <  2:2509799/150‑1 (MQ=255)

       ccATCAAGAATGCCTTCCGGGATTTCATTGCCTTCGCTGTCCCAAATCACAGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCa                                                                                                                                                >  2:4140677/1‑150 (MQ=255)

            aaGAATGCCTTCCGGGATTTCATTGCCTTCGCTGTCCCAAATCACAGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAg                                                                                                                                           >  2:695170/1‑150 (MQ=255)

            aaGAATGCCTTCCGGGATTTCATTGCCTTCGCTGTCCCAAATCACAGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAg                                                                                                                                           >  2:5064106/1‑150 (MQ=255)

                aTGCCTTCCGGGATTTCATTGCCTTCGCTGTCCCAAATCACAGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTa                                                                                                                                       >  1:4331058/1‑150 (MQ=255)

                  gCCTTCCGGGATTTCATTGCCTTCGCTGTCCCAAATCACAGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACg                                                                                                                                     >  2:607593/1‑150 (MQ=255)

                      tCCGGGATTTCATTGCCTTCGCTGTCCCAAATCACAGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACg                                                                                                                                 >  2:3974540/1‑150 (MQ=255)

                       ccGGGATTTCATTGCCTTCGCTGTCCCAAATCACAGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGa                                                                                                                                >  1:3103971/1‑150 (MQ=255)

                        cGGGATTTCATTGCCTTCGCTGTCCCAAATCACAGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGAt                                                                                                                               >  1:3924297/1‑150 (MQ=255)

                         gggATTTCATTGCCTTCGCTGTCCCAAATCACAGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGAtt                                                                                                                              <  1:1100747/150‑1 (MQ=255)

                          ggATTTCATTGCCTTCGCTGTCCCAAATCACAGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACTATCGTCATTCAGTTTACGCACGAttt                                                                                                                             <  1:3777548/150‑1 (MQ=255)

                          ggATTTCATTGCCTTCGCTGTCCCAAATCACAGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGAttt                                                                                                                             <  2:2724033/150‑1 (MQ=255)

                           gATTTCATTGCCTTCGCTGTCCCAAATCACAGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTg                                                                                                                            <  1:1138864/150‑1 (MQ=255)

                            aTTTCATTGCCTTCGCTGTCCCAAATCACAGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACTATCGTCATTCAGTTTACGCACGATTTGc                                                                                                                           <  1:4790482/150‑1 (MQ=255)

                             tttCATTGCCTTCGCTGTCCCAAATCACAGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGcc                                                                                                                          >  1:2336341/1‑150 (MQ=255)

                             tttCATTGCCTTCGCTGTCCCAAATCACAGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGcc                                                                                                                          <  1:1817953/150‑1 (MQ=255)

                               tCATTGCCTTCGCTGTCCCAAATCACAGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGa                                                                                                                        >  1:2715528/1‑150 (MQ=255)

                                  ttGCCTTCGCTGTCCCAAATCACAGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCg                                                                                                                     <  2:4546435/150‑1 (MQ=255)

                                    gCCTTCGCTGTCCCAAATCACAGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGtt                                                                                                                   >  1:1192172/1‑150 (MQ=255)

                                       ttCGCTGTCCCAAATCACAGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGttttt                                                                                                                <  1:686660/150‑1 (MQ=255)

                                        tCGCTGTCCCAAATCACAGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTtc                                                                                                               <  1:3574512/150‑1 (MQ=255)

                                         cGCTGTCCCAAATCACAGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTCATGACGATCGTCATTCAGTTTACGCACGATGTGCCGACCGTTTTtct                                                                                                              >  2:3227337/1‑150 (MQ=255)

                                         cGCTGTCCCAAATCACAGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTtct                                                                                                              >  2:351719/1‑150 (MQ=255)

                                          gCTGTCCCAAATCACAGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTtctc                                                                                                             <  2:305842/150‑1 (MQ=255)

                                           cTGTCCCAAATCACAGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTcc                                                                                                            >  2:4268537/1‑150 (MQ=255)

                                           cTGTCCCAAATCACAGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTcc                                                                                                            >  2:4942565/1‑150 (MQ=255)

                                            tGTCCCAAATCACAGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCa                                                                                                           >  1:230644/1‑150 (MQ=255)

                                                ccAAATCACAGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCAtttt                                                                                                       >  1:2425792/1‑150 (MQ=255)

                                                 cAAATCACAGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTtc                                                                                                      <  1:2019893/150‑1 (MQ=255)

                                                  aaaTCACAGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGCGCGTCAGTCAGTTTACGCACGATTTGCCGACCGTTTTGCTACaggtgct                                                                                                     >  2:3112052/1‑144 (MQ=255)

                                                  aaaTCACAGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTtct                                                                                                     <  1:1007634/150‑1 (MQ=255)

                                                  aaaTCACAGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTtct                                                                                                     >  2:3112185/1‑150 (MQ=255)

                                                   aaTCACAGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTtctc                                                                                                    <  2:1004104/150‑1 (MQ=255)

                                                    aTCACAGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTctct                                                                                                   >  1:3516252/1‑150 (MQ=255)

                                                        caGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGAACGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCa                                                                                               <  2:3810710/150‑1 (MQ=255)

                                                         aGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCTGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCag                                                                                              >  2:4391394/1‑150 (MQ=255)

                                                          ggAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCaga                                                                                             >  2:1453794/1‑150 (MQ=255)

                                                            aaTGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGt                                                                                           >  1:2500056/1‑150 (MQ=255)

                                                               ggTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTccc                                                                                        <  1:1478139/150‑1 (MQ=255)

                                                                gTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCACTTTTCTCTTGCAGAGTcccc                                                                                       >  2:4978303/1‑150 (MQ=255)

                                                                   atcaAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACTATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGc                                                                                    >  1:818766/1‑150 (MQ=255)

                                                                    tcaAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCa                                                                                   >  1:126748/1‑150 (MQ=255)

                                                                      aaaTGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATc                                                                                 >  1:624155/1‑150 (MQ=255)

                                                                       aaTGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCa                                                                                >  1:389115/1‑150 (MQ=255)

                                                                       aaTGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCa                                                                                >  2:3961428/1‑150 (MQ=255)

                                                                           aCCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGcc                                                                            >  1:1575516/1‑150 (MQ=255)

                                                                                cGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACTATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCcagca                                                                       >  2:4712202/1‑150 (MQ=255)

                                                                                  ttGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAgg                                                                     <  2:1325916/150‑1 (MQ=255)

                                                                                   tGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGt                                                                    >  1:2031347/1‑150 (MQ=255)

                                                                                       gATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCGGGTTGCg                                                                <  1:3837728/150‑1 (MQ=255)

                                                                                       gATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGTTGCg                                                                >  2:2329545/1‑150 (MQ=255)

                                                                                           aaCAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGTTGCGGTAc                                                            >  2:4098304/1‑150 (MQ=255)

                                                                                           aaCAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGACCAGCAGGTTGCGGTAc                                                            <  1:1379700/150‑1 (MQ=255)

                                                                                             cagcagCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGTTGCGGTAcag                                                          >  1:3158794/1‑150 (MQ=255)

                                                                                              agcagcGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGTTGCGGTAcagc                                                         >  1:4894390/1‑150 (MQ=255)

                                                                                              agcagcGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGTTGCGGTAcagc                                                         >  2:553486/1‑150 (MQ=255)

                                                                                                    gagcgTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGTTGCGGTACAGCAGGAtt                                                   <  1:1964390/150‑1 (MQ=255)

                                                                                                    gagcgTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGTTGCGGTACAGCAGGAtt                                                   <  1:3447404/150‑1 (MQ=255)

                                                                                                     agcgTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGTTGCGGTACAGCAGGAttt                                                  <  1:2218226/150‑1 (MQ=255)

                                                                                                     agcgTCCGTGCAGAGAAATGTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGTTGCGGTACAGCAGGAttt                                                  <  2:5164282/150‑1 (MQ=255)

                                                                                                         tCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGTTGCGGTACAGCAGGATTTTATc                                              >  2:4619843/1‑150 (MQ=255)

                                                                                                         tCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGTTGCGGTACAGCAGGATTTTATc                                              <  1:3695429/150‑1 (MQ=255)

                                                                                                          ccGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGTTGCGGTACAGCAGGATTTTATCt                                             <  1:5115801/150‑1 (MQ=255)

                                                                                                          ccGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGTTGCGGTACAGCAGGATTTTATCt                                             >  2:281562/1‑150 (MQ=255)

                                                                                                          ccGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGGTGCGGTACAGCAGGATTTTATCt                                             >  2:4568837/1‑150 (MQ=255)

                                                                                                           cGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGTTGCGGTACAGCAGGATTTTATCtt                                            >  2:20415/1‑150 (MQ=255)

                                                                                                            gTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATTACGATCGTCATTCAGATTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGTTGCGTTACAGCAGGATTTTATCTTc                                           <  2:4740690/150‑1 (MQ=255)

                                                                                                            gTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATGTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGTTGCGGTACAGCAGGATTTTATCTTc                                           <  1:1864110/150‑1 (MQ=255)

                                                                                                             tGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGTTGCGGTCCAGCAGGATTTTATCTTcc                                          >  2:1573518/1‑150 (MQ=255)

                                                                                                                agagAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGTTGCGGTACAGCAGGATTTTATCTTCCGCa                                       >  2:755877/1‑150 (MQ=255)

                                                                                                                agagAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGTTGCGGTACAGCAGGATTTTATCTTCCGCa                                       <  1:4184219/150‑1 (MQ=255)

                                                                                                                 gagaAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGTTGCGGTACAGCAGGATTTTATATTCCGCAt                                      >  2:2944069/1‑150 (MQ=255)

                                                                                                                    aaaTTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGTTGCGGTACAGCAGGATTTTATCTTCCGCATCaa                                   <  2:129752/150‑1 (MQ=255)

                                                                                                                    aaaTTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGTTGCGGTACAGCAGGATTTTATCTTCCGCATCaa                                   <  2:4215470/150‑1 (MQ=255)

                                                                                                                        ttCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGTTGCGGTACAGCAGGATTTTATCTTCCGCATCAAccgc                               >  1:1666170/1‑150 (MQ=255)

                                                                                                                        ttCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGTTGCGGTACAGCAGGATTTTATCTTCCGCATCAAccgc                               <  1:3241789/150‑1 (MQ=255)

                                                                                                                         tCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGTTGCGGTACAGCAGGATTTTATCTTCCGCATCAAccgcc                              >  2:207161/1‑150 (MQ=255)

                                                                                                                            gagCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGTTGCGGTACAGCAGGATTTTATCTTCCGCATCAACCGCcgcg                           <  1:2788566/150‑1 (MQ=255)

                                                                                                                                cATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGACCAGCAGGTTGCGGTACAGCAGGATTTTATCTTCCGCATCAACCGCCGCGACCg                       >  1:4580557/1‑150 (MQ=255)

                                                                                                                                 ataggcggcggTTTAATGACTATCGAGAGTAAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTAGCAGATTCCCCTGCATAAGGCCCAGCAGGATGCGGTACAGAATGATTTTATCTTCCGCATCAACCGACGCGAACGa                      >  2:1814905/4‑150 (MQ=255)

                                                                                                                                 aTCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGTTGCGGTACAGCAGGATTTTATCTTCCGCATCAACCGCCGCGACCGa                      >  2:29876/1‑150 (MQ=255)

                                                                                                                                 aTCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGTTGCGGTACAGCAGGATTTTATCTTCCGCATCAACCGCCGCGACCGa                      >  2:4131021/1‑150 (MQ=255)

                                                                                                                                  tCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGTTGCGGTACAGCAGGATTTTATCTTCCGCATCAACCGCCGCGACCGaa                     >  2:4033465/1‑150 (MQ=255)

                                                                                                                                   cggcggcggTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGTTGCGGTACAGCAGGATTTTATCTTCCGCATCAACCGCCGCGACCGAAt                    <  1:3018899/150‑1 (MQ=255)

                                                                                                                                   cggcggcggTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGTTGCGGTACAGCAGGATTTTATCTTCCGCATCAACCGCCGCGACCGAAt                    >  2:2912493/1‑150 (MQ=255)

                                                                                                                                       ggcggTGTAATGACTATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGTTGCGGTACAGCAGGATTTTATCTTCCGCATCAACCGCCGCGACCGAATCTTc                >  1:676222/1‑150 (MQ=255)

                                                                                                                                       ggcggTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGTTGCGGTACAGCAGGATTTTATCTTCCGCATCAACCGCCGCGACCGAATCTTc                >  2:3277010/1‑150 (MQ=255)

                                                                                                                                         cggTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGTTGCGGTACAGCAGGATTTTATCTTCCGCATCAACCGCCGCGACCGAATCTTCGc              <  1:3270670/150‑1 (MQ=255)

                                                                                                                                          ggTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGTTGCGGTACAGCAGGATTTTATCTTCCGCATCAACCGCCGCGACCGAATCTTCGCa             >  2:3339694/1‑150 (MQ=255)

                                                                                                                                             gtAATGACTATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGTTGCGGTACAGCAGGATTTTATCTTCCGCATCAACCGCCGCGACCGAATCTTCGCAGTc          >  2:4841443/1‑150 (MQ=255)

                                                                                                                                             gtAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGTTGCGGTACAGCAGGATTTTATCTTCCGCATCAACCGCCGCGACCGAATCTTCGCAGTc          <  1:2921668/150‑1 (MQ=255)

                                                                                                                                              tAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGTTGCGGTACAGCAGGATTTTATCTTCCGCATCAACCGCCGCGACCGAATCTTCGCAGTCg         <  2:236935/150‑1 (MQ=255)

                                                                                                                                               aaTGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGTTGCGGTACAGCAGGATTTTATCTTCCGCATCAACCGCCGCGACCGAATCTTCGCAGTCga        >  1:4235504/1‑150 (MQ=255)

                                                                                                                                                aTGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCAATTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGTTGCGGTACAGCAGGATTTTATCTTCCGCATCAACCGCCGCGACCGAATCTTCGCAGTCgag       <  2:1477133/150‑1 (MQ=255)

                                                                                                                                                    cGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGTTGCGGTACAGCAGGATTTTATCTTCCGCATCAACCGCCGCGACCGAATCTTCGCAGTCGAGAATg   <  1:4068858/150‑1 (MQ=255)

                                                                                                                                                     gATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGTTGCGGTACAGCAGGATTTTATCTTCCGCATCAACCGCCGCGACCGAATCTTCGCAGTCGAGAATgg  <  2:3934663/150‑1 (MQ=255)

CATGACGCCATCAAGAATGCCTTCCGGGATTTCATTGCCTTCGCTGTCCCAAATCACAGGAATGGTCATCAAATGACCCACGTTGCGGATAAACAGCAGCGAGCGTCCGTGCAGAGAAATTTCAGAGCCATCGGCGGCGGTGTAATGACGATCGTCATTCAGTTTACGCACGATTTGCCGACCGTTTTTCTCCATTTTCTCTTGCAGAGTCCCCTGCATCAGGCCCAGCAGGTTGCGGTACAGCAGGATTTTATCTTCCGCATCAACCGCCGCGACCGAATCTTCGCAGTCGAGAATGG  >  CP009273/3116067‑3116365

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