BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	4,341,023	A→T	100%	Y153N (TAC→AAC)	dcuS ←	sensory histidine kinase in two‑component regulatory system with DcuR, regulator of anaerobic fumarate respiration

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	4,341,023	0	A	T	99.6%	829.4 / ‑5.8	252	Y153N (TAC→AAC)	dcuS	sensory histidine kinase in two‑component regulatory system with DcuR, regulator of anaerobic fumarate respiration
Reads supporting (aligned to +/- strand): ref base A (1/1); new base T (136/114); total (137/115)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 252 total aligned reads displayed.

GAATGCAGGTGCCAATCAGTCCGACCAGCATGCCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTAGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACATTTTCTTCGCCATTCAGCGCTTTAAGGATGTCATCACCTTTAAATGGCTGACCAATACGCTGGGCTTCAGGATGCGAGTAGCGAAGACTT  >  CP009273/4340879‑4341168

gAATGCAGGTGCCAATCAGTCCGACCAGCATGCCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTTGATg                                                                                                                                               >  1:524120/1‑149 (MQ=255)

  aTGCAGGTGCCAATCAGTCCGACCAGCATGCCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTTGATggg                                                                                                                                             >  1:703942/1‑149 (MQ=255)

          gCCAATCAGTCCGACCAGCATGCCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTTGATGGGGGTAAATa                                                                                                                                     <  2:891543/149‑1 (MQ=255)

           ccAATCAGTCCGACCAGCATGCCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTTGATGGGGGTAAATAc                                                                                                                                    <  2:3540648/149‑1 (MQ=255)

           ccAATCAGTCCGACCAGCATGCCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTTGATGGGGGTAAATAc                                                                                                                                    >  1:2985743/1‑149 (MQ=255)

              aTCAGTCCGACCAGCATGCCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTTGATGGGGGTAAATAcgcg                                                                                                                                 >  1:3584933/1‑149 (MQ=255)

                cAGTCCGACCAGCATGCCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTTGATGGGGGTAAATACGCGTa                                                                                                                               >  2:521540/1‑149 (MQ=255)

                     cGACCAGCATGCCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTTGATGGGGGTAAATACGCGTAAAGcc                                                                                                                          <  2:3296806/149‑1 (MQ=255)

                      gACCAGCATGCCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTTGATGGGGGTAAATACGCGTAAAGCCt                                                                                                                         >  1:163125/1‑149 (MQ=255)

                       aCCAGCTTGCCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTTGATGGGGGTAAATACGCGTAAAGCCTg                                                                                                                        <  1:649239/149‑1 (MQ=255)

                        ccAGCATGCCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTTGATGGGGGTAAATACGCGTAAAGCCTgc                                                                                                                       >  2:2099634/1‑149 (MQ=255)

                           gCATGCCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTTGATGGGGGTAAATACGCGTAAAGCCTGCGcc                                                                                                                    >  1:3498365/1‑149 (MQ=255)

                           gCATGCCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTTGATGGGGGTAAATACGCGTAAAGCCTGCGcc                                                                                                                    <  1:4099809/149‑1 (MQ=255)

                           gCATGCCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTTGATGGGGGTAAATACGCGTAAAGCCTGCGcc                                                                                                                    <  1:1729994/149‑1 (MQ=255)

                                 cAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTTGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGaaaa                                                                                                              <  2:1899244/149‑1 (MQ=255)

                                  aaaTAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTTGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAAc                                                                                                             <  1:402866/149‑1 (MQ=255)

                                        taCCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTTGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGACAACCGCGAt                                                                                                       >  1:3229796/1‑149 (MQ=255)

                                        taCCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTTGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGAt                                                                                                       >  1:203526/1‑149 (MQ=255)

                                          ccGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTTGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTg                                                                                                     >  2:1085362/1‑149 (MQ=255)

                                           cGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTTGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGa                                                                                                    <  2:2110599/149‑1 (MQ=255)

                                            gACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTTGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGAt                                                                                                   <  2:729678/149‑1 (MQ=255)

                                            gACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTTGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGAt                                                                                                   <  2:3886832/149‑1 (MQ=255)

                                            gACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTTGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGAt                                                                                                   >  2:1546657/1‑149 (MQ=255)

                                             aCCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTTGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATa                                                                                                  >  2:912636/1‑149 (MQ=255)

                                             aCCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAAGTTGTTTATGATTTTCATCGTTGATGTGGGTAAATACGCGTAAAGCCTGCGCCAGAAACACGCGATTGATa                                                                                                  >  2:3760174/1‑149 (MQ=255)

                                              ccAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTTGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGAGTGATAg                                                                                                 >  2:896633/1‑149 (MQ=255)

                                                   tAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTTGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACa                                                                                            <  2:3603335/149‑1 (MQ=255)

                                                    aaggCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTTGTTGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACAt                                                                                           <  1:3783280/146‑1 (MQ=255)

                                                     aTGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTTGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACAtt                                                                                          >  2:254804/1‑149 (MQ=255)

                                                      tGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCTTTTTGTTTTTGTTTTTCATCGTTGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACAttt                                                                                         >  1:1997056/1‑149 (MQ=255)

                                                      tGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTTGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACAttt                                                                                         >  1:747704/1‑149 (MQ=255)

                                                      tGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTTGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACAttt                                                                                         >  2:2439673/1‑149 (MQ=255)

                                                       gCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTTGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACAtttt                                                                                        >  2:3128259/1‑149 (MQ=255)

                                                         tCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTTGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACATTttct                                                                                      <  1:1113651/149‑1 (MQ=255)

                                                          ccAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTTGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACATTttctt                                                                                     >  1:1003243/1‑149 (MQ=255)

                                                          ccAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTTGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACATTttctt                                                                                     >  1:1002977/1‑149 (MQ=255)

                                                            aGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTTGATGTGGGTAAATACGCGTAAAGCCTGCTCCAGAAAACCGCGATTGATAGCGACATTTTCTTCg                                                                                   >  2:75909/1‑149 (MQ=255)

                                                             gCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTTGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACATTTTCTTCGc                                                                                  >  2:1806159/1‑149 (MQ=255)

                                                              cGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTTGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACATTTTCTTCGcc                                                                                 <  1:956822/149‑1 (MQ=255)

                                                                aCTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTTGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACATTTTCTTCGCCAt                                                                               <  2:3091709/149‑1 (MQ=255)

                                                                    tCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTTGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACATTTTCTTCGCCATTCAg                                                                           >  1:1777767/1‑149 (MQ=255)

                                                                    tCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTTGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACATTTTCTTCGCCATTCAg                                                                           >  2:1628862/1‑149 (MQ=255)

                                                                    tCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTTGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACATTTTCTTCGCCATTCAg                                                                           >  1:668660/1‑149 (MQ=255)

                                                                     cATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTTGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACATTTTCTTCGCCATTCAgc                                                                          <  1:2353678/149‑1 (MQ=255)

                                                                        tGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTTGATGTGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACATTTTCTTCGCCATTCAGCGCt                                                                       >  1:4012379/1‑149 (MQ=255)

                                                                          aTCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTTGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACATTTTCTTCGCCATTCAGCGCttt                                                                     <  1:540603/149‑1 (MQ=255)

                                                                           tCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTTGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACATTTTCTTCGCCATTCAGCGCTTTa                                                                    <  2:781146/149‑1 (MQ=255)

                                                                             tGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTTGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACATTTTCTTCGCCATTCAGCGCTTTAAg                                                                  <  2:2344158/149‑1 (MQ=255)

                                                                              gTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTTGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACATTTTCTTCGCCATTCAGCGCTTTAAgg                                                                 >  1:2481731/1‑149 (MQ=255)

                                                                                   gTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTTGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACATTTTCTTCGCCATTCAGCGCTTTAAGGATGtc                                                            >  2:2740933/1‑149 (MQ=255)

                                                                                     cacaCGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTTGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACATTTTCTTCGCCATTCAGCGCTTTAAGGATGtcat                                                          <  2:3875606/149‑1 (MQ=255)

                                                                                       cacGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTTGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACATTTTCTTCGCCATTCAGCGCTTTAAGGATGtcatca                                                        >  1:1527668/1‑149 (MQ=255)

                                                                                          ggCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTTGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACATTTTCTTCGCCATTCAGCGCTTTAAGGATGTCATCACCt                                                     <  1:1194070/149‑1 (MQ=255)

                                                                                           gCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTTGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACATTTTCTTCGCCATTCAGCGCTTTAAGGATGTCATCACCtt                                                    >  2:1401344/1‑149 (MQ=255)

                                                                                            cTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTTGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACATTTTCTTCGCCATTCAGCGCTTTAAGGATGTCATCACCttt                                                   >  1:1832504/1‑149 (MQ=255)

                                                                                             ttAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTTGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACATTTTCTTCGCCATTCAGCGCTTTAAGGATGTCATCACCTTTa                                                  <  1:1821624/149‑1 (MQ=255)

                                                                                             ttAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTTGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACATTTTCTTCGCCATTCAGCGCTTTAAGGATGTCATCACCTTTa                                                  <  1:2788741/149‑1 (MQ=255)

                                                                                               aaCTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTTGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACATTTTCTTCGCCATTCAGCGCTTTAAGGATGTCATCACCTTTaaa                                                >  2:2956665/1‑149 (MQ=255)

                                                                                                 cTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTTGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACATTTTCTTCGCCATTCAGCGCTTTAAGGATGTCATCACCTTTAAATg                                              >  2:225787/1‑149 (MQ=255)

                                                                                                     aGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTTGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACATTTTCTTCGCCATTCAGCGCTTTAAGGATGTCATCACCTTTAAATGGCTg                                          >  2:2398877/1‑149 (MQ=255)

                                                                                                     aGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTTGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACATTTTCTTCGCCATTCAGCGCTTTAAGGATGTCATCACCTTTAAATGGCTg                                          >  1:3904487/1‑149 (MQ=255)

                                                                                                      ggCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTTGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACATTTTCTTCGCCATTCAGCGCTTTAAGGATGTCATCACCTTTAAATGGCTGa                                         <  1:4074718/149‑1 (MQ=255)

                                                                                                        ccGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTTGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACATTTTCTTCGCCATTCAGCGCTTTAAGGATGTCATCACCTTTAAATGGCTGAcc                                       <  1:2099634/149‑1 (MQ=255)

                                                                                                         cGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTTGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACATTTTCTTCGCCATTCAGCGCTTTAAGGATGTCATCACCTTTAAATGGCTGACCa                                      >  2:3985635/1‑149 (MQ=255)

                                                                                                          gATCGCCACCACGCCAATTTGTTTATGATTCTCATCGTTGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACATTTTCTTCGCCATTCAGCGCTTTAAGGATGTCATCACCTTTAAATGGCTGACCaa                                     <  2:1134806/149‑1 (MQ=255)

                                                                                                           aTCGCCACGACGCCAATTTGTTTATGATTTTCATCGTTGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACATTTTCTTCGCCATTCAGCGCTTTAAGGATGTCATCACCTTTAAATGGCTGACCAAt                                    >  2:734135/1‑149 (MQ=255)

                                                                                                           aTCGCCACCACGCCAATTTGTTTATGATTTTCATCGTTGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACATTTTCTTCGCCATTCAGCGCTTTAAGGATGTCATCACCTTTAAATGGCTGACCAAt                                    >  1:1353152/1‑149 (MQ=255)

                                                                                                            tCGCCACCACGCCAATTTGTTTATGATTTTCATCGTTGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACATTTTCTTCGCCATTCAGCGCTTTAAGGATGTCATCACCTTTAAATGGCTGACCAATa                                   >  2:2550589/1‑149 (MQ=255)

                                                                                                            tCGCCACCACGCCAATTTGTTTATGATTTTCATCGTTGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACATTTTCTTCGCCATTCAGCGCTTTAAGGATGTCATCACCTTTAAATGGCTGACCAATa                                   <  1:3373947/149‑1 (MQ=255)

                                                                                                            tCGCCACCACGCCAATTTGTTTATGATTTTCATCGTTGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACATTTTCTTCGCCATTCAGCGCTTTAAGGATGTCATCACCTTTAAATGGCTGACCAATa                                   <  1:570317/149‑1 (MQ=255)

                                                                                                            tCGCCACCACGCCAATTTGTTTATGATTTTCATCGTTGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACATTTTCTTCGCCATTCAGCGCTTTAAGGATGTCATCACCTTTAAATGGCTGACCAATa                                   >  1:2503849/1‑149 (MQ=255)

                                                                                                             cGCCACCACGCCAATTTGTTTATGATTTTCATCGTTGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACATTTTCTTCGCCATTCAGCGCTTTAAGGATGTCATCACCTTTAAATGGCTGACCAATAc                                  <  2:488383/149‑1 (MQ=255)

                                                                                                             cGCCACCACGCCAATTTGTGTATGATTTTCATCGTAGATGGGGGTAAATCCGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACATTTTCTTCGCCATTCAGCGCTTTAAGGATGTCATCACCTTTAAATGGCTGACCAATAc                                  <  2:2657734/149‑1 (MQ=255)

                                                                                                              gCCACCACGCCAATTTGTTTATGATTTTCATCGTTGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACATTTTCTTCGCCATTCAGCGCTTTAAGGATGTCATCACCTTTAAATGTCTGACCAATACg                                 >  2:828167/1‑149 (MQ=255)

                                                                                                               ccaccaCGCCAATTTGTTTATGATTTTCATCGTTGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACATTTTCTTCGCCATTCAGCGCTTTAAGGATGTCATCACCTTTAAATGGCTGACCAATACGc                                <  1:254804/149‑1 (MQ=255)

                                                                                                                 accacGCCAATTTGTTTATGATTTTCATCGTTGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACATTTTCTTCGCCATTCAGCGCTTTAAGGATGTCATCACCTTTAAATGGCTGACCAATACGCTg                              <  1:2887910/149‑1 (MQ=255)

                                                                                                                 accacGCCAATTTGTTTATGATTTTCATCGTTGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACATTTTCTTCGCCATTCAGCGCTTTAAGGATGTCATCACCTTTAAATGGCTGACCAATACGCTg                              >  1:3430860/1‑149 (MQ=255)

                                                                                                                    tcGCCAATTTTTTTATGATTTTCATCGTTGATTGGTGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACATTTTCTTCGCCATTCAGCGCTTTAAGGATGTCATCACCTTTAAATGGCTGACCAATACGCTGGGc                           <  2:3167204/148‑1 (MQ=255)

                                                                                                                      gCCAATTTGTTTATGATTTTCATCGTTGATGGGGGTAAATACGCGTAAAGCCTGCTCCAGAAAACCGCGATTGATAGCGACATTTTCTTCGCCATTCAGCGCTTTAAGGATGTCATCACCTTTAAATGGCTGACCAATACGCTGGGCtt                         >  1:2235288/1‑149 (MQ=255)

                                                                                                                        cAATTTGTTTATGATTTTCATCGTTGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACATTTTCTTCGCCATTCAGCGCTTTAAGGATGTCATCACCTTTAAATGGCTGACCAATACGCTGGGCTTCa                       >  2:249451/1‑149 (MQ=255)

                                                                                                                          aTTTGTTTATGATTTTCATCGTTGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACATTTTCTTCGCCATTCAGCGCTTTAAGGATGTCATCACCTTTAAATGGCTGACCAATACGCTGGGCTTCAgg                     <  2:3484395/149‑1 (MQ=255)

                                                                                                                           tttGTTTATGATTTTCATCGTTGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACATTTTCTTCGCCATTCTGCGCTTTAAGGATGTCATCACCTTTAAATGGCTGACCAATACGCTGGGCTTCAGGa                    >  2:277777/1‑149 (MQ=255)

                                                                                                                           tttGTTTATGATTTTCATCGTTGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACATTTTCTTCGCCATTCAGCGCTTTAAGGATGTCATCACCTTTAAATGGCTGACCAATACGCTGGGCTTCAGGa                    >  2:3477334/1‑149 (MQ=255)

                                                                                                                           tttGTTTATGATTTTCATCGTTGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACATTTTCTTCGCCATTCAGCGCTTTAAGGATGTCATCACCTTTAAATGGCTGACCAATACGCTGGGCTTCAGGa                    >  2:1429693/1‑149 (MQ=255)

                                                                                                                              gTTTATGATTTTCATCGTTGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACATTTTCTTCGCCATTCAGCGCTTTAAGGATGTCATCACCTTTAAATGGCTGACCAATACGCTGGGCTTCAGGATGc                 >  2:3694961/1‑149 (MQ=255)

                                                                                                                                ttATGATTTTCATCGTTGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACATTTTCTTCGCCATTCAGCGCTTTAAGGATGTCATCACCTTTAAATGGCTGACCAATACGCTGGGCTTCAGGATGCGa               <  2:345497/149‑1 (MQ=255)

                                                                                                                                ttATGATTTTCATCGTTGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACATTTTCTTCGCCATTCAGCGCTTTAAGGATGTCATCACCTTTAAATGGCTGACCAATACGCTGGGCTTCAGGATGCGa               >  2:406972/1‑149 (MQ=255)

                                                                                                                                  aTGATTTTCATCGTTGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACATTTTCTTCGCCATTCAGCGCTTTAAGGATGTCATCACCTTTAAATGGCTGACCAATACGCTGGGCTTCAGGATGCGAGt             >  2:1502137/1‑149 (MQ=255)

                                                                                                                                   tGATTTTCATCGTTGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACATTTTCTTCGCCATTCAGCGCTTTAAGGATGTCATCACCTTTAAATGGCTGACCAATACGCTGGGCTTCAGGATGCGAGTa            <  1:3336473/149‑1 (MQ=255)

                                                                                                                                   tGATTTTCATCGTTGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACATTTTCTTCGCCATTCAGCGCTTTAAGGATGTCATCACCTTTAAATGGCTGACCAATACGCTGGGCTTCAGGATGCGAGTa            <  1:2857030/149‑1 (MQ=255)

                                                                                                                                      ttttCATCGTTGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACATTTTCTTCGCCATTCAGCGCTTTAAGGATGTCATCACCTTTAAATGGCTGACCAATACGCTGGGCTTCAGGATGCGAGTAGCg         >  2:471105/1‑149 (MQ=255)

                                                                                                                                       tttCATCGTTGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACATTTTCTTCGCCATTCAGCGCTTTAAGGATGTCATCACCTTTAAATGGCTGACCAATACGCTGGGCTTCAGGATGCGAGTAGCGa        <  1:2125233/149‑1 (MQ=255)

                                                                                                                                        ttCATCGTTGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACATTTTCTTCGCCATTCAGCGCTTTAAGGATGTCATCACCTTTAAATGGCTGACCAATACGCTGGGCTTCAGGATGCGAGTAGCGaa       <  1:3277907/149‑1 (MQ=255)

                                                                                                                                          cATCGTTGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACATTTTCTTCGCCATTCAGCGCTTTAAGGATGTCATCACCTTTAAATGGCTGACCAATACGCTGGGCTTCAGGATGCGAGTAGCGAAGa     <  1:2439673/149‑1 (MQ=255)

                                                                                                                                          cATCGTTGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACATTTTCTTCGCCATTCAGCGCTTTAAGGATGTCATCACCTTTAAATGGCTGACCAATACGCTGGGCTTCAGGATGCGAGTAGCGAAGa     <  1:1247213/149‑1 (MQ=255)

                                                                                                                                           aTCGTTGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACATTTTCTTCGCCATTCAGCGCTTTAAGGATGTCATCACCTTTAAATGGCTGACCAATACGCTGGGCTTCAGGATGCGAGTAGCGAAGAc    >  1:372100/1‑149 (MQ=255)

                                                                                                                                           aTCGTTGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACATTTTCTTCGCCATTCAGCGCTTTAAGGATGTCATCACCTTTAAATGGCTGACCAATACGCTGGGCTTCAGGATGCGAGTAGCGAAGAc    >  1:372142/1‑149 (MQ=255)

                                                                                                                                           aTCGTTGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACATTTTCTTCGCCATTCAGCGCTTTAAGGATGTCATCACCTTTAAATGGCTGACCAATACGCTGGGCTTCAGGATGCGAGTAGCGAAGAc    >  2:43205/1‑149 (MQ=255)

                                                                                                                                             cGTTGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACATTTTCTTCGCCATTCAGCGCTTTAAGGATGTCATCACCTTTAAATGGCTGACCAATACGCTGGGCTTCAGGATGCGAGTAGCGAAGACtt  >  2:676612/1‑149 (MQ=255)

                                                                                                                                             cGTTGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACATTTTCTTCGCCATTCAGCGCTTTAAGGATGTCATCACCTTTAAATGGCTGACCAATACGCTGGGCTTCAGGATGCGAGTAGCGAAGACtt  <  2:3498365/149‑1 (MQ=255)

GAATGCAGGTGCCAATCAGTCCGACCAGCATGCCAAATAATACCGACCAGATAATGCTCCAGCGACTGTCATTGATCTGTTGGGTCACACGGCTTAACTCAAGGCCGATCGCCACCACGCCAATTTGTTTATGATTTTCATCGTAGATGGGGGTAAATACGCGTAAAGCCTGCGCCAGAAAACCGCGATTGATAGCGACATTTTCTTCGCCATTCAGCGCTTTAAGGATGTCATCACCTTTAAATGGCTGACCAATACGCTGGGCTTCAGGATGCGAGTAGCGAAGACTT  >  CP009273/4340879‑4341168

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