BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	1,932,296	G→T	58.1%	D131Y (GAT→TAT)	pykA →	pyruvate kinase II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	1,932,296	0	G	T	58.1%	45.1 / inf	251	D131Y (GAT→TAT)	pykA	pyruvate kinase II
Reads supporting (aligned to +/- strand): ref base G (58/47); new base T (77/69); total (135/116)
Fisher's exact test for biased strand distribution p-value = 7.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 255 total aligned reads displayed.

TTAAAGAAGGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGTTTGTCGGCTGAAGCGCTGACCGAAAAAGACAAAGCA  >  CP009273/1932147‑1932445

ttAAAGAAGGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACg                                                                                                                                                       >  2:3745265/1‑150 (MQ=255)

  aaaGAAGGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTAt                                                                                                                                                     <  1:1164572/150‑1 (MQ=255)

   aagaagGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATg                                                                                                                                                    >  2:1969173/1‑150 (MQ=255)

     gaagGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGt                                                                                                                                                  >  2:791346/1‑150 (MQ=255)

     gaagGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGt                                                                                                                                                  >  1:535450/1‑150 (MQ=255)

      aagGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTc                                                                                                                                                 <  1:2890016/150‑1 (MQ=255)

        ggCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTcgc                                                                                                                                               >  2:959961/1‑150 (MQ=255)

           aaaGTTTTCCTCACTATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTc                                                                                                                                            <  1:829667/150‑1 (MQ=255)

            aaGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTcc                                                                                                                                           >  2:1865665/1‑150 (MQ=255)

            aaGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTcc                                                                                                                                           >  1:4913783/1‑150 (MQ=255)

               ttttCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGt                                                                                                                                        >  2:3250423/1‑150 (MQ=255)

               ttttCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGt                                                                                                                                        <  2:3395770/150‑1 (MQ=255)

               ttttCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGt                                                                                                                                        <  2:1799367/150‑1 (MQ=255)

                tttCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGtt                                                                                                                                       >  2:2668352/1‑150 (MQ=255)

                 ttCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTa                                                                                                                                      <  2:2780168/150‑1 (MQ=255)

                 ttCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTa                                                                                                                                      <  2:507515/150‑1 (MQ=255)

                       aaTATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTa                                                                                                                                >  2:4881131/1‑150 (MQ=255)

                        atatTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTAc                                                                                                                               >  1:1972933/1‑150 (MQ=255)

                        atatTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTAc                                                                                                                               >  1:1102515/1‑150 (MQ=255)

                         tatTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACt                                                                                                                              <  1:2481272/150‑1 (MQ=255)

                          atTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTg                                                                                                                             >  1:2394311/1‑150 (MQ=255)

                           ttGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTgg                                                                                                                            >  2:1546913/1‑150 (MQ=255)

                           ttGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTgg                                                                                                                            >  2:4563244/1‑150 (MQ=255)

                            tGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGa                                                                                                                           >  1:3082026/1‑150 (MQ=255)

                              gggATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAg                                                                                                                         >  2:4239147/1‑150 (MQ=255)

                              gggATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAg                                                                                                                         >  2:84388/1‑150 (MQ=255)

                                gATAAATTACTGCTCGACGCCAACCTGGTTAAAGGAGAAGGAGACAAAGAAAAAGTAGGTATCGACTACAAAGGCATGCCTGCTGACTTCGTGCCTGGTGAAATCCTGATGCTGTACTATGGTAGCGTCCAGTTAAAAGTACTGGAAGtt                                                                                                                       >  1:2752981/1‑150 (MQ=255)

                                 aTAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTc                                                                                                                      >  2:3175678/1‑150 (MQ=255)

                                 aTAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTc                                                                                                                      <  2:1266189/150‑1 (MQ=255)

                                   aaaTTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAg                                                                                                                    <  1:5155581/150‑1 (MQ=255)

                                   aaaTTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAg                                                                                                                    >  1:4951824/1‑150 (MQ=255)

                                     aTTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAggg                                                                                                                  >  1:447366/1‑150 (MQ=255)

                                         cTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATg                                                                                                              <  1:1174136/150‑1 (MQ=255)

                                         cTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATg                                                                                                              >  2:2144165/1‑150 (MQ=255)

                                           gCTCGACGCCAACCTGGGTAAAGGTGTAGGCGACAATGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCTTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGaa                                                                                                            <  1:3605527/150‑1 (MQ=255)

                                               gACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAgtg                                                                                                        >  2:1209040/1‑150 (MQ=255)

                                               gACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAgtg                                                                                                        >  1:1597137/1‑150 (MQ=255)

                                                 cGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGtt                                                                                                      <  2:4598275/150‑1 (MQ=255)

                                                 cGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGtt                                                                                                      <  1:644295/150‑1 (MQ=255)

                                                   ccAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCa                                                                                                    <  1:3497732/150‑1 (MQ=255)

                                                       ccTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGATGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGTGCATGAAAGAGTTCACCGa                                                                                                >  1:894410/1‑150 (MQ=255)

                                                        cTGGGTAATGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGaa                                                                                               >  1:467583/1‑150 (MQ=255)

                                                         tGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAg                                                                                              <  2:3710275/150‑1 (MQ=255)

                                                           ggTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTc                                                                                            <  1:3994908/150‑1 (MQ=255)

                                                           ggTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTc                                                                                            >  2:1380836/1‑150 (MQ=255)

                                                            gTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCa                                                                                           >  2:2782821/1‑150 (MQ=255)

                                                              aaaGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCAcc                                                                                         <  1:1865665/150‑1 (MQ=255)

                                                               aaGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAATTCACCg                                                                                        >  2:1400894/1‑150 (MQ=255)

                                                                   tGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCgg                                                                                    >  2:3907304/1‑150 (MQ=255)

                                                                     aaGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCggtg                                                                                  >  2:2278793/1‑150 (MQ=255)

                                                                     aaGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCggtg                                                                                  <  2:4628144/150‑1 (MQ=255)

                                                                      aGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCggtgg                                                                                 >  2:576270/1‑150 (MQ=255)

                                                                          gACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTccc                                                                             >  1:3941507/1‑150 (MQ=255)

                                                                            cAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCct                                                                           >  1:923223/1‑150 (MQ=255)

                                                                              aaGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCAAGTTAAAAGTACTGTAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCctct                                                                         <  2:2951552/150‑1 (MQ=255)

                                                                               aGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCtctc                                                                        >  2:3858295/1‑150 (MQ=255)

                                                                                 aaaaaGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCAGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCca                                                                      >  1:4290314/1‑150 (MQ=255)

                                                                                  aaaaGTCGGTATCGACTACAAAGTCCTGCATGCTGACGTCGTGCCTGGTGACATCCTGCTGCTTGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCcat                                                                     <  2:4421173/150‑2 (MQ=255)

                                                                                  aaaaGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCcaa                                                                     >  2:5238574/1‑150 (MQ=255)

                                                                                    aaGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCcaaca                                                                   <  1:3201120/150‑1 (MQ=255)

                                                                                      gTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAAt                                                                 >  2:4070310/1‑150 (MQ=255)

                                                                                         ggTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATaaa                                                              <  2:1129196/150‑1 (MQ=255)

                                                                                          gTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACCATAAAg                                                             >  2:4164599/1‑150 (MQ=255)

                                                                                           tATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAgg                                                            >  1:2446672/1‑150 (MQ=255)

                                                                                           tATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAgg                                                            <  2:4760362/150‑1 (MQ=255)

                                                                                           tATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAgg                                                            <  1:3317282/150‑1 (MQ=255)

                                                                                             tCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTa                                                          <  1:1691737/150‑1 (MQ=255)

                                                                                               gACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATc                                                        >  1:1318992/1‑150 (MQ=255)

                                                                                                aCTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATca                                                       <  1:3070995/150‑1 (MQ=255)

                                                                                                aCTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATca                                                       <  2:1082278/150‑1 (MQ=255)

                                                                                                aCTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATca                                                       <  1:29363/150‑1 (MQ=255)

                                                                                                 cTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATcaa                                                      >  1:1739988/1‑150 (MQ=255)

                                                                                                   aCAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATcaaca                                                    >  2:1474191/1‑150 (MQ=255)

                                                                                                     aaaGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACaaa                                                  >  1:2630649/1‑150 (MQ=255)

                                                                                                      aaGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAAc                                                 <  2:152269/150‑1 (MQ=255)

                                                                                                       aGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACt                                                <  1:3797301/150‑1 (MQ=255)

                                                                                                        ggCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACtt                                               >  2:3556010/1‑150 (MQ=255)

                                                                                                        ggCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACtt                                               <  1:4153065/150‑1 (MQ=255)

                                                                                                           ctgcctgcTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTggc                                            <  1:2356331/150‑1 (MQ=255)

                                                                                                               ctgctgACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAATGTATCAACAAACATGgcggcg                                        >  1:1936064/1‑150 (MQ=255)

                                                                                                               ctgctgACGTCGTGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGgcggcg                                        >  1:625982/1‑150 (MQ=255)

                                                                                                                       gtcgtGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGTTTGTCg                                <  1:4378868/150‑1 (MQ=255)

                                                                                                                        tcgtGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGTTTGTCgg                               <  2:2495314/150‑1 (MQ=255)

                                                                                                                         cgtGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGTTTGTCGGc                              >  1:2552593/1‑150 (MQ=255)

                                                                                                                           tGCCTGGTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGTTTGTCGGCTg                            <  2:523273/150‑1 (MQ=255)

                                                                                                                             ccTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGTTTGTCGGCTGaa                          >  2:4557023/1‑150 (MQ=255)

                                                                                                                                tgTGACATCCTGATGCTGGACTATGGTCGCTTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGTAGTCACCGTCGGTGGTCCCCTCTCCAACAATCAAGGTATCAACAAACTTGGCGGCGGTTTGACGGCTGAAgcg                       >  2:4794980/2‑150 (MQ=255)

                                                                                                                                 gTGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGTTTGTCGGCTGAAgcgc                      >  1:4989138/1‑150 (MQ=255)

                                                                                                                                  tGACATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGTTTGTCGGCTGAAGCGCt                     >  2:2974326/1‑150 (MQ=255)

                                                                                                                                     cATCCTGCTGCTGGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCACTCCAACAAAAAAGGTATCAACAAACATGGCGGCGGTTTTTCGGCTGAAGCGCTGAc                  >  1:2309450/1‑150 (MQ=255)

                                                                                                                                         ctgctgctgGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGTTTGTCGGCTGAAGCGCTGACCGaa              >  1:2124724/1‑150 (MQ=255)

                                                                                                                                          tgctgctgGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGTTTGTCGGCTGAAGCGCTGACCGaaa             >  1:1053345/1‑150 (MQ=255)

                                                                                                                                          tgctgctgGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGTTTGTCGGCTGAAGCGCTGACCGaaa             >  1:464329/1‑150 (MQ=255)

                                                                                                                                               ctgGACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGTTTGTCGGCTGAAGCGCTGACCGAAAAAGAc        <  2:447443/150‑1 (MQ=255)

                                                                                                                                               ctgGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGTTTGTCGGCTGAAGCGCTGACCGAAAAAGAc        <  1:643173/150‑1 (MQ=255)

                                                                                                                                               ctgGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGTTTGTCGGCTGAAGCGCTGACCGAAAAAGAc        <  1:2668352/150‑1 (MQ=255)

                                                                                                                                                tgGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGTTTGTCGGCTGAAGCGCTGACCGAAAAAGACa       <  1:2548574/150‑1 (MQ=255)

                                                                                                                                                  gACTATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGTTTGTCGGCTGAAGCGCTGACCGAAAAAGACaaa     <  2:63117/150‑1 (MQ=255)

                                                                                                                                                     gATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGTTTGTCGGCTGAAGCGCTGACCGAAAAAGACAAAGCa  <  2:5097979/150‑1 (MQ=255)

                                                                                                                                                     gATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGTTTGTCGGCTGAAGCGCTGACCGAAAAAGACAAAGCa  >  2:3452218/1‑150 (MQ=255)

TTAAAGAAGGCAAAGTTTTCCTCAATATTGGGGATAAATTCCTGCTCGACGCCAACCTGGGTAAAGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGCGGTTTGTCGGCTGAAGCGCTGACCGAAAAAGACAAAGCA  >  CP009273/1932147‑1932445

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