BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	853,320	A→T	7.5%	W400R (TGG→AGG)	ybiU ←	DUF1479 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	853,320	0	A	T	7.5%	1031.7 / 53.8	293	W400R (TGG→AGG)	ybiU	DUF1479 family protein
Reads supporting (aligned to +/- strand): ref base A (137/134); new base T (11/11); total (148/145)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

Only 100 of 298 total aligned reads displayed.

TTAACAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCTTTTTCCAGTGCGGCCTTCACCTTGTGCGCGTAGGCAAGATTTTTCTCGCACATCGGTGCCGCAGGAATGTACATCACGTTGCCCCAACCTTGTTGATTTTCAA  >  CP009273/853171‑853469

ttAACAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCa                                                                                                                                                       >  1:4721987/1‑150 (MQ=255)

     aaaGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAgtttg                                                                                                                                                  >  1:1285088/1‑150 (MQ=255)

       aGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAgtttgtt                                                                                                                                                >  1:1829272/1‑150 (MQ=255)

        gcgcATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAgtttgttt                                                                                                                                               <  1:4840313/150‑1 (MQ=255)

         cgcATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTc                                                                                                                                              <  1:1646166/150‑1 (MQ=255)

           cATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCAt                                                                                                                                            >  2:1525365/1‑150 (MQ=255)

             tttGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCCATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAg                                                                                                                                          >  2:4427270/1‑150 (MQ=255)

             tttGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAg                                                                                                                                          <  1:4410525/150‑1 (MQ=255)

              ttGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGt                                                                                                                                         >  1:794764/1‑150 (MQ=255)

              ttGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAAACCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGt                                                                                                                                         <  1:1952101/150‑1 (MQ=255)

                 tCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCt                                                                                                                                      >  2:3330844/1‑150 (MQ=255)

                  cAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTc                                                                                                                                     <  1:152392/150‑1 (MQ=255)

                   aGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTcg                                                                                                                                    >  1:2103182/1‑150 (MQ=255)

                   aGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTcg                                                                                                                                    >  1:2379433/1‑150 (MQ=255)

                        tCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCgg                                                                                                                               >  1:3810188/1‑150 (MQ=255)

                         cTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCggg                                                                                                                              >  2:1891681/1‑150 (MQ=255)

                           gAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGaa                                                                                                                            <  1:2292972/150‑1 (MQ=255)

                           gAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGaa                                                                                                                            <  1:2293274/150‑1 (MQ=255)

                            aGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAg                                                                                                                           <  2:623799/150‑1 (MQ=255)

                              tAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTc                                                                                                                         <  1:1809811/150‑1 (MQ=255)

                                 cGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGcc                                                                                                                      >  2:2987329/1‑150 (MQ=255)

                                 cGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGcc                                                                                                                      <  1:1009870/150‑1 (MQ=255)

                                  gTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGccc                                                                                                                     >  1:1292997/1‑150 (MQ=255)

                                    cGAATCGATTTTCTTTATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCTACGCTCGCTTACCGTGAATGTTGAGGTCGGAAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGTAAGTCGCCCgg                                                                                                                   <  1:3093861/150‑1 (MQ=255)

                                    cGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAAAATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCTGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCgg                                                                                                                   >  1:1335472/1‑150 (MQ=255)

                                      aaTCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCg                                                                                                                 >  2:4857857/1‑150 (MQ=255)

                                          gATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCTGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATgc                                                                                                             >  1:435755/1‑150 (MQ=255)

                                           aTTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATgcg                                                                                                            <  2:2586028/150‑1 (MQ=255)

                                              ttCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCt                                                                                                         <  2:4814791/150‑1 (MQ=255)

                                                cTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCttt                                                                                                       <  2:642152/150‑1 (MQ=255)

                                                  tAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCttttt                                                                                                     >  2:2847134/1‑150 (MQ=255)

                                                   aaTCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCTTTTTc                                                                                                    >  1:3951521/1‑150 (MQ=255)

                                                          caTCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCTTTTTCCAGTGCg                                                                                             >  2:3843612/1‑150 (MQ=255)

                                                           aTCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCTTTTTCCAGTGCgg                                                                                            <  2:738833/150‑1 (MQ=255)

                                                            tCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCTTTTTCCAGTGCGGc                                                                                           <  2:418280/150‑1 (MQ=255)

                                                              aaTTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCTTTTTCCAGTGCGGCCt                                                                                         >  2:2677126/1‑150 (MQ=255)

                                                              aaTTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCTTTTTCCAGTGCGGCCt                                                                                         <  2:904300/150‑1 (MQ=255)

                                                              aaTTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGCGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCTGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCTTTTTCCAGTGCGGCCt                                                                                         >  2:107637/1‑150 (MQ=255)

                                                               aTTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCTTTTTCCAGTGCGGCCtt                                                                                        <  1:4041949/150‑1 (MQ=255)

                                                                  cACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCTTTTTCCAGTGCGGCCTTCAc                                                                                     >  2:2507440/1‑150 (MQ=255)

                                                                  cACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCTTTTTCCAGTGCGGCCTTCAc                                                                                     >  1:2989583/1‑150 (MQ=255)

                                                                   aCTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCTTTTTCGAGTGCGGCCTTCAcc                                                                                    <  1:3271472/150‑1 (MQ=255)

                                                                       ccATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCTTTTTCCAGTGCGGCCTTCACCTtgt                                                                                >  1:2848553/1‑150 (MQ=255)

                                                                       ccATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCTTTTTCCAGTGCGGCCTTCACCTtgt                                                                                >  2:3953972/1‑150 (MQ=255)

                                                                         aTTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCTTTTTCCAGTGCGGCCTTCACCTTGTgc                                                                              >  1:4197233/1‑150 (MQ=255)

                                                                          ttAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCTTTTTCCAGTGCGGCCTTCACCTTGTgcg                                                                             >  1:1029041/1‑150 (MQ=255)

                                                                            aaTTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCTGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCTTTTTCCAGTGCGGCCTTCACCTTGTGcgcg                                                                           >  1:2912551/1‑150 (MQ=255)

                                                                            aaTTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCTTTTTCCAGTGCGGCCTTCACCTTGTGcgcg                                                                           >  2:2728402/1‑150 (MQ=255)

                                                                             aTTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCTTTTTCCAGTGCGGCCTTCACCTTGTGCGCGt                                                                          >  2:2906436/1‑150 (MQ=255)

                                                                              ttCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCTTTTTCCAGTGCGGCCTTCACCTTGTGCGCGTa                                                                         <  1:36261/150‑1 (MQ=255)

                                                                                  cAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCTTTTTCCAGTGCGGCCTTCACCTTGTGCGCGTAGGCa                                                                     >  2:2401572/1‑150 (MQ=255)

                                                                                  cAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCTTTTTCCAGTGCGGCCTTCACCTTGTGCGCGTAGGCa                                                                     >  1:2714267/1‑150 (MQ=255)

                                                                                   aaaCATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCTTTTTCCAGTGCGGCCTTCACCTTGTGCGCGTAGGCaa                                                                    >  1:3413568/1‑150 (MQ=255)

                                                                                   aaaCATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCTTTTTCCAGTGCGGCCTTCACCTTGTGCGCGTAGGCaa                                                                    <  2:2092014/150‑1 (MQ=255)

                                                                                     aCATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCTTTTTCCAGTGCGGCCTTCACCTTGTGCGCGTAGGCAAGa                                                                  >  1:1165245/1‑150 (MQ=255)

                                                                                      cATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCTTTTTCCAGTGCGGCCTTCACCTTGTGCGCGTAGGCACGAt                                                                 >  1:4686713/1‑150 (MQ=255)

                                                                                       aTCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCTTTTTCCAGTGCGGCCTTCACCTTGTGCGCGTAGGCAAGAtt                                                                >  1:3011076/1‑150 (MQ=255)

                                                                                        tCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCTTTTTCCAGTGCGGCCTTCACCTTGTGCGCGTAGGCAAGAttt                                                               >  2:1125555/1‑150 (MQ=255)

                                                                                        tCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCTTTTTCCAGTGCGGCCTTCACCTTGTGCGCGTAGGCAAGAttt                                                               >  2:4787063/1‑150 (MQ=255)

                                                                                        tCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCTTTTTCCAGTGCGGCCTTCACCTTGTGCGCGTAGGCAAGAttt                                                               >  2:1125445/1‑150 (MQ=255)

                                                                                        tCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCTTTTTCCAGTGCGGCCTTCACCTTGTGCGCGTAGGCAAGAttt                                                               <  2:3258283/150‑1 (MQ=255)

                                                                                         ccATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCTTTTTCCAGTGCGGCCTTCACCTTGTGCGCGTAGGCAAGAtttt                                                              <  2:1455348/150‑1 (MQ=255)

                                                                                           aTGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCTTTTTCCAGTGCGGCCTTCACCTTGTGCGCGTAGGCAAGATTTTtc                                                            <  1:1027402/150‑1 (MQ=255)

                                                                                           aTGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCTTTTTCCAGTGCGGCCTTCACCTTGTGCGCGTAGGCAAGATTTTtc                                                            >  1:3541377/1‑150 (MQ=255)

                                                                                             gCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCTTTTTCCAGTGCGGCCTTCACCTTGTGCGCGTAGGCAAGATTTTtctc                                                          <  1:1466869/150‑1 (MQ=255)

                                                                                                cAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCTTTTTCCAGTGCGGCCTTCACCTTGTGCGCGTAGGCAAGATTTTTCTCGca                                                       <  2:2323613/150‑1 (MQ=255)

                                                                                                cAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCTTTTTCCAGTGCGGCCTTCACCTTGTGCGCGTAGGCAAGATTTTTCTCGca                                                       <  2:3450100/150‑1 (MQ=255)

                                                                                                cAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCTTTTTCCAGTGCGGCCTTCACCTTGTGCGCGTAGGCAAGATTTTTCTCGca                                                       <  2:3449275/150‑1 (MQ=255)

                                                                                                cAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCTTTTTCCAGTGCGGCCTTCACCTTGTGCGCGTAGGCAAGATTTTTCTCGca                                                       <  2:4395251/150‑1 (MQ=255)

                                                                                                 aaCGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCTTTTTCCAGTGCGGCCTTCACCTTGTGCGCGTAGGCAAGATTTTTCTCGcac                                                      >  2:1753174/1‑150 (MQ=255)

                                                                                                    gctcgctTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCTTTTTCCAGTGCGGCCTTCACCTTGTGCGCGTAGGCAAGATTTTTCTCGCACATc                                                   <  2:794764/150‑1 (MQ=255)

                                                                                                    gctcgctTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCTTTTTCCAGTGCGGCCTTCACCTTGTGCGCGTAGGCAAGATTTTTCTCGCACATc                                                   <  2:3574112/150‑1 (MQ=255)

                                                                                                         ctTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCTTTTTCCAGTGCGGCCTTCACCTTGTGCGCGTAGGCAAGATTTTTCTCGCACATCGGTGc                                              <  2:2848919/150‑1 (MQ=255)

                                                                                                         ctTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCTTTTTCCAGTGCGGCCTTCACCTTGTGCGCGTAGGCAAGATTTTTCTCGCACATCGGTGc                                              <  2:3871472/150‑1 (MQ=255)

                                                                                                           tACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCTTTTTCCAGTGCGGCCTTCACCTTGTGCGCGTAGGCAAGATTTTTCTCGCACATCGGTGCCg                                            >  2:329030/1‑150 (MQ=255)

                                                                                                             ccGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCTTTTTCCAGTGCGGCCTTCACCTTGTGCGCGTAGGCAAGATTTTTCTCGCACATCGGTGCCGCa                                          <  1:2611964/150‑1 (MQ=255)

                                                                                                                   aTGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCTTTTTCCAGTGCGGCCTTCACCTTGTGCGCGTAGGCAAGATTTTTCTCGCACATCGGTGCCGCAGGAATg                                    >  1:4305018/1‑150 (MQ=255)

                                                                                                                        gAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCTTTTTCCAGTGCGGCCTTCACCTTGTGCGCGTAGGCAAGATTTTTCTCGCACATCGGTGCCGCAGGAATGTACAt                               >  1:1375725/1‑150 (MQ=255)

                                                                                                                         aGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCTTTTTCCAGTGCGGCCTTCACCTTGTGCGCGTAGGCAAGATTTTTCTCGCACATCGGTGCCGCAGGAATGTACATc                              <  2:2877367/150‑1 (MQ=255)

                                                                                                                          ggTCGGCAAGCGTAAAGCGTCCTTCCCTGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCTTTTTCCAGTGCGGCCTTCACCTTGTGCGCGTAGGCAAGATTTTTCTCGCACATCGGTGCCGCAGGAATGTACATCa                             <  2:2634854/150‑1 (MQ=255)

                                                                                                                          ggTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCTTTTTCCAGTGCGGCCTTCACCTTGTGCGCGTAGGCAAGATTTTTCTCGCACATCGGTGCCGCAGGAATGTACATCa                             >  1:1874654/1‑150 (MQ=255)

                                                                                                                             cGGCAAGCGTAAAGCGTCCTTCCCTGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCTTTTTCCAGTGCGGCCTTCACCTTGTGCGCGTAGGCAAGATTTTTCTCGCACATCGGTGCCGCAGGAATGTACATCACGt                          <  2:4633684/150‑1 (MQ=255)

                                                                                                                              ggCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCTTTTTCCAGTGCGGCCTTCACCTTGTGCGCGTAGGCAAGATTTTTCTCGCACATCGGTGCCGCAGGAATGTACATCACGtt                         >  1:1061476/1‑150 (MQ=255)

                                                                                                                                   gCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCTTTTTCCAGTGCGGCCTTCACCTTGTGCGCGTAGGCAAGATTTTTCTCGCACATCGGTGCCGCAGGAATGTACATCACGTTGcccc                    <  1:553347/150‑1 (MQ=255)

                                                                                                                                    cGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCTTTTTCCAGTGCGGCCTTCACCTTGTGCGCGTAGGCAAGATTTTTCTCGCACATCGGTGCCGCAGGAATGTACATCACGTTGCCCCa                   <  1:254993/150‑1 (MQ=255)

                                                                                                                                     gTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGTGCCTTTTTCCAGTGCGGCCTTCACCTTGTGCGCGTAGGCAAGATTTTTCTCGCACATCGGTGCCGCAGGAATGTACATCACGTTGCCCCaa                  >  2:4572207/1‑150 (MQ=255)

                                                                                                                                        aaGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCTTTTTCCAGTGCGGCCTTCACCTTGTGCGCGTAGGCAAGATTTTTCTCGCACATCGGTGCCGCAGGAATGTACATCACGTTGCCCCAACCt               <  2:2856131/150‑1 (MQ=255)

                                                                                                                                         aGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCTTTTTCCAGTGCGGCCTTCACCTTGTGCGCGTAGGCAAGATTTTTCTCGCACATCGGTGCCGCAGGAATGTACATCACGTTGCCCCAACCtt              <  1:966544/150‑1 (MQ=255)

                                                                                                                                         aGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCTTTTTCCAGTGCGGCCTTCACCTTGTGCGCGTAGGCAAGATTTTTCTCGCACATCGGTGCCGCAGGAATGTACATCACGTTGCCCCAACCtt              >  2:1667534/1‑150 (MQ=255)

                                                                                                                                            gTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCTTTTTCCAGTGCGGCCTTCACCTTGTGCGCGTAGGCAAGATTTTTCTCGCACATCGGTGCCGCAGGAATGTACATCACGTTGCCCCAACCttgtt           >  2:2173467/1‑150 (MQ=255)

                                                                                                                                             tCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCTTTTTCCAGTGCGGCCTTCACCTTGTGCGCGTAGGCAAGATTTTTCTCGCACATCGGTGCCGCAGGAATGTACATCACGTTGCCCCAACCttgttg          <  2:905162/150‑1 (MQ=255)

                                                                                                                                             tCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCTTTTTCCAGTGCGGCCTTCACCTTGTGCGCGTAGGCAAGATTTTTCTCGCACATCGGTGCCGCAGGAATGTACATCACGTTGCCCCAACCttgttg          >  2:4745201/1‑150 (MQ=255)

                                                                                                                                             tCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCTTTTTCCAGTGCGGCCTTCACCTTGTGCGCGTAGGCAAGATTTTTCTCGCACATCGGTGCCGCAGGAATGTACATCACGTTGCCCCAACCttgttg          <  2:3970373/150‑1 (MQ=255)

                                                                                                                                               cTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCTTTTTCCAGTGCGGCCTTCACCTTGTGCGCGTAGGCAAGATTTTTCTCGCACATCGGTGCCGCAGGAATGTACATCACGTTGCCCCAACCTTGTTGAt        >  1:2018212/1‑150 (MQ=255)

                                                                                                                                                 tCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCTTTTTCCAGTGCGGCCTTCACCTTGTGCGCGTAGGCAAGATTTTTCTCGCACATCGGTGCCGCAGGAATGTACATCACGTTGCCCCAACCTTGTTGAttt      <  1:2209263/150‑1 (MQ=255)

                                                                                                                                                 tCCCAGTTAGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCTTTTTCCAGTGCGGCCTTCACCTTGTGCGCGTAGGCAAGATTTTTCTCGCACATCGGTGCCGCAGGAATGTACATCACGTTGCCCCAACCTTGTTGAttt      <  1:1717149/150‑1 (MQ=255)

                                                                                                                                                  cccAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCTTTTTCCAGTGCGGCCTTCACCTTGTGCGCGTAGGCAAGATTTTTCTCGCACATCGGTGCCGCAGGAATGTACATCACGTTGCCCCAACCTTGTTGAtttt     <  1:626750/150‑1 (MQ=255)

                                                                                                                                                     aGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCTTTTTCCAGTGCGGCCTTCACCTTGTGCGCGTAGGCAAGATTTTTCTCGCACATCGGTGCCGCAGGAATGTACATCACGTTGCCCCAACCTTGTTGATTTTCaa  >  2:1416707/1‑150 (MQ=255)

                                                                                                                                                     aGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCTTTTTCCAGTGCGGCCTTCACCTTGTGCGCGTAGGCAAGATTTTTCTCGCACATCGGTGCCGCAGGAATGTACATCACGTTGCCCCAACCTTGTTGATTTTCaa  >  2:4222/1‑150 (MQ=255)

                                                                                                                                                     aGTGTGATTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCTTATTCCAGAGCGGCCTTCACCTTGTGAGCGTCGGCAAGATTTATCTCGCACATCGGTGCCGCAGGAAACAACATCACGTTGCCCAAACCTAGTTGATTTTCaa  >  2:3726/1‑150 (MQ=255)

TTAACAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCTTTTTCCAGTGCGGCCTTCACCTTGTGCGCGTAGGCAAGATTTTTCTCGCACATCGGTGCCGCAGGAATGTACATCACGTTGCCCCAACCTTGTTGATTTTCAA  >  CP009273/853171‑853469

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