BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	3,512,321	A→T	13.0%	S35T (TCC→ACC)	aroK ←	shikimate kinase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	3,512,321	0	A	T	13.0%	866.5 / 92.8	300	S35T (TCC→ACC)	aroK	shikimate kinase I
Reads supporting (aligned to +/- strand): ref base A (120/141); new base T (19/20); total (139/161)
Fisher's exact test for biased strand distribution p-value = 8.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 306 total aligned reads displayed.

AGAGCCGCCGCCAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCA‑TTTTTTCGGTACTACTAAGACTATTCGTTAATGATAAACCCGCTTCGC  >  CP009273/3512172‑3512470

  agCCGCCGCCAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAAt                                                                                                                                                      <  1:2152735/150‑1 (MQ=255)

   gccgccgccAGTAGCCAGCACAAGACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATc                                                                                                                                                     <  2:4507940/150‑1 (MQ=255)

      gccgccAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTa                                                                                                                                                  <  2:224484/150‑1 (MQ=255)

          ccAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAAt                                                                                                                                              >  1:335870/1‑150 (MQ=255)

           cAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAAtt                                                                                                                                             <  1:657532/150‑1 (MQ=255)

           cAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAAtt                                                                                                                                             <  1:3243209/150‑1 (MQ=255)

                 ccAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCatat                                                                                                                                       <  2:841442/150‑1 (MQ=255)

                 ccAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCatat                                                                                                                                       <  1:4322412/150‑1 (MQ=255)

                  cAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATAtt                                                                                                                                      >  2:4277331/1‑150 (MQ=255)

                   aGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTg                                                                                                                                     <  1:1506892/150‑1 (MQ=255)

                   aGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCACCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTg                                                                                                                                     <  1:3320414/150‑1 (MQ=255)

                     cacaATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAg                                                                                                                                   <  1:2367613/150‑1 (MQ=255)

                     cacaATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAg                                                                                                                                   <  2:1500451/150‑1 (MQ=255)

                      acaATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGTATCGTAAAATTCCATATTGAgt                                                                                                                                  <  1:3681650/150‑1 (MQ=255)

                      acaATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAgt                                                                                                                                  <  1:3751507/150‑1 (MQ=255)

                        aaTACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAgttg                                                                                                                                >  2:4635944/1‑150 (MQ=255)

                          tACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAgttgtt                                                                                                                              <  2:2982540/150‑1 (MQ=255)

                          tACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAgttgtt                                                                                                                              <  1:2095989/150‑1 (MQ=255)

                           aCCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCAATATTGAGttgttg                                                                                                                             >  2:258633/1‑150 (MQ=255)

                           aCCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCCTTCGTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGTATCGTAAAATTCCATATTGAGttgttg                                                                                                                             <  2:2643943/150‑1 (MQ=255)

                            cccTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGa                                                                                                                            >  2:487292/1‑150 (MQ=255)

                             ccTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGTATCGTAAAATTCCATATTGAGTTGTTGAg                                                                                                                           >  1:1047016/1‑150 (MQ=255)

                             ccTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAg                                                                                                                           <  2:3784554/150‑1 (MQ=255)

                               tGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCt                                                                                                                         >  2:4715016/1‑150 (MQ=255)

                                gTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTa                                                                                                                        <  2:2796684/150‑1 (MQ=255)

                                 tttCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTATTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTaa                                                                                                                       >  2:3842926/1‑150 (MQ=255)

                                   tctcGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACt                                                                                                                     >  2:4224889/1‑150 (MQ=255)

                                    ctcGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTg                                                                                                                    >  2:1807915/1‑150 (MQ=255)

                                    ctcGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTg                                                                                                                    <  1:1793869/150‑1 (MQ=255)

                                     tcGGTCAACTCATTGATGACCTTTTCTTCGCGAGCGCGGAAGCCTTCGTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTgg                                                                                                                   <  1:949982/150‑1 (MQ=255)

                                      cGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGTATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGgc                                                                                                                  <  2:4806565/150‑1 (MQ=255)

                                      cGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGgc                                                                                                                  >  2:1721076/1‑150 (MQ=255)

                                      cGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGgc                                                                                                                  <  2:4327572/150‑1 (MQ=255)

                                          cAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGccc                                                                                                              <  1:2889178/150‑1 (MQ=255)

                                          cAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGccc                                                                                                              <  1:113484/150‑1 (MQ=255)

                                             cTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAAt                                                                                                           <  1:2480984/150‑1 (MQ=255)

                                                    atgaCCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCtt                                                                                                    <  1:473574/150‑1 (MQ=255)

                                                      gaCCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCtttt                                                                                                  <  2:1285548/150‑1 (MQ=255)

                                                       aCCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCttttt                                                                                                 <  2:2913184/150‑1 (MQ=255)

                                                       aCCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCG‑TTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTc                                                                                                >  1:4687946/1‑150 (MQ=255)

                                                        ccTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTc                                                                                                >  1:3961381/1‑150 (MQ=255)

                                                         cTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTcc                                                                                               <  2:3526150/150‑1 (MQ=255)

                                                         cTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTcc                                                                                               <  2:1635481/150‑1 (MQ=255)

                                                           tttCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCgg                                                                                             <  1:1115652/150‑1 (MQ=255)

                                                           tttCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCgg                                                                                             <  1:536000/150‑1 (MQ=255)

                                                              cttcGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGTATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCAc                                                                                          >  2:3373141/1‑150 (MQ=255)

                                                               ttcGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCAcc                                                                                         >  1:611040/1‑150 (MQ=255)

                                                                  gcgATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCAt                                                                                      <  1:487292/150‑1 (MQ=255)

                                                                   cgATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATa                                                                                     <  2:3206607/150‑1 (MQ=255)

                                                                   cgATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATa                                                                                     >  1:589472/1‑150 (MQ=255)

                                                                     aTCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCACTAGTGCTTTTTCCGGCACCCATAgg                                                                                   >  1:3071894/1‑150 (MQ=255)

                                                                       cgcgGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGTATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGcc                                                                                 <  2:3602412/150‑1 (MQ=255)

                                                                       cgcgGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGcc                                                                                 >  1:811163/1‑150 (MQ=255)

                                                                       cgcgGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGcc                                                                                 <  1:2024975/150‑1 (MQ=255)

                                                                       cgcgGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGcc                                                                                 >  2:1024679/1‑150 (MQ=255)

                                                                            aaGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAAcc                                                                            >  2:2106241/1‑150 (MQ=255)

                                                                            aaGCCTTCTTCGCCTTCTAAATCGAAAAACCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTTATCGGAATCGTAAAATTCAATATTGAGTTTTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCTGCACCCATAGGCCCAAcc                                                                            >  1:4911565/1‑150 (MQ=255)

                                                                             aGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCa                                                                           <  1:270009/150‑1 (MQ=255)

                                                                              gCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAg                                                                          <  1:2271803/150‑1 (MQ=255)

                                                                               ccTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGa                                                                         >  1:1888765/1‑150 (MQ=255)

                                                                                cttcttCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGaa                                                                        >  2:1971383/1‑150 (MQ=255)

                                                                                cttcttCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGaa                                                                        >  1:408074/1‑150 (MQ=255)

                                                                                cttcttCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGaa                                                                        >  2:2382261/1‑150 (MQ=255)

                                                                                  tcttcGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAg                                                                      <  1:13086/150‑1 (MQ=255)

                                                                                   cttcGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGa                                                                     >  1:1406513/1‑150 (MQ=255)

                                                                                    ttcGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGTATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGat                                                                    >  2:3582007/1‑150 (MQ=255)

                                                                                    ttcGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGat                                                                    >  2:2267501/1‑150 (MQ=255)

                                                                                     tcGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGata                                                                   <  1:4982891/150‑1 (MQ=255)

                                                                                      cGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGatat                                                                  <  1:1559858/150‑1 (MQ=255)

                                                                                       gCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATAtt                                                                 <  2:532303/150‑1 (MQ=255)

                                                                                           tCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGt                                                             >  2:4284126/1‑150 (MQ=255)

                                                                                            cTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGAGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGtt                                                            >  2:2302256/1‑150 (MQ=255)

                                                                                                 tCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTctct                                                       <  1:387587/150‑1 (MQ=255)

                                                                                                   gAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGc                                                     >  1:2398403/1‑150 (MQ=255)

                                                                                                    aaaaCCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGcc                                                    <  2:1732660/150‑1 (MQ=255)

                                                                                                      aaCCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGTATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCA‑t                                                 <  2:1413847/150‑1 (MQ=255)

                                                                                                       aCCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGTATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCA‑tt                                                <  1:2936689/150‑1 (MQ=255)

                                                                                                        cccAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCA‑ttt                                               >  2:3145441/1‑150 (MQ=255)

                                                                                                          cAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCA‑ttttt                                             <  1:4169200/150‑1 (MQ=255)

                                                                                                          cAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCA‑ttttt                                             >  1:4206801/1‑150 (MQ=255)

                                                                                                           aGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCA‑tttttt                                            >  2:1440970/1‑150 (MQ=255)

                                                                                                               cacaTCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCA‑TTTTTTCGGt                                        <  1:2993723/150‑1 (MQ=255)

                                                                                                                acaTCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCA‑TTTTTTCGGta                                       >  1:4635558/1‑150 (MQ=255)

                                                                                                                 caTCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCA‑TTTTTTCGGtac                                      <  1:1972785/150‑1 (MQ=255)

                                                                                                                  aTCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCA‑TTTTTTCGGtact                                     <  1:4489600/150‑1 (MQ=255)

                                                                                                                  aTCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGGTTCTCTGCCA‑TTTTTTCGGtact                                     >  1:2602095/1‑150 (MQ=255)

                                                                                                                   tCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGTATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCA‑TTTTTTCGGtacta                                    >  2:4107494/1‑150 (MQ=255)

                                                                                                                   tCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGTATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCA‑TTTTTTCGGtacta                                    >  2:4161898/1‑150 (MQ=255)

                                                                                                                     aGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCA‑TTTTTTCGGTactact                                  >  2:511418/1‑150 (MQ=255)

                                                                                                                      gCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTTAGCTAACTGGCGCCCAATAGTTCTTTTTCCGGCACCCATAGGCCCAAACAGAAAGATATTGCGTTTCTCTGCCA‑TTTTTTAGGTActacta                                 >  1:3503749/1‑150 (MQ=255)

                                                                                                                       cTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCA‑TTTTTTCGGTACTACTaa                                <  1:4715010/150‑1 (MQ=255)

                                                                                                                              tCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCA‑TTTTTTCGGTACTACTAAGACTAtt                         >  2:3722633/1‑150 (MQ=255)

                                                                                                                               cGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCATCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCA‑TTTTTTCGGTACTACTAAGACTATTc                        <  2:3985084/150‑1 (MQ=255)

                                                                                                                                gTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTTCGGTACTACTAAGACTATTc                        <  2:2041284/150‑1 (MQ=255)

                                                                                                                                         tctcTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCA‑TTTTTTCGGTACTAATAAGACTATTCGTTAATGATa              >  2:8088/1‑150 (MQ=255)

                                                                                                                                           tctTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGTGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCA‑TTTTTTCGGTACTACTAAGACTATTCGTTAATGATaaa            <  2:2487621/150‑1 (MQ=255)

                                                                                                                                            ctTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCA‑TTTTTTCGGTACTACTAAGACTATTAGTTAATGATAAAc           >  2:1165882/1‑150 (MQ=255)

                                                                                                                                              tGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCA‑TTTTTTCGGTACTACTAAGACTATTCGTTAATGATAAAccc         <  2:687284/150‑1 (MQ=255)

                                                                                                                                              tGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCA‑TTTTTTCGGTACTACTAAGACTATTCGTTAATGATAAAccc         >  2:5046088/1‑150 (MQ=255)

                                                                                                                                                   ggAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCA‑TTTTTTCGGTACTACTAAGACTATTCGTTAATGATAAACCCGCTTc    <  2:875387/150‑1 (MQ=255)

AGAGCCGCCGCCAGTAGCCAGCACAATACCCTGTTTCTCGGTCAACTCATTGATGACCTTTTCTTCGCGATCGCGGAAGCCTTCTTCGCCTTCTAAATCGAAAACCCAGCCCACATCAGCTCCGGTTCGTTTCTCAATCTCTTGATCGGAATCGTAAAATTCCATATTGAGTTGTTGAGCTAACTGGCGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCA‑TTTTTTCGGTACTACTAAGACTATTCGTTAATGATAAACCCGCTTCGC  >  CP009273/3512172‑3512470

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