BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	4,297,485	G→C	45.9%	G44G (GGC→GGG)	alsK ←	D‑allose kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	4,297,485	0	G	C	45.9%	109.7 / inf	283	G44G (GGC→GGG)	alsK	D‑allose kinase
Reads supporting (aligned to +/- strand): ref base G (87/66); new base C (67/63); total (154/129)
Fisher's exact test for biased strand distribution p-value = 4.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

Only 100 of 286 total aligned reads displayed.

CGCCGCTGTTAACGGCAGGTTAGGCGTAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGC‑TTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGCTGTTTTTGCATGAGGCTGTACCTCGCTT  >  CP009273/4297338‑4297633

  ccgcTGTTAACGGCAGGTTAGGCGTAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCgg                                                                                                                                                   <  1:4370907/150‑1 (MQ=255)

   cgcTGTTAACGGCAGGTTAGGCGTAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCggg                                                                                                                                                  >  1:3572831/1‑150 (MQ=255)

    gcTGTTAACGGCAGGTTAGGCGTAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGa                                                                                                                                                 <  2:2755768/150‑1 (MQ=255)

         tAACGGCAGGTTAGGCGTAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAAt                                                                                                                                            <  2:2611915/150‑1 (MQ=255)

         tAACGGCAGGTTAGGCGTAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAAt                                                                                                                                            >  2:1415738/1‑150 (MQ=255)

            cGGCAGGTTAGGCGTAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGAc                                                                                                                                         <  1:239058/150‑1 (MQ=255)

            cGGCAGGTTAGGCGTAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGAc                                                                                                                                         >  2:1105094/1‑150 (MQ=255)

             ggCAGGTTAGGCGTAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACt                                                                                                                                        >  1:4526133/1‑150 (MQ=255)

              gCAGGTTAGGCGTAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACtt                                                                                                                                       >  2:4330302/1‑150 (MQ=255)

              gCAGGTTAGGCGTAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACtt                                                                                                                                       <  2:4941411/150‑1 (MQ=255)

               cAGGTTAGGCGTAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTc                                                                                                                                      >  2:522408/1‑150 (MQ=255)

               cAGGTTAGGCGTAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTc                                                                                                                                      >  2:522847/1‑150 (MQ=255)

                tggTTAGGCGTAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATTACTTCt                                                                                                                                     <  1:2794624/149‑1 (MQ=255)

                aGGTTAGGCGTAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCt                                                                                                                                     <  1:3217987/150‑1 (MQ=255)

                aGGTTAGGCGTAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCt                                                                                                                                     >  2:281007/1‑150 (MQ=255)

                  gTTAGGCGTAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGc                                                                                                                                   <  1:4211876/150‑1 (MQ=255)

                        cGTAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCg                                                                                                                             <  2:2418461/150‑1 (MQ=255)

                         gTAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGc                                                                                                                            >  1:3155465/1‑150 (MQ=255)

                         gTAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGc                                                                                                                            >  1:1034074/1‑150 (MQ=255)

                         gTAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGc                                                                                                                            >  1:2550226/1‑150 (MQ=255)

                           aGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGC‑tt                                                                                                                         <  2:1947300/150‑1 (MQ=255)

                           aGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGC‑tt                                                                                                                         <  2:4594779/150‑1 (MQ=255)

                            gAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGC‑ttt                                                                                                                        >  1:5074612/1‑150 (MQ=255)

                                    ggTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGC‑TTTTTTTCGCa                                                                                                                <  2:1748191/150‑1 (MQ=255)

                                     gTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGC‑TTTTTTTCGCAg                                                                                                               <  1:2074649/150‑1 (MQ=255)

                                     gTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGC‑TTTTTTTCGCAg                                                                                                               >  1:4780604/1‑150 (MQ=255)

                                      tGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGC‑TTTTTTTCGCAgt                                                                                                              >  1:4096227/1‑150 (MQ=255)

                                      tGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGC‑TTTTTTTCGCAgt                                                                                                              >  2:285447/1‑150 (MQ=255)

                                      tGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGC‑TTTTTTTCGCAgt                                                                                                              <  1:4296879/150‑1 (MQ=255)

                                       gCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGC‑TTTTTTTCGCAgtg                                                                                                             >  2:850524/1‑150 (MQ=255)

                                         gTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGC‑TTTTTTTCGCAGTGTa                                                                                                           <  1:1734338/150‑1 (MQ=255)

                                         gTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGC‑TTTTTTTCGCAGTGTa                                                                                                           <  2:4627884/150‑1 (MQ=255)

                                             aTCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGC‑TTTTTTTCGCAGTGTAGCGt                                                                                                       >  1:2630439/1‑150 (MQ=255)

                                              tCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGC‑TTTTTTTCGCAGTGTAGCGtt                                                                                                      >  2:4465898/1‑150 (MQ=255)

                                               cTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGC‑TTTTTTTCGCAGTGTAGCGttt                                                                                                     <  1:646756/150‑1 (MQ=255)

                                                    cTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGC‑TTTTTTTCGCAGTGTAGCGTTTCACCt                                                                                                <  2:1569633/150‑1 (MQ=255)

                                                      gACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGC‑TTTTTTTCGCAGTGTAGCGTTTCACCTTc                                                                                              >  1:4592232/1‑150 (MQ=255)

                                                        ccAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGC‑TTTTTTTCGCAGTGTAGCGTTTCACCTTctg                                                                                            <  2:939830/150‑1 (MQ=255)

                                                        ccAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGC‑TTTTTTTCGCAGTGTAGCGTTTCACCTTctg                                                                                            >  2:306440/1‑150 (MQ=255)

                                                          aGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGC‑TTTTTTTCGCAGTGTAGCGTTTCACCTTctgct                                                                                          >  1:1641499/1‑150 (MQ=255)

                                                              ccGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGC‑TTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTcc                                                                                      <  2:4153100/150‑1 (MQ=255)

                                                              ccGGAAATCCCATCACCAGACCAGGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGC‑TTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTcc                                                                                      <  2:3108967/150‑1 (MQ=255)

                                                               cGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGC‑TTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCg                                                                                     <  1:1996722/150‑1 (MQ=255)

                                                                ggAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGC‑TTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGc                                                                                    <  2:898798/150‑1 (MQ=255)

                                                                ggAAATCCCATCACCAGACCAGGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCCGGCCGGGAGCAATGACTTCTGCGGTCCGC‑TTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGc                                                                                    <  2:666819/150‑1 (MQ=255)

                                                                  aaaTCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGC‑TTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAg                                                                                  <  2:3217063/150‑1 (MQ=255)

                                                                  aaaTCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGC‑TTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAg                                                                                  >  2:3412600/1‑150 (MQ=255)

                                                                  aaaTCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGC‑TTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAg                                                                                  >  1:56735/1‑150 (MQ=255)

                                                                  aaaTCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGC‑TTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAg                                                                                  >  1:56835/1‑150 (MQ=255)

                                                                    aTCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGC‑TTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGAc                                                                                >  2:4206541/1‑150 (MQ=255)

                                                                     tCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGC‑TTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACa                                                                               <  1:5007985/150‑1 (MQ=255)

                                                                       ccATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGC‑TTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACaaa                                                                             >  2:2721663/1‑150 (MQ=255)

                                                                        cATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGC‑TTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACaaaa                                                                            >  2:4799802/1‑150 (MQ=255)

                                                                        cATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGC‑TTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACaaaa                                                                            <  2:2309980/150‑1 (MQ=255)

                                                                          tCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGC‑TTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGc                                                                          <  1:4811462/150‑1 (MQ=255)

                                                                              cAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGC‑TTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGat                                                                      <  1:2997491/150‑1 (MQ=255)

                                                                              cAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGC‑TTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGat                                                                      <  2:80353/150‑1 (MQ=255)

                                                                                 aCCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGC‑TTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATg                                                                   >  1:159244/1‑150 (MQ=255)

                                                                                   cATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGC‑TTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCg                                                                 >  1:916544/1‑150 (MQ=255)

                                                                                    aTGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGC‑TTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGt                                                                <  2:646309/150‑1 (MQ=255)

                                                                                      gACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGC‑TTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTg                                                              >  1:270125/1‑150 (MQ=255)

                                                                                          gcgagcgTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGC‑TTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGcccc                                                          >  2:1718121/1‑150 (MQ=255)

                                                                                           cgagcgTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGC‑TTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGccccc                                                         <  2:3908406/150‑1 (MQ=255)

                                                                                             agcgTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGC‑TTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCat                                                       >  1:1792726/1‑150 (MQ=255)

                                                                                             agcgTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGC‑TTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCat                                                       >  1:3011795/1‑150 (MQ=255)

                                                                                                  tAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGC‑TTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCa                                                  >  2:4704701/1‑150 (MQ=255)

                                                                                                   aaaGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGC‑TTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCAc                                                 <  1:1883773/150‑1 (MQ=255)

                                                                                                   aaaGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGC‑TTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCAc                                                 <  1:758781/150‑1 (MQ=255)

                                                                                                    aaGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGC‑TTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACg                                                >  1:5108186/1‑150 (MQ=255)

                                                                                                          cTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGCTTTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGc                                           >  1:3960936/1‑150 (MQ=255)

                                                                                                           tGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGC‑TTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTa                                         <  2:3946978/150‑1 (MQ=255)

                                                                                                           tGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGC‑TTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTa                                         >  2:50833/1‑150 (MQ=255)

                                                                                                            gAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGC‑TTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTac                                        >  2:2182888/1‑150 (MQ=255)

                                                                                                              gTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGC‑TTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTacga                                      <  1:1160298/150‑1 (MQ=255)

                                                                                                              gTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGC‑TTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTacga                                      <  2:1034074/150‑1 (MQ=255)

                                                                                                               ttGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGC‑TTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTacgac                                     >  1:1194374/1‑150 (MQ=255)

                                                                                                                 gCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGC‑TTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGt                                   >  2:3518935/1‑150 (MQ=255)

                                                                                                                   tCGTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGC‑TTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTa                                 >  2:3664969/1‑150 (MQ=255)

                                                                                                                     gTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGC‑TTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATg                               >  1:862821/1‑150 (MQ=255)

                                                                                                                     gTCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGC‑TTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATg                               <  1:2350101/150‑1 (MQ=255)

                                                                                                                      tCAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGC‑TTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGc                              >  1:4995401/1‑150 (MQ=255)

                                                                                                                       cAATCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGC‑TTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGCt                             >  1:383091/1‑150 (MQ=255)

                                                                                                                        aaTCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGC‑TTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGCTg                            >  1:1184941/1‑150 (MQ=255)

                                                                                                                          tCATTTCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGC‑TTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGCTGtt                          >  1:1000813/1‑150 (MQ=255)

                                                                                                                             tttCGCCGATACCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGC‑TTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGCGGttttt                       >  1:3683567/1‑150 (MQ=255)

                                                                                                                               tCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGC‑TTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGCTGTTTTTGc                     >  2:3634445/1‑150 (MQ=255)

                                                                                                                                 gCCGATACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGC‑TTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGCTGTTTTTGCAt                   >  1:1560620/1‑150 (MQ=255)

                                                                                                                                  ccGATACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGC‑TTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGCTGTTTTTGCATg                  >  2:1657191/1‑150 (MQ=255)

                                                                                                                                    gATACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGC‑TTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGCTGTTTTTGCATGAg                >  1:2302193/1‑150 (MQ=255)

                                                                                                                                       aCCCGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGC‑TTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGCTGTTTTTGCATGAGGCt             <  1:1233429/150‑1 (MQ=255)

                                                                                                                                       aCCCGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGC‑TTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGCTGTTTTTGCATGAGGCt             <  2:4831843/150‑1 (MQ=255)

                                                                                                                                         ccGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGC‑TTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGCTGTTTTTGCATGAGGCTGt           >  2:3035597/1‑150 (MQ=255)

                                                                                                                                         ccGACACCAGCCCGGGAGCAATGACTTCTGCGGTCCGC‑TTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGCTGTTTTTGCATGAGGCTGt           >  1:3133365/1‑150 (MQ=255)

                                                                                                                                           gACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGC‑TTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGCTGTTTTTGCATGAGGCTGTAc         <  1:4727012/150‑1 (MQ=255)

                                                                                                                                             cccccggCCGGGAGCAATGACTTCTGCGGGCCGC‑TTTTTTGCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGCTGTTTTTGCATGAGGCTGTACCt       <  2:2964418/145‑1 (MQ=255)

                                                                                                                                              acCAGGCCGGGAGCAATGACTTCTGCGGTCCGC‑TTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGCTGTTTTTGCATGAGGCTGTACCTc      >  2:3126789/1‑150 (MQ=255)

                                                                                                                                              acCAGCCCGGGAGCAATGACTTCTGCGGTCCGC‑TTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGCTGTTTTTGCATGAGGCTGTCCCTc      >  1:2141556/1‑150 (MQ=255)

                                                                                                                                              acCAGCCCGGGAGCAATGACTTCTGCGGTCCGC‑TTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGCTGTTTTTGCATGAGGCTGTACCTc      >  1:2141582/1‑150 (MQ=255)

                                                                                                                                              acCAGCCCGGGAGCAATGACTTCTGCGGTCCGC‑TTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGCTGTTTTTGCATGAGGCTGTACCTc      >  1:3566466/1‑150 (MQ=255)

                                                                                                                                               ccAGGCCGGGAGCAATGACTTCTGCGGTCCGC‑TTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGCTGTTTTTGCATGAGGCTGTACCTCg     <  2:4278558/150‑1 (MQ=255)

CGCCGCTGTTAACGGCAGGTTAGGCGTAGAAATAATGGTGCGTTTATCTTTACTGACCAGCGCCGGAAATCCCATCACCAGACCATGACAGCGAGCGTTAAAGCGCCTGAGTTGCTCGTCAATCATTTCGCCGATACCCGACACCAGGCCGGGAGCAATGACTTCTGCGGTCCGC‑TTTTTTTCGCAGTGTAGCGTTTCACCTTCTGCTGTCCGCAGACAAAAGCGGATATGCGTTGCCCCCATATCCACGCCCGCTACGACGTTATGCTGTTTTTGCATGAGGCTGTACCTCGCTT  >  CP009273/4297338‑4297633

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