BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	4,402,665	T→G	6.0%	G154G (GGT→GGG)	yjfM →	DUF1190 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	4,402,665	0	T	G	6.0%	812.8 / 7.6	292	G154G (GGT→GGG)	yjfM	DUF1190 family protein
Reads supporting (aligned to +/- strand): ref base T (122/145); new base G (21/4); total (143/149)
Fisher's exact test for biased strand distribution p-value = 2.45e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 302 total aligned reads displayed.

AATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCC‑GTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGTTCTTC  >  CP009273/4402516‑4402814

aaTATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCC‑GTTTGTTTATAACAGCGGCGGt                                                                                                                                                       <  2:2524775/150‑1 (MQ=255)

  tatGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCC‑GTTTGTTTATAACAGCGGCGGTTc                                                                                                                                                     <  1:1284654/150‑1 (MQ=255)

  tatGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCC‑GTTTGTTTATAACAGCGGCGGTTc                                                                                                                                                     <  2:541885/150‑1 (MQ=255)

   atGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCC‑GTTTGTTTATAACAGCGGCGGTTcc                                                                                                                                                    <  1:4439979/150‑1 (MQ=255)

        tCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTGGGGATTTTTTTTAATCCGGGGTATTGGTAAAGATCGCGATGAGCC‑GGTTGTTTATAAAACCGGCGGGTACTCttt                                                                                                                                               >  2:4783224/1‑150 (MQ=255)

         cagcagAACTGTCAGTCTAAGTACTAAAATTGCTAATACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCC‑GTTTGTTTATAACAGCGGCGGTTCCTCtttt                                                                                                                                              <  1:2256886/150‑1 (MQ=255)

          agcagAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGGTTTTTTTTAAGCCGGGTTTTTAGTAAAGATAGCGATGAGCC‑GTTTGTTTATAACAGCGGCGGTTCCGCTGTTg                                                                                                                                             >  2:2871893/1‑150 (MQ=255)

          agcagAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCC‑GTTTGTTTATAACAGCGGCGGTTCCTCTTTTg                                                                                                                                             <  2:4976519/150‑1 (MQ=255)

             agAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTTTTAATCCGGGTTATTAGTTAATATAGCGATGAGCC‑GTTTTTTTATAAAAGAGGCGGTTGGACTTTTGCtt                                                                                                                                          >  2:1464146/1‑150 (MQ=255)

             agAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCC‑GTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCtt                                                                                                                                          <  2:1187703/150‑1 (MQ=255)

             agAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATAGCGATGTGCC‑GTTTGTTTATAAAAGCGGCGGTTGCTCTTTTGcgt                                                                                                                                          >  2:691172/1‑148 (MQ=255)

              gAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTTTTAATCCGGGTTATTCGTAAAATTAGCGATTATCC‑GTTTGTTTATAAAAGCGGCGGTTACTCTTTTGCTTc                                                                                                                                         >  2:3322219/1‑150 (MQ=255)

                 cTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCC‑GTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTcgcg                                                                                                                                      <  1:3445351/150‑1 (MQ=255)

                  tGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTTTTAAGCCGGGTTATTAGTAAAGATAGCGATGAGCC‑GTTTGTTTATAACAGCGGCGGTTCCTCATTTGATTAgcgc                                                                                                                                     >  2:2341560/1‑150 (MQ=255)

                    tCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCC‑GTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGccc                                                                                                                                   <  2:3107781/150‑1 (MQ=255)

                    tCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTGTGGGTTTTTTTTAATCCGGGGTTTTAGGTAATATAGCGATGAGCC‑GGTTGTTTTTAACAGCGGCGGTTGCTCTTTTGCTTAGCGccc                                                                                                                                   >  1:2804202/1‑150 (MQ=255)

                     cAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCC‑GTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCa                                                                                                                                  <  2:4433764/150‑1 (MQ=255)

                     cAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCC‑GTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCa                                                                                                                                  <  2:556617/150‑1 (MQ=255)

                      aGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCC‑GTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAg                                                                                                                                 <  1:4349006/150‑1 (MQ=255)

                      aGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCC‑GTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAg                                                                                                                                 <  2:1579234/150‑1 (MQ=255)

                          tAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCC‑GTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTg                                                                                                                             <  1:4542529/150‑1 (MQ=255)

                          tAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCC‑GTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTg                                                                                                                             <  1:4250324/150‑1 (MQ=255)

                          tAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTAGTAAATATCGCGATTAGCC‑GTTTGTTTATAACAGCGGCGGGTCCTCTTTTGCTTCGCGCCCAggcag                                                                                                                             >  2:2882069/1‑146 (MQ=255)

                            aGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTATGGATTTTTTTTAATCCGGGTTATTAGTAAATATAGCGATGATCC‑GTTTGTTTATAAAAGCGGCGGTTCCTCATTTTCTTCGCGCCCAGTCTGgc                                                                                                                           >  2:229690/1‑150 (MQ=255)

                              tACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTGTGGATTTTTTTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCC‑GTTTGGTTTTAACAGCGGCGGTTCCTCTTTTTGTTAGCGCCCAGTCTGgcgc                                                                                                                         >  2:5061013/1‑150 (MQ=255)

                              tACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTCGGGGTTTTTGTTAAGCCCGGGTATTAGTAAAGAGAGCCAGGAGCCGGTTTGTTT‑TAAACGCCGCCGGTCCACGTTTGCGTCGCGCCCCCTCGGgcgc                                                                                                                         >  2:591498/1‑150 (MQ=255)

                               aCGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCC‑GTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCgca                                                                                                                        <  1:889348/150‑1 (MQ=255)

                                cGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCC‑GTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCgcag                                                                                                                       <  1:2708916/150‑1 (MQ=255)

                                cGAAAAGTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCC‑GTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCgcag                                                                                                                       <  1:108780/150‑1 (MQ=255)

                                      ttGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCC‑CTTTGTTTATAACAGCGGCGGTTCCTCTTTTGGGTAGCGCCCCGGCTGGCGCAGCACTTc                                                                                                                 >  2:2728550/1‑150 (MQ=255)

                                        gCTACTACGACAATGTTGAACAGAGTTTGATCCCGGTCGTTTCTGTATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCC‑GTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTg                                                                                                               <  1:459035/150‑1 (MQ=255)

                                        gCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAATCCGGGTTATTCGTTAAGATAGCGATGAGCC‑GTTTGTTTATAAAAGCGGCGGTTGCTATTTTTCTTAGCGCCCAGTCTGGCGCAGCACTTcag                                                                                                               >  2:2573973/1‑148 (MQ=255)

                                        gCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCC‑GTTTGATTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTg                                                                                                               <  2:5105213/150‑1 (MQ=255)

                                          tactacGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTGTTGTTTTTTTTTAATCCGGGTTTTTGGTAAATATCGCGAAGAGCC‑GTTTTTTTTTAAAAGCGGCGGTTGCTCTTTTTCTTAGCGCCCAGTGTGGCGCAGCACATATGGt                                                                                                             >  2:1516286/1‑150 (MQ=255)

                                          tactacGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCC‑GTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGt                                                                                                             <  1:2978488/150‑1 (MQ=255)

                                               cGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCC‑GTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTa                                                                                                        <  2:2947464/150‑1 (MQ=255)

                                                gACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGGTATTAGTAAAGATCGCGATGAGCC‑GTTTGTTTATAACAGCGGCGGTTCCTCTTTTTCTTCGCGCCCAGTCTGGCGCAGCACTTATGGTGATTAc                                                                                                       >  2:4277395/1‑150 (MQ=255)

                                                  cAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCC‑GTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTc                                                                                                     <  2:460445/150‑1 (MQ=255)

                                                       ttGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCC‑GTTTGTTTATAACAGCGGCGGTTCCACTTTTGCTTCGCGCCCAGTCTGGCGCCGCACATCTGGTGATTACTCCTGgc                                                                                                >  2:725680/1‑150 (MQ=255)

                                                          aaCAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCC‑GTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCt                                                                                             <  2:3129844/150‑1 (MQ=255)

                                                           aCAGAGTTGGATCCCGGTCGTTTCTGGATTTTTTTTAAGCCGGGTTATTCGTAAATATAGCGATGAGCC‑GGTTGTTTATAACAGCGGCGGGTGCTCTTTTGCTTCGCGCCCACTCTGGCGCAGCACTTATGGGGATTACACCAGgcgcgc                                                                                            >  2:1646994/1‑148 (MQ=255)

                                                            cAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCC‑GTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTcc                                                                                           <  1:4067948/150‑1 (MQ=255)

                                                            cAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCC‑GTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTcc                                                                                           <  2:1916843/150‑1 (MQ=255)

                                                                gTTGGATCCCGGTCGTGTGTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGTGCC‑GTTTGTTTATAAAAGCGGCGGTTCCAATTTTTGTTGGCGCCCACTCTGGCGCAGCACATCTGGTGGTTACTCCTGGCGCTCCgggg                                                                                       >  2:4984113/1‑148 (MQ=255)

                                                                   ggATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCC‑GTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTg                                                                                    <  1:8357/150‑1 (MQ=255)

                                                                    gATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCC‑GTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTgg                                                                                   <  1:930714/150‑1 (MQ=255)

                                                                    gATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCC‑GTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTgg                                                                                   <  2:3251925/150‑1 (MQ=255)

                                                                      tCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCC‑GTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCa                                                                                 <  2:3668084/150‑1 (MQ=255)

                                                                       cccGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCC‑GTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCaa                                                                                <  2:643022/150‑1 (MQ=255)

                                                                            tCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCC‑GTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAag                                                                           >  2:754300/1‑150 (MQ=255)

                                                                             cGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCC‑GTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAaga                                                                          >  2:1143176/1‑150 (MQ=255)

                                                                              gTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCC‑GTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAagag                                                                         >  2:1774998/1‑150 (MQ=255)

                                                                                 tCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCC‑GTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCt                                                                      >  2:2756321/1‑150 (MQ=255)

                                                                                     gATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCC‑GTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTAct                                                                  >  2:2068485/1‑150 (MQ=255)

                                                                                      aTTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCC‑GTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTActc                                                                 >  2:3650359/1‑150 (MQ=255)

                                                                                      aTTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCC‑GTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTActc                                                                 <  1:4115001/150‑1 (MQ=255)

                                                                                       tttttGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCC‑GTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTActct                                                                <  2:3590896/150‑1 (MQ=255)

                                                                                        ttttGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCC‑GTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCtt                                                               >  2:3410093/1‑150 (MQ=255)

                                                                                        ttttGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCC‑GTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCtt                                                               >  2:832601/1‑150 (MQ=255)

                                                                                             ttAAGCCGGGTTATTCGTAAAGATCGCGATGAGCC‑GTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGggc                                                          >  2:3825553/1‑150 (MQ=255)

                                                                                              tAAGCCGGGTTATTCGTAAAGATCGCGATGAGCC‑GTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGggcg                                                         <  2:4580910/150‑1 (MQ=255)

                                                                                               aaGCCGGGTTATTCGTAAAGATCGCGATGAGCC‑GTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGggcgg                                                        >  1:3765160/1‑150 (MQ=255)

                                                                                               aaGCCGGGTTATTCGTAAAGATCGCGATGAGCC‑GTTTGTTTATAACAGCGGCGGGTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGggcgg                                                        <  2:3238170/150‑1 (MQ=255)

                                                                                                  ccGGGTTATTCGTAAAGATCGCGATGAGCC‑GTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCtt                                                     >  1:172820/1‑150 (MQ=255)

                                                                                                  ccGGGTTATTCGTAAAGATCGCGATGAGCC‑GTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCtt                                                     <  1:1839194/150‑1 (MQ=255)

                                                                                                   cGGGTTATTCGTAAAGATCGCGATGAGCC‑GTTTGTTTCTAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCttt                                                    <  2:4042562/150‑1 (MQ=255)

                                                                                                    gggTTATTCGTAAAGATCGCTATGAGCC‑GTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTa                                                   <  1:3520677/150‑1 (MQ=255)

                                                                                                        tagtCGTAAAGATCGCGATGAGCC‑GTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTaccnc                                               <  1:2832305/147‑1 (MQ=255)

                                                                                                            cGTAAAGATCGCGATGAGCC‑GTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCaaa                                           >  1:2107391/1‑150 (MQ=255)

                                                                                                            cGTAAAGATCGCGATGAGCC‑GTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCaaa                                           >  1:553494/1‑150 (MQ=255)

                                                                                                            cGTAAAGATCGCGATGAGCC‑GTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCaaa                                           >  1:580383/1‑150 (MQ=255)

                                                                                                                aaGATCGCGATGAGCC‑GTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAg                                       <  1:3927050/150‑1 (MQ=255)

                                                                                                                 aGATCGCGATGAGCC‑GTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGc                                      >  1:2377080/1‑150 (MQ=255)

                                                                                                                     cgcgATGAGCC‑GTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCt                                  >  2:4005690/1‑150 (MQ=255)

                                                                                                                     cgcgATGAGCC‑GTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCt                                  <  2:2403886/150‑1 (MQ=255)

                                                                                                                       cgATGAGCC‑GTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTAc                                <  2:4873782/150‑1 (MQ=255)

                                                                                                                        gATGAGCC‑GTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTAcc                               <  1:2156411/150‑1 (MQ=255)

                                                                                                                        gATGAGCC‑GTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTAcc                               >  2:3786066/1‑150 (MQ=255)

                                                                                                                         aTGAGCC‑GTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCg                              >  2:1895737/1‑150 (MQ=255)

                                                                                                                           gAGCC‑GTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGtt                            >  2:1303578/1‑150 (MQ=255)

                                                                                                                            aGCC‑GTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGttt                           <  2:2348682/150‑1 (MQ=255)

                                                                                                                               c‑GTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTtctc                        <  2:4914530/150‑1 (MQ=255)

                                                                                                                               c‑GTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTtctc                        >  1:2882392/1‑150 (MQ=255)

                                                                                                                                     gtttATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCgcggc                   <  1:2655167/150‑1 (MQ=255)

                                                                                                                                      tttATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCgcggcg                  >  1:1472536/1‑150 (MQ=255)

                                                                                                                                         ataACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCt               >  2:1283413/1‑150 (MQ=255)

                                                                                                                                         ataACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCt               >  1:974707/1‑150 (MQ=255)

                                                                                                                                          taACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTa              <  1:1606452/150‑1 (MQ=255)

                                                                                                                                           aaCAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTAt             <  1:691172/150‑1 (MQ=255)

                                                                                                                                            aCAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATg            >  2:658486/1‑150 (MQ=255)

                                                                                                                                            aCAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATg            <  1:3019639/150‑1 (MQ=255)

                                                                                                                                              aGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGt          >  1:4479105/1‑150 (MQ=255)

                                                                                                                                                cggcggTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCg        >  1:4811263/1‑150 (MQ=255)

                                                                                                                                                cggcggTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCg        <  1:1143176/150‑1 (MQ=255)

                                                                                                                                                 ggcggTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGt       >  2:4685166/1‑150 (MQ=255)

                                                                                                                                                 ggcggTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGt       >  2:2131616/1‑150 (MQ=255)

                                                                                                                                                 ggcggGTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTctt       <  1:1519560/150‑3 (MQ=255)

                                                                                                                                                   cggTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGttc     <  2:1975546/150‑1 (MQ=255)

                                                                                                                                                   cggTTACTCTTTTTCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCATATAATCGTCTACCGTTTCTCGCTGCGGCTATGGTCGttc     <  1:3320935/150‑1 (MQ=255)

                                                                                                                                                     gTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGttctt   >  1:949415/1‑150 (MQ=255)

AATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGATTTTTGTTAAGCCGGGTTATTCGTAAAGATCGCGATGAGCC‑GTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGAGTCTTACTCTTCGGGCGGCTTTACCACCAAAAAAGCGTCTACCGTTTCTCGCGGCGGCTATGGTCGTTCTTC  >  CP009273/4402516‑4402814

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 11 ≤` `ATCG/ATCG` `< 25 ≤` `ATCG/ATCG` `< 37 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`