BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	1,872,520	G→T	95.0%	P183P (CCG→CCT)	yeaP →	diguanylate cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	1,872,520	0	G	T	95.0%	887.5 / 33.5	278	P183P (CCG→CCT)	yeaP	diguanylate cyclase
Reads supporting (aligned to +/- strand): ref base G (8/6); new base T (115/149); total (123/155)
Fisher's exact test for biased strand distribution p-value = 4.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 278 total aligned reads displayed.

TGGAGTGAACGCGCGGAACAGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCGAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCT  >  CP009273/1872374‑1872666

  gAGTGAACGCGCGGAACAGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATcg                                                                                                                                               <  2:1537385/150‑1 (MQ=255)

  gAGTGAACGCGCGGAACAGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATcg                                                                                                                                               <  1:2089929/150‑1 (MQ=255)

    gTGAACGCGCGGAACAGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATcggc                                                                                                                                             <  1:3035026/150‑1 (MQ=255)

    gTGAACGCGCGGAACAGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATcggc                                                                                                                                             >  1:1577701/1‑150 (MQ=255)

     tGAACGCGCGGAACAGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATcggcg                                                                                                                                            >  2:18848/1‑150 (MQ=255)

     tGAACGCGCGGAACAGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATcggcg                                                                                                                                            <  1:2813980/150‑1 (MQ=255)

     tGAACGCGCGGAACAGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATcggcg                                                                                                                                            <  1:3464987/150‑1 (MQ=255)

     tGAACGCGCGGAACAGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATcggcg                                                                                                                                            <  1:771821/150‑1 (MQ=255)

        aCGCGCGGAACAGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCggcgggc                                                                                                                                         <  1:4361260/150‑1 (MQ=255)

           cgcgGAACAGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCTCTGGTTGAACAGCTGCGCGAAGCCAATTCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGAt                                                                                                                                      <  1:4422237/150‑1 (MQ=255)

            gcgGAACAGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGAtt                                                                                                                                     >  2:3230215/1‑150 (MQ=255)

             cgGAACAGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGAttt                                                                                                                                    <  1:2341865/150‑1 (MQ=255)

               gAACAGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGTTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAAACGGGCTACCTAATCGGCGGGCGAttttt                                                                                                                                  >  1:980864/1‑150 (MQ=255)

               gAACAGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGAttttt                                                                                                                                  >  1:4322573/1‑150 (MQ=255)

                aaCAGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTg                                                                                                                                 >  2:537491/1‑150 (MQ=255)

                aaCAGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTg                                                                                                                                 >  2:3286817/1‑150 (MQ=255)

                 aCAGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGa                                                                                                                                <  2:864919/150‑1 (MQ=255)

                  cAGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGaa                                                                                                                               >  2:1037314/1‑150 (MQ=255)

                   aGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGaaa                                                                                                                              >  2:3054963/1‑150 (MQ=255)

                    ggTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGaaaa                                                                                                                             <  1:2732754/150‑1 (MQ=255)

                    ggTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGaaaa                                                                                                                             <  2:2367951/150‑1 (MQ=255)

                     gTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTGAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAAt                                                                                                                            <  1:479103/150‑1 (MQ=255)

                          aCAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGATAATCTgac                                                                                                                       <  1:168039/150‑1 (MQ=255)

                               tATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGacac                                                                                                                  >  2:719893/1‑150 (MQ=255)

                                aTTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACt                                                                                                                 <  2:1848423/150‑1 (MQ=255)

                                aTTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACt                                                                                                                 >  2:2055676/1‑150 (MQ=255)

                                aTTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACt                                                                                                                 <  1:4262629/150‑1 (MQ=255)

                                   cGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGtt                                                                                                              <  2:4372640/150‑1 (MQ=255)

                                    gCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGCTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGttt                                                                                                             <  2:3900689/150‑1 (MQ=255)

                                         aCTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCt                                                                                                        <  1:3565367/150‑1 (MQ=255)

                                         aCTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCt                                                                                                        <  1:2274917/150‑1 (MQ=255)

                                           tGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCg                                                                                                      <  1:159361/150‑1 (MQ=255)

                                             aTTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGcc                                                                                                    >  1:3668942/1‑150 (MQ=255)

                                                   caATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCAt                                                                                              <  2:1024658/150‑1 (MQ=255)

                                                    aaTATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATc                                                                                             >  1:923630/1‑150 (MQ=255)

                                                     atatatTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCt                                                                                            >  2:3160301/1‑150 (MQ=255)

                                                     atatatTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCt                                                                                            <  1:2845544/150‑1 (MQ=255)

                                                     atatatTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACAGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCt                                                                                            <  1:1006868/150‑1 (MQ=255)

                                                          ttCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAAcc                                                                                       <  2:3854165/150‑1 (MQ=255)

                                                          ttCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAAcc                                                                                       <  1:1400149/150‑1 (MQ=255)

                                                            cAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCAt                                                                                     <  2:2508846/150‑1 (MQ=255)

                                                             aaaaaGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATa                                                                                    <  1:2870292/150‑1 (MQ=255)

                                                               aaaGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAg                                                                                  >  2:1954942/1‑150 (MQ=255)

                                                               aaaGAGGCACTGGATGAACAGCAGCGCGAAGCCAATGCTGCGCTGATTGCGAAAACGTATACCGACACGTTAACCGGGCTACCTAATCTCCGGGCGATTTTTGAAAATCTGACGACACAGTTTTCCCTCGCCCGGCATCTTAACCATAAg                                                                                  >  2:3522589/1‑150 (MQ=255)

                                                                  gagGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATa                                                                               >  1:3673586/1‑150 (MQ=255)

                                                                       aCTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGAt                                                                          <  2:2275351/150‑1 (MQ=255)

                                                                         tGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATcg                                                                        <  1:2210897/150‑1 (MQ=255)

                                                                         tGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATcg                                                                        >  1:4525956/1‑150 (MQ=255)

                                                                         tGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATcg                                                                        >  1:1874121/1‑150 (MQ=255)

                                                                            ttGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGt                                                                     >  1:2448783/1‑150 (MQ=255)

                                                                            ttGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATAGCGt                                                                     >  1:3956331/1‑150 (MQ=255)

                                                                              gAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGttt                                                                   <  2:1209638/150‑1 (MQ=255)

                                                                              gAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGttt                                                                   >  1:3497281/1‑150 (MQ=255)

                                                                              gAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGttt                                                                   >  1:1711075/1‑150 (MQ=255)

                                                                                aCAGCTGCGCGAAGCCAATTCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTAt                                                                 >  2:1240856/1‑150 (MQ=255)

                                                                                aCAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTAt                                                                 >  1:1150343/1‑150 (MQ=255)

                                                                                 cAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATc                                                                <  1:4092678/150‑1 (MQ=255)

                                                                                     tGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATc                                                            >  2:1800723/1‑150 (MQ=255)

                                                                                     tGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCGAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATc                                                            >  1:991696/1‑150 (MQ=255)

                                                                                      gcgcGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCt                                                           <  2:1254889/150‑1 (MQ=255)

                                                                                      gcgcGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCt                                                           <  2:2879299/150‑1 (MQ=255)

                                                                                       cgcgAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTg                                                          <  1:4341184/150‑1 (MQ=255)

                                                                                         cgAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGa                                                        >  2:4092711/1‑150 (MQ=255)

                                                                                          gAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGAt                                                       <  2:2791736/150‑1 (MQ=255)

                                                                                           aaGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATa                                                      >  1:3206341/1‑150 (MQ=255)

                                                                                               cAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACtt                                                  <  1:820107/150‑1 (MQ=255)

                                                                                               cAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACtt                                                  <  1:1225824/150‑1 (MQ=255)

                                                                                                aaTGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTc                                                 >  2:4448031/1‑150 (MQ=255)

                                                                                                  tgctgcGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCaa                                               >  2:3375448/1‑150 (MQ=255)

                                                                                                   gctgcGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCaaa                                              <  2:2177975/150‑1 (MQ=255)

                                                                                                     tgcGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAAtt                                            >  2:1077620/1‑150 (MQ=255)

                                                                                                      gcgcTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCGAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTa                                           <  2:2567117/150‑1 (MQ=255)

                                                                                                         cTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATc                                        >  2:1393663/1‑150 (MQ=255)

                                                                                                         cTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATc                                        <  2:3641139/150‑1 (MQ=255)

                                                                                                          tGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCa                                       <  2:1697682/150‑1 (MQ=255)

                                                                                                            aTTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAAt                                     <  2:2467061/150‑1 (MQ=255)

                                                                                                            aTTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAAt                                     <  1:2569645/150‑1 (MQ=255)

                                                                                                             ttGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCATATTAATCAATg                                    <  2:4134727/150‑1 (MQ=255)

                                                                                                             ttGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATg                                    <  2:937432/150‑1 (MQ=255)

                                                                                                              tGCGCAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGa                                   <  1:2322366/150‑1 (MQ=255)

                                                                                                               gcgcAATCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGAt                                  <  1:741120/150‑1 (MQ=255)

                                                                                                                   aaTCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGtt                              >  2:4229369/1‑150 (MQ=255)

                                                                                                                   aaTCGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGtt                              >  1:4115662/1‑150 (MQ=255)

                                                                                                                      cGTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTg                           <  1:56305/150‑1 (MQ=255)

                                                                                                                       gTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTgg                          >  1:1515644/1‑150 (MQ=255)

                                                                                                                       gTATACCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTgg                          >  2:2543336/1‑150 (MQ=255)

                                                                                                                        tataCCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGt                         >  2:1630764/1‑150 (MQ=255)

                                                                                                                         ataCCGACTCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTc                        >  2:1386080/1‑150 (MQ=255)

                                                                                                                              gACTCGTTAACCGGGCTACCTAATCGGCGTGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGTCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCATATTAATCAATGATCGTTTTGGTCataat                   <  1:2038740/150‑1 (MQ=255)

                                                                                                                                 tCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGt                <  1:10513/150‑1 (MQ=255)

                                                                                                                                 tCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGt                <  1:1928273/150‑1 (MQ=255)

                                                                                                                                 tCGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGt                <  1:4162515/150‑1 (MQ=255)

                                                                                                                                  cGTTAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTg               <  1:3778574/150‑1 (MQ=255)

                                                                                                                                    ttAACCGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGc             <  2:762025/150‑1 (MQ=255)

                                                                                                                                        ccGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATc         <  2:2893728/150‑1 (MQ=255)

                                                                                                                                         cGGGCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCt        <  2:3447568/150‑1 (MQ=255)

                                                                                                                                            gCTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGtt     >  1:4387360/1‑150 (MQ=255)

                                                                                                                                             cTACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGttt    >  1:2141854/1‑150 (MQ=255)

                                                                                                                                              tACCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTtc   <  2:2953307/150‑1 (MQ=255)

                                                                                                                                               aCCTAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTtct  >  1:1762959/1‑150 (MQ=255)

TGGAGTGAACGCGCGGAACAGGTTTTACAGTTATTCGCCGGACTGATTGCACAATATATTCAAAAAGAGGCACTGGTTGAACAGCTGCGCGAAGCCAATGCTGCGCTGATTGCGCAATCGTATACCGACTCGTTAACCGGGCTACCGAATCGGCGGGCGATTTTTGAAAATCTGACGACACTGTTTTCCCTCGCCCGGCATCTTAACCATAAGATAATGATCGCGTTTATCGATCTGGATAACTTCAAATTAATCAATGATCGTTTTGGTCATAATAGTGGCGATCTGTTTCT  >  CP009273/1872374‑1872666

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