BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	4,276,940	G→T	56.2%	F120L (TTC→TTA)	acs ←	acetyl‑CoA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	4,276,940	0	G	T	56.2%	4.6 / inf	336	F120L (TTC→TTA)	acs	acetyl‑CoA synthetase
Reads supporting (aligned to +/- strand): ref base G (88/59); new base T (81/108); total (169/167)
Fisher's exact test for biased strand distribution p-value = 2.13e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 338 total aligned reads displayed.

CTTCCGGCGAGAAGCCGCCGAAAATCACCGAATGCACCGCGCCAATGCGGGCGCAGGCCAGCATCGCAACCGCGGCTTCCGGCACCATCGGCATATAAATCGCCACCACATCACCTTTTTTAATGCCCAGCTCGAGCAGGGTATTGGCGAAGCGGCAGACGTCGCGGTGCAGCTCTTTATAGCTGATATGTTTGCTCTGGCTGGCGTCGTCGCCTTCCCAGATGATGGCGGTACGATCGCCGTTTTCTTGCAGATGGCGGTCAAGGCAGTTTGCCGCCAGATTCAGCGTGCCGTCCTC  >  CP009273/4276792‑4277089

cTTCCGGCGAGAAGCCGCCGAAAATCACCGAATGCACCGCGCCAATGCGGGCGCAGGCCAGCATCGCAACCGCGGCTTCCGGCACCATCGGCATATAAATCGCCACCACATCACCTTTTTTAATGCCCAGCTCGAGCAGGGTATTGGCGa                                                                                                                                                      >  1:413054/1‑150 (MQ=255)

cTTCCGGCGAGAAGCCGCCGAAAATCACCGAATGCACCGCGCCAATGCGGGCGCAGGCCAGCATCGCAACCGCGGCTTCCGGCACCATCGGCATATAAATCGCCACCACATCACCTTTTTTAATGCCCAGCTCGAGCAGGGTATTGGCGa                                                                                                                                                      >  1:4051507/1‑150 (MQ=255)

 ttCCGGCGAGAAGCCGCCGAAAATCACCGAATGCACCGCGCCAATGCGGGCGCAGGCCAGCATCGCAACCGCGGCTTCCGGCACCATCGGCATATAAATCGCCACCACATCACCTTTTTTAATGCCCAGCTCGAGCAGGGTATTGGCTaa                                                                                                                                                     <  2:2623167/150‑1 (MQ=255)

 ttCCGGCGAGAAGCCGCCGAAAATCACCGAATGCACCGCGCCAATGCGGGCGCAGGCCAGCATCGCAACCGCGGCTTCCGGCACCATCGGCATATAAATCGCCACCACATCACCTTTTTTAATGCCCAGCTCGAGCAGGGTATTGGCTaa                                                                                                                                                     >  2:2661314/1‑150 (MQ=255)

 ttCCGGCGAGAAGCCGCCGAAAATCACCGAATGCACCGCGCCAATGCGGGCGCAGGCCAGCATCGCAACCGCGGCTTCCGGCACCATCGGCATATAAATCGCCACCACATCACCTTTTTTAATGCCCAGCTCGAGCAGGGTATTGGCTaa                                                                                                                                                     >  1:3333181/1‑150 (MQ=255)

       cGAGAAGCCGCCGAAAATCACCGAATGCACCGCGCCAATGCGGGCGCAGGCCAGCATCGCAACCGCGGCTTCCGGCACCATCGGCATATAAATCGCCACCACATCACCTTTTTTAATGCCCAGCTCGAGCAGGGTATTGGCTAAGCGGCa                                                                                                                                               >  2:3552157/1‑150 (MQ=255)

             gccgccGAAAATCACCGAATGCACCGCGCCAATGCGGGCGCAGGCCAGCATCGCAACCGCGGCTTCCGGCACCATCGGCATATAAATCGCCACCACATCACCTTTTTTAATGCCCAGCTCGAGCAGGGTATTGGCTAAGCGGCAGACGTc                                                                                                                                         <  1:4129519/150‑1 (MQ=255)

             gccgccGAAAATCACCGAATGCACCGCGCCAATGCGGGCGCAGGCCAGCATCGCAACCGCGGCTTCCGGCACCATCGGCATATAAATCGCCACCACATCACCATTTTTAATGCCCAGCTCGAGCAGGGTATTGGCTAAGCGGCAGACGTc                                                                                                                                         <  2:3387780/150‑1 (MQ=255)

              ccgccgAAAATCACCGAATGCACCGCGCCAATGCGGGCGCAGGCCAGCATCGCAACCGCGGCTTCCGGCACCATCGGCATATAAATCGCCACCACATCACCTTTTTTAATGCCCAGCTCGAGCAGGGTATTGGCGAAGCGGCAGACGTcg                                                                                                                                        >  2:3732953/1‑150 (MQ=255)

               cgccgAAAATCACCGAATGCACCGCGCCAATGCGGGCGCAGGCCAGCATCGCAACCTCGGCTTCCGGCACCATCGGCATATAAATCGCCACCACATCACCTTTTTTAATGCCCAGCTCGAGCAGGGTATTGGCGAAGCGGCAGACGTcgc                                                                                                                                       >  2:3357421/1‑150 (MQ=255)

               cgccgAAAATCACCGAATGCACCGCGCCAATGCGGGCGCAGGCCAGCATCGCAACCGCGGCTTCCGGCACCATCGGCATATAAATCGCCACCACATCACCTTTTTTAATGCCCAGCTCGAGCAGGGTATTGGCGAAGCGGCAGACGTcgc                                                                                                                                       >  2:4552767/1‑150 (MQ=255)

                gccgAAAATCACCGAATGCACCGCGCCAATGCGGGCGCAGGCCAGCATCGCAACCGCGGCTTCCGGCACCATCGGCATATAAATCGCCACCACATCACCTTTTTTAATGCCCAGCTCGAGCAGGGTATTGGCGAAGCGGCAGACGTcgcg                                                                                                                                      <  2:3954908/150‑1 (MQ=255)

                   gAAAATCACCGAATGCACCGCGCCAATGCGGGCGCAGGCCAGCATCGCAACCGCGGCTTCCGGCACCATCGGCATATAAATCGCCACCACATCACCTTTTTTAATGCCCAGCTCGAGCAGGGTATTGGCTAAGCGGCAGACGTCGCGGTg                                                                                                                                   <  1:575099/150‑1 (MQ=255)

                   gAAAATCACCGAATGCACCGCGCCAATGCGGGCGCAGGCCAGCATCGCAACCGCGGCTTCCGGCACCATCGGCATATAAATCGCCACCACATCACCTTTTTTAATGCCCAGCTCGAGCAGGGTATTGGCGAAGCGGCAGACGTCGCGGTg                                                                                                                                   >  2:3416690/1‑150 (MQ=255)

                    aaaaTCACCGAATGCACCGCGCCAATGCGGGCGCAGGCCAGCATCGCAACCGCGGCTTCCGGCACCATCGGCATATAAATCGCCACCACATCACCTTTTTTAATGCCCAGCTCGAGCAGGGTATTGGCTAAGCGGCAGACGTCGCGGTGc                                                                                                                                  >  2:4338464/1‑150 (MQ=255)

                     aaaTCACCGAATGCACCGCGCCAATGCGGGCGCAGGCCAGCATCGCAACCGCGGCTTCCGGCACCATCGGCATATAAATCGCCACCACATCACCTTTTTTAATGCCCAGCTCGAGCAGGGTATTGGCTAAGCGGCAGACGTCGCGGTGCa                                                                                                                                 <  1:1797935/150‑1 (MQ=255)

                     aaaTCACCGAATGCACCGCGCCAATGCGGGCGCAGGCCAGCATCGCAACCGCGGCTTCCGGCACCATCGGCATATAAATCGCCACCACATCACCTTTTTTAATGCCCAGCTCGAGCAGGGTATTGGCTAAGCGGCAGACGTCGCGGTGCa                                                                                                                                 >  2:736690/1‑150 (MQ=255)

                      aaTCACCGAATGCACCGCGCCAATGCGGGCGCAGGCCAGCATCGCAACCGCGGCTTCCGGCACCATCGGCATATAAATCGCCACCACATCACCTTTTTTAATGCCCAGCTCGAGCAGGGTATTGGCGAAGCGGCAGACGTCGCGGTGCAg                                                                                                                                >  2:648033/1‑150 (MQ=255)

                       aTCACCGAATGCACCGCGCCAATGCGGGCGCAGGCCAGCATCGCAACCGCGGCTTCCGGCACCATCGGCATATAAATCGCCACCACATCACCTTTTTTAATGCCCAGCTCGAGCAGGGTATTGGCTAAGCGGCAGACGTCGCGGTGCAGc                                                                                                                               <  2:275998/150‑1 (MQ=255)

                       aTCACCGAATGCACCGCGCCAATGCGGGCGCAGGCCAGCATCGCAACCGCGGCTTCCGGCACCATCGGCATATAAATCGCCACCACATCACCTTTTTTAATGCCCAGCTCGAGCAGGGTATTGGCGAAGCGGCAGACGTCGCGGTGCAGc                                                                                                                               >  1:1910003/1‑150 (MQ=255)

                         cACCGAATGCACCGCGCCAATGCGGGCGCAGGCCAGCATCGCAACCGCGGCTTCCGGCACCATCGGCATATAAATCGCCACCACATCACCTTTTTTAATGCCCAGCTCGAGCAGGGTATTGGCTAAGCGGCAGACGTCGCGGTGCAGctc                                                                                                                             >  1:2128606/1‑150 (MQ=255)

                            cGAATGCACCGCGCCAATGCGGGCGCAGGCCAGCATCGCAACCGCGGCTTCCGGCACCATCGGCATATAAATCGCCACCACATCACCTTTTTTAATGCCCATCTCGAGCAGGGTATTGGCTAAGAGGCAGAAGTCGCGGTGCAGCTCttt                                                                                                                          <  2:360894/150‑1 (MQ=255)

                              aaTGCACCGCGCCAATGCGGGCGCAGGCCAGCATCGCAACCGCGGCTTCCGGCACCATCGGCATATAAATCGCCACCACATCACCTTTTTTAATGCCCAGCTCGAGCAGGGTATTGGCTAAGCGGCAGACGTCGCGGTGCAGCTCTTtat                                                                                                                        <  2:3070913/150‑1 (MQ=255)

                                tGCACCGCGCCAATGCGGGCGCAGGCCAGCATCGCAACCGCGGCTTCCGGCACCATCGGCATATAAATCGCCACCACATCACCTTTTTTAATGCCCAGCTCGAGCAGGGTATTGGCTAAGCGGCAGACGTCGCGGTGCAGCTCTTTATAg                                                                                                                      <  2:4087196/150‑1 (MQ=255)

                                tGCACCGCGCCAATGCGGGCGCAGGCCAGCATCGCAACCGCGGCTTCCGGCACCATCGGCATATAAATCGCCACCACATCACCTTTTTTAATGCCCAGCTCGAGCAGGGTATTGGCTAAGCGGCAGACGTCGCGGTGCAGCTCTTTATAg                                                                                                                      >  2:2262520/1‑150 (MQ=255)

                                tGCACCGCGCCAATGCGGGCGCAGGCCAGCATCGCAACCGCGGCTTCCGGCACCATCGGCATATAAATCGCCACCACATCACCTTTTTTAATGCCCAGCTCGAGCAGGGTATTGGCGAAGCGGCAGACGTCGCGGTGCAGCTCTTTATAg                                                                                                                      <  1:356792/150‑1 (MQ=255)

                                 gCACCGCGCCAATGCGGGCGCAGGCCAGCATCGCAACCGCGGCTTCCGGCACCATCGGCATATAAATCGCCACCACATCACCTTTTTTAATGCCCAGCTCGAGCAGGGTATTGGCTAAGCGGCAGACGTCGCGGTGCAGCTCTTTATAGc                                                                                                                     <  1:3353551/150‑1 (MQ=255)

                                     cgcgCCAATGCGGGCGCAGGCCAGCATCGCAACCGCGGCTTCCGGCACCATCGGCATATAAATCGCCACCACATCACCTTTTTTAATGCCCAGCTCGAGCAGGGTATTGGCTAAGCGGCAGACGTCGCGGTGCAGCTCTTTATAGCTGat                                                                                                                 <  2:2978254/150‑1 (MQ=255)

                                     cgcgCCAATGCGGGCGCAGGCCAGCATCGCAACCGCGGCTTCCGGCACCATCGGCATATAAATCGCCACCACATCACCTTTTTTAATGCCCAGCTCGAGCAGGGTATTGGCTAAGCGGCAGACGTCGCGGTGCAGCTCTTTATAGCTGat                                                                                                                 >  1:1474391/1‑150 (MQ=255)

                                     cgcgCCAATGCGGGCGCAGGCCAGCATCGCAACCGCGGCTTCCGGCACCATCGGCATATAAATCGCCACCACATCACCTTTTTTAATGCCCAGCTCGAGCAGGGTATTGGCTAAGCGGCAGACGTCGCGGTGCAGCTCTTTATAGCTGat                                                                                                                 >  1:3799514/1‑150 (MQ=255)

                                        gcCAATGCGGGCGCAGGCCAGCATCGCAACCGCGGCTTCCGGCACCATCGGCATATAAATCGCCACCACATCACCTTTTTTAATGCCCAGCTCGAGCAGGGTATTGGCTAAGCGGCAGACGTCGCGGTGCAGCTCTTTATAGCTGATATg                                                                                                              <  1:320017/150‑1 (MQ=255)

                                         ccAATGCGGGCGCAGGCCAGCATCGCAACCGCGGCTTCCGGCACCATCGGCATATAACTCGCCACCACATCACCTTTGTTAATGCCCAGCTCGAGCAGGGTATTGGCGAAGCGGCAGACGTCGCGGTGCAGCTCTTTATAGCTGAGATgg                                                                                                             >  1:174448/1‑149 (MQ=255)

                                          cAATGCGGGCGCAGGCCAGCATCGCAACCGCGGCTTCCGGCACCATCGGCATATAAATCGCCACCACATCACCTTTTTTAATGCCCAGCTCGAGCAGGGTATTGGCTAAGCGGCAGACGTCGCGGTGCAGCTCTTTATAGCTGATATGtt                                                                                                            <  1:2400149/150‑1 (MQ=255)

                                           aaTGCGGGCGCAGGCCAGCATCGCAACCGCGGCTTCCGGCACCATCGGCATATAAATCGCCACCACATCACCTTTTTTAATGCCCAGCTCGAGCAGGGTATTGGCTAAGCGGCAGACGTCGCGGTGCAGCTCTTTATAGCTGATATGttt                                                                                                           >  2:3459736/1‑150 (MQ=255)

                                            aTGCGGGCGCAGGCCAGCATCGCAACCGCGGCTTCCGGCACCATCGGCATATAAATCGCCACCACATCACCTTTTTTAATGCCCAGCTCGAGCAGGGTATTGGCTAAGCGGCAGACGTCGCGGTGCAGCTCTTTATAGCTGATATGTTTg                                                                                                          <  2:2350663/150‑1 (MQ=255)

                                              gCGGGCGCAGGCCAGCATCGCAACCGCGGCTTCCGGCACCATCGGCATATAAATCGCCACCACATCACCTTTTTTAATGCCCAGCTCGAGCAGGGTATTGGCTAAGCGGCAGACGTCGCGGTGCAGCTCTTTATAGCTGATATGTTTGct                                                                                                        <  2:1563181/150‑1 (MQ=255)

                                                gggCGCAGGCCAGCATCGCAACCGCGGCTTCCGGCACCATCGGCATATAAATCGCCACCACATCACCTTTTTTAATGCCCAGCTCGAGCAGGGTATTGGCTAAGCGGCAGACGTCGCGGTGCAGCTCTTTATAGCTGATATGTTTGctct                                                                                                      >  1:1581026/1‑150 (MQ=255)

                                                   cgcAGGCCAGCATCGCAACCGCGGCTTCCGGCACCATCGGCATATAAATCGCCACCACATCACCTTTTTTAATGCCCAGCTCGAGCAGGGTATTGGCTAAGCGGCAGACGTCGCGGTGCAGCTCTTTATAGCTGATATGTTTGCTctggc                                                                                                   >  2:3110466/1‑150 (MQ=255)

                                                    gcAGGCCAGCATCGCAACCGCGGCTTCCGGCACCATCGGCATATAAATCGCCACCACATCACCTTTTTTAATGCCCAGCTCGAGCAGGGTATTGGCTAAGCGGCAGACGTCGCGGTGCAGCTCTTTATAGCTGATATGTTTGCTctggct                                                                                                  <  1:2534212/150‑1 (MQ=255)

                                                       ggCCAGCATCGCAACCGCGGCTTCCGGCACCATCGGCATATAAATCGCCACCACATCACCTTTTTTAATGCCCAGCTCGAGCAGGGTATTGGCGAAGCGGCAGACGTCGCGGTGCAGCTCTTTATAGCTGATATGTTTGCTCtggctggc                                                                                               <  1:3610694/150‑1 (MQ=255)

                                                           aGCATCGCAACCGCGGCTTCCGGCACCATCGGCATATAAATCGCCACCACATCACCTTTTTTAATGCCCAGCTCGAGCAGGGTATTGGCTAAGCGGCAGACGTCGCGGTGCAGCTCTTTATAGCTGATATGTTTGCTCTGGCTGGcgtcg                                                                                           <  1:3465542/150‑1 (MQ=255)

                                                               tCGCAACCGCGGCTTCCGGCACCATCGGCATATAAATCGCCACCACATCACCTTTTTTAATGCCCAGCTCGAGCAGGGTATTGGCGAAGCGGCAGACGTCGCGGTGCAGCTCTTTATAGCTGATATGTTTGCTCTGGCTGGCGTCGTCGc                                                                                       >  1:1920695/1‑150 (MQ=255)

                                                                  cAACCGCGGCTTCCGGCACCATCGGCATATAAATCGCCACCACATCACCTTTTTTAATGCCCAGCTCGAGCAGGGTATTGGCGAAGCGGCAGACGTCGCGGTGCAGCTCTTTATAGCTGATATGTTTGCTCTGGCTGGCGTCGTCGCCtt                                                                                    <  2:4622218/150‑1 (MQ=255)

                                                                   aaCCGCGGCTTCCGGCACCATCGGCATATAAATCGCCACCACATCACCTTTTTTAATGCCCAGCTCGAGCAGGGTATTGGCTAAGCGGCAGACGTCGCGGTGCAGCTCTTTATAGCTGATATGTTTGCTCTGGCTGGCGTCGTCGCCTTc                                                                                   <  1:1100126/150‑1 (MQ=255)

                                                                      cgcgGCTTCCGGCACCATCGGCATATAAATCGCCACCACATCACCTTTTTTAATGCCCAGCTCGAGCAGGGTATTGGCTAAGCGGCAGACGTCGCGGTGCAGCTCTTTATAGCTGATATGTTTGCTCTGGCTGGCGTCGTCGCCTTCCCa                                                                                <  2:4613323/150‑1 (MQ=255)

                                                                         ggCTTCCGGCACCATCGGCATATAAATCGCCACCACATCACCTTTTTTAATGCCCAGCTCGAGCAGGGTATTGGCGAAGCGGCAGACGTCGCGGTGCAGCTCTTTATAGCTGATATGTTTGCTCTGGCTGGCGTCGTCGCCTTCCCAgat                                                                             <  1:1193472/150‑1 (MQ=255)

                                                                               cGGCACCATCGGCATATAAATCGCCACCACATCACCTTTTTTAATGCCCAGCTCGAGCAGGGTATTGGCGAAGCGGCAGACGTCGCGGTGCAGCTCTTTATAGCTGATATGTTTGCTCTGGCTGGCGTCGTCGCCTTCCCAGATGATGGc                                                                       >  1:1532502/1‑150 (MQ=255)

                                                                               cGGCACCATCGGCATATAAATCGCCACCACATCACCTTTTTTAATGCCCAGCTCGAGCAGGGTATTGGCGAAGCGGCAGACGTCGCGGTGCAGCTCTTTATAGCTGATATGTTTGCTCTGGCTGGCGTCGTCGCCTTCCCAGATGATGGc                                                                       <  1:560249/150‑1 (MQ=255)

                                                                               cGGCACCATCGGCATATAAATCGCCACCACATCACCTTTTTTAATGCCCAGCTCGAGCAGGGTATTGGCGAAGCGGCAGACGTCGCGGTGCAGCTCTTTATAGCTGATATGTTTGCTCTGGCTGGCGTCGTCGCCTTCCCAGATGATGGc                                                                       <  1:954672/150‑1 (MQ=255)

                                                                                ggCACCATCGGCATATAAATCGCCACCACATCACCTTTTTTAATGCCCAGCTCGAGCAGGGTATTGGCTAAGCGGCAGACGTCGCGGTGCAGCTCTTTATAGCTGATATGTTTGCTCTGGCTGGCGTCGTCGCCTTCCCAGATGATGGCg                                                                      >  1:3015346/1‑150 (MQ=255)

                                                                                 gCACCATCGGCATAGAAATCGCCACCACATCACCTTTTTTAATGCCCAGCTCGAGCAGGGTATTGGCGAAGCGGCAGACGTCGCGGTGCAGCTCTTTATAGCTGATATGTTTGCTCTGGCTGGCGTCGTCGCCTTCCCAGATGATGGCgg                                                                     >  2:3761409/1‑150 (MQ=255)

                                                                                  cACCATCGGCATATAAATCGCCACCACATCACCTTTTTTAATGCCCAGCTCGAGCAGGGTATTGGCTAAGCGGCAGACGTCGCGGTGCAGCTCTTTATAGCTGATATGTTTGCTCTGGCTGGCGTCGTCGCCTTCCCAGATGATGGCGGt                                                                    <  1:740363/150‑1 (MQ=255)

                                                                                   aCCATCGGCATATAAATCGCCACCACATCACCTTTTTTAATGCCCAGCTCGAGCAGGGTATTGGCTAAGCGGCAGACGTCGCGGTGCAGCTCTTTATAGCTGATATGTTTGCTCTGGCTGGCGTCGTCGCCTTCCCAGATGATGGCGGTa                                                                   <  2:3793275/150‑1 (MQ=255)

                                                                                   aCCATCGGCATATAAATCGCCACCACATCACCTTTTTTAATGCCCAGCTCGAGCAGGGTATTGGCTAAGCGGCAGACGTCGCGGTGCAGCTCTTTATAGCTGATATGTTTGCTCTGGCTGGCGTCGTCGCCTTCCCAGATGATGGCGGTa                                                                   >  1:3660316/1‑150 (MQ=255)

                                                                                   aCCATCGGCATATAAATCGCCACCACATCACCTTTTTTAATGCCCAGCTCGAGCAGGGTATTGGCTAAGCGGCAGACGTCGCGGTGCAGCTCTTTATAGCTGATATGTTTGCTCTGGCTGGCGTCGTCGCCTTCCCAGATGATGGCGGTa                                                                   <  1:390639/150‑1 (MQ=255)

                                                                                     cATCGGCATATAAATCGCCACCACATCACCTTTTTTAATGCCCAGCTCGAGCAGGGTATTGGCGAAGCGGCAGACGTCGCGGTGCAGCTCTTTATAGCTGATATGTTTGCTCTGGCTGGCGTCGTCGCCTTCCCAGATGATGGCGGTACg                                                                 <  2:1845198/150‑1 (MQ=255)

                                                                                     cATCGGCATATAAATCGCCACCACATCACCTTTTTTAATGCCCAGCTCGAGCAGGGTATTGGCGAAGCGGCAGACGTCGCGGTGCAGCTCTGTATAGCTGATATGTTTGCTCTGGCTGGCGTCGTCGCCTTCCCAGATGATGGCGGTACg                                                                 >  2:4350923/1‑150 (MQ=255)

                                                                                      aTCGGCATATAAATCGCCACCACATCACCTTTTTTAATGCCCAGCTCGAGCAGGGTATTGGCTAAGCGGCAGACGTCGCGGTGCAGCTCTTTATAGCTGATATGTTTGCTCTGGCTGGCGTCGTCGCCTTCCCAGATGATGGCGGTACGa                                                                >  2:175308/1‑150 (MQ=255)

                                                                                        cGGCATATAAATCGCCACCACATCACCTTTTTTAATGCCCAGCTCGAGCAGGGTATTGGCTAAGCGGCAGACGTCGCGGTGCAGCTCTTTATAGCTGATATGTTTGCTCTGGCTGGCGTCGTCGCCTTCCCAGATGATGGCGGTACGATc                                                              <  1:2785021/150‑1 (MQ=255)

                                                                                        cGGCATATAAATCGCCACCACATCACCTTTTTTAATGCCCAGCTCGAGCAGGGTATTGGCTAAGCGGCAGACGTCGCGGTGCAGCTCTTTATAGCTGATATGTTTGCTCTGGCTGGCGTCGTCGCCTTCCCAGATGATGGCGGTACGATc                                                              <  2:1993752/150‑1 (MQ=255)

                                                                                         ggCATATAAATCGCCACCACATCACCTTTTTTAATGCCCAGCTCGAGCAGGGTATTGGCTAAGCGGCAGACGTCGCGGTGCAGCTCTTTATAGCTGATATGTTTGCTCTGGCTGGCGTCGTCGCCTTCCCAGATGATGGCGGTACGATCg                                                             <  1:2728572/150‑1 (MQ=255)

                                                                                         ggCATATAAATCGCCACCACATCACCTTTTTTAATGCCCAGCTCGAGCAGGGTATTGGCTAAGCGGCAGACGTCGCGGTGCAGCTCTTTATAGCTGATATGTTTGCTCTGGCTGGCGTCGTCGCCTTCCCAGATGATGGCGGTACGATCg                                                             <  1:3407363/150‑1 (MQ=255)

                                                                                          gCATATAAATCGCCACCACATCACCTTTTTTAATGCCCAGCTCGAGCAGGGTATTGGCTAAGCGGCAGACGTCGCGGTGCAGCTCTTTATAGCTGATATGTTTGCTCTGGCTGGCGTCGTCGCCTTCCCAGATGATGGCGGTACGATCGc                                                            <  1:2592043/150‑1 (MQ=255)

                                                                                             tataAATCGCCACCACATCACCTTTTTTAATGCCCAGCTCGAGCAGGGTATTGGCGAAGCGGCAGACGTCGCGGTGCAGCTCTTTATAGCTGATATGTTTGCTCTGGCTGGCGTCGTCGCCTTCCCAGATGATGGCGGTACGATCGCCGt                                                         <  1:465492/150‑1 (MQ=255)

                                                                                              ataAATCGCCACCACATCACCTTTTTTAATGCCCAGCTCGAGCAGGGTATTGGCTAAGCGGCAGACGTCGCGGTGCAGCTCTTTATAGCTGATATGTTTGCTCTGGCTGGCGTCGTCGCCTTCCCAGATGATGGCGGTACGATCGCCgct                                                        <  2:3434110/150‑3 (MQ=255)

                                                                                                aaaTCGCCACCACATCACCTTTTTTAATGCCCAGCTCGAGCAGGGTATTGGCGAAGCGGCAGACGTCGCGGTGCAGCTCTTTATAGCTGATATGTTTGCTCTGGCTGGCGTCGTCGCCTTCCCAGATGATGGCGGTACGATCGCCGtttt                                                      >  2:927634/1‑150 (MQ=255)

                                                                                                 aaTCGCCACCACATCACCTTTTTTAATGCCCAGCTCGAGCAGGGTATTGGCTAAGCGGCAGACGTCGCGGTGCAGCTCTTTATAGCTGATATGTTTGCTCTGGCTGGCGTCGTCGCCTTCCCAGATGATGGCGGTACGATCGCCGTTTTc                                                     >  2:4470834/1‑150 (MQ=255)

                                                                                                  aTCGCCACCACATCACCTTTTTTAATGCCCAGCTCGAGCAGGGTATTGGCTAAGCGGCAGACGTCGCGGTGCAGCTCTTTATAGCTGATATGTTTGCTCTGGCTGGCGTCGTCGCCTTCCCAGATGATGGCGGTACGATCGCCGTTTTCt                                                    <  1:3977185/150‑1 (MQ=255)

                                                                                                  aTCGCCACCACATCACCTTTTTTAATGCCCAGCTCGAGCAGGGTATTGGCGAAGCGGCAGACGTCGCGGTGCAGCTCTTTATAGCTGATATGTTTGCTCTGGCTGGCGTCGTCGCCTTCCCAGATGATGGCGGTACGATCGCCGTTTTCt                                                    <  2:2566531/150‑1 (MQ=255)

                                                                                                   tCGCCACCACATCACCTTTTTTAATGCCCAGCTCGAGCAGGGTATTGGCTAAGCGGCAGACGTCGCGGTGCAGCTCTTTATAGCTGATATGTTTGCTCTGGCTGGCGTCGTCGCCTTCCCAGATGATGGCGGTACGATCGCCGTTTTCtt                                                   <  2:4253980/150‑1 (MQ=255)

                                                                                                    cGCCACCACATCACCTTTTTTAATGCCCAGCTCGAGCAGGGTATTGGCTAAGCGGCAGACGTCGCGGTGCAGCTCTTTATAGCTGATATGTTTGCTCTGGCTGGCGTCGTCGCCTTCCCAGATGATGGCGGTACGATCGCCGTTTTCTTg                                                  <  1:2855101/150‑1 (MQ=255)

                                                                                                         ccacATCACCTTTTTTAATGCCCAGCTCGAGCAGGGTATTGGCTAAGCGGCAGACGTCGCGGTGCAGCTCTTTATAGCTGATATGTTTGCTCTGGCTGGCGTCGTCGCCTTCCCAGATGATGGCGGTACGATCGCCGTTTTCTTGCAGAt                                             >  2:2815283/1‑150 (MQ=255)

                                                                                                         ccacATCACCTTTTTTAATGCCCAGCTCGAGCAGGGTATTGGCGAAGCGGCAGACGTCGCGGTGCAGCTCTTTATAGCTGATATGTTTGCTCTGGCTGGCGTCGTCGCCTTCCCAGATGATGGCGGTACGATCGCCGTTTTCTTGCAGAt                                             <  1:2621242/150‑1 (MQ=255)

                                                                                                           acaTCACCTTTTTTAATGCCCAGCTCGAGCAGGGTATTGGCTAAGCGGCAGACGTCGCGGTGCAGCTCTTTATAGCTGATATGTTTGCTCTGGCTGGCGTCGTCGCCTTCCCAGATGATGGCGGTACGATCGCCGTTTTCTTGCAGATgg                                           <  1:1777831/150‑1 (MQ=255)

                                                                                                            caTCACCTTTTTTAATGCCCAGCTCGAGCAGGGTATTGGCTAAGCGGCAGACGTCGCGGTGCAGCTCTTTATAGCTGATATGTTTGCTCTGGCTGGCGTCGTCGCCTTCCCAGATGATGGCGGTACGATCGCCGTTTTCTTGCAGATGGc                                          >  1:2020129/1‑150 (MQ=255)

                                                                                                            caTCACCTTTTTTAATGCCCAGCTCGAGCAGGGTATTGGCGAAGCGGCAGACGTCGCGGTGCAGCTCTTTATAGCTGATATGTTTGCTCTGGCTGGCGTCGTCGCCTTCCCAGATGATGGCGGTACGATCGCCGTTTTCTTGCAGATGGc                                          >  1:4237752/1‑150 (MQ=255)

                                                                                                             aTCACCTTTTTTAATGCCCAGCTCGAGCAGGGTATTGGCGAAGCGGCAGGCGTCGCGGTGCAGCTCTTTATAGCTGATATGTTTGCTCTGGCTGGCGTCGTCGCCTTCCCAGATGATGGCGGTACGATCGCCGTTTTCTTGCAGATGGCg                                         >  1:246965/1‑150 (MQ=255)

                                                                                                               cACCTTTTTTAATGCCCAGCTCGAGCAGGGTATTGGCTAAGCGGCAGACGTCGCGGTGCAGCTCTTTATAACTGATATGTTTGCTCTGGCTGGCGTCGTCGCCTTCCCAGATGATGGCGGTACGATCGCCGTTTTCTTGCAGATGGCGGt                                       >  2:1116422/1‑150 (MQ=255)

                                                                                                                aCCTTTTTTAATGCCCAGCTCGAGCAGGGTATTGGCTAAGCGGCAGACGTCGCGGTGCAGCTCTTTATAGCTGATATGTTTGCTCTGGCTGGCGTCGTCGCCTTCCCAGATGATGGCGGTACGATCGCCGTTTTCTTGCAGATGGCGGTc                                      >  2:4501982/1‑150 (MQ=255)

                                                                                                                      tttATTTCCCAGCTCGAGCAGGGTATTGGCTCTGCGGCAGACGTCGCGGTGCAGCTCTTTATAGCTGATATGTTTGCTCTGGCTGGCGTCGTCGCCTTCCCCGATGATGGCGGTACGATCGCCGTTTTCTTGCAGATGTCGGTCAAGGCa                                <  2:233043/150‑1 (MQ=255)

                                                                                                                      tttAATGCCCAGCTCGAGCAGGGTATTGGCGAAGCGGCAGACGTCGCGGTGCAGCTCTTTATAGCTGATATGTTTGCTCTGGCTGGCGTCGTCGCCTTCCCAGATGATGGCGGTACGATCGCCGTTTTCTTGCAGATGGCGGTCAAGGCa                                >  2:3341110/1‑150 (MQ=255)

                                                                                                                          aTGCCCAGCTCGAGCAGGGTCTTGGCTAAGCGGCAGACGTCGCGGTGCAGCTCTTTATAGCTGATATGTTTGCTCTGGCTGGCGTCGTCGCCTTCCCAGATGATGGCGGTACGATCGCCGTTTTCTTGCAGATGGCGGTCAAGGCAGttt                            >  1:4601216/1‑150 (MQ=255)

                                                                                                                          aTGCCCAGCTCGAGCAGGGTATTGGCTAAGCGGCAGACGTCGCGGTGCAGCTCTTTATAGCTGATATGTTTGCTCTGGCTGGCGTCGTCGCCTTCCCAGATGATGGCGGTACGATCGCCGTTTTCTTGCAGATGGCGGTCAAGGCAGttt                            <  2:882355/150‑1 (MQ=255)

                                                                                                                           tGCCCAGCTCGAGCAGGGTATTGGCTAAGCGGCAGACGTCGCGGTGCAGCTCTTTATAGCTGATATGTTTGCTCTGGCTGGCGTCGTCGCCTTCCCAGATGATGGCGGTACGATCGCCGTTTTCTTGCAGATGGCGGTCAAGGCAGTTTg                           >  1:1667310/1‑150 (MQ=255)

                                                                                                                            gCCCAGCTCGAGCAGGGTATTGGCTAAGCGGCAGACGTCGCGGTGCAGCTCTTTATAGCTGATATGTTTGCTCTGGCTGGCGTCGTCGCCTTCCCAGATGATGGCGGTACGATCGCCGTTTTCTTGCAGATGGCGGTCAAGGCAGTTTgc                          <  1:468840/150‑1 (MQ=255)

                                                                                                                                aGCTCGAGCAGGGTATTGGCGAAGCGGCAGACGTCGCGGTGCAGCTCTTTATAGCTGATATGTTTGCTCTGGCTGGCGTCGTCGCCTTCCCAGATGATGGCGGTACGATCGCCGTTTTCTTGCAGATGGCGGTCAAGGCAGTTTgccgcc                      >  1:4457797/1‑150 (MQ=255)

                                                                                                                                  cTCGAGCAGGGTATTGGCTAAGCGGCAGACGTCGCGGTGCAGCTCTTTATAGCTGATATGTTTGCTCTGGCTGGCGTCGTCGCCTTCCCAGATGATGGCGGTACGATCGCCGTTTTCTTGCAGATGGCGGTCAAGGCAGTTTGCCGCCAg                    <  2:1970219/150‑1 (MQ=255)

                                                                                                                                  cTCGAGCAGGGTATTGGCTAAGCGGCAGACGTCGCGGTGCAGCTCTTTATAGCTGATATGTTTGCTCTGGCTGGCGTCGTCGCCTTCCCAGATGATGGCGGTACGATCGCCGTTTTCTTGCAGATGGCGGTCAAGGCAGTTTGCCGCCAg                    >  2:3282923/1‑150 (MQ=255)

                                                                                                                                   tCGAGCAGGGTATTGGCGAAGCGGCAGACGTCGCGGTGCAGCTCTTTATAGCTGATATGTTTGCTCTGGCTGGCGTCGTCGCCTTCCCAGATGATGGCGGTACGATCGCCGTTTTCTTGCAGATGGCGGTCAAGGCAGTTTGCCGCCAGa                   >  1:232764/1‑150 (MQ=255)

                                                                                                                                    cGAGCAGGGTATTGGCGAAGCGGCAGACGTCGCGGTGCAGCTCTTTATAGCTGATATGTTTGCTCTGGCTGGCGTCGTCGCCTTCCCAGATGATGGCGGTACGATCGCCGTTTTCTTGCAGATGGCGGTCAAGGCAGTTTGCCGCCAGAt                  <  1:3357145/150‑1 (MQ=255)

                                                                                                                                    cGAGCAGGGTATTGGCGAAGCGGCAGACGTCGCGGTGCAGCTCTTTATAGCTGATATGTTTGCTCTGGCTGGCGTCGTCGCCTTCCCAGATGATGGCGGTACGATCGCCGTTTTCTTGCAGATGGCGGTCAAGGCAGTTTGCCGCCAGAt                  <  1:4552767/150‑1 (MQ=255)

                                                                                                                                        cAGGGTATTGGCGAAGCGGCAGACGTCGCGGTGCAGCTCTTTATAGCTGATATGTTTGCTCTGGCTGGCGTCGTCGCCTTCCCAGATGATGGCGGTACGATCGCCGTTTTCTTGCAGATGGCGGTCAAGGCAGTTTGCCGCCAGATTCAg              >  2:760388/1‑150 (MQ=255)

                                                                                                                                         aGGGTATTGGCTAAGCGGCAGACGTCGCGGTGCAGCTCTTTATAGCTGATATGTTTGCTCTGGCTGGCGTCGTCGCCTTCCCAGATGATGGCGGTACGATCGCCGTTTTCTTGCAGATGGCGGTCAAGGCAGTTTGCCGCCAGATTCAGc             <  2:131649/150‑1 (MQ=255)

                                                                                                                                         aGGGTATTGGCGAAGCGGCAGACGTCGCGGTGCAGCTCTTTATAGCTGATATGTTTGCTCTGGCTGGCGTCGTCGCCTTCCCAGATGATGGCGGTACGATCGCCGTTTTCTTGCAGATGGCGGTCAAGGCAGTTTGCCGCCAGATTCAGc             >  1:39213/1‑150 (MQ=255)

                                                                                                                                            gTATTGGCTAAGCGGCAGACGTCGCGGTGCAGCTCTTTATAGCTGATATGTTTGCTCTGGCTGGCGTCGTCGCCTTCCCAGATGATGGCGGTACGATCGCCGTTTTCTTGCAGATGGCGGTCAAGGCAGTTTGCCGCCAGATTCAGCGTg          >  2:3989689/1‑150 (MQ=255)

                                                                                                                                              aTTGGCTAAGCGGCAGACGTCGCGGTGCAGCTCTTTATAGCTGATATGTTTGCTCTGGCTGGCGTCGTCGCCTTCCCAGATGATGGCGGTACGATCGCCGTTTTCTTGCAGATGGCGGTCAAGGCAGTTTGCCGCCAGATTCAGCGTGcc        <  1:2428033/150‑1 (MQ=255)

                                                                                                                                              aTTGGCGAAGCGGCAGACGTCGCGGTGCAGCTCTTTATAGCTGATATGTTTGCTCTGGCTGGCGTCGTCGCCTTCCCAGATGATGGCGGTACGATCGCCGTTTTCTTGCAGATGGCGGTCAAGGCAGTTTGCCGCCAGATTCAGCGTGcc        >  1:863521/1‑150 (MQ=255)

                                                                                                                                              aTTGGCGAAGCGGCAGACGTCGCGGTGCAGCTCTTTATAGCTGATATGTTTGCTCTGGCTGGCGTCGTCGCCTTCCCAGATGATGGCGGTACGATCGCCGTTTTCTTGCAGATGGCGGTCAAGGCAGTTTGCCGCCAGATTCAGCGTGcc        <  2:1392424/150‑1 (MQ=255)

                                                                                                                                                    gAAGCGGCAGACGTCGCGGTGCAGCTCTTTATAGCTGATATGTTTGCTCTGGCTGGCGTCGTCGCCTTCCCAGATGATGGCGGTACGATCGCCGTTTTCTTGCAGATGGCGGTCAAGGCAGTTTGCCGCCAGATTCAGCGTGCCGTCCTc  <  1:1017263/150‑1 (MQ=255)

                                                                                                                                                    gAAGCGGCAGACGTCGCGGTGCAGCTCTTTATAGCTGATATGTTTGCTCTGGCTGGCGTCGTCGCCTTCCCAGATGATGGCGGTACGATCGCCGTTTTCTTGCAGATGGCGGTCAAGGCAGTTTGCCGCCAGATTCAGCGTGCCGTCCTc  >  1:3720862/1‑150 (MQ=255)

CTTCCGGCGAGAAGCCGCCGAAAATCACCGAATGCACCGCGCCAATGCGGGCGCAGGCCAGCATCGCAACCGCGGCTTCCGGCACCATCGGCATATAAATCGCCACCACATCACCTTTTTTAATGCCCAGCTCGAGCAGGGTATTGGCGAAGCGGCAGACGTCGCGGTGCAGCTCTTTATAGCTGATATGTTTGCTCTGGCTGGCGTCGTCGCCTTCCCAGATGATGGCGGTACGATCGCCGTTTTCTTGCAGATGGCGGTCAAGGCAGTTTGCCGCCAGATTCAGCGTGCCGTCCTC  >  CP009273/4276792‑4277089

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