BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	4,357,095	A→C	47.2%	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	47.2%	79.8 / inf	310	V350V (GTT→GTG)	aspA	aspartate ammonia‑lyase
Reads supporting (aligned to +/- strand): ref base A (88/76); new base C (73/73); total (161/149)
Fisher's exact test for biased strand distribution p-value = 5.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.60e-01

Only 100 of 313 total aligned reads displayed.

ACTTCTTTGTTAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGCCTGCAGTTCCGGCAGGTTGATCTCGTTCAGGCCGGCACGTGGGCCTGAAGAGAG  >  CP009273/4356946‑4357244

  ttcttTGTTAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAAc                                                                                                                                                     <  1:2388687/150‑1 (MQ=255)

  ttcttTGTTAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAAc                                                                                                                                                     >  1:1536083/1‑150 (MQ=255)

   tcttTGTTAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTcaca                                                                                                                                                    <  1:3006905/150‑1 (MQ=255)

    cttTGTTAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAg                                                                                                                                                   >  1:3147192/1‑150 (MQ=255)

      ttGTTAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAgtg                                                                                                                                                 <  1:1946374/150‑1 (MQ=255)

       tGTTAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAgtgg                                                                                                                                                >  1:1259143/1‑150 (MQ=255)

         ttAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAgtggtg                                                                                                                                              <  1:4120584/150‑1 (MQ=255)

         ttAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAgtggtg                                                                                                                                              <  2:1566447/150‑1 (MQ=255)

          tAGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGgtgt                                                                                                                                             >  1:2394891/1‑150 (MQ=255)

           aGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATACCGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTc                                                                                                                                            <  1:3990254/150‑1 (MQ=255)

           aGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTc                                                                                                                                            >  2:4365360/1‑150 (MQ=255)

           aGCAGTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTc                                                                                                                                            <  2:1093781/150‑1 (MQ=255)

               gTGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTa                                                                                                                                        <  2:3735166/150‑1 (MQ=255)

                tGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTAc                                                                                                                                       <  1:953726/150‑1 (MQ=255)

                tGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTAc                                                                                                                                       <  1:3680447/150‑1 (MQ=255)

                tGATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATATGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTAc                                                                                                                                       >  1:1174444/1‑150 (MQ=255)

                 gATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTAcc                                                                                                                                      <  1:3411687/150‑1 (MQ=255)

                 gATGCCGTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTAcc                                                                                                                                      <  1:3381370/150‑1 (MQ=255)

                       gTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGAc                                                                                                                                <  1:2159929/150‑1 (MQ=255)

                       gTTAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGAc                                                                                                                                >  2:1410662/1‑150 (MQ=255)

                         tAATGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGACtt                                                                                                                              >  2:3195854/1‑150 (MQ=255)

                          aaTGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGACttt                                                                                                                             <  2:4116582/150‑1 (MQ=255)

                          aaTGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGACttt                                                                                                                             <  2:4178016/150‑1 (MQ=255)

                          aaTGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGACttt                                                                                                                             <  2:4382994/150‑1 (MQ=255)

                          aaTGCATTTTTCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACttt                                                                                                                             <  2:1690940/150‑1 (MQ=255)

                                 ttttCCAGCATGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGACTTTGAAGCAt                                                                                                                      <  2:4537149/150‑1 (MQ=255)

                                    tCCAGCAGGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGACTTTGAAGCATAcc                                                                                                                   >  1:1886296/1‑150 (MQ=255)

                                       agcagGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTgg                                                                                                                >  2:704023/1‑150 (MQ=255)

                                         cagGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGtt                                                                                                              <  1:1042868/150‑1 (MQ=255)

                                         cagGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGtt                                                                                                              <  1:2332442/150‑1 (MQ=255)

                                          agGTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTa                                                                                                             >  2:3038770/1‑150 (MQ=255)

                                           ggTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTaa                                                                                                            >  2:1408783/1‑150 (MQ=255)

                                           ggTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTaa                                                                                                            >  1:392281/1‑150 (MQ=255)

                                            gTTGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAAc                                                                                                           >  1:4586891/1‑150 (MQ=255)

                                             ttGTAGCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAAcc                                                                                                          <  2:2481854/150‑1 (MQ=255)

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

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

                                                  gCAAGCGTTGGTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTc                                                                                                     >  2:3341489/1‑150 (MQ=255)

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

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

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

                                                           ggTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGaacaac                                                                                            <  2:2394528/150‑1 (MQ=255)

                                                           ggTCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGaacaac                                                                                            >  1:1690115/1‑150 (MQ=255)

                                                             tCAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCg                                                                                          <  2:484965/150‑1 (MQ=255)

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

                                                              cAGAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCgg                                                                                         >  2:4046421/1‑150 (MQ=255)

                                                                gAATGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGt                                                                                       >  1:2783535/1‑150 (MQ=255)

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

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

                                                                  aTGTGAACGGATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGttt                                                                                     <  1:3176038/150‑1 (MQ=255)

                                                                          ggATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGc                                                                             <  2:813494/150‑1 (MQ=255)

                                                                          ggATTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGc                                                                             >  1:3556277/1‑150 (MQ=255)

                                                                            aTTCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTg                                                                           >  1:967361/1‑150 (MQ=255)

                                                                             ttCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTgg                                                                          <  1:2555881/150‑1 (MQ=255)

                                                                             ttCGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTgg                                                                          >  1:3365476/1‑150 (MQ=255)

                                                                               cGAACATGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCa                                                                        >  2:2069911/1‑150 (MQ=255)

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

                                                                                    aTGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCatgatg                                                                   <  2:3866480/150‑1 (MQ=255)

                                                                                    aTGGCCTGGCCAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCatgatg                                                                   <  2:2522577/150‑1 (MQ=255)

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

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

                                                                                     tggcctggccAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATgg                                                                  <  2:1927281/150‑1 (MQ=255)

                                                                                     tggcctggccAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATgg                                                                  >  2:340168/1‑150 (MQ=255)

                                                                                     tggcctggccAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATgg                                                                  >  2:2998618/1‑150 (MQ=255)

                                                                                       gcctggccAATGACCGGCTCCATAACGTTCAACTGCAGCGGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGaa                                                                <  1:668720/150‑1 (MQ=255)

                                                                                        cctggccAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAag                                                               <  1:1378758/150‑1 (MQ=255)

                                                                                         ctggccAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAaga                                                              >  1:2899319/1‑150 (MQ=255)

                                                                                          tggccAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGTTGGAagag                                                             <  1:1891514/150‑1 (MQ=255)

                                                                                           ggccAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGc                                                            >  2:1436025/1‑150 (MQ=255)

                                                                                           ggccAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGCCATTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGc                                                            <  1:4162838/150‑1 (MQ=255)

                                                                                            gccAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGcc                                                           <  2:2600453/150‑1 (MQ=255)

                                                                                             ccAATGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGccc                                                          <  1:4365433/150‑1 (MQ=255)

                                                                                               aaTGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCga                                                        <  2:1526406/150‑2 (MQ=255)

                                                                                                 tGACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGCCt                                                      >  1:4356903/1‑150 (MQ=255)

                                                                                                  gACCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACGTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGCCTg                                                     <  2:1499091/150‑1 (MQ=255)

                                                                                                   aCCGGCTCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGCCTGc                                                    <  1:1778464/150‑1 (MQ=255)

                                                                                                         tCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGCCTGCAGTTcc                                              >  1:1062742/1‑150 (MQ=255)

                                                                                                         tCCATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGCCTGCAGTTcc                                              <  2:1862768/150‑1 (MQ=255)

                                                                                                          ccATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGCCTGCAGTTCCg                                             <  1:861706/150‑1 (MQ=255)

                                                                                                          ccATAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGCCTGCAGTTCCg                                             <  2:1957285/150‑1 (MQ=255)

                                                                                                             tAACGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGCCTGCAGTTCCGGCa                                          <  1:2465545/150‑1 (MQ=255)

                                                                                                              aaCGTTCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGCCTGCAGTTCCGGCAg                                         >  2:4326543/1‑150 (MQ=255)

                                                                                                                   tCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGCCTGCAGTTCCGGCAGGTTGa                                    <  1:4222484/150‑1 (MQ=255)

                                                                                                                   tCAACTGCAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGCCTGCAGTTCCGGCAGGTTGa                                    <  1:2292765/150‑1 (MQ=255)

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

                                                                                                                     aaCTGCAGCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGCCTGCAGTTCCGGCAGGTTGAtc                                  >  2:2146088/1‑150 (MQ=255)

                                                                                                                          cAGCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGCCTGCAGTTCCGGCAGGTTGATCTCGtt                             >  2:2585081/1‑150 (MQ=255)

                                                                                                                           aNCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGCCTGCAGTTCCGGCAGGTTGATCTCGTTc                            >  2:3254732/1‑150 (MQ=255)

                                                                                                                            gCTGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGCCTGCAGTTCCGGCAGGTTGATCTCGTTCa                           <  1:2408034/150‑1 (MQ=255)

                                                                                                                            gCTGACCTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGCCTGCAGTTCCGGCAGGTTGATCTCGTTCa                           <  2:2323498/150‑1 (MQ=255)

                                                                                                                              tGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCTGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGCCTGCAGTTCCGGCAGGTTGATCTCGTTCAgg                         >  1:4035806/1‑150 (MQ=255)

                                                                                                                              tGACCTGCTTCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGCCTGCAGTTCCGGCAGGTTGATCTCGTTCAgg                         <  2:130449/150‑1 (MQ=255)

                                                                                                                                 ccTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGCCTGCAGTTCCGGCAGGTTGATCTCGTTCAGGCCg                      >  2:1193189/1‑150 (MQ=255)

                                                                                                                                 ccTGCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGCCTGCAGTTCCGGCAGGTTGATCTCGTTCAGGCCg                      >  1:576639/1‑150 (MQ=255)

                                                                                                                                    gCTTCTGCTGCCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGCCTGCAGTTCCGGCAGGTTGATCTCGTTCAGGCCGGCa                   >  1:3026857/1‑150 (MQ=255)

                                                                                                                                      ttCTGCTGCCATGGTCACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGCCTGCAGTTCCGGCAGGTTGATCTCGTTCAGGCCGGCACg                 <  1:4142903/150‑1 (MQ=255)

                                                                                                                                        ctgctgCCATGGTCACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGCCTGCAGTTCCGGCAGGTTGATCTCGTTCAGGCCGGCACGTg               <  1:1205176/150‑1 (MQ=255)

                                                                                                                                            tgcCATGGTAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGCCGGCAGTACCGGAAGTTGGTACTCGTAAAGCCCGCCAGGTggcca           >  1:3760733/1‑147 (MQ=255)

                                                                                                                                                    tAACAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGCCTGCAGTTCCGGCAGGTTGATCTCGTTCAGGCCGGCACGTGGGCCTGAagaga   >  2:2702626/1‑150 (MQ=255)

                                                                                                                                                     aaCAGTGGTGTCGTTACCGATGACTTTGAAGCATACCTGGTTAACCACTTCCGGAACAACCGGGTTTACTTTAGCTGGCATGATGGAAGAGCCCGCCTGCAGTTCCGGCAGGTTGATCTCGTTCAGGCCGGCACGTGGGCCTGAagagag  <  1:2617837/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: `‑`