BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	1,606,883	A→C	6.8%	F209C (TTT→TGT)	glsB ←	glutaminase 2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	1,606,883	0	A	C	6.8%	603.3 / 5.4	223	F209C (TTT→TGT)	glsB	glutaminase 2
Reads supporting (aligned to +/- strand): ref base A (108/94); new base C (2/19); total (110/113)
Fisher's exact test for biased strand distribution p-value = 1.16e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 227 total aligned reads displayed.

TCGCCGGTAGCCCCACCCGCCAGGCAAACTCCCCCGCGTTCTGGTACATACCACTGGTCGCCATCAGCGCGTTAATTTGCCGCGCCTGCATTGGCGTCACCACTGGTTCATCAATATGAATAGCTTTCCCCTGATTAGCCAGAAAGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTCGCGCGGAATGTTCAAATTCGGA  >  CP009273/1606734‑1607031

 cGCCGGTAGCCCCACCCGCCAGGCAAACTCCCCCGCGTTCTGGTACATACCACTGGTCGCCATCAGCGCGTTAATTTGCCGCGCCTGCATTGGCGTCACCACTGGTTCATCAATATGAATAGCTTTCCCCTGATTAGCCAGAAAGACaaa                                                                                                                                                     <  1:2481898/150‑1 (MQ=255)

  gCCGGTAGCCCCACCCGCCAGGCAAACTCCCCCGCGTTCTGGTACATACCACTGGTCGCCATCAGCGCGTTAATTTGCCGCGCCTGCATTGGCGTCACCAGTGGTTCATCAATATGAATAGCTTTCCCCTGATTAGCCAGAAAGACAAAc                                                                                                                                                    <  2:5087650/150‑1 (MQ=255)

     ggTAGCCCCACCCGCCAGGCAAACTCCCCCGCGTTCTGGTACATACCACTGGTCGCCATCAGCGCGTTAATTTGCCGCGCCTGCATTGGCGTCACCACTGGTTCATCAATATGAATAGCTTTCCCCTGATTAGCCAGAAAGACAAACGTc                                                                                                                                                 <  1:3487907/150‑1 (MQ=255)

     ggTAGCCCCACCCGCCAGGCAAACTCCCCCGCGTTCTGGTACATACCACTGGTCGCCATCAGCGCGTTAATTTGCCGCGCCTGCATTGGCGTCACCACTGGTTCATCAATATGAATAGCTTTCCCCTGATTAGCCAGAAAGACAAACGTc                                                                                                                                                 >  1:5090599/1‑150 (MQ=255)

     ggTAGCCCCACCCGCCAGGCAAACTCCCCCGCGTTCTGGTACATACCACTGGTCGCCATCAGCGCGTTAATTTGCCGCGCCTGCATTGGCGTCACCACTGGTTCATCAATATGAATAGCTTTCCCCTGATTAGCCAGAAAGACAAACGTc                                                                                                                                                 <  2:2921267/150‑1 (MQ=255)

          ccccACCCGCCAGGCAAACTCCCCCGCGTTCTGGTACATACCACTGGTCGCCATCAGCGCGTTAATTTGCCGCGCCTGCATTGGCGTCACCACTGGTTCATCAATATGAATAGCTTTCCACTGATTAGCCAGAAAGACAAACGTCCGGGc                                                                                                                                            >  1:3951707/1‑150 (MQ=255)

           cccACCCGCCCGGCAAACTCCCCGGCGTTTTGGTACATACCACTGGTCGCCATCAGCGCGTTATTTTGCCGCGCCTGCTTTGGCGTCACCACTGGTTCATCAATATTAATCGCTTTCCCCTTATTAGCAATAAAGACAAACGTCCGGGcc                                                                                                                                           <  2:1606422/150‑1 (MQ=255)

           cccACCCGCCAGGCAAACTCCCCCGCGTTCTGGTACATACCACTGGTCGCCATCAGCGCGTTAATTTGCCGCGCCTGCATTGGCGTCACCACTGGTTCATCAATATGAATAGCTTTCCCCTGATTAGCCAGAAAGACAAACGTCCGGGcc                                                                                                                                           >  1:1843733/1‑150 (MQ=255)

            ccACCCGCCAGGCAAACTCCCCCGCGTTCTGGTACATACCACTGGTCGCCATCAGCGCGTTAATTTGCCGCGCCTGCATTGGCGTCACCACTGGTTCATCAATATGAATAGCTTTCCCCTGATTAGCCAGAAAGACAAACGTCCGGGCCa                                                                                                                                          >  2:1864852/1‑150 (MQ=255)

              aCCCGCCAGGCAAACTCCCCCGCGTTCTGGTACATACCACTGGTCGCCATCAGCGCGTTAATTTGCCGCGCCTGCATTGGCGTCACCACTGGTTCATCAATATGAATAGCTTTCCCCTGATTCGCCAGAACGACAAACGTCCGGGCCAGc                                                                                                                                        >  2:4996684/1‑150 (MQ=255)

              aCCCGCCAGGCAAACTCCCCCGCGTTCTGGTACATACCACTGGTCGCCATCAGCGCGTTAATTTGCCGCGCCTGCATTGGCGTCACCACTGGTTCATCAATATGAATAGCTTTCCCCTGATTAGCCAGAAAGACAAACGTCCGGGCCAGc                                                                                                                                        <  2:4627569/150‑1 (MQ=255)

                ccGCCAGGCAAACTCCCCCGCGTTCTGGTACATACCACTGGTCGCCATCAGCGCGTTAATTTGCCGCGCCTGCATTGGCGTCACCACTGGTTCATCAATATGAATAGCTTTCCCCTGATTAGCCAGAAAGACAAACGTCCGGGCCAGctc                                                                                                                                      <  2:4538697/150‑1 (MQ=255)

                 cGCCAGGCAAACTCCCCCGCGTTCTGGTACATACCACTGGTCGCCATCAGCGCGTTAATTTGCCGCGCCTGCATTGGCGTCACCACTGGTTCATCAATATGAATAGCTTTCCCCTGAGTAGCCAGAAAGACAAACGTCCGGGCCAGctct                                                                                                                                     >  2:76597/1‑150 (MQ=255)

                  gCCAGGCAAACTCCCCCGCGTTCTGGTACATACCACTGGTCGCCATCAGCGCGTTAATTTGCCGCGCCTGCATTGGCGTCACCACTGGTTCATCAATATCAATAGCTTTCCCCTGATTAGCCAGAAAGACAAACGTCCGGGCCAGCTCTa                                                                                                                                    <  2:3441430/150‑1 (MQ=255)

                    cAGGCAAACTCCCCCGCGTTCTGGTACATACCACTGGTCGCCATCAGCGCGTTAATTTGCCGCGCCTGCATTGGCGTCACCACTGGTTCATCAATATGAATAGCTTTCCCCTGATTAGCCAGAAAGACAAACGTCCGGGCCAGCTCTaca                                                                                                                                  >  1:651697/1‑150 (MQ=255)

                    cAGGCAAACTCCCCCGCGTTCTGGTACATACCACTGGTCGCCATCAGCGCGTTAATTTGCCGCGCCTGCATTGGCGTCACCACTGGTTCATAATTATTAATAGCTTTCCCCTGATTAGCCAGAAAGACAAACGTCCGGGCCAGCTCTaca                                                                                                                                  <  2:174156/150‑1 (MQ=255)

                     aGGCAAACTCCCCCGCGTTCTGGTACATACCACTGGTCGCCATCAGCGCGTTAATTTGCCGCGCCTGCATTGGCGTCACCACTGGTTCATCAATATGAATAGCTTTCCCCTGATTAGCCAGAAAGACAAACGTCCGGGCCAGCTCTacac                                                                                                                                 >  2:491895/1‑150 (MQ=255)

                      ggCAAACTCCCCCGCGTTCTGGTACATACCACTGGTCGCCATCAGCGCGTTAATTTGCCGCGCCTGCATTGGCGTCACCACTGGTTCATCAATATGAATAGCTTTCCCCTGATTAGCCAGAAAGACAAACGTCCGGGCCAGCTCTAcaca                                                                                                                                >  2:249760/1‑150 (MQ=255)

                      ggCAAACTCCCCCGCGTTCTGGTACATACCACTGGTCGCCATCAGCGCGTTAATTTGCCGCGCCTGCATTGGCGTCACCACTGGTTCATCAATATGAATAGCTTTCCCCTGATTAGCCAGAAAGACAAACGTCCGGGCCAGCTCTAcaca                                                                                                                                >  1:4648891/1‑150 (MQ=255)

                      ggCAAACTCCCCCGCGTTCTGGTACATACCACTGGTCGCCATCAGCGCGTTAATTTGCCGCGCCTGCATTGGCGTCACCACTGGTTCATCAATATGAATAGCTTTCCCCTGATTAGCCAGAAAGACAAACGTCCGGGCCAGCTCTAcaca                                                                                                                                >  2:3161644/1‑150 (MQ=255)

                       gCAAACTCCCCCGCGTTCTGGTACATACCACGGGTCGCCATCTGCGCGTTAATTGGCGGCGCCTGCATGGGCGTCACCACTGGTTCATCAATATTATTAGCTTTCCCCTGATTAGCCATAAAGACAAACGTCCGGGCCAGCTCTACACAg                                                                                                                               <  1:2886734/150‑1 (MQ=255)

                         aaaCTCCCCCGCGTTCTGGTACATACCACTGGTCGCCATCAGCGCGTTAATTTGCCGCGCCTGCATTGGCGTCACCACTGGTTCATCAATATGAATAGCTTTCCCCTGATTAGCCAGAAAGACAAACGTCCGGGCCAGCTCTACACAGCt                                                                                                                             >  2:201913/1‑150 (MQ=255)

                         aaaCTCCCCCGCGTTCTGGTACATACCACTGGTCGCCATCAGCGCGTTAATTTGCCGCGCCTGCATTGGCGTCACCACTGGTTCATCAATATGAATAGCTTTCCCCTGATTAGCCAGAAAGACAAACGTCCGGGCCAGCTCTACACAGCt                                                                                                                             >  1:3977873/1‑150 (MQ=255)

                            cTCCCCCGCGTTCTGGTACATACCACTGGTCGCCATCAGCGCGTTAATTTGCCGCGCCTGCATTGGCGTCACCACTGGTTCATCAATATGAATAGCTTTCCCCTGATTAGCCAGAAAGACAAACGTCCGGGCCAGCTCTACACAGCTCAt                                                                                                                          >  1:3440918/1‑150 (MQ=255)

                              cccccGCGTTCTGGTACATACCACTGGTCGCCATCAGCGCGTTAATTTGCCGCGCCTGCATTGGCGTCACCACTGGTTCATCAATATGAATAGCTTTCCCCTGATTAGCCAGAAAGACAAACGTCCGGGCCAGCTCTACACAGCTCAttt                                                                                                                        >  1:1228263/1‑150 (MQ=255)

                                cccGCGTTCTGGTACATACCACTGGTCGCCATCAGCGCGTTAATTTGCCGCGCCTGCATTGGCGTCACCACTGGTTCATCAATATGAATAGCTTTCCCCTGATTAGCCAGAAAGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTc                                                                                                                      >  2:2329452/1‑150 (MQ=255)

                                cccGCGTTCTGGTACATACCACTGGTCGCCATCAGCGCGTTAATTTGCCGCGCCTGCATTGGCGTCACCACTGGTTCATCAATATGAATAGCTTTCCCCTGATTAGCCAGAAAGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTc                                                                                                                      >  1:964263/1‑150 (MQ=255)

                                    cgTTCTGGTACATACCACTGGTCGCCATCAGCGCGTTAATTTGCCGCGCCTGCATTGGCGTCACCACTGGTTCATCAATATGAATAGCTTTCCCCTGATTAGCCAGAAAGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCagag                                                                                                                  >  1:3213850/1‑150 (MQ=255)

                                    cgTTCTGGTACATACCACTGGTCGCCATCAGCGCGTTAATTTGCCGCGCCTGCATTGGCGTCACCACTGGTTCATCAATATGAATAGCTTTCCCCTGATTAGCCAGAAAGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCagag                                                                                                                  >  2:2587091/1‑150 (MQ=255)

                                     gTTCTGGTACATACCACTGGTCGCCATCAGCGCGTTAATTTGCCGCGCCTGCATTGGCGTCACCACTGGTTCATCAATATGAATAGCTTTCCCCTGATTAGCCAGAAAGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAgc                                                                                                                 >  2:1881870/1‑150 (MQ=255)

                                      ttCTGGTACATACCACTGGTCGCCATCAGCGCGTTAATTTGCCGCGCCTGCATTGGCGTCACCACTGGTTCATCAATATGAATAGCTTTCCCCTGATTAGCCAGAAAGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAgcg                                                                                                                >  2:1200652/1‑150 (MQ=255)

                                       tCTGGTACATACCCCGGGTCGCCATCAGCGCGTTTATTTGCCGCGCCTGCATTGGCGTCACCACTGGTTCAACAATATGAATAGCTTTCCCCTGTTTAGCAATAAAGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAgcgc                                                                                                               <  2:3053255/150‑1 (MQ=255)

                                         tGGTACATACCACTGGTCGCCATCAGCGCGTTAATTTGCCGCGCCTGCATTGGCGTCACCACTGGTTCATCAATATGAATAGCTTTCCCCTGATTAGCCAGAAAGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCACAg                                                                                                             <  2:4747608/150‑1 (MQ=255)

                                           gTACATACCACTGGTCGCCATCAGCGCGTTAATTTGCCGCGCCTGCATTGGCGTCACCACTGGTTCATCAATATGAATAGCTTTCCCCTGATTAGCCAGAAAGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTa                                                                                                           >  2:1732550/1‑150 (MQ=255)

                                            tACATCCCAGTGGTCGCCATCTGCGCTTTATTTTGCCGCGCTTGCTTTGGCTTCACCACTGTTTAATCAATATTATTAGCTTTCCCCTTTTTAGCCATAAAGACACACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTaa                                                                                                          <  1:1427767/150‑1 (MQ=255)

                                              cATACCACTGGCCGCAATCAGCGCGTTAATTTGCCGCGCCTGATTTGGCTTCACCACTGGTTCATAAATATTAATAGCTTCCCCCTTATTAGCCACAAACACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATg                                                                                                        <  1:1382664/150‑1 (MQ=255)

                                               aTACCACTGGTCGCCATCAGCGCGTTAATTTGCCGCGCCTGCATTGGCGTCACCACTGGTTCATCAATATGAATAGCTTTCCCCTGATTAGCCAGAAAGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGa                                                                                                       >  2:5095645/1‑150 (MQ=255)

                                                 aCCACTGGTCGCCATCAGCGCGTTAATTTGCCGCGCCTGCATTGGCGTCACCACTGGTTCATCAATATGAATAGCTTTCCCCTGATTAGCCAGAAAGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGaaa                                                                                                     >  2:696518/1‑150 (MQ=255)

                                                  ccccTGGTCGCCATCAGCGCGTTATTTTGCCGCGCCTGCATGGGCGTCACCACTGGTTCATCAATATGAATAGCTTTCCCCTGATTAGCCAGAAAGACACACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAag                                                                                                    <  1:2918436/147‑1 (MQ=255)

                                                   caggggTCGCCATCTGCGGTTTATTTTGCCGCGCCTGCTTGGGCGTCACCACTGGTTAATCAATATGAATAGCTTCCCCCTTTTTCGCCAGAAAGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAagt                                                                                                   <  1:49563/146‑1 (MQ=255)

                                                   cACTGGTCGCCATCAGCGCTTTAATTTGCCGCGCCTGCATTGGCTTCACCACTGTTTCATCAATATTAATAGCTTTCCCCTGATTACCCATAAATACACACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAagt                                                                                                   <  1:49421/150‑1 (MQ=255)

                                                     cTGGTCGCCATCAGCGCGTTATTTTGCCGCGCCTGCATGGGCGTCACCACTGGTTCATCAATATGAATAGCTTTCCCCTGATTAGCCATAAAGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAagtag                                                                                                 <  1:2550586/150‑1 (MQ=255)

                                                      gggTCCCCACCAGCCCCTTCATTTGCCGCGCCTGCATGGGCGTCACCACTGGTTCAACAAGCCGAATAGCTTTCCCCTGAATAGCCAGAAAGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAagtagt                                                                                                <  2:3840088/149‑1 (MQ=255)

                                                         tCGCCATCAGCGCGTTAATTTGCCGCGCCTGCATTGGCGTCACCACTGGTTCATCAATATGAATAGCGTTCCCCTGATTAGCCAGAAAGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGtttt                                                                                             <  2:3237589/150‑1 (MQ=255)

                                                           gCCATCAGCGCGTTAATTTGCCGCGCCTGCTTTGGCGTCACCACTGGTTCATCAATATGAATAGCTTTCCCCTGATTAGACATAAAGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGc                                                                                           <  2:3406134/150‑1 (MQ=255)

                                                                ccgcgcTTTATTTTGCCGCGCCTGCTTGGGCGTCACCACTGTTTCATAAATATTAATAGCTGTCCCCTGATTAGCCACAAAGACACACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAAc                                                                                      <  1:2289637/148‑1 (MQ=255)

                                                                cAGCGCGTTAATTTGCCGCGCCTGCATTGGCGTCACCACTGGTTCATCAATATGAATAGCTTTCCCCTGATTAGCCAGAAAGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAAc                                                                                      >  1:423282/1‑150 (MQ=255)

                                                                cAGCGCGTTAATTTGCCGCGCCTGCATTGGCGTCACCACTGGTTCATCAATATGAATAGCTTTCCCCTGATTAGCCAGAAAGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAAc                                                                                      >  2:4299378/1‑150 (MQ=255)

                                                                 cgcgcGTTAAGTTGCCGCGCCTGCAGTGGCGTCACCACTGGTTCATCCAGATGAATAGCTTTCCCCGGATTAGCCAGAAAGACACACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACg                                                                                     <  2:4045811/149‑1 (MQ=255)

                                                                    gcgTTATTTGGCCGCGCCTGCATTGGCGTCACCACTGGTTCATCAATATGAATAGCTTTCCCCTGATTAGCCAGAAAGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGtt                                                                                  <  1:4049013/150‑1 (MQ=255)

                                                                       ttAATTGGCCGCGCCTGCTTTGGCGTCACCACTGGTTCATCAATATTAATCGCTTTCCCCTGATTAGCCAGAAACACACACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTc                                                                               <  1:164425/150‑1 (MQ=255)

                                                                         aaTTTGCCGCGCCTGCATTGGCGTCACCACTGGTTCATCAATATGAATAGCTTTCCCCTGTTTAGCCAGAAAGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCAc                                                                             <  1:3356883/150‑1 (MQ=255)

                                                                         aaTTTGCCGCGCCTGCATTGGCGTCACCACTGGTTCATCAATATGAATAGCTTTCCCCTGATTAGCCAGAAAGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCAc                                                                             >  1:4170828/1‑150 (MQ=255)

                                                                            ttGCCGCGCCTGCATTGGCGTCACCACTGGTTCATCAATATCAATAGCTTTCCCCTGATTAGCCATAAAGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTc                                                                          <  2:3513392/150‑1 (MQ=255)

                                                                            ttGCCGCGCCTGCATGGGCGTCACCACTGGTTCATCAATATCAATAGCTTTCCCCTCATTAGCCAGAAAGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTc                                                                          <  2:3293739/150‑1 (MQ=255)

                                                                               ccGCGCCTGCATTGGCGTCACCACTGGTTCATCAATATGAATAGCTTTCCCCTGATTAGCCAGAAAGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCatg                                                                       >  1:5157471/1‑150 (MQ=255)

                                                                                  cgcCTGCATGGGCGTCACCACTGGTTCATCAATAGGAATAGCTTTCCCCTGATTAGCCAGAAAGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCatgatg                                                                    <  1:279507/150‑1 (MQ=255)

                                                                                   gcCTGCATTGGCTTCACCACTGTTTCATCAATATTAATCGCTTTCCCCTGATTACCCAGAAAGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATgg                                                                   <  1:1929108/150‑1 (MQ=255)

                                                                                   gcCTGCATTGGCGTCACCACTGGTTCATCAATATGAATAGCTTTCCCCTGATTAGCCAGAAAGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATgg                                                                   >  2:3363678/1‑150 (MQ=255)

                                                                                   gcCGGCATTGGCGTCACCACTGTTTCATCAATATGAATCGCGTTCCCCTTATTAGCCCTAAAGACACACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATgg                                                                   <  1:1928900/150‑1 (MQ=255)

                                                                                    ccTGCATTGGCGTCACCACTGGTTCATAAATATGAATAGCTTTCCCCTGATTAGCCAGAAAGACAATCGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTTCAGAACGGTTGTCACGTCATGATGGa                                                                  >  1:4199601/1‑150 (MQ=255)

                                                                                       gCATTGGCGTCACCACTGGTTCATCAATATGAATAGCTTTCCCCTGATTAGCCAGAAAGACACACGTCCGGGCCAGCTCTACACAGCTCATCTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAAt                                                               <  1:2012879/150‑1 (MQ=255)

                                                                                        cATTGGCGTCACCACTGGTTCATCAATATGAATAGCTTTCCCCTGATTAGCCAGAAAGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAAtt                                                              >  1:3966122/1‑150 (MQ=255)

                                                                                         aTTGGCGTCACCACTGGTTCATCAATATGAATCGCTTTCCCCTGATTAGCCAGAAAGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTg                                                             <  1:2590573/150‑1 (MQ=255)

                                                                                              cGTCACCACTGGTTCATCAATATGAATAGCTTTCCCCTGATTAGCCAGAAAGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCaaa                                                        >  2:1386228/1‑150 (MQ=255)

                                                                                               gTCACCACTGGTTCATCAATATTATTCGCTTTCCCGTGATTAGCCATAAACACACACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAAc                                                       <  1:3478996/150‑1 (MQ=255)

                                                                                                 caccacTGGTTCCTCAATCTGAATAGCTTTCCCCTGATTAGCCAGAAAGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGa                                                     <  2:1473856/150‑1 (MQ=255)

                                                                                                   ccacTGGTTCATCAATATGAATAGCTTTCCCCTGATTAGCCAGAAAGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACt                                                   >  1:4166283/1‑150 (MQ=255)

                                                                                                   ccacTGGTTCATCAATATGAATAGCTTTCCCCTGATTAGCCAGAAAGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACt                                                   >  1:3532908/1‑150 (MQ=255)

                                                                                                      cTGTTTCATCACTATGAATAGCTTTCCCCTGATTAGCCAGAAAGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTtca                                                <  2:3715278/150‑1 (MQ=255)

                                                                                                      cTGGTTCATCAATATGAATAGCTTTCCCCTGATTAGCCAGAAAGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTtca                                                <  2:4743846/150‑1 (MQ=255)

                                                                                                       tGGTTCATCAATATGAATAGCTTTCCCCTGATTAGCCAGAAAGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTtcat                                               <  2:3214786/150‑1 (MQ=255)

                                                                                                        ggTTCATCAATATGAATCGCTTTCCCCTGATTAGCCATAAAGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTtcatc                                              <  2:563706/150‑1 (MQ=255)

                                                                                                         gTTCATCAATATGAATAGCTTTCCCCTGATTAGCCAGAAAGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTtcatca                                             <  2:2688412/150‑1 (MQ=255)

                                                                                                         gTTCATCAATATGAATAGCTTTCCCCTGATTAGCCAGAAAGACAAACGTCCGGGCCAGCTCTACACAGCTCATATTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGGTGTCACGTAATGATGTAAATTGCCAAACGACTtaatca                                             >  1:1469460/1‑150 (MQ=255)

                                                                                                           tcatcaATATGAATAGCTTTCCCCTGATTAGCCATAAAGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGc                                           <  2:2454812/150‑1 (MQ=255)

                                                                                                           tcatcaATATGAATAGCTTTCCCCTGATTAGCCAGAAAGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGc                                           >  2:3003576/1‑150 (MQ=255)

                                                                                                            catcaATATGAATAGCTTTCCCCTGATTAGCCAGAAAGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCCGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGcc                                          >  1:3502475/1‑150 (MQ=255)

                                                                                                             atcaATATGAATAGCTTTCCCCTGATTAGCCAGAAAGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGCCa                                         >  1:2004331/1‑150 (MQ=255)

                                                                                                                   atGAATAGCTTTCCCCTGATTAGCCAGAAAGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCCTCAGCCCGGCGAt                                   >  2:1592337/1‑150 (MQ=255)

                                                                                                                   atGAATAGCTTTCCCCTGATTAGCCAGAAAGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGCCAGGCGAt                                   <  2:4322401/150‑1 (MQ=255)

                                                                                                                    ttaaTAGCTTCCCCATTATTACCAACAATCACAAACGTCCGGGCCAGCTCTAAAAAGCTCATTTTCATAGCGAAGTATTGAAAGTAGTTTTGCAGAACGTTTTTCACGTCATGTTGTAATTTGCAAAACGACTTCATCAGCCAGGCGATa                                  <  1:1573753/148‑1 (MQ=255)

                                                                                                                    tGATTAGCTTTCCCCTGATTAGCCAGAAAGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGCCAGGCGATa                                  <  2:2616417/150‑1 (MQ=255)

                                                                                                                    tGAATAGCTTTCCCCTGATTAGCCAGAAAGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGCCAGGCGATa                                  <  2:1129228/150‑1 (MQ=255)

                                                                                                                     gactAGCTTTCCCCAGAATAGCAATAAAGACAAAAGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGAGCAGTAATGAAAGTAGATTTGCAGAAAGGTTGTCACGTCATGATGGACATTGCCAAACGACTTCATCAGCCAGGCGATAg                                 >  2:3388334/4‑150 (MQ=255)

                                                                                                                             ttCCCCTGATTAGCCAGAAAGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGCCAGGCGATAGCCGCAttt                         >  2:132394/1‑150 (MQ=255)

                                                                                                                              tcccaagATTAGCCAGAATGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGAAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTc                        <  2:3600965/144‑1 (MQ=255)

                                                                                                                              tCCCCTGATTAGCCAGAAAGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTc                        <  2:1172869/150‑1 (MQ=255)

                                                                                                                              cccctttaTTCGCCATAAAGACACACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTc                        <  1:4882733/143‑1 (MQ=255)

                                                                                                                              cccccTGATTAGCCAGAAAGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTc                        <  2:1271186/149‑1 (MQ=255)

                                                                                                                               ccccTGATTAGCCAGAAAGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTcg                       >  2:1347342/1‑150 (MQ=255)

                                                                                                                               ccccTGATTAGCCAGAAAGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTcg                       >  1:1960245/1‑150 (MQ=255)

                                                                                                                                        aGCCAGAAAGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATTAAAGTAGTTTTGCAGAACGGTTGTCACGTCAAGATGGAAATTGCCAAACGACTTCATCAGCCAGGCGATATCCGCAATTCGCGCGGAATg              >  1:192083/1‑150 (MQ=255)

                                                                                                                                           cAGAAAGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTCGCGCGGAATGTTc           >  1:4429708/1‑150 (MQ=255)

                                                                                                                                            aGAAAGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTCGCGCGGAATGTTCa          >  1:3394389/1‑150 (MQ=255)

                                                                                                                                             gAAAGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTCGCGCGGAATGTTCaa         <  2:715792/150‑1 (MQ=255)

                                                                                                                                               aaGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTCGCGCGGAATGTTCatat       <  1:1732623/150‑1 (MQ=255)

                                                                                                                                               aaGAAAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTCGCGCGGAATGTTCAAAt       <  1:1260032/150‑1 (MQ=255)

                                                                                                                                                aGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTCGCGCGGAATGTTCAAAtt      >  1:3741267/1‑150 (MQ=255)

                                                                                                                                                  aCAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATAGCCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTCGCGCGGAATGTTCAAATTCg    >  2:2129752/1‑150 (MQ=255)

TCGCCGGTAGCCCCACCCGCCAGGCAAACTCCCCCGCGTTCTGGTACATACCACTGGTCGCCATCAGCGCGTTAATTTGCCGCGCCTGCATTGGCGTCACCACTGGTTCATCAATATGAATAGCTTTCCCCTGATTAGCCAGAAAGACAAACGTCCGGGCCAGCTCTACACAGCTCATTTTCAGAGCGCAGTAATGAAAGTAGTTTTGCAGAACGGTTGTCACGTCATGATGGAAATTGCCAAACGACTTCATCAGCCAGGCGATAGCCGCATTTCGCGCGGAATGTTCAAATTCGGA  >  CP009273/1606734‑1607031

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