BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	1,484,221:1	+T	100%	pseudogene (1/252 nt) pseudogene (750/750 nt)	gapC ← gapC ←	pseudogene, GAP dehydrogenase; glyceraldehyde‑3‑phosphate dehydrogenase (second fragment) pseudogene, GAP dehydrogenase; glyceraldehyde‑3‑phosphate dehydrogenase (second fragment)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	1,484,221	1	.	T	99.1%	412.5 / ‑3.1	106	pseudogene (1/252 nt) pseudogene (750/750 nt)	gapC gapC	pseudogene, GAP dehydrogenase; glyceraldehyde‑3‑phosphate dehydrogenase (second fragment) pseudogene, GAP dehydrogenase; glyceraldehyde‑3‑phosphate dehydrogenase (second fragment)
Reads supporting (aligned to +/- strand): ref base . (1/0); new base T (62/43); total (63/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.35e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

Only 100 of 116 total aligned reads displayed.

GGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCAC‑TTTTTTCCGAGAATCGATACCAGTTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTGGGGAATGATATTT  >  CP009273/1484072‑1484365

ggcggTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCAc                                                                                                                                                   >  1:778327/1‑150 (MQ=255)

  cggTAATTTCCGTTTGCGTGGCATCAAGCGCCGAACGGAGATGGCTGCCAATAATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAATGCGTTAGTCACGTCTTCGGCAGTCAC‑tt                                                                                                                                                <  2:366752/150‑1 (MQ=255)

  cggTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCAC‑tt                                                                                                                                                >  2:768599/1‑150 (MQ=255)

        tttCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCAC‑tttttttc                                                                                                                                          >  2:1006046/1‑148 (MQ=255)

         ttCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCAC‑tttttttcc                                                                                                                                         >  2:924934/1‑147 (MQ=255)

         ttCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCAC‑tttttttcc                                                                                                                                         >  1:715079/1‑147 (MQ=255)

           ccGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCAC‑tttttttcc                                                                                                                                         >  1:287874/1‑145 (MQ=255)

            cGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCAC‑tttttttcc                                                                                                                                         >  1:644352/1‑144 (MQ=255)

                tGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTTCCGAGAAtg                                                                                                                                   >  1:12203/1‑149 (MQ=255)

                  cGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTTCCGAGAAtgg                                                                                                                                  <  1:274404/149‑3 (MQ=255)

                  cGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTTCCGAGAATGGa                                                                                                                                 >  2:910653/1‑150 (MQ=255)

                     ggCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTTCCGAGAATGGacac                                                                                                                              >  1:528452/1‑150 (MQ=255)

                      gCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTTCCGAGAATGGacac                                                                                                                              >  1:578114/1‑149 (MQ=255)

                      gCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTTCCGAGAATGGACAcc                                                                                                                             <  1:55626/150‑1 (MQ=255)

                       cATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTTCCGAGAATGGACACCa                                                                                                                            <  1:487447/150‑1 (MQ=255)

                         tCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTTCCGAGAATGGACACCAgc                                                                                                                          >  2:888596/1‑149 (MQ=255)

                         tCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTTCCGAGAATGGACACCAgc                                                                                                                          <  1:56936/150‑2 (MQ=255)

                           aaaCACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTTCCGAGAATGGACACCAgct                                                                                                                         >  2:439429/1‑147 (MQ=255)

                            aaCACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTTCCGAGAATGGACACCAgct                                                                                                                         >  2:53615/1‑146 (MQ=255)

                            aaCACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTTCCGAGAATGGACACCAGCTCa                                                                                                                       <  2:469555/150‑1 (MQ=255)

                             acacCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTTCCGAGAATGGACACCAGCTCAg                                                                                                                      >  2:1013172/1‑150 (MQ=255)

                              cacCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTTCCGAGAATGGACACCAGCTCAGt                                                                                                                     >  1:184227/1‑150 (MQ=255)

                              cacCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTTCCGAGAATGGACACCAGCTCAGt                                                                                                                     <  1:858531/150‑1 (MQ=255)

                                    accgaaATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTTCCGAGAATGGACACCAGCTCAGTGACCGa                                                                                                               <  2:168570/150‑1 (MQ=255)

                                      cgaaATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTTCCGAGAATGGACACCAGCTCAGTGACCGAAc                                                                                                             <  1:971534/150‑1 (MQ=255)

                                       gaaATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTTCCGAGAATGGACACCAGCTCAGTGACCGAAcc                                                                                                            <  1:756350/150‑1 (MQ=255)

                                       gaaATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTTCCGAGAATGGACACCAGCTCAGTGACCGAAcc                                                                                                            >  1:668963/1‑150 (MQ=255)

                                         aaTGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTTCCGAGAATGGACACCAGCTCAGTGACCGAACCTg                                                                                                          >  1:136899/1‑150 (MQ=255)

                                          aTGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTTCCGAGAATGGACACCAGCTCAGTGACCGAACCTGt                                                                                                         >  2:518064/1‑150 (MQ=255)

                                                 ccAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTTCCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCAcc                                                                                                  >  2:441871/1‑150 (MQ=255)

                                                   aaTGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTTCCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCgg                                                                                                <  2:364378/150‑1 (MQ=255)

                                                    aTGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTTCCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCgg                                                                                                >  1:517465/1‑149 (MQ=255)

                                                    aTGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTTCCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGc                                                                                               <  1:791332/150‑1 (MQ=255)

                                                     tGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTTCCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCa                                                                                              >  2:214799/1‑150 (MQ=255)

                                                     tGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTTCCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCa                                                                                              >  1:1011210/1‑150 (MQ=255)

                                                       atatCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTTCCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCAcg                                                                                            >  2:651096/1‑150 (MQ=255)

                                                       atatCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTTCCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCAcg                                                                                            <  2:12203/150‑1 (MQ=255)

                                                       atatCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTTCCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCAc                                                                                             <  1:367206/149‑1 (MQ=255)

                                                         atCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGGGGGTGCTTGTTTAAGTGCGTTATTCCCCTCTTCGGCAGTCACTTTTTTTCCGAGAATGGACACCAGCTCAGTGACCCAACCTGTTTTCACCGGCCcgcg                                                                                          >  2:42127/1‑150 (MQ=255)

                                                          tCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTTCCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCAcgcg                                                                                          <  2:371079/149‑1 (MQ=255)

                                                          tCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTTCCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCAcgcg                                                                                          <  2:633608/149‑1 (MQ=255)

                                                            ggAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTTCCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTg                                                                                       >  1:803878/1‑150 (MQ=255)

                                                               agagACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTTCCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTgcg                                                                                     >  1:3975/1‑149 (MQ=255)

                                                                gagaCTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTTCCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCa                                                                                   >  2:431638/1‑150 (MQ=255)

                                                                gagaCTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTTCCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCa                                                                                   >  2:475946/1‑150 (MQ=255)

                                                                gagaCTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTTCCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCa                                                                                   >  2:252885/1‑150 (MQ=255)

                                                                       tttCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTTCCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCtt                                                                            <  2:427074/150‑1 (MQ=255)

                                                                         tcttcATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTTCCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCttt                                                                           <  2:968767/149‑1 (MQ=255)

                                                                             cATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTTCCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAg                                                                        <  1:214799/148‑1 (MQ=255)

                                                                             cATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTTCCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAg                                                                        <  2:460629/148‑1 (MQ=255)

                                                                             cATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTTCCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGt                                                                       <  1:79056/149‑1 (MQ=255)

                                                                               tCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTTCCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTg                                                                    >  1:849362/1‑150 (MQ=255)

                                                                                 ggTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTTCCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGcc                                                                  <  1:1021624/150‑1 (MQ=255)

                                                                                  gTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTTCCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCg                                                                 <  1:439429/150‑1 (MQ=255)

                                                                                    ataACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTTCCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCt                                                               >  2:834158/1‑150 (MQ=255)

                                                                                    ataACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTTCCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCt                                                               >  1:226872/1‑150 (MQ=255)

                                                                                     taACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTTCCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCt                                                               <  2:383657/149‑1 (MQ=255)

                                                                                       aCCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTTCCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAg                                                            >  2:882714/1‑150 (MQ=255)

                                                                                       aCCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTTCCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAg                                                            >  1:733809/1‑150 (MQ=255)

                                                                                       aCCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTTCCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAg                                                            >  2:768848/1‑150 (MQ=255)

                                                                                       aCCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTTCCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTCGCCGCTCAg                                                            >  2:786005/1‑150 (MQ=255)

                                                                                          aaaTGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTTCCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTc                                                         >  1:593050/1‑150 (MQ=255)

                                                                                             tGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTTCCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCgg                                                      <  2:61225/150‑1 (MQ=255)

                                                                                             tGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTTCCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCAgg                                                      <  1:568462/150‑1 (MQ=255)

                                                                                              gACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTTCCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCgg                                                      >  2:316710/1‑149 (MQ=255)

                                                                                                   gTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTTCCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCAc                                                >  2:685451/1‑150 (MQ=255)

                                                                                                   gTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTTCCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCAc                                                >  1:41043/1‑150 (MQ=255)

                                                                                                    ttATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTTCCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCAcc                                               <  2:1011210/150‑1 (MQ=255)

                                                                                                      aTTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTTCCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCAc                                                >  2:205837/1‑147 (MQ=255)

                                                                                                         ggtggtTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTTCCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGAc                                           >  1:1019838/1‑149 (MQ=255)

                                                                                                          gtggtTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTTCCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGAcc                                          <  2:644352/149‑1 (MQ=255)

                                                                                                          gtggtTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTTCCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGAcc                                          >  1:798481/1‑149 (MQ=255)

                                                                                                             gtTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTTCCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATg                                      <  1:668758/150‑1 (MQ=255)

                                                                                                              ttGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTTCCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATgg                                     >  2:491810/1‑150 (MQ=255)

                                                                                                                gCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTTCCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCt                                   >  2:143997/1‑150 (MQ=255)

                                                                                                                 cTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTTCCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCtt                                  <  1:837415/150‑1 (MQ=255)

                                                                                                                 cTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTTCCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCtt                                  <  2:287874/150‑1 (MQ=255)

                                                                                                                  ttGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTTCCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCttt                                 >  2:216322/1‑150 (MQ=255)

                                                                                                                  ttGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTTCCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCttt                                 >  2:581315/1‑150 (MQ=255)

                                                                                                                   tGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTTCCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCtttt                                >  1:154403/1‑150 (MQ=255)

                                                                                                                    gTTTAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTTCCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCttt                                 <  2:398216/148‑1 (MQ=255)

                                                                                                                     tttAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTTCCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTgc                              <  1:888596/150‑1 (MQ=255)

                                                                                                                     tttAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTTCCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTg                               <  1:611591/149‑1 (MQ=255)

                                                                                                                       tAAGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTTCCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTgccg                            >  2:719559/1‑150 (MQ=255)

                                                                                                                        aaGTGCGTTATTCACCTCTTCGGCAGTCACTTTTTTTCCGAGAATGGACACAAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTgccgc                           >  2:751755/1‑150 (MQ=255)

                                                                                                                        aaGTGCGTTATTCACCTCTTCGGCAGTCAC‑TTTTTTCCGAGAATCGATACCAGTTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTgccgcc                          >  1:346774/1‑150 (MQ=255)

                                                                                                                            gCGTTATTCACCTCTTCGGCAGTCACTTTTTTTCCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGccccc                       >  1:345964/1‑150 (MQ=255)

                                                                                                                             cGTTATTCACCTCTTCGGCAGTCACTTTTTTTCCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGccccc                       <  1:12382/149‑1 (MQ=255)

                                                                                                                              gTTATTCACCTCTTCGGCAGTCACTTTTTTTCCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCg                      >  1:897057/1‑149 (MQ=255)

                                                                                                                                    cACCTCTTCGGCAGTCACTTTTTTTCCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAgtgt                <  1:205597/149‑1 (MQ=255)

                                                                                                                                    cACCTCTTCGGCAGTCACTTTTTTTCCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGtgtg               >  2:970588/1‑150 (MQ=255)

                                                                                                                                     aCCTCTTCGGCAGTCACTTTTTTTCCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGtgtg               <  1:1006046/149‑1 (MQ=255)

                                                                                                                                       ctctTCGGCAGTCACTTTTTTTCCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTgggg            <  1:768848/150‑1 (MQ=255)

                                                                                                                                       ctctTCGGCAGTCACTTTTTTTCCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTggg             >  1:1026614/1‑149 (MQ=255)

                                                                                                                                       ctctTCGGCAGTCACTTTTTTTCCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTggg             >  1:471587/1‑149 (MQ=255)

                                                                                                                                        tctTCGGCAGTCACTTTTTTTCCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTgggg            >  2:27911/1‑149 (MQ=255)

                                                                                                                                        tctTCGGCAGTCACTTTTTTTCCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTgggg            >  1:202997/1‑149 (MQ=255)

                                                                                                                                            cGGCAGTCACTTTTTTTCCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTGGGGAATGa       >  2:41322/1‑150 (MQ=255)

                                                                                                                                             ggCAGTCACTTTTTTTCCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTGGGGAATGat      <  2:799571/150‑1 (MQ=255)

                                                                                                                                               cAGTCACTTTTTTTCCGAGAATGGACACCAGCTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTGGGGAATGatat    >  1:694918/1‑150 (MQ=255)

GGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGGCAGTCAC‑TTTTTTCCGAGAATCGATACCAGTTCAGTGACCGAACCTGTTTTCACCGGCACGCGTTGCGCATGACCTTTCAGTTTGCCGCTCAGTTCCGGGATCACCAGACCAATGGCTTTTGCCGCCCCCGTAGTGTGGGGAATGATATTT  >  CP009273/1484072‑1484365

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