BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	4,062,424	G→T	65.3%	I71I (ATC→ATA)	yihQ ←	alpha‑glucosidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	4,062,424	0	G	T	65.3%	313.6 / inf	323	I71I (ATC→ATA)	yihQ	alpha‑glucosidase
Reads supporting (aligned to +/- strand): ref base G (44/68); new base T (75/135); total (119/204)
Fisher's exact test for biased strand distribution p-value = 5.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 326 total aligned reads displayed.

AGGCGCAGCCAGATACGGTTGTGGTTAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTATGGGTTAAAATAAGACGTTGTTGAAAATGTAG  >  CP009273/4062275‑4062573

   cgcAGCCAGATACGGTTGTGGTTAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATg                                                                                                                                                    <  1:2961599/150‑1 (MQ=255)

     cAGCCAGATACGGTTGTGGTTAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACGATGGc                                                                                                                                                  <  2:5227752/150‑1 (MQ=255)

         cAGATACGGTTGTGGTTAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCg                                                                                                                                              <  2:5143056/150‑1 (MQ=255)

         cAGATACGGTTGTGGTTAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACGATGGCGTCg                                                                                                                                              >  1:1735727/1‑150 (MQ=255)

           gATACGGTTGTGGTTAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGt                                                                                                                                            <  1:2298852/150‑1 (MQ=255)

            aTACGGTTGTGGTTAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCTTGGCGCTAATGTCAGAACCACGGCTGAAAAGAATTAACCAACCATCCTGTGACTGGCTGACTATGTCGTCGGTa                                                                                                                                           >  2:2716226/1‑150 (MQ=255)

            aTACGGTTGTGGTTAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGTTGACTGGCTGACGATGGCGTCGGTa                                                                                                                                           >  2:2546477/1‑150 (MQ=255)

            aTACGGTTGTGGTTAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTa                                                                                                                                           <  1:204773/150‑1 (MQ=255)

               cGGTTGTGGTTAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAgc                                                                                                                                        >  2:3843421/1‑150 (MQ=255)

                ggTTGTGGTTAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAgcg                                                                                                                                       >  1:3271441/1‑150 (MQ=255)

                ggTTGTGGTTAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAgcg                                                                                                                                       >  1:3318604/1‑150 (MQ=255)

                 gTTGTGGTTAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACGATGGCGTCGGTAAgcgc                                                                                                                                      >  1:1257045/1‑150 (MQ=255)

                  ttGTGGTTAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACGATGGCGTCGGTAAGCGCa                                                                                                                                     <  1:2514414/150‑1 (MQ=255)

                      ggTTAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAAttt                                                                                                                                 >  2:4770864/1‑150 (MQ=255)

                       gTTAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAAtttt                                                                                                                                <  2:2855128/150‑1 (MQ=255)

                        ttAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTtc                                                                                                                               >  1:4537490/1‑150 (MQ=255)

                          aaGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTtctc                                                                                                                             >  1:1717687/1‑150 (MQ=255)

                           aGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCACCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATTTCAGAACCACGGATTAAATGAATTAACCAACCATCCGGTGACTGGATGACTATGGCTTCTGTAAGCGCAATTTTCTcc                                                                                                                            >  2:1585449/1‑150 (MQ=255)

                            ggTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCt                                                                                                                           <  2:1799499/150‑1 (MQ=255)

                             gTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTg                                                                                                                          >  1:4901177/1‑150 (MQ=255)

                                gTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAg                                                                                                                       <  2:2171498/150‑1 (MQ=255)

                                  cGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGtt                                                                                                                     <  1:1861800/150‑1 (MQ=255)

                                  cGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGtt                                                                                                                     <  1:443357/150‑1 (MQ=255)

                                   gTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGttt                                                                                                                    <  1:4383471/150‑1 (MQ=255)

                                   gTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGttt                                                                                                                    <  1:3391312/150‑1 (MQ=255)

                                   gTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGttt                                                                                                                    <  1:1032242/150‑1 (MQ=255)

                                    ttttGTATTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGATAATGTCAGAACCACGTCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTa                                                                                                                   <  1:1324377/150‑1 (MQ=255)

                                     tttGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTAt                                                                                                                  <  1:707041/150‑1 (MQ=255)

                                         tAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCttt                                                                                                              <  2:5076529/150‑1 (MQ=255)

                                         tAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCttt                                                                                                              >  1:3445395/1‑150 (MQ=255)

                                         tAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCttt                                                                                                              <  1:2817914/150‑1 (MQ=255)

                                           gTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTTAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTaa                                                                                                            <  2:5193500/150‑1 (MQ=255)

                                           gTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTTACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTaa                                                                                                            >  1:4633815/1‑150 (MQ=255)

                                            ttCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGGAGTTTATCTTTAAt                                                                                                           >  2:2637432/1‑150 (MQ=255)

                                                 gCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGa                                                                                                      <  2:3622718/150‑1 (MQ=255)

                                                   aaTAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGaaa                                                                                                    >  1:1619132/1‑150 (MQ=255)

                                                   aaTAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGaaa                                                                                                    >  1:4214322/1‑150 (MQ=255)

                                                   aaTAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGaaa                                                                                                    >  1:2647947/1‑150 (MQ=255)

                                                    aTAAACGCACCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAAt                                                                                                   <  1:2724936/150‑1 (MQ=255)

                                                       aaCGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTAc                                                                                                >  1:2571033/1‑150 (MQ=255)

                                                        aCGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTAcc                                                                                               >  1:3814331/1‑150 (MQ=255)

                                                        aCGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTAcc                                                                                               <  1:4488849/150‑1 (MQ=255)

                                                         cGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACcg                                                                                              >  1:4735276/1‑150 (MQ=255)

                                                          gCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACcgc                                                                                             <  2:1968293/150‑1 (MQ=255)

                                                                 aTCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACat                                                                                      <  2:2405469/150‑1 (MQ=255)

                                                                    gtcgGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATc                                                                                   <  2:3739067/150‑1 (MQ=255)

                                                                    gtcgGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATc                                                                                   <  2:4011412/150‑1 (MQ=255)

                                                                     tcgGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCg                                                                                  <  1:2235960/150‑1 (MQ=255)

                                                                         cAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGatat                                                                              <  1:190854/150‑1 (MQ=255)

                                                                         cAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGatat                                                                              >  2:2277291/1‑150 (MQ=255)

                                                                         cAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGatat                                                                              >  1:329872/1‑150 (MQ=255)

                                                                           gagaTATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATcc                                                                            <  1:3635485/150‑1 (MQ=255)

                                                                                atTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCtat                                                                       <  2:3034738/150‑1 (MQ=255)

                                                                                 ttCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCtata                                                                      <  1:501104/150‑1 (MQ=255)

                                                                                  tCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATAc                                                                     <  1:4015661/150‑1 (MQ=255)

                                                                                    aGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCt                                                                   <  1:2187110/150‑1 (MQ=255)

                                                                                     gCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTg                                                                  <  1:1392844/150‑1 (MQ=255)

                                                                                      cGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGa                                                                 <  1:3082746/150‑1 (MQ=255)

                                                                                        tGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGc                                                               <  1:3728510/150‑1 (MQ=255)

                                                                                         ggCGCTAATGTCAGATCCACGGCTGAAATGAAGTAACCAACCATCCGGTGACTGGATGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGcc                                                              <  2:3155304/150‑1 (MQ=255)

                                                                                            tatAATGAAAGAACCACGGCTGAAATGAATTAAACAACCATCCGTTTACTGGATTACGAATGAGTAGTTAAGAGCAATTTTCTAGTGTAGTTTATCTTTAATTATGAAATAACAGAGGAACAAAACGATATCCTCTAAAACTGAGACaaa                                                           >  2:4560382/3‑149 (MQ=255)

                                                                                            gcTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAAt                                                           <  2:4047069/150‑1 (MQ=255)

                                                                                             cTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATc                                                          <  2:2647878/150‑1 (MQ=255)

                                                                                             cTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATc                                                          <  2:4214322/150‑1 (MQ=255)

                                                                                                aTGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCAt                                                       <  1:3065177/150‑1 (MQ=255)

                                                                                                   tCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATaaa                                                    <  2:4329804/150‑1 (MQ=255)

                                                                                                    cAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAAc                                                   <  2:2318372/150‑1 (MQ=255)

                                                                                                      gAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACaa                                                 >  1:3250900/1‑150 (MQ=255)

                                                                                                       aaCCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAg                                                >  1:1768922/1‑150 (MQ=255)

                                                                                                       aaCCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAg                                                >  2:1000781/1‑150 (MQ=255)

                                                                                                       aaCCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAg                                                <  1:4663374/150‑1 (MQ=255)

                                                                                                       aaCCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAg                                                >  2:2752944/1‑150 (MQ=255)

                                                                                                        aCCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAgg                                               >  2:748642/1‑150 (MQ=255)

                                                                                                        aCCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAgg                                               <  1:2745103/150‑1 (MQ=255)

                                                                                                        aCCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAgg                                               >  2:1077764/1‑150 (MQ=255)

                                                                                                         ccACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGa                                              <  1:110738/150‑1 (MQ=255)

                                                                                                         ccACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGa                                              >  1:4617834/1‑150 (MQ=255)

                                                                                                            cGGCTGAAATGAATTACCCAACCATCCGGTGTCTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTa                                           <  2:3392734/150‑1 (MQ=255)

                                                                                                            cGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTa                                           <  1:83683/150‑1 (MQ=255)

                                                                                                             ggCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTAt                                          >  2:4989389/1‑150 (MQ=255)

                                                                                                             ggCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTAt                                          >  2:2086231/1‑150 (MQ=255)

                                                                                                                tGAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCtt                                       <  2:3855740/150‑1 (MQ=255)

                                                                                                                 gAAATGAATTAACCAACCATCCGGTGACTGGCTGACTTTGGCGTCGGTTAGCGCAATTTTCTCCTGTAGTTTAACTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCttt                                      <  2:4369156/150‑1 (MQ=255)

                                                                                                                 gAAATGAATTAACCAACCATCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCttt                                      <  2:621233/150‑1 (MQ=255)

                                                                                                                     tGAATTAACCAACCATCCGGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTa                                  >  1:5083874/1‑150 (MQ=255)

                                                                                                                      gAATTAACCAACCATCCGGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTAt                                 >  1:572694/1‑150 (MQ=255)

                                                                                                                        aTTAACCAACCATCCGGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTATgg                               >  2:3380245/1‑150 (MQ=255)

                                                                                                                         ttAACCAACCCTCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTATggg                              <  2:855794/150‑1 (MQ=255)

                                                                                                                          tAACCAACCATCCGGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTATGGGt                             <  1:5218735/150‑1 (MQ=255)

                                                                                                                            accaaccaTCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTATGGGTTa                           >  1:4820365/1‑150 (MQ=255)

                                                                                                                              caaccaTCCGGTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTATGGGTTaaa                         >  1:117264/1‑150 (MQ=255)

                                                                                                                                   aTCCGGTGACTGGCTGACTATGTCGACTGTAAGCGCAATTTTCACCTCTAGTTTATCTTAAATGCTGAAATTACAGAGGAAAATATCGATATCCGCTATACCTGAGCCAAACCATAAACAATGATTATCTTTGCTCTGGGTTAAAATAAg                    >  2:4664111/1‑150 (MQ=255)

                                                                                                                                     ccGGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTATGGGTTAAAATAAGAc                  <  2:881906/150‑1 (MQ=255)

                                                                                                                                       ggTGACTGGCTGACTCTGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTATGGGTTAAAATAAGACgt                <  2:5109897/150‑1 (MQ=255)

                                                                                                                                        gTGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTATGGGTTAAAATAAGACgtt               <  2:3933446/150‑1 (MQ=255)

                                                                                                                                         tGACTGGCTGACTATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTATGGGTTAAAATAAGACgttg              <  1:2433271/150‑1 (MQ=255)

                                                                                                                                           aCTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTATGGGTTAAAATAAGACgttttt            >  1:3871208/1‑150 (MQ=255)

                                                                                                                                           aCTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTATGGGTTAAAATAAGACgttgtt            >  1:1299316/1‑150 (MQ=255)

                                                                                                                                           aCTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTATGGGTTAAAATAAGACgttgtt            <  2:3250900/150‑1 (MQ=255)

                                                                                                                                                     gATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTATGGGTTAAAATAAGACGTTGTTGAAAATGTAg  >  1:4121639/1‑150 (MQ=255)

AGGCGCAGCCAGATACGGTTGTGGTTAAGGTTGTCGTTTTGTAGTTCCAGCAATAAACGCCCCTGATCGTCGGCAGAGATATTCAGCGTGGCGCTAATGTCAGAACCACGGCTGAAATGAATTAACCAACCATCCGGTGACTGGCTGACGATGGCGTCGGTAAGCGCAATTTTCTCCTGTAGTTTATCTTTAATGCTGAAATTACCGCGGAACATATCGATATCCGCTATACCTGAGCCAATCCATAAACAAGGATTATCTTTGCTATGGGTTAAAATAAGACGTTGTTGAAAATGTAG  >  CP009273/4062275‑4062573

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