BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	3,969,710	G→T	94.5%	G25C (GGT→TGT)	wzxE →	O‑antigen translocase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	3,969,710	0	G	T	94.5%	920.2 / 41.0	293	G25C (GGT→TGT)	wzxE	O‑antigen translocase
Reads supporting (aligned to +/- strand): ref base G (8/8); new base T (141/135); total (150/143)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01

Only 100 of 295 total aligned reads displayed.

TGCCGCTGTTCTACAACCTGTCGCCCGTCAATCAGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAA  >  CP009273/3969563‑3969859

tGCCGCTGTTCTACAACCTGTCGCCCGTCAATCAGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGt                                                                                                                                                     <  1:253454/150‑1 (MQ=255)

 gCCGCTGTTCTACAACCTGTCGCCCGTCAATCAGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTa                                                                                                                                                    <  1:3208034/150‑1 (MQ=255)

  ccGCTGTTCTACAACCTGTCGCCCGTCAATCAGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTaa                                                                                                                                                   >  2:709980/1‑150 (MQ=255)

     cTGTTCTACAACCTGTCGCCCGTCAATCAGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGtt                                                                                                                                                <  2:3387512/150‑1 (MQ=255)

     cTGTTCTACAACCTGTCGCCCGTCAATCAGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGtt                                                                                                                                                <  1:4004229/150‑1 (MQ=255)

      tGTTCTACAACCTGTCGCCCGTCAATCAGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTg                                                                                                                                               <  2:2072194/150‑1 (MQ=255)

        ttCTACAACCTGTCGCCCGTCAATCAGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCt                                                                                                                                             <  2:2899890/150‑1 (MQ=255)

           tACAACCTGTCGCCCGTCAATCAGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGc                                                                                                                                          <  1:3612382/150‑1 (MQ=255)

           tACAACCTGTCGCCCGTCAATCAGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGc                                                                                                                                          >  2:2973904/1‑150 (MQ=255)

           tACAACCTGTCGCCCGTCAATCAGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGc                                                                                                                                          >  2:3200493/1‑150 (MQ=255)

            aCAACCTGTCGCCCGTCAATCAGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCg                                                                                                                                         >  1:3376157/1‑150 (MQ=255)

              aaCCTGTCGCCCGTCAATCAGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGgt                                                                                                                                       <  1:2078683/150‑1 (MQ=255)

                   gTCGCCCGTCAATCAGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCAt                                                                                                                                  <  2:2687302/150‑1 (MQ=255)

                   gTCGCCCGTCAATCAGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCAt                                                                                                                                  <  1:760049/150‑1 (MQ=255)

                    tCGCCCGTCAATCAGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCAtt                                                                                                                                 <  2:1818964/150‑1 (MQ=255)

                       cccGTCAATCAGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTgg                                                                                                                              >  2:1058155/1‑150 (MQ=255)

                       cccGTCAATCAGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTACTGGCGGTGTCATTTgg                                                                                                                              >  2:2483956/1‑150 (MQ=255)

                        ccGTCAATCAGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGt                                                                                                                             >  1:4197812/1‑150 (MQ=255)

                        ccGTCAATCAGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTAAACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGAGGTGTCATTTGGt                                                                                                                             >  1:878768/1‑150 (MQ=255)

                        ccGTCAATCAGCGTACGGTAATTGAGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGTTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGt                                                                                                                             <  1:2484694/150‑1 (MQ=255)

                          gTCAATCAGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTcc                                                                                                                           >  2:3980620/1‑150 (MQ=255)

                           tCAATCAGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCcg                                                                                                                          <  2:2839080/150‑1 (MQ=255)

                           tCAATCAGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCcg                                                                                                                          <  1:559914/150‑1 (MQ=255)

                            cAATCAGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCcgg                                                                                                                         >  2:2229863/1‑150 (MQ=255)

                              aTCAGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCcggcg                                                                                                                       <  1:3703079/150‑1 (MQ=255)

                               tCAGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCcggcgg                                                                                                                      >  1:2504477/1‑150 (MQ=255)

                                cAGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCggg                                                                                                                     <  1:776294/150‑1 (MQ=255)

                                       cGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTggg                                                                                                              <  2:2804048/150‑1 (MQ=255)

                                       cGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTggg                                                                                                              <  2:4417183/150‑1 (MQ=255)

                                       cGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTggg                                                                                                              >  2:3353565/1‑150 (MQ=255)

                                         gTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCt                                                                                                            <  1:3358780/150‑1 (MQ=255)

                                          tAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTg                                                                                                           >  2:780795/1‑150 (MQ=255)

                                            aTTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTggc                                                                                                         >  1:2111475/1‑150 (MQ=255)

                                            aTTGCGACTTTGTTGAACTACTTATCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGATGCTGGCGGAGTCATTAGGTCCGGCGGGGCTTGGGCTggc                                                                                                         >  1:2646830/1‑150 (MQ=255)

                                             ttGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTggcg                                                                                                        >  2:1500354/1‑150 (MQ=255)

                                             ttGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTggcg                                                                                                        <  2:4297770/150‑1 (MQ=255)

                                              tGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTggcgg                                                                                                       >  1:1038954/1‑150 (MQ=255)

                                              tGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTggcgg                                                                                                       >  1:3820343/1‑150 (MQ=255)

                                              tGCGACTGTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTggcgt                                                                                                       >  2:303095/1‑149 (MQ=255)

                                                cGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCa                                                                                                     >  2:327645/1‑150 (MQ=255)

                                                      tgttgAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTc                                                                                               >  1:4047582/1‑150 (MQ=255)

                                                            actactTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAg                                                                                         <  2:3764417/150‑1 (MQ=255)

                                                             ctactTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGAGTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGt                                                                                        >  2:4455366/1‑150 (MQ=255)

                                                              tactTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGtt                                                                                       <  1:1650884/150‑1 (MQ=255)

                                                               actTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTg                                                                                      >  1:3409624/1‑150 (MQ=255)

                                                                ctTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGa                                                                                     >  2:2368025/1‑150 (MQ=255)

                                                                  tttCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGAtt                                                                                   >  1:3822442/1‑150 (MQ=255)

                                                                    tCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTAc                                                                                 <  1:4213847/150‑1 (MQ=255)

                                                                     ccTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTAcc                                                                                >  1:364285/1‑150 (MQ=255)

                                                                       tGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGt                                                                              >  1:923621/1‑150 (MQ=255)

                                                                             gTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCgg                                                                        >  2:3864948/1‑150 (MQ=255)

                                                                              tCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGc                                                                       >  2:4603588/1‑150 (MQ=255)

                                                                                 ttGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTg                                                                    >  2:440195/1‑150 (MQ=255)

                                                                                    gCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCtt                                                                 <  2:4078597/150‑1 (MQ=255)

                                                                                    gCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCtt                                                                 <  2:1400332/150‑1 (MQ=255)

                                                                                      aaaaGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGc                                                               >  1:2126452/1‑150 (MQ=255)

                                                                                           cGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCgggg                                                          <  1:791545/150‑1 (MQ=255)

                                                                                           cGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCgggg                                                          >  1:63531/1‑150 (MQ=255)

                                                                                           cGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGATTGGGCTGGCGGAAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCgggg                                                          >  2:1758299/1‑150 (MQ=255)

                                                                                            gTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGc                                                         <  1:299288/150‑1 (MQ=255)

                                                                                             tCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCt                                                        <  2:3600955/150‑1 (MQ=255)

                                                                                             tCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATCTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCt                                                        <  1:2402060/150‑1 (MQ=255)

                                                                                              ccTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTg                                                       >  2:596901/1‑150 (MQ=255)

                                                                                                 tgtgGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCa                                                    <  2:4143595/150‑1 (MQ=255)

                                                                                                 tgtgGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCa                                                    <  1:2784222/150‑1 (MQ=255)

                                                                                                  gtgGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCAt                                                   >  2:218058/1‑150 (MQ=255)

                                                                                                     gACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCtt                                                >  2:3121252/1‑150 (MQ=255)

                                                                                                     gACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCtt                                                <  2:4262368/150‑1 (MQ=255)

                                                                                                     gACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCtt                                                <  2:1805024/150‑1 (MQ=255)

                                                                                                      aCGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCttt                                               >  2:2415166/1‑150 (MQ=255)

                                                                                                       cggcggCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTa                                              >  1:4318492/1‑150 (MQ=255)

                                                                                                           ggcCAGTACACTGGTCAAGATTGGTGCCGGNTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACgg                                          >  1:4284352/1‑150 (MQ=255)

                                                                                                           ggcCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACgg                                          >  2:3475238/1‑150 (MQ=255)

                                                                                                            gcCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGgt                                         <  1:306123/150‑1 (MQ=255)

                                                                                                               aGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTa                                      <  1:96925/150‑1 (MQ=255)

                                                                                                               aGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTa                                      <  1:2193791/150‑1 (MQ=255)

                                                                                                               aGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTa                                      <  1:3132970/150‑1 (MQ=255)

                                                                                                               aGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTa                                      <  1:4486552/150‑1 (MQ=255)

                                                                                                                 tACACTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAAc                                    >  1:608119/1‑150 (MQ=255)

                                                                                                                  acacTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAAcc                                   >  1:2345971/1‑150 (MQ=255)

                                                                                                                    acTGGTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCaa                                 <  2:1264662/150‑1 (MQ=255)

                                                                                                                       ggTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATa                              <  1:2630948/150‑1 (MQ=255)

                                                                                                                       ggTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATa                              <  1:2848329/150‑1 (MQ=255)

                                                                                                                       ggTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATa                              <  1:4106310/150‑1 (MQ=255)

                                                                                                                        gTCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATAc                             <  2:1568278/150‑1 (MQ=255)

                                                                                                                         tCAAGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACg                            <  2:1115777/150‑1 (MQ=255)

                                                                                                                           aaGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGtt                          >  2:4522344/1‑150 (MQ=255)

                                                                                                                            aGATTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTg                         <  1:1480461/150‑1 (MQ=255)

                                                                                                                              aTTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGcc                       >  1:1069698/1‑150 (MQ=255)

                                                                                                                              aTTGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGcc                       >  1:1321737/1‑150 (MQ=255)

                                                                                                                              aTTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGcc                       >  1:192916/1‑150 (MQ=255)

                                                                                                                               ttGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGccc                      >  2:1719786/1‑150 (MQ=255)

                                                                                                                                tGGTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCa                     <  2:2727982/150‑1 (MQ=255)

                                                                                                                                 ggTGCCGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAg                    <  2:3545307/150‑1 (MQ=255)

                                                                                                                                     ccGGGTTACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTAcc                <  1:1176882/150‑1 (MQ=255)

                                                                                                                                          ttACTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGAt           <  1:2719454/150‑1 (MQ=255)

                                                                                                                                            aCTGGTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATaa         >  2:4205746/1‑150 (MQ=255)

                                                                                                                                               ggTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATcc      <  1:4027767/150‑1 (MQ=255)

                                                                                                                                                gTCTGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCg     >  1:1005097/1‑150 (MQ=255)

                                                                                                                                                   ggTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCaa  >  2:4224021/1‑150 (MQ=255)

TGCCGCTGTTCTACAACCTGTCGCCCGTCAATCAGCGTACGGTAATTGCGACTTTGTTGAACTACTTTTCCTGATATGTCGTTGGCAAAAGCGTCCTTGTGGACGGCGGCCAGTACACTGGTCAAGATTGGTGCCGGGTTACTGGTCGGTAAGTTGCTGGCGGTGTCATTTGGTCCGGCGGGGCTTGGGCTGGCGGCAAATTTCCGCCAGTTGATTACCGTGCTCGGCGTGCTTGCCGGGGCTGGCATCTTTAACGGTGTAACCAAATACGTTGCCCAGTACCATGATAATCCGCAA  >  CP009273/3969563‑3969859

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