BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	4,357,095	A→C	49.5%	V350V (GTT→GTG)	aspA ←	aspartate ammonia‑lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	4,357,095	0	A	C	49.5%	102.7 / inf	267	V350V (GTT→GTG)	aspA	aspartate ammonia‑lyase
Reads supporting (aligned to +/- strand): ref base A (92/43); new base C (64/68); total (156/111)
Fisher's exact test for biased strand distribution p-value = 1.26e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.05e-01

Only 100 of 275 total aligned reads displayed.

ACTTCTTTGTTAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGCCTGCAGTTCCGGCAGGTTGATCTCGTTCAGGCCGGCACGTGGGCCTGAAGAGAG  >  CP009273/4356946‑4357244

acTTCTTTGTTAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTa                                                                                                                                                       >  2:4123079/1‑150 (MQ=255)

 cttcttTGTTAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTaa                                                                                                                                                      <  2:2205297/150‑1 (MQ=255)

 cttcttTGTTAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTaa                                                                                                                                                      >  2:1122431/1‑150 (MQ=255)

  ttcttTGTTAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCAc                                                                                                                                                     >  2:3217517/1‑150 (MQ=255)

  ctcttTGTTAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAAc                                                                                                                                                     >  1:2181473/2‑150 (MQ=255)

    cttTGTTAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAg                                                                                                                                                   >  1:3973306/1‑150 (MQ=255)

         ttAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAgtggtg                                                                                                                                              >  2:2874450/1‑150 (MQ=255)

          tAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGgtgt                                                                                                                                             >  2:2594951/1‑150 (MQ=255)

           aGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTc                                                                                                                                            >  2:1468593/1‑150 (MQ=255)

               gTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGTCCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTa                                                                                                                                        <  1:2920849/150‑1 (MQ=255)

                tGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTAc                                                                                                                                       >  2:1096370/1‑150 (MQ=255)

                tGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTAc                                                                                                                                       >  1:1232067/1‑150 (MQ=255)

                tGATGCCGTTAATGCA‑TTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTAcc                                                                                                                                      >  1:1696488/1‑150 (MQ=255)

                 gATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTAcc                                                                                                                                      >  1:1933312/1‑150 (MQ=255)

                  aTGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCg                                                                                                                                     <  2:2858193/150‑1 (MQ=255)

                    gCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGAt                                                                                                                                   >  2:1594262/1‑150 (MQ=255)

                    gCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGAt                                                                                                                                   >  1:50261/1‑150 (MQ=255)

                    gCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGAt                                                                                                                                   >  2:4164386/1‑150 (MQ=255)

                     ccGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATg                                                                                                                                  <  2:791349/150‑1 (MQ=255)

                     ccGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATg                                                                                                                                  <  1:2347867/150‑1 (MQ=255)

                      cGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGa                                                                                                                                 >  2:3797944/1‑150 (MQ=255)

                       gTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGAc                                                                                                                                >  1:3159023/1‑150 (MQ=255)

                       cttAATGCATTTTTCCAGCAGGTTGAAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCAAGGTAACAGAGGTGTCGATACCCATGAc                                                                                                                                >  2:3542674/2‑150 (MQ=255)

                         tAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACtt                                                                                                                              >  1:3071602/1‑150 (MQ=255)

                           aTGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGACTTTg                                                                                                                            >  2:3451525/1‑150 (MQ=255)

                           aTGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGACTTTg                                                                                                                            <  2:2581825/150‑1 (MQ=255)

                            tGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGACTTTGa                                                                                                                           <  1:3582674/150‑1 (MQ=255)

                             gCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGACTTTGaa                                                                                                                          <  1:514558/150‑1 (MQ=255)

                              cATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAg                                                                                                                         >  1:688957/1‑150 (MQ=255)

                               aTTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGc                                                                                                                        >  1:3760381/1‑150 (MQ=255)

                                tttttCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCa                                                                                                                       >  1:2451982/1‑150 (MQ=255)

                                     ccAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCt                                                                                                                  <  1:2814525/150‑1 (MQ=255)

                                         cagGTTTTAGATAGTGTTGTTCTGAATGTGAACGTATTCTATCTTGGCCTTGCGAATGACGGTCTACATAACTTTAAACTTCATCTGAGCTGATTCTGCTTACATGGTCACAGTGGTGTCGTTACCGATGAATTTGAAGCATACATTGtt                                                                                                              <  2:4043129/150‑1 (MQ=255)

                                               gTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCAc                                                                                                        >  2:1350264/1‑150 (MQ=255)

                                                tAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACt                                                                                                       >  1:434753/1‑150 (MQ=255)

                                                   cAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTcc                                                                                                    >  2:3084221/1‑150 (MQ=255)

                                                    aaGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCg                                                                                                   <  2:481117/150‑1 (MQ=255)

                                                    aaGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCg                                                                                                   >  1:2026099/1‑150 (MQ=255)

                                                      gCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGa                                                                                                 <  2:3690388/150‑1 (MQ=255)

                                                        gTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGaac                                                                                               >  2:1244838/1‑150 (MQ=255)

                                                          tGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGaacat                                                                                             <  2:1040578/150‑2 (MQ=255)

                                                          tGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGaacaa                                                                                             <  1:2298951/150‑1 (MQ=255)

                                                           ggTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGaacaac                                                                                            >  1:1585498/1‑150 (MQ=255)

                                                            gTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAAcc                                                                                           >  2:3552992/1‑150 (MQ=255)

                                                            gTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAAcc                                                                                           >  2:1997014/1‑150 (MQ=255)

                                                              cAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCgg                                                                                         <  2:4006838/150‑1 (MQ=255)

                                                               aGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCggg                                                                                        >  1:3982065/1‑150 (MQ=255)

                                                               aGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCggg                                                                                        <  1:4119706/150‑1 (MQ=255)

                                                                 aaTGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGtt                                                                                      >  1:1352663/1‑150 (MQ=255)

                                                                    gtgAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTaa                                                                                   <  2:486529/150‑2 (MQ=255)

                                                                     tgAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACt                                                                                  <  1:2236603/150‑1 (MQ=255)

                                                                     ngAACGGANTCGAACANGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGNAGCATACCTGGTTAACCACTTCCGGAANAACCGGGTTTACt                                                                                  >  2:732697/2‑150 (MQ=255)

                                                                      gAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACtt                                                                                 <  1:1079905/150‑1 (MQ=255)

                                                                        aCGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTa                                                                               >  2:2395009/1‑150 (MQ=255)

                                                                           gATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCt                                                                            <  1:808220/150‑1 (MQ=255)

                                                                             ttCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTgg                                                                          >  1:700348/1‑150 (MQ=255)

                                                                              tCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGc                                                                         <  2:3178383/150‑1 (MQ=255)

                                                                                gAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCat                                                                       >  2:1761242/1‑150 (MQ=255)

                                                                                 aaCATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCatg                                                                      >  2:3977896/1‑150 (MQ=255)

                                                                                 aaCATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCatg                                                                      >  2:2426419/1‑150 (MQ=255)

                                                                                  aCATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGACTTTGAAGAATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCatga                                                                     >  1:3501033/1‑150 (MQ=255)

                                                                                   cATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCatgat                                                                    <  2:2767937/150‑1 (MQ=255)

                                                                                    aTGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCatgatg                                                                   >  2:1978493/1‑150 (MQ=255)

                                                                                    aTGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCatgatg                                                                   >  2:1111346/1‑150 (MQ=255)

                                                                                        cctggccAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAag                                                               <  2:2814637/150‑1 (MQ=255)

                                                                                         ctggccAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAaga                                                              <  2:3319997/150‑1 (MQ=255)

                                                                                           ggccAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGCAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGc                                                            <  2:2106613/150‑1 (MQ=255)

                                                                                              cAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCg                                                         >  2:3800590/1‑150 (MQ=255)

                                                                                                aTGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGcc                                                       >  1:303488/1‑150 (MQ=255)

                                                                                                aTGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGcc                                                       <  2:428614/150‑1 (MQ=255)

                                                                                                aTGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGcc                                                       <  2:207755/150‑1 (MQ=255)

                                                                                                aTGACCGGCTCCATAAAGTTCAACTGCAGCTTACCTGCATCTGCTGCCATCGTCAAAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACATCCGGAACAACCGGGTTAACTTTAGCAGGCATGATGGAACAGCCCGcc                                                       >  1:2817387/1‑150 (MQ=255)

                                                                                                 tGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGCCt                                                      <  2:3877325/150‑1 (MQ=255)

                                                                                                  gACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGCCTg                                                     >  1:3388048/1‑150 (MQ=255)

                                                                                                     cGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGCCTGCAg                                                  >  1:853918/1‑150 (MQ=255)

                                                                                                      ggCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGCCTGCAGt                                                 >  2:1803025/1‑150 (MQ=255)

                                                                                                      ggCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGCCTGCAGt                                                 <  1:1634546/150‑1 (MQ=255)

                                                                                                         tcaaTAACGTTCATCTGCAGCAGAACTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACGTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGCCTGCAGTTcc                                              <  2:1640112/147‑1 (MQ=255)

                                                                                                           cATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGCCTGCAGTTCCgg                                            >  2:151858/1‑150 (MQ=255)

                                                                                                              aaCGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGCCTGCAGTTCCGGCAg                                         >  1:3131251/1‑150 (MQ=255)

                                                                                                                  ttCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGCCTGCAGTTCCGGCAGGTTg                                     <  1:1997030/150‑1 (MQ=255)

                                                                                                                  ttCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGCCTGCAGTTCCGGCAGGTTg                                     <  1:1411524/150‑1 (MQ=255)

                                                                                                                   tCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGCCTGCAGTTCCGGCAGGTTGa                                    <  1:962435/150‑1 (MQ=255)

                                                                                                                        tGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGCCTGCAGTTCCGGCAGGTTGATCTCg                               >  1:1247873/1‑150 (MQ=255)

                                                                                                                        tGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGCCTGCAGTTCCGGCAGGTTGATCTCg                               >  1:3578751/1‑150 (MQ=255)

                                                                                                                          cAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGCCTGCAGTTCCGGCAGGTTGATCTCGtt                             >  1:1655832/1‑150 (MQ=255)

                                                                                                                              tGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGCCTGCAGTTCCGGCAGGTTGATCTCGTTCAgg                         >  1:1782939/1‑150 (MQ=255)

                                                                                                                                 ccTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGCCTGCAGTTCCGGCAGGTTGATCTCGTTCAGGCCg                      >  2:612652/1‑150 (MQ=255)

                                                                                                                                  cTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGCCTGCAGTTCCGGCAGGTTGATCTCGTTCAGGCCgg                     >  1:1363850/1‑150 (MQ=255)

                                                                                                                                     cTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGCCTGCAGTTCCGGCAGGTTGATCTCGTTCAGGCCGGCAc                  <  1:3224665/150‑1 (MQ=255)

                                                                                                                                         tgctgcCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCAGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGCCTGCAGTTCCGGCAGGTTGATCTCGTTCAGGCCGGCACGTgg              >  2:1844854/1‑150 (MQ=255)

                                                                                                                                          gctgcCATGGTAACAGTGATGTCGGTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGCCTGCAGTTCCGGCAGGTTGATCTCGTTCAGGCCGGCACGTggg             <  2:2638624/150‑1 (MQ=255)

                                                                                                                                             gcCATGGTCACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGCCTGCAGTTCCGGCAGGTTGATCTCGTTCAGGCCGGCACGTGGGCCt          >  2:2689016/1‑150 (MQ=255)

                                                                                                                                             gcCATGGTCACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGCCTGCAGTTCCGGCAGGTTGATCTCGTTCAGGCCGGCACGTGGGCCt          <  2:2352791/150‑1 (MQ=255)

                                                                                                                                               cATGGTCACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGCCTGCAGTTCCGGCAGGTTGATCTCGTTCAGGCCGGCACGTGGGCCTGa        >  2:563604/1‑150 (MQ=255)

                                                                                                                                                 tGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGCCTGCAGTTCCGGCAGGTTGATCTCGTTCAGGCCGGCACGTGGGCCTGAag      <  1:255777/150‑1 (MQ=255)

                                                                                                                                                 tGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGCCTGCAGTTCCGGCAGGTTGATCTCGTTCAGGCCGGCACGTGGGCCTGAag      <  1:255808/150‑1 (MQ=255)

                                                                                                                                                 tGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGCCTGCAGTTCCGGCAGGTTGATCTCGTTCAGGCCGGCACGTGGGCCTGAag      <  1:3652337/150‑1 (MQ=255)

                                                                                                                                                  ggTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGCCTGCAGTTCCGGCAGGTTGATCTCGTTCAGGCCGGCACGTGGGCCTGAaga     <  2:3908094/150‑1 (MQ=255)

                                                                                                                                                     aaCAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGCCTGCAGTTCCGGCAGGTTGATCTCGTTCAGGCCGGCACGTGGGCCTGAagagag  <  1:2110526/150‑1 (MQ=255)

ACTTCTTTGTTAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGCCTGCAGTTCCGGCAGGTTGATCTCGTTCAGGCCGGCACGTGGGCCTGAAGAGAG  >  CP009273/4356946‑4357244

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