BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	1,689,180	T→G	45.5%	E383D (GAA→GAC)	uidA ←	beta‑D‑glucuronidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	1,689,180	0	T	G	45.5%	9.8 / inf	285	E383D (GAA→GAC)	uidA	beta‑D‑glucuronidase
Reads supporting (aligned to +/- strand): ref base T (68/87); new base G (66/63); total (135/150)
Fisher's exact test for biased strand distribution p-value = 2.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.15e-01

Only 100 of 292 total aligned reads displayed.

CGCGTTGCTTCCGCCAGTGGCGCGAAATATTCCCGTGCACCTTGCGGACGGGTATCCGGTTCGTTGGCAATACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTTTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGACAGCAGCAGTTTCATCAATCACCACGATGCCATGTTCATCTGCCCAGTCGAGCATCTCTTCAGC  >  CP009273/1689031‑1689329

cgcgTTGCTTCCGCCAGTGGCGCGAAATATTCCCGTGCACCTTGCGGACGGGTATCCGGTTCGTTGGCAATACTCCACATCACCACGCTTGTGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGtt                                                                                                                                                       <  1:1541167/150‑1 (MQ=255)

    ttGCTTCCGCCAGTGGCGCGAAATATTCCCGTGCACCTTGCGGACGGGTATCCGGTTCGTTGGCAATACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTccc                                                                                                                                                   >  2:3970531/1‑150 (MQ=255)

        ttCCGCCAGTGGCGCGAAATATTCCCGTGCACCTTGCGGACGGGTATCCGGTTCGTTGGCAATACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTTTCCCCGtt                                                                                                                                               >  2:788038/1‑150 (MQ=255)

           cGCCAGTGGCGCGAAATATTCCCGTGCACCTTGCGGACGGGTATCCGGTTCGTTGGCAATACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGAc                                                                                                                                            >  2:1015688/1‑150 (MQ=255)

           cGCCAGTGGCGCGAAATAGTCCCGTGCACCTTGCGGACGGGTATCCGGTTCGTTGGCAATACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGAc                                                                                                                                            <  2:1305424/150‑1 (MQ=255)

               aGTGGCGCGAAATATTCCCGTGCACCTTGCGGACGGGTATCCGGTTCGTTGGCAATACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTTTCCCCGTTGACTGcc                                                                                                                                        >  2:2609409/1‑150 (MQ=255)

               aGTGGCGCGAAATATTCCCGTGCACCTTGCGGACGGGTATCCGGTTCGTTGGCAATACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTTTCCCCGTTGACTGcc                                                                                                                                        >  2:4149529/1‑150 (MQ=255)

                 tGGCGCGAAATATTCCCGTGCACCTTGCGGACGGGTATCCGGTTCGTTGGCAATACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCctc                                                                                                                                      >  1:1495491/1‑150 (MQ=255)

                       gAAATATTCCCGTGCACCTTGCGGACGGGTATCCGGTTCGTTGGCAATACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTTTCCCCGTTGACTGCCTCTTCGCt                                                                                                                                >  2:3185962/1‑150 (MQ=255)

                          atatTCCCGTGCACCTTGCGGACGGGTATCCGGTTCGTTGGCAATACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTTTCCCCGTTGACTGCCTCTTCGCTGTa                                                                                                                             >  1:1371798/1‑150 (MQ=255)

                            atTCCCGTGCACCTTGCGGACGGGTATCCGGTTCGTTGGCAATACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTTTCCCCGTTGACTGCCTCTTCGCTGTACa                                                                                                                           <  2:3680804/150‑1 (MQ=255)

                               cccGTGCACCTTGCGGACGGGTATCCGGTTCGTTGGCAATACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTTTCCCCGTTGACTGCCTCTTCGCTGTACAGtt                                                                                                                        >  2:1718290/1‑150 (MQ=255)

                                 cGTGCACATTGAGGACGTGTATCCTGTTAGTTGGCTATACTCCACATCACCAAGCTTGTGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCATGTAAGTGCGCTTGCTGATTGTACCCGTTGACTGCCTCTTCGCTGTACAGTTCt                                                                                                                      <  1:2714713/150‑1 (MQ=255)

                                  gTGCACCTTGCGGACGGGTATCCGGTTCGTTGGCAATACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCtt                                                                                                                     <  1:3797308/150‑1 (MQ=255)

                                     cACCTTGCGGACGGGTATCCGGTTCGTTGGCAATACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTTTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCg                                                                                                                  <  1:1196093/150‑1 (MQ=255)

                                     cACCTTGCGGACGGGTATCCGGTTCGTTGGCAATACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTTTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCg                                                                                                                  <  2:2189271/150‑1 (MQ=255)

                                       ccTTGCGGACGGGTATCCGGTTCGTTGGCAATACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGc                                                                                                                >  2:1556249/1‑150 (MQ=255)

                                       ccTTGCGGACGGGTATCCGGTTCGTTGGCAATACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGc                                                                                                                >  1:751597/1‑150 (MQ=255)

                                         ttGCGGACGGGTATCCGGTTCGTTGGCAATACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTTTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCtt                                                                                                              <  1:3146392/150‑1 (MQ=255)

                                          tGCGGACGGGTATCCGGTTCGTTGGCAATACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTTTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCttg                                                                                                             <  2:1371798/150‑1 (MQ=255)

                                            cGGACGGGTATCCGGTTCGTTGGCAATACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTTTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCttgtt                                                                                                           >  2:2374356/1‑150 (MQ=255)

                                               aCGGGTATCGGGTTCGTTGGCAATACTCCACATCTCCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGcc                                                                                                        <  1:3899353/150‑1 (MQ=255)

                                               aCGGGTATCCGGTTCGTTGGCAATACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGcc                                                                                                        <  1:1714148/150‑1 (MQ=255)

                                                    tATCCGGTTCGTTGGCAATACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTTTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCtt                                                                                                   <  2:4227093/150‑1 (MQ=255)

                                                    tATCCGGTTCGTTGGCAATACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCtt                                                                                                   >  1:529585/1‑150 (MQ=255)

                                                     aTCCGGTTCGTTGGCAATACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTTTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTc                                                                                                  <  1:1811860/150‑1 (MQ=255)

                                                     aTCCGGTTCGTTGGCAATACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTc                                                                                                  <  1:4099508/150‑1 (MQ=255)

                                                     aTCCGGTTCGTTGGCAATACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTc                                                                                                  <  1:4330808/150‑1 (MQ=255)

                                                        cGGTTCGTTGGCAATACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGaa                                                                                               >  2:2221894/1‑150 (MQ=255)

                                                        cGGTTCGTTGGCAATACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGaa                                                                                               <  2:1417985/150‑1 (MQ=255)

                                                          gTTCGTTGGCAATACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAAc                                                                                             <  2:2281346/150‑1 (MQ=255)

                                                          gTTCGTTGGCAATACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAAc                                                                                             <  2:71883/150‑1 (MQ=255)

                                                            tCGTTGGCTATACTCCAGCTCACCTCGCTTGGGTGGTTTTTGTATCGCGCTTTCAGCTGTTTTATCGCGTGTAAGTGCGCTTGCTGAGTGTCGCCGTTGACTGCCTCTTGGCTGTACAGTTCTTTCGGCTTGTTGCCCGGTTGGatagca                                                                                           <  2:2462656/150‑6 (MQ=255)

                                                              gTTGGCAATACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTTTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAAt                                                                                         <  1:2922525/150‑1 (MQ=255)

                                                              gTTGGCAATACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTTTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAAt                                                                                         <  1:1946795/150‑1 (MQ=255)

                                                               ttGGCAATACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATg                                                                                        <  1:372155/150‑1 (MQ=255)

                                                                tGGCAATACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTTTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGc                                                                                       <  2:1341871/150‑1 (MQ=255)

                                                                 ggCAATACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTTTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGcc                                                                                      <  2:2764516/150‑1 (MQ=255)

                                                                  gaaaTACACCACATCACCACGCTTGGGTGGTTTTATTCACGCGCTATCAGCTCATTAATCGCCTGTAAGTGCGCTTGCTGAGTATCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTATCGGCTTGTTGCCCGCTTCGAAACCAATGCCt                                                                                     >  2:3626497/3‑150 (MQ=255)

                                                                  gCAATACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTTTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCt                                                                                     >  2:1434071/1‑150 (MQ=255)

                                                                    aaTACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTTTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTaa                                                                                   <  1:4083282/150‑1 (MQ=255)

                                                                     aTACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTaaa                                                                                  >  1:2662250/1‑150 (MQ=255)

                                                                        cTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAagag                                                                               >  2:180316/1‑150 (MQ=255)

                                                                          ccACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTTTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAagagag                                                                             >  1:3994154/1‑150 (MQ=255)

                                                                          ccACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAagagag                                                                             <  2:2562233/150‑1 (MQ=255)

                                                                          ccACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAagagag                                                                             <  1:682349/150‑1 (MQ=255)

                                                                            acaTCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTTTCCCCGTTGACTGCCTCTTCGCTGTTCAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGt                                                                           <  1:3081922/150‑1 (MQ=255)

                                                                              aTCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTTTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTa                                                                         >  1:1980397/1‑150 (MQ=255)

                                                                              aTCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTTTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTa                                                                         <  1:2351761/150‑1 (MQ=255)

                                                                              aTCACCACGCTTGGATGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTTTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTa                                                                         >  2:1675329/1‑150 (MQ=255)

                                                                               tCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTaa                                                                        >  1:949659/1‑150 (MQ=255)

                                                                                caccacGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTaaa                                                                       >  1:4238837/1‑150 (MQ=255)

                                                                                caccacGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTaaa                                                                       <  1:3553227/150‑1 (MQ=255)

                                                                                 accacGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTTTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAg                                                                      >  2:1850666/1‑150 (MQ=255)

                                                                                 accacGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAg                                                                      <  1:2568676/150‑1 (MQ=255)

                                                                                      gCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGAc                                                                 <  1:3152308/150‑1 (MQ=255)

                                                                                      gCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGAc                                                                 <  1:807579/150‑1 (MQ=255)

                                                                                       cTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTTTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGAca                                                                <  2:3273140/150‑1 (MQ=255)

                                                                                       cTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTTTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGAca                                                                >  1:2300576/1‑150 (MQ=255)

                                                                                       cTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTTTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGAca                                                                <  2:1941232/150‑1 (MQ=255)

                                                                                       cTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGAca                                                                <  1:1258786/150‑1 (MQ=255)

                                                                                          gggTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTTTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGAcagca                                                             <  1:3240422/150‑1 (MQ=255)

                                                                                             tggtTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTTTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGTTTAAATCCGAcagcagca                                                          >  2:2763052/1‑150 (MQ=255)

                                                                                              ggtTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGAcagcagcag                                                         >  2:4086949/1‑150 (MQ=255)

                                                                                                 ttttGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTTTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGACAGCAGCAGttt                                                      <  2:899857/150‑1 (MQ=255)

                                                                                                 ttttGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTTTCCCCGTTGACTGCCTCCTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGACAGCAGCAGttt                                                      <  2:447221/150‑1 (MQ=255)

                                                                                                   ttGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGACAGCAGCAGTTtca                                                    <  1:2323246/150‑1 (MQ=255)

                                                                                                    tGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTTTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGACAGCAGCAGTTtcat                                                   >  2:1616833/1‑150 (MQ=255)

                                                                                                    tGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGACTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGACAGCAGCAGTTtcat                                                   <  2:3736838/150‑1 (MQ=255)

                                                                                                    agTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTACCCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGACAGCAGCAGTAtcaa                                                   >  2:3491928/2‑149 (MQ=255)

                                                                                                      tCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGACAGCAGCAGTTtcatca                                                 >  1:1634488/1‑150 (MQ=255)

                                                                                                      tCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGACAGCAGCAGTTtcatca                                                 >  1:1880258/1‑150 (MQ=255)

                                                                                                      tCACGCGCTATCAGCTATTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGACAGCAGCAGTTtcatca                                                 <  2:339385/150‑1 (MQ=255)

                                                                                                          gcgcTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAACCCAATGCCTAAAGAGAGGTTAAAGCCGACAGCAGCAGTTTCatcaatca                                             <  1:1015688/150‑1 (MQ=255)

                                                                                                          gcgcTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGACAGCAGCAGTTTCatcaatca                                             <  1:1015272/150‑1 (MQ=255)

                                                                                                            gcTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTTTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGACAGCAGCAGTTTCATCAATcacc                                           >  1:4356944/1‑150 (MQ=255)

                                                                                                              tATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTTTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGACAGCAGCAGTTTCATCAATcaccac                                         <  2:1751167/150‑1 (MQ=255)

                                                                                                                 cAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGACAGCAGCAGTTTCATCAATCACCACGAt                                      >  1:3191267/1‑150 (MQ=255)

                                                                                                                   gCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGACAGCAGCAGTTTCATCAATCACCACGATGc                                    >  1:1928691/1‑150 (MQ=255)

                                                                                                                      ctTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGACAGCAGCAGTTTCATCAATCACCACGATGCCAt                                 >  2:4354481/1‑150 (MQ=255)

                                                                                                                            tCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGACAGCAGCAGTTTCATCAATCACCACGATGCCATGTTCAt                           <  2:75771/150‑1 (MQ=255)

                                                                                                                            tCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGACAGCAGCAGTTTCATCAATCACCACGATGCCATGTTCAt                           <  2:1495491/150‑1 (MQ=255)

                                                                                                                             cGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGACAGCAGCAGTTTCATCAATCACCACGATGCCATGTTCATc                          >  2:2800854/1‑150 (MQ=255)

                                                                                                                                 tGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGACAGCAGCAGTTTCATCAATCACCACGATGCCATGTTCATCTGcc                      >  2:3017753/1‑150 (MQ=255)

                                                                                                                                   tAAGTGCGCTTGCTGAGTTTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGACAGCAGCAGTTTCATCAATCACCACGATGCCATGTTCATCTGCCCa                    >  1:859166/1‑150 (MQ=255)

                                                                                                                                    aaGTGCGCTTGCTGAGTTTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGACAGCAGCAGTTTCATCAATCACCACGATGCCATGTTCATCTGCCCAg                   >  2:4013237/1‑150 (MQ=255)

                                                                                                                                      gTGCGCTTGCTGAGTTTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGACAGCAGCAGTTTCATCAATCACCACGATGCCATGTTCATCTGCCCAGTc                 >  2:508256/1‑150 (MQ=255)

                                                                                                                                       tGCGCTTGCTGAGTTTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGACAGCAGCAGTTTCATCAATCACCACGATGCCATGTTCATCTGCCCAGTCg                <  1:3272159/150‑1 (MQ=255)

                                                                                                                                             tGCTGAGTTTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGACAGCAGCAGTTTCATCAATCACCACGATGCCATGTTCATCTGCCCAGTCGAGCAtc          >  2:1122057/1‑150 (MQ=255)

                                                                                                                                               cTGAGTTTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGACAGCAGCAGTTTCATCAATCACCACGATGCCATGTTCATCTGCCCAGTCGAGCAtctc        <  2:3746061/150‑1 (MQ=255)

                                                                                                                                               cTGAGTTTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGACAGCAGCAGTTTCATCAATCACCACGATGCCATGTTCATCTGCCCAGTCGAGCAtctc        >  2:1040552/1‑150 (MQ=255)

                                                                                                                                               cTGAGTTTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGACAGCAGCAGTTTCATCAATCACCACGATGCCATGTTCATCTGCCCAGTCGAGCAtctc        >  2:1897528/1‑150 (MQ=255)

                                                                                                                                               cTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGACAGCAGCAGTTTCATCAATCACCACGATGCCATGTTCATCTGCCCAGTCGAGCAtctc        <  1:3789088/150‑1 (MQ=255)

                                                                                                                                               cTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGACAGCAGCAGTTTCATCAATCACCACGATGCCATGTTCATCTGCCCAGTCGAGCAtctc        <  1:1601564/150‑1 (MQ=255)

                                                                                                                                                 gAGTTTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGACAGCAGCAGTTTCATCAATCACCACGATGCCATGTTCATCTGCCCAGTCGAGCATCTCtt      >  1:2119820/1‑150 (MQ=255)

                                                                                                                                                 gAGTTTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGACAGCAGCAGTTTCATCAATCACCACGATGCCATGTTCATCTGCCCAGTCGAGCATCTCtt      >  1:4306573/1‑150 (MQ=255)

                                                                                                                                                  aGTTTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGACAGCAGCAGTTTCATCAATCACCACGATGCCATGTTCATCTGCCCAGTCGAGCATCTCTTc     >  1:2059455/1‑150 (MQ=255)

                                                                                                                                                  aGTTTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGACAGCAGCAGTTTCATCAATCACCACGATGCCATGTTCATCTGCCCAGTCGAGCATCTCTTc     >  1:3395275/1‑150 (MQ=255)

                                                                                                                                                   gTTTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGACAGCAGCAGTTTCATCAATCACCACGATGCCATGTTCATCTGCCCAGTCGAGCATCTCTTCa    <  2:487958/150‑1 (MQ=255)

                                                                                                                                                     ttCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGACAGCAGCAGTTTCATCAATCACCACGATGCCATGTTCATCTGCCCAGTCGAGCATCTCTTCAGc  >  1:3223521/1‑150 (MQ=255)

CGCGTTGCTTCCGCCAGTGGCGCGAAATATTCCCGTGCACCTTGCGGACGGGTATCCGGTTCGTTGGCAATACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTTTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGACAGCAGCAGTTTCATCAATCACCACGATGCCATGTTCATCTGCCCAGTCGAGCATCTCTTCAGC  >  CP009273/1689031‑1689329

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