BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	1,282,578	T→G	13.1%	noncoding (115/171 nt) R19R (AGG→CGG)	rttR ← tpr ←	rtT sRNA, processed from tyrT transcript protamine‑like protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	1,282,578	0	T	G	13.1%	268.1 / 2.9	127	noncoding (115/171 nt) R19R (AGG→CGG)	rttR tpr	rtT sRNA, processed from tyrT transcript protamine‑like protein
Reads supporting (aligned to +/- strand): ref base T (32/69); new base G (17/2); total (55/72)
Fisher's exact test for biased strand distribution p-value = 2.91e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 128 total aligned reads displayed.

GCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCT‑CATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTT  >  CP009273/1282429‑1282726

gCCGAACCTTAGTCGTAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCt                                                                                                                                                       <  2:333430/150‑1 (MQ=255)

 ccGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTg                                                                                                                                                      >  2:1219567/1‑150 (MQ=255)

  cGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGc                                                                                                                                                     >  2:3400462/1‑150 (MQ=255)

  cGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCAGCGCGCTTCGCCCTTCGGGTCGTTGCCTgg                                                                                                                                                     >  2:3855903/1‑149 (MQ=255)

     aCCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTGGGCTGCGGc                                                                                                                                                  >  1:3062991/1‑150 (MQ=255)

             cGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGctct                                                                                                                                          <  1:3802551/150‑1 (MQ=255)

              gAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCCTCTGCGCTGCGCGCGTCGCCCTGCGGGTCGTTGGCGGGGGCAACGcacta                                                                                                                                         >  2:3051638/1‑149 (MQ=14)

              gAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTAAGCGCGGAGAGCTTCGCCCTACGGGTCGTGGGCTGTGGCCAAGctaga                                                                                                                                         >  2:4035505/1‑147 (MQ=14)

                 gCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGAGCTGCGCGCTTCGCCCATCGGGTCGTTGCGGGCGGCCACGccagctcg                                                                                                                                      >  1:3057222/1‑143 (MQ=14)

                    tctcATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTg                                                                                                                                   <  1:3146017/150‑1 (MQ=37)

                     ctcATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACGCACTCTCGCTgg                                                                                                                                  >  1:2829678/1‑150 (MQ=255)

                     ctcATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGGGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTgg                                                                                                                                  >  2:940065/1‑150 (MQ=255)

                      tcATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAAAGCGCAATCGCAGgc                                                                                                                                 >  1:3231084/1‑150 (MQ=14)

                      tcATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGACGTCAGGGCCTTTTGAAAGTGATGGTGGTGGGGGTAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGgc                                                                                                                                 <  2:1684988/150‑1 (MQ=255)

                           cTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGATTTCAGGGCCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCg                                                                                                                            <  2:1056126/150‑1 (MQ=255)

                            ttCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGAGCTTCGCCCTTCGGGTCGTGTCAGGACCCAACGCCCGCTCGCTGgcggacga                                                                                                                           >  2:1671995/1‑145 (MQ=14)

                             tCCCCGGATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGTCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAg                                                                                                                          <  2:2621326/150‑1 (MQ=255)

                              ccccGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGTATTACTCAGCGCTGCGCGCTTCGACCTTCGGGTCGTTTCCGGCTGCAACACTCTATAGCAGGCCCACGaat                                                                                                                         >  2:477924/1‑148 (MQ=14)

                              ccccGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCCCTTCGCCCGTCGGGTCCTTTCCGGCCGCACCGCGCTCACGCAGGCGCTCGAGt                                                                                                                         >  2:2057924/1‑150 (MQ=14)

                                 cGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCACCCTTCGGGTCGTTGCCTGAGGCAACGCTCTCTCGCTTGCGCCCGAATCGa                                                                                                                      >  2:1381449/1‑150 (MQ=14)

                                  gCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGaa                                                                                                                     <  1:3400462/150‑1 (MQ=255)

                                    aTGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAAcc                                                                                                                   <  1:2456384/150‑1 (MQ=255)

                                      gggCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCtt                                                                                                                 <  2:2867936/150‑1 (MQ=255)

                                         cAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCCACGCCCACACGCAGGGCCGCGAGTCGAAcccgggg                                                                                                              >  2:69801/1‑145 (MQ=14)

                                          aGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGTCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTc                                                                                                             <  2:970651/150‑1 (MQ=255)

                                           gAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCACTCTCGCTGTCGCTCGAGTCGAACCTTGGTCg                                                                                                            >  2:1784908/1‑150 (MQ=255)

                                             atatTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTAAGAGCTGCGCGCTTCGCCCTTCGGGTCGTTTCGTGAGGCAACGCGCTCTCGCTGGCCCAAGCGTCGAACCTTGGTAGaa                                                                                                          >  2:76639/1‑150 (MQ=14)

                                                tttGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCt                                                                                                       <  2:1613537/150‑1 (MQ=255)

                                                 ttGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTGGCCGGCGGCAAACCACACACGCCGGCGCGCGAGTCGAACCAAGGTCGaaacat                                                                                                      >  1:2455707/1‑146 (MQ=14)

                                                  tGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCCGCGGCAACGCACACACGCAGGCCCCCCAGTCGCACCCAGGTCGaaaacac                                                                                                     >  1:3545060/1‑145 (MQ=14)

                                                   gATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCTCATCGCCCTTCGGGACGTTTCCAGCGGCAACCCTCTCTCGCTGGCGCACGAGTCGAACCGCGGACGAATCTtca                                                                                                    >  2:930578/1‑149 (MQ=14)

                                                    aTTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTTCCGGCCGCAACCCACACTCGCGGGCGCACCAGTCGAACCATGGTAGAACCTtcaa                                                                                                   >  1:3945557/1‑148 (MQ=14)

                                                     ttGCGGATTCGATTGAGATTTCAGGGTCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCT‑Ca                                                                                                 <  1:62002/150‑1 (MQ=255)

                                                      tGCGGCTTCGATTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCT‑CAt                                                                                                <  2:1083265/150‑1 (MQ=255)

                                                         ggATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCACTCTCGCTGGCGCTCGACTAGAAACTTGGTAGAAGCATAT‑CAtacg                                                                                             >  2:1325659/1‑147 (MQ=14)

                                                             tCGCTTGAGAGTTCTGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGACCCTTGGTCGAAGCTTCT‑CATCCTTccc                                                                                         <  1:714932/150‑1 (MQ=255)

                                                              cgattGAGACTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCT‑CATCCTTcccc                                                                                        <  1:2054227/147‑1 (MQ=255)

                                                               gCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTAATCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCCACCCTCTCTAGCGGGCGCTCGATTCGTAACTTGGTCGAATCATCT‑CATCCATCCCCg                                                                                       >  2:3786361/1‑150 (MQ=14)

                                                                cTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGTCGGCAACGCACTCTCGCTGGCGCTCGAGTAGAACCTTGGGCGAAGCTTCT‑CATACTTCCCCGc                                                                                      >  2:4298429/1‑150 (MQ=14)

                                                                 ttGAGAGTTCAGGGACTTTTGAAAGTTATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCCACGCTCGCTCGCGTGCGCACCAGGCGAACCATTTGCGAAGCTTCG‑CAACCGTCCCCGcc                                                                                     >  2:3507195/1‑149 (MQ=14)

                                                                  tGAGAGTTCTGGGACTTTGGTAAGTGTTGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCGTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTAGAACCTTGGTCGAAGCTTCT‑CATCCTTCCCCGCAt                                                                                    <  1:1234129/150‑1 (MQ=14)

                                                                  atcgagTTCCGGGTCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCT‑CATCCTTCCCCGCAt                                                                                    <  1:196204/147‑1 (MQ=255)

                                                                   gagaGTTCAGGGACTTTTGTAAGTGATGGTGGTGGGGGATGGTTTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCT‑CATCCTTCCCCGCATg                                                                                   <  1:2155527/150‑1 (MQ=255)

                                                                   gagaGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCCGCGGCAACGCTCACACGCTGGCGCACGAGTCGAACCATGGGCGAAGATTCT‑CAACCATCCCCCCagg                                                                                   >  2:3883190/1‑148 (MQ=14)

                                                                     tagTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCACTCTCGCTGGCGCTCGAGTCGAACCGTGGTCGAAGCTTCT‑AATCCATCCCCCCATggg                                                                                 >  2:1767984/2‑150 (MQ=14)

                                                                      cgTTCCGGGGCTTTTGAAATGGATGGTGGTGGGGGTAGGATTACTCTGGGCTGCGCGCTTCGCCCTTCTGGTCGTTGCCTGCGGCCACGCTCTCTCGCTGGCGCTCGAGTCGATCCTTGGTCGAAGCTTCT‑CATCCTTCCCCGCATGGGc                                                                                <  1:3834255/149‑1 (MQ=11)

                                                                        ttCCGGGACTATGGTATGTGATGGTGGTGGGGGTTTGTTTAGTCTGGGGTGCGGGCGTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCT‑CATCCTTCCCCGCATGGGCAg                                                                              <  1:3335798/150‑1 (MQ=11)

                                                                         tCAGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCT‑CATCCTTCCCCGCATGGGCAGa                                                                             <  1:2857833/150‑1 (MQ=37)

                                                                         tCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCT‑CATCCTTCCCCGCATGGGCAGa                                                                             <  2:3684037/150‑1 (MQ=255)

                                                                         tCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCT‑CATCCTTCCCCGCATGGGCAGa                                                                             <  2:3683973/150‑1 (MQ=255)

                                                                         tCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCT‑CATCCTTCCCCGCATGGGCAGa                                                                             <  2:3683873/150‑1 (MQ=255)

                                                                           aGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGAGTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCT‑CATCCTTCCCCGCATGGGCAGAat                                                                           <  1:2479658/150‑1 (MQ=255)

                                                                            gggTCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCT‑CATCCTTCCCCGCATGGGCAGAata                                                                          <  1:1699327/150‑1 (MQ=37)

                                                                            gggACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTGTCCGGGGGCAACGCTCACACGCTGGCGCGCGAGTAGTAACGTGGTCGCAGCTTCT‑AAACCTTACCCCCAAGGGCCGaaaa                                                                          >  2:3815917/1‑148 (MQ=14)

                                                                              gACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCT‑CATCCTTCCCCGCATGGGCAGAATAtt                                                                        <  1:463410/150‑1 (MQ=38)

                                                                                   ttGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCT‑CATCCTTCCCCGCATGGGCAGAATATTTAAtt                                                                   <  1:1498720/150‑1 (MQ=39)

                                                                                         gTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCT‑CATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGAt                                                             <  1:393481/150‑1 (MQ=38)

                                                                                          cgATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCCAGACGAACCTTTGTCGAAGCTTCT‑CATACTTCCCCCCATGGGCCGAATATTTTATTTCGGAtt                                                            >  2:3081402/2‑150 (MQ=11)

                                                                                              ggtggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCT‑CATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGtt                                                        <  1:2181912/150‑1 (MQ=37)

                                                                                                tggtggGGGAAGGATTTCTAAGCGATGCTGTCTTCGAGGTTCGGTTCGGTGTGAGCGGCACAGCTCTCTCGCTGTCGCTCGAGTCGAACCTTGGTCGAAGCTTCT‑CATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTgg                                                      <  1:831854/150‑1 (MQ=11)

                                                                                                tggtggGGGAAGGATTACTCAGCGCTGCTCTCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCT‑CATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTgg                                                      <  1:831793/150‑1 (MQ=255)

                                                                                                  gtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGAGGCAAAGCTCTCGCGCTGGCGCTCGAGTCGAACCTTTGTGGAAGCGTCT‑CATACTTAACCCCATGGGCAGAATAATTAATTGAGGGTtgggtgggc                                                    >  2:498234/1‑149 (MQ=11)

                                                                                                     ggggAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCT‑CATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGt                                                 <  1:3535668/150‑1 (MQ=255)

                                                                                                      gggAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCT‑CATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGtt                                                <  1:4096809/150‑1 (MQ=255)

                                                                                                      gggAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCT‑CATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGtt                                                <  1:3198644/150‑1 (MQ=255)

                                                                                                       ggAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCT‑CATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTc                                               <  1:2456117/150‑1 (MQ=255)

                                                                                                       ggAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACGCACTCTCGCTGGCGCTCGAGTAGCACCCTGGTCGAAGCTTATCCAACATACCCCCCA‑GGGCAGAATAATTAAATGCGGATTCGGTGGGAAGGTc                                               >  2:3544009/1‑150 (MQ=14)

                                                                                                        gAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTGGCCTGGGGCAACGCTCTCTCGCTGGCGCGCGTGTCGAACCTGGGTCGAAGCGTCT‑CATACTTGCCCGCATGGGGAGAAAATTTAATTGGGGATTGGTTGGGAAGttaa                                              >  2:2351308/1‑148 (MQ=14)

                                                                                                         atggATTTCTCAGCGCTGCGCGCGTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCT‑CATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAg                                             <  1:3666810/148‑1 (MQ=255)

                                                                                                         aaGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCT‑CATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAg                                             <  1:735586/150‑1 (MQ=255)

                                                                                                         aaGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCT‑CATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAg                                             <  1:3029009/150‑1 (MQ=255)

                                                                                                          aGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCT‑CATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAgg                                            <  1:940065/150‑1 (MQ=255)

                                                                                                           ggATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCT‑CATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAggg                                           <  2:3187337/150‑1 (MQ=255)

                                                                                                            gATTTCTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCT‑CATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGa                                          <  1:380285/150‑1 (MQ=255)

                                                                                                            gATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCT‑CATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGa                                          <  1:511577/150‑1 (MQ=255)

                                                                                                             aTTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCCAAGCTCTCTCGGTGGCGCGCGCGTCGAACCGGGGTCGCAGCGTGT‑AATGCTTACCCGCCGGGGGAGCATAATTAATTGCGGATTGGTTGGGGAGttgggggga                                         >  2:2331262/1‑143 (MQ=2)

                                                                                                             aTTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCACTCTCGCTGGCGCGCGAGTCGAAACATGGTCGTAGCATCT‑CATACTTACCCCCAAGGGCAGCAAATTTAATTGCGGGTTCGTTTGGTAGTTCAGGGAc                                         >  2:3569546/1‑150 (MQ=14)

                                                                                                             aTTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAAAGCTCTCTCGCTGGCGCTCGGGTCGAACCTTGGTGGAAGCTTCA‑CAAACTTACCCGCCGGGGCCGAATAATTAAATTGGGATTAGTTGGGGAGTTGAgggtc                                         >  2:2331101/1‑148 (MQ=14)

                                                                                                                 cTCAGCGCTGCTCGCTTCTCCATTCGGGTCGTTGCCTGCGTCTAGGCTCTCTCGCTGTCGCTCGTGTAGATCCTTTTTCGAAGCTTCT‑CATCCTTCCCGGCTTGGGCAGAATTTTTAATTGCGGATTCGTTGGGAAGTTCAGGGACtttt                                     <  1:2112510/150‑1 (MQ=14)

                                                                                                                      cgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGGGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTTGTCGAAGCTTCT‑CAACCTTACCCGCATTGGGCGAAAATTTAAATTCGGATTCGTGGTGACGGTAAGGGTCTTTTGaaat                                >  2:2501902/1‑149 (MQ=14)

                                                                                                                        cTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCT‑CATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTg                              <  1:1604861/150‑1 (MQ=255)

                                                                                                                        cTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCGCTCTCGCGGGCGCGCGAGGCGAACCTTTGTCGAAGCGTCT‑CAACCTTCCCCGCCGGGGCAGAATATTTTATTTGGGCTTGGGTGGGTAGTTAAGGGGCTTTTGTTAGTg                              >  2:4310267/1‑150 (MQ=2)

                                                                                                                             cgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCT‑CATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGAtggt                         <  1:614478/150‑1 (MQ=255)

                                                                                                                                  cGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCT‑CATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGAtggtggtgg                    <  2:3368338/150‑1 (MQ=255)

                                                                                                                                   gCCCTTCGGGTCGTTGCCAGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCATGGTCGAATCTTCT‑CAACCATCCCCGCCAGGGCAGAATATTTAATTGCGGAGTCGGTGGGTAGTTCAGGGGCCTTTGTAAGTGATGGTGGTggg                   >  2:2888056/1‑150 (MQ=14)

                                                                                                                                    cccTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCTCTCGAGTCGTACCTTGGTCGAAGCTTCT‑CATCCTTCACCGCTTGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTgggg                  <  1:1826543/150‑1 (MQ=255)

                                                                                                                                     ccTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCT‑CATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTggggg                 <  2:3511758/150‑1 (MQ=255)

                                                                                                                                       ttCGGGTCGTTGCCTGCGGCAACCCACTATCGCGGGCGCGCGAGTCGAACCGTGGTCGCAGCTTCT‑GAAACATACCCCCCTGGGCAGAAAATTTTATGGCGGGTTGGTTGGGGAGTTGCGGGACTGTTGGAAG‑GGGGGTgggggggggaa              >  2:4266941/1‑147 (MQ=2)

                                                                                                                                        tCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCT‑CATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAg              <  1:428946/150‑1 (MQ=255)

                                                                                                                                        tCGGGTCGGTGCCGGCGGCAACGCTCACTCGCTGTGGCTCGAGTCGATCACTGGTGGCAGCTTGT‑CAGACGTCCCCGCATGGGCAGAATATTTTATTGCGGATTCGTTGGTCAGTTCAGGGTCTTTTGAAAGTGATGGTGGTGGGGGAAg              <  1:3632519/150‑1 (MQ=14)

                                                                                                                                             tCGTTGCCTTCGGAAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCT‑CATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTc         <  2:3210824/150‑1 (MQ=255)

                                                                                                                                             tCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCT‑CATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTc         <  1:3051638/150‑1 (MQ=255)

                                                                                                                                             tCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCT‑CATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTc         <  2:3210349/150‑1 (MQ=255)

                                                                                                                                             tCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTTGTCGAAGCTTCT‑CATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTc         <  1:440527/150‑1 (MQ=255)

                                                                                                                                              cGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCT‑CATACTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCg        >  2:1919187/1‑150 (MQ=255)

                                                                                                                                              cGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTAT‑CATCCTTCCCCGCATGGGCAGAATATTTAATTGGGGATTCGTTGGGTAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGGAGGAttag        >  1:1071587/1‑148 (MQ=255)

                                                                                                                                              cGTGGCCGGCGGCACCGCTCTCTCGCGGGCGCACGAGTCGAACCCTGGTCGAAGCTTCC‑CATCCTTCCCCGCATGGGCAGAATACTTTCTTGCGGATTCGTTGGGCAGTTCAGGGCCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCg        <  2:1784728/150‑1 (MQ=14)

                                                                                                                                               gTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCT‑CATCCTTCCCCGCATGGGCAGAATTTTTAATTGCGGATTCGTTGGGTAGTTGAGGGAGTTTTGAAAGTGTTGGTGGTGGGGGAAGGATTCgt       >  2:2587196/1‑149 (MQ=17)

                                                                                                                                               gTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCT‑CATCCTTCCCCGCATGGGCAGAATCTGTAATTGCGGATTCGTTGGGACGTTCAGGGACTTTTGAAAGTGAGGGGGGTGGGGGAAGGATTCGa       >  2:1272580/1‑150 (MQ=17)

                                                                                                                                                    cTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCT‑CATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCtt  <  1:4220317/150‑1 (MQ=255)

GCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCT‑CATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTT  >  CP009273/1282429‑1282726

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