BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,287,122	T→G	7.5%	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
*	NC_000913	1,287,122	0	T	G	7.5%	1435.6 / 45.4	622	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 (268/262); new base G (54/4); total (350/272)
Fisher's exact test for biased strand distribution p-value = 3.59e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 636 total aligned reads displayed.

GCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTC  >  NC_000913/1286973‑1287271

gCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCt                                                                                                                                                       >  1:4256364/1‑150 (MQ=255)

gCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCt                                                                                                                                                       <  1:6231052/150‑1 (MQ=255)

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

 ccGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTACGGGTCGATTCCTg                                                                                                                                                      >  2:484673/1‑150 (MQ=255)

 ccGAACCTTAGTCGAAGCTTCTCATCCGTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTg                                                                                                                                                      <  2:9043445/150‑1 (MQ=255)

   gAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCg                                                                                                                                                    <  2:12854961/150‑1 (MQ=255)

           gTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCGGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGct                                                                                                                                            >  2:6142950/1‑150 (MQ=255)

            tCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTTCCAGCGGCAACGcac                                                                                                                                           >  1:11265740/1‑148 (MQ=255)

              gAAGCTTCTCATCCTTCCCCGCATGGGCAGACTATTTGATTGCGGATTCGCTTGAGCGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGCTTACAAAGCTCTGAGCGCTTGACCCATCGGGACGTGGCCAGACGCAACGcacgc                                                                                                                                         >  1:10039250/1‑146 (MQ=14)

              gAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACGcactc                                                                                                                                         >  1:12762339/1‑150 (MQ=255)

               aaGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGctctct                                                                                                                                        >  1:13968307/1‑150 (MQ=255)

                     ctcATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCGGCGCGCTTCGCCCTTCGGGTCGTTGCCAGCGGCCACCCTATCACGCggg                                                                                                                                  >  1:3202930/1‑150 (MQ=14)

                      tcATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCACCGCTCTCTCGCTGgc                                                                                                                                 >  1:2049021/1‑150 (MQ=255)

                      tcATCCTTCCCCGCATGGGCAGAATATTTGAGTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCCACGCTCTCTCGCTGgc                                                                                                                                 >  2:5250864/1‑150 (MQ=255)

                       cATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGAGGCAAAGCGCTCTCGCGGgcc                                                                                                                                >  2:14563847/1‑149 (MQ=255)

                        aTCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGgcgc                                                                                                                               >  1:10166645/1‑150 (MQ=255)

                         tCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGGCGGCGGCAACGCTCACAAGCTGGCGCt                                                                                                                              >  1:6656590/1‑150 (MQ=14)

                         tCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCt                                                                                                                              >  1:1224582/1‑150 (MQ=255)

                            ttCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTTCCTGAGGCAAAGCTCTCTCGCTGgcgcgaga                                                                                                                           >  1:5033118/1‑146 (MQ=255)

                            ttCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCGCTCTCGCTGGCGCGCGa                                                                                                                           >  2:6687845/1‑150 (MQ=255)

                              ccccGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGt                                                                                                                         <  1:3561360/150‑1 (MQ=255)

                               cacGCATGTTCAGAATATTTGATTGCGGATTTGCTTGAGATTTCCGGGTCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTc                                                                                                                        <  1:12675742/148‑1 (MQ=11)

                                ccGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCGCGAGTCg                                                                                                                       >  1:6400746/1‑150 (MQ=255)

                                 cGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGGGGCAACGGTCTCTCGCTTGCGCTCGAGTCGa                                                                                                                      >  2:8544064/1‑150 (MQ=255)

                                 cGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCACCGCTGCGAGCTTCGCCCTTCGGGTCGTTGCATGCCGCAACGCTCTCTCTCTGGCGCGCGAGACGa                                                                                                                      >  2:7760325/1‑150 (MQ=14)

                                  gCATGGGCAGAATATTTGATTGCGGATTCGATTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGaa                                                                                                                     <  2:13034908/150‑1 (MQ=255)

                                       ggCAGAATATTTGATTGCGGATTCGTTTGAGAATTCAGGGTCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTg                                                                                                                <  1:4683464/150‑1 (MQ=12)

                                        gCAGAATATATGATTGCGGATTCGCTAGAGAGTTCAGGGACTTTTCAAAGTGATGGTGGTGGGGGACGGAATACTCACCGCTCCGCCCATCGCCCATCGGGACGTAGCCTCAGGCAACGCTCTCTCGCTGGCCCTCGAGTCGAAccaagg                                                                                                               >  2:10651839/1‑146 (MQ=18)

                                              tatTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAg                                                                                                         <  1:9045294/150‑1 (MQ=255)

                                              tatTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAg                                                                                                         <  1:5331025/150‑1 (MQ=255)

                                                  tGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTtc                                                                                                     <  2:3553072/150‑1 (MQ=255)

                                                     ttGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCACTCACGCTGGCGCTCGAGTCGAACCATGGTCGAAGCTTATCa                                                                                                  >  1:11472565/1‑150 (MQ=255)

                                                      tGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACACTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCAt                                                                                                 >  2:10407945/1‑150 (MQ=255)

                                                          gATTCGTTTGAGATTTCAGGGTCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCtt                                                                                             <  2:14401173/150‑1 (MQ=255)

                                                            ttCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGCGTCGAACCATGGTCGAAGCTTCTCATccatac                                                                                           >  1:1843293/1‑146 (MQ=255)

                                                                cTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCaa                                                                                                                                                >  1:13421440/1‑93 (MQ=255)

                                                                 ttGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCAGCGGCAACCCTCACTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCAAACTTCCCCCCa                                                                                      >  1:1097179/1‑150 (MQ=14)

                                                                      agTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCAGCGGCAAAGCTCTCACGCAGGCGCTCCAGTCGAACCTTGGTCGAAGCATATCAAACTTACCCGCATGGGc                                                                                 >  1:10340550/1‑150 (MQ=14)

                                                                       gTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCa                                                                                <  2:1927553/150‑1 (MQ=255)

                                                                       gTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAATCGAACCTTGGTCGAAGCTTCTCAACCTTACCCCCATGGGCa                                                                                >  1:15330537/1‑150 (MQ=255)

                                                                        ttCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAg                                                                               <  1:9915742/150‑1 (MQ=255)

                                                                         tCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTg                                                                                                                                                      <  1:8050132/78‑1 (MQ=39)

                                                                          cAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGaa                                                                             <  1:6870657/150‑1 (MQ=40)

                                                                           aGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAat                                                                            >  2:824022/1‑150 (MQ=40)

                                                                           aGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACGCACACAAGCTGGCGCACCAGTCGCACCAGGGTCGAAGCATATCAACCATACCCCCCAGGGCCGAat                                                                            >  1:3897987/1‑150 (MQ=14)

                                                                           aGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCCGCCACCCAATCACCCAGGCGCGCCACTCGAACCATGGTCGAAACATCTCAAAAATACCCGCAAGGGCAGAat                                                                            >  1:11792646/1‑150 (MQ=9)

                                                                            gggACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAata                                                                           <  1:5250864/150‑1 (MQ=39)

                                                                            gggACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCCAACCTTGGTCGAAGCTTCTCAACCTTCCCCCCAAGGGCAGaaaa                                                                           >  2:4413313/1‑148 (MQ=255)

                                                                            gggACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCACCAGTCGAACCGTGGTCGAAGCTTCTCAACCGTCCCCGCATGGGCAGAata                                                                           >  2:3721048/1‑150 (MQ=14)

                                                                            gggACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCGCTCTCGCTGGCGCGCGAGGCGAACCTTGGTCGGAGCGTCTCAACCTTACCCGCAGGGGGCGAata                                                                           >  1:13183883/1‑150 (MQ=11)

                                                                             ggACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCAACCTTCCCCGAATAGGCACAaaat                                                                          >  2:14029577/1‑150 (MQ=14)

                                                                                 ttttGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCAACCTTCCCCGCATGGGCAGAATATTTaa                                                                      >  2:6499243/1‑150 (MQ=38)

                                                                                  tttGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCAGCGGCAAAGCTCTCTAGCATGCGCTCGAGTCGCACCTTGGTAGTAGCATCTAAAACATGCCCCCAAGGGGAGAAAAATTaaa                                                                     >  1:339087/1‑149 (MQ=2)

                                                                                      aaaGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCg                                                                 <  1:96371/150‑1 (MQ=39)

                                                                                           gATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTc                                                            <  2:9193746/150‑1 (MQ=37)

                                                                                           gATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTc                                                            <  1:2648828/150‑1 (MQ=37)

                                                                                           gATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTc                                                            <  2:1565149/150‑1 (MQ=37)

                                                                                               gtggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTTCCTGCGGCAACGCACTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTg                                                        >  2:10871214/1‑150 (MQ=255)

                                                                                               gtggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTTGCGCTCGACTCGAACCTTGGTCGAAGCTTCTCAACCTTCCCCGCAGGGGGAGAATAATTAAATGCGGAttagatg                                                        >  2:3299395/1‑145 (MQ=11)

                                                                                               gtggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTTGTCGAAGCTTCTCATACATACCAACATGGGCATAATAGTAAATTGAGGagtagatg                                                        >  2:8743267/1‑143 (MQ=11)

                                                                                               gtggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTg                                                        <  1:14839113/150‑1 (MQ=38)

                                                                                                   ttggggTAGGATTTCTCTGCGCTGCGCGCGTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGaa                                                    <  2:9726626/148‑1 (MQ=255)

                                                                                                   ttggggAAGGATTACTCAGCGCTGCGCGATTCGCCATTCGGGTCGTTGCCAGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGaa                                                    <  1:1942540/148‑1 (MQ=255)

                                                                                                   ggggggTTGGATTCCTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGGTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGaa                                                    <  1:13677682/149‑1 (MQ=255)

                                                                                                     ggggAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGt                                                  >  2:14929655/1‑150 (MQ=255)

                                                                                                       ggAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTc                                                                                            >  2:14172732/1‑106 (MQ=32)

                                                                                                       ggAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCAGCGGCAACGCACTCTCGCTGGCGCGCGAGTCGAACCATGGTCGAAGCTTCACAACCTTACCCGCAGGGGCAGAATATTTAATTGCGGATTAGTTGGGAAGTTc                                                >  2:12115858/1‑150 (MQ=14)

                                                                                                         taGGATTACTCAGCTCTGCGCGCTTCGCCCTTCGGGCCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAg                                              <  2:12386023/149‑1 (MQ=255)

                                                                                                         aaGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAg                                              <  1:14547353/150‑1 (MQ=255)

                                                                                                            gATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCGTGGTCGAAGCTTCTCATCCTTCCCCGCCTGGGCCGCATATTTAATTGCGGCTTAGTTGGGGAGTTGAgggc                                           >  2:2845032/1‑149 (MQ=14)

                                                                                                              ttACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACt                                         <  2:11226663/150‑1 (MQ=255)

                                                                                                                  tCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTg                                     >  2:2293987/1‑150 (MQ=255)

                                                                                                                    aGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGaa                                   >  2:13811370/1‑150 (MQ=255)

                                                                                                                     gcgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCATAATATTTAATTGCGGATTCGTTGGGAAGTTCAGTGACTTTTGaaa                                  <  1:8320790/150‑1 (MQ=255)

                                                                                                                     gcgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGaaa                                  <  1:9362568/150‑1 (MQ=255)

                                                                                                                        cTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACGCTCTCGCGCTGGCGCGCGAGTCGAACCGTGGTCGAAGCTTCTCATCCTTCCCCGCAGGGGCAGCATATTTAATTGCGGGTTCGGTGGGGAGTTCAGGGACTTTGGAGAGTg                               >  2:1529671/1‑150 (MQ=14)

                                                                                                                         tGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTAGCTGGCGCTCGAGTCGAACCGTGGTCGAAGCTTCTCAAACATCCCCGCAGGGGCAGAATAGTTAATTTCGGATTGGTTGGGAAGTTGAGGGAGGTTTGAAAGTgg                              >  2:11765728/1‑149 (MQ=17)

                                                                                                                         tGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCGCTCTCGCTGGCGCTCGAGTCGGACCGTGGTCGAAGCTTCTCATCCTTCCCCGCCTGGGCCGGCTATTTAATTGCGGCTTCGGTGGGGAGTTCGGGGACGTTTGACAGTgg                              >  2:1780924/1‑149 (MQ=14)

                                                                                                                             cgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGCCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGAtggt                          <  2:13623891/150‑1 (MQ=255)

                                                                                                                                ttCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGAtggtggt                       <  1:15062861/150‑1 (MQ=255)

                                                                                                                                ttCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCCGAATATTTAATTGCGGATTGGTTGGGTCGTTCATGGACTTTTGGAAGTGAtggtggt                       >  2:747482/1‑150 (MQ=14)

                                                                                                                                  cGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGAtggtggtgg                     <  1:1201571/150‑1 (MQ=255)

                                                                                                                                    cccTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTTGCGCACGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTCCAGGGACTTTTGAAAGTGATGGTGGTgggg                   >  2:11459555/1‑150 (MQ=255)

                                                                                                                                       ttCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGaa                <  1:14151537/150‑1 (MQ=255)

                                                                                                                                       ttCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATACATACCCGCATTGGCAGAATATTTAATTGCGGATTAGTTGGGAAGTTAAGGGACTTTTGAAAGTGATGGGGGTGGGGGac                >  2:13426961/1‑149 (MQ=17)

                                                                                                                                           gagcGTTGCCGCAGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAg               <  1:12420582/144‑1 (MQ=255)

                                                                                                                                            agcGTTGCCGCAGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGa                                    >  1:1060234/3‑125 (MQ=255)

                                                                                                                                             tCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAATCTTCTCAACCTTCACCGCATGGGCAGAATATTAAATTGCGGAATCGGTTTGAAGATCAGGGTCTTGTGAAAGTGATGGTGGTTGGGGAAGGAGTc          >  2:14739970/1‑150 (MQ=14)

                                                                                                                                             tCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTc          <  1:12524442/150‑1 (MQ=255)

                                                                                                                                              cGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACtttt                                      <  1:776690/121‑1 (MQ=255)

                                                                                                                                               ggtGCCTGCGGCACCGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGa        <  1:7464988/148‑1 (MQ=255)

                                                                                                                                               gTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTcccc                                                                                         >  2:4052827/1‑69 (MQ=35)

                                                                                                                                               gTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATAATTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGACAGTGATGGTGGTGGGGGAAGGATTCGa        >  2:468148/1‑150 (MQ=255)

                                                                                                                                                ttGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGaa       >  2:7050413/1‑150 (MQ=255)

                                                                                                                                                ttGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGaa       >  2:5732015/1‑150 (MQ=255)

                                                                                                                                                ttGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCCCCTGGGCCGAATATTTAATTGCGGATTCGTTGGGTAGTTCAGGGACGTTTGGAAGTGGTGGTGGTGGGGGGAGGATTCgca       >  1:2346562/1‑148 (MQ=17)

                                                                                                                                                ttGCCTGCGGCAACGCTCTCTCGCTGCCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGaa       <  1:575140/150‑1 (MQ=255)

                                                                                                                                                 tGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAAc      <  1:8530879/150‑1 (MQ=255)

                                                                                                                                                     tGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGACGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTc  >  2:8010198/1‑150 (MQ=255)

                                                                                                                                                     tGCGCCAACGCTCTCTCGCTGTCGATCGAGTCGAACCTTGGTATAAGCTTCTCATCCTTCCCAGCATGGGAACAATATTTAATTGCGGATTCTTTTGTAACTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTc  <  1:8851052/150‑1 (MQ=14)

GCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTC  >  NC_000913/1286973‑1287271

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