BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,287,122	T→G	9.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
*	NC_000913	1,287,122	0	T	G	9.1%	774.9 / 31.0	351	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 (162/132); new base G (36/2); total (216/135)
Fisher's exact test for biased strand distribution p-value = 6.26e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 356 total aligned reads displayed.

CCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTC  >  NC_000913/1286979‑1287271

  ttAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCTCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAAc                                                                                                                                               >  2:4682378/1‑150 (MQ=255)

     gTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACACAGCGCTGAGAGCTTCGCCCGACGGGTAGTTTCCGGCCGCAAAGcg                                                                                                                                            >  1:1696803/1‑149 (MQ=14)

      tCGAAGCTTCTCATCCATCCCCGCATGGGCCGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCATCGGGTCGTTGCCTGCGGCCACGcgc                                                                                                                                           >  1:2466704/1‑148 (MQ=255)

        gAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGctctc                                                                                                                                         >  2:916219/1‑150 (MQ=255)

             ttCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCt                                                                                                                                    <  1:9997117/150‑1 (MQ=255)

                tcATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCTCCCTTCGGGTCGTTTCCTGCGGCAACGCTCTCTCGCTGgc                                                                                                                                 >  2:6098387/1‑150 (MQ=255)

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

                    ccTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTc                                                                                                                             <  1:362379/150‑1 (MQ=255)

                     cTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACGCACTCTCGCTGGcgcgcg                                                                                                                            >  2:2342830/1‑150 (MQ=255)

                      ttCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGa                                                                                                                           <  1:9954344/150‑1 (MQ=255)

                      ttCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCACTCGCTGGCGCTCGa                                                                                                                           >  1:1333020/1‑150 (MQ=255)

                      ttCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAAAGCTCTCACGCTGGCGCACGa                                                                                                                           >  2:10389350/1‑150 (MQ=255)

                      ttCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGAGGCAACGCACTCTCGCTGGCGCGCGa                                                                                                                           >  1:9053348/1‑150 (MQ=255)

                          ccGCTTGGGCTGAATATTTGATTGCGGATTCGCTTGAGTGTTCAGGGTCTTTTGAATGTGATGGTGGTGGGGGTAGGATTTCTCAGCGCTGCGCGCTTCTCCCTTCGGGTCGTTGCCTGCGTCTTCGCTCTCTCGCTGGCGCTCGAGTCg                                                                                                                       <  1:8128089/150‑1 (MQ=14)

                           cGCATGGGCAGCATAGTAGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGa                                                                                                                      <  2:8812560/150‑1 (MQ=255)

                           cGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCGCTCTCGCTGGCGCTCGAGTCGa                                                                                                                      >  1:8851069/1‑150 (MQ=255)

                           cGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAAAGCTCTAACGCTGGCGCACGAGTAGa                                                                                                                      >  1:2552807/1‑150 (MQ=14)

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

                             cATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTACGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAAc                                                                                                                    <  1:5517426/150‑1 (MQ=255)

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

                                 ggCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCGCTCGAGCTGGCGCTCGAGTCGAACCATg                                                                                                                >  1:9732039/1‑150 (MQ=255)

                                  gCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCATCGCCCTTCGGGTCGTTGCCTGAGGCAACGCTCTCTCGCTTGCGCTCGAATCGAACCttag                                                                                                               >  2:10153295/1‑148 (MQ=255)

                                  gCAGAATACTCGATTGCGGAGTCGATTGAGATTTCAGGGTCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTgg                                                                                                               <  2:4009334/150‑1 (MQ=11)

                                   cAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGt                                                                                                              >  2:9100741/1‑150 (MQ=255)

                                   cAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCATCGCCCTTCGGGTCGTTGCCTGCGGCAACCCTCTCACGCTTTCGCTCGACTCGAAccagata                                                                                                              >  2:84279/1‑145 (MQ=14)

                                      aaTATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGa                                                                                                           <  2:2845394/150‑1 (MQ=255)

                                       atagTTGTTGGCGGTTTCGCTTGAGAGTTCAGGGACTTTTGTAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAAGCTTGGTCGaa                                                                                                          <  1:1880530/150‑1 (MQ=14)

                                        tatTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCATGGTCGAAg                                                                                                         >  1:2725800/1‑150 (MQ=255)

                                          tttGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCTCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTATCTCTGGCGCTCGAGTCGACCCTTGGTCGAAGCt                                                                                                       >  2:5605609/1‑150 (MQ=255)

                                           ttGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTCCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCtt                                                                                                      <  1:3094299/150‑1 (MQ=255)

                                            tGAGTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTtc                                                                                                     <  2:6574366/150‑1 (MQ=255)

                                             gtttGGGGATTCGTTTGAGATTTCAGGGTCTTTTGATAGTGTTGGTGGTGGGGGAAGGATTACTCAGCGGTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTtct                                                                                                    <  1:8216279/148‑1 (MQ=11)

                                                  cGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATcc                                                                                               >  2:8285395/1‑150 (MQ=255)

                                                   ggATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGAATACTAAGCGCTGCGCGCTGCGCCCTTCGGGTCGTTGCCGGCGGCAACGCTCTCTCGCTCGCGCTCGAGACGAACCGTGGTCGAAGCGTCGCATccg                                                                                              >  2:7720604/1‑149 (MQ=14)

                                                    gATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCtt                                                                                             >  2:7186652/1‑150 (MQ=255)

                                                      ttCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTcc                                                                                           >  2:10846078/1‑150 (MQ=255)

                                                      ttCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCGCACTCGCTGGCGCACGAGGCGAACCATGGTCGAAGCTTCTCaaacatac                                                                                           >  1:9125388/1‑143 (MQ=14)

                                                         gCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCAGAGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCg                                                                                        >  2:11143413/1‑150 (MQ=255)

                                                          cTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGCACCATGGTCGAAGCTTATCATCCATCCCCGc                                                                                       >  1:6872015/1‑150 (MQ=255)

                                                           ttGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAATCTTCTCAACCAACCCCGCa                                                                                      >  1:10003639/1‑150 (MQ=255)

                                                           ttGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCATGGTCGAAGCTTCTCATCCATacccccc                                                                                      >  1:7068517/1‑147 (MQ=14)

                                                           ttGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACGCTCTCTCGCGGGCGCTCGAGTCGCACCGTGGTCGAAGCGGCTCATCCTTCCCCGcc                                                                                      >  2:3450571/1‑149 (MQ=14)

                                                           ttGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTAAGCGCTTCGCGCTTCGCCCTTCGGGTCGTTGCCTGACGCAACGCGCTCTCGCTAGCGCAAGAGTCGAACCTAGGTCGAATCTTCTAATCCTTCCCCGct                                                                                      >  2:6479310/1‑149 (MQ=14)

                                                            tGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCAGCGGCAACCCACTCACCCAGGCGCACCAGTCGAACCAGGGTCGAAGCATCCCAACCAAccccccca                                                                                     >  1:6733717/1‑146 (MQ=14)

                                                             gagaGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCAGCGGCAACGCACTCACGCGGGCGCTCGCGTCGAACCGTGGTCGAAGCCTCTCAACCATcccccccag                                                                                    >  2:4978843/1‑145 (MQ=14)

                                                               gagTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCATCGCCCTTCGGGTCGTTGCCTGCGGCAACGAGCTATCGCTTGCGCAAGAGTCGAACCTTGGTCGAAGCTACACATCCTTCCCCGCAAggg                                                                                  >  2:1841644/1‑150 (MQ=14)

                                                                 gTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCATCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTAGATTCGAAACTTGGACGAATCATCACATCCTTCCCCGCATGGGCa                                                                                >  2:4646786/1‑150 (MQ=14)

                                                                  ttCAGGGACTTTTGAAAGTGATGGTGTTGGTGGAAGGATTAATCAGCGCTGCGCTCATCGCCCTTCGGGGCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTAGAACCATGGTCGAAGGTTCTCAACCTTCCCCGCagcgtgcg                                                                               >  2:6516401/1‑143 (MQ=14)

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

                                                                     aGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCATCGCTGCGCTCTTCGCCCTTCGGGACGTTGCCTGCGGCAACGCTCTCTCGCTGGCCCGCGAGTCGAACCTTGGTAGAATCTTCTAATACTTCCCCGCATTGGCCGAat                                                                            >  2:4220842/1‑150 (MQ=14)

                                                                     aGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGGCGGCGGCAACGCTCTCTCGCGGGCGCACGAGTCGGACCTTGGTCGGAGCTTCACATCCATCCCCCCCTGGGccgcaa                                                                            >  1:5032934/1‑145 (MQ=14)

                                                                     aGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGTGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAat                                                                            >  2:1464743/1‑150 (MQ=255)

                                                                      gggACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTTGTCGAAGCTTATCAACCATCCCCGCATGGGCAGAata                                                                           >  2:5495710/1‑150 (MQ=255)

                                                                        gACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGACGTTGCCATCGGCAACCCTCTCACGCATGCGCTCGAGTCGAACCGTGGTCGAAGCTTATCAACCATACCCGCATGGGCCGaaaaat                                                                         >  1:9627678/1‑146 (MQ=11)

                                                                           ttttGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTaa                                                                      <  1:3226977/150‑1 (MQ=38)

                                                                           ttttGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCATCCCCGCATGGGCAGAATATTTaa                                                                      >  2:7477266/1‑150 (MQ=38)

                                                                           ttttGAAAGTGATGGTGGTGGGGGAAGGATTACTAACCGCTGCGCGCTTCGCACTTCGGGTCGATTCCAGCGTCAACGCACACTCCCTTGCCCTCAAGAAGAACCGGGGTCCGAGCTTCTAAACCTTCCCCGCATGGGCAGAATATTTaa                                                                      >  2:10701371/1‑150 (MQ=2)

                                                                             ttGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAAtt                                                                    >  2:7380579/1‑150 (MQ=39)

                                                                              tGAAAGTGATGGTGGTGGGGGAAGGATTACTCCGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTg                                                                   <  1:7516913/150‑1 (MQ=255)

                                                                               gAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGc                                                                  <  1:8395795/150‑1 (MQ=39)

                                                                                    tGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGACTGCCGCAACCCTCACTAGCAGGCGCTCGAGTCGAAACATGGTCCAACCTTCAAAACCATACCCCCAAGGGCAGAAAAATTAATTGAGGAtt                                                             >  2:3070054/1‑150 (MQ=2)

                                                                                      aTGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTCATTGCGGATTCg                                                           <  1:3898290/150‑1 (MQ=21)

                                                                                        ggtggtggGGGAAGGATTACTCCGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTCAATTGCGGATTCGtt                                                         <  1:4388789/150‑1 (MQ=255)

                                                                                        ggtggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGtt                                                         <  1:9539707/150‑1 (MQ=37)

                                                                                        ggtggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTCAATTGCGGATTCGtt                                                         <  1:605864/150‑1 (MQ=37)

                                                                                        ggtggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCCCACACGCTGGCGCTCGAGTCGCACCGTGGTCGAAGCGTATAAACCATCCCCCCCAGGGCAGAATATTTAAATGCGGATTCGtt                                                         >  1:7854159/1‑150 (MQ=11)

                                                                                            gtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGt                                                                                                                         <  1:6824018/82‑1 (MQ=33)

                                                                                            gtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGa                                                     <  2:8066478/150‑1 (MQ=40)

                                                                                             tggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGaa                                                    <  1:4700911/150‑1 (MQ=40)

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

                                                                                                           cTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACtttt                                      <  1:5371593/150‑1 (MQ=255)

                                                                                                            tCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGAt                                                              <  1:583767/125‑1 (MQ=38)

                                                                                                            tCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGAt                                                              >  2:583767/1‑125 (MQ=38)

                                                                                                            tCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCGCGCTGGCGCGCGAGCCGAACCGGGGTCGAAGCTTCTCATCCTTCCCCGCCTGGGCCGAAGATTTGCTTGCGGTTTCGCTGGGCAGTTCAGGGCCGTTTg                                     >  2:2512213/1‑150 (MQ=11)

                                                                                                             cAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGa                                    <  2:2725800/150‑1 (MQ=255)

                                                                                                             cAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCAACCTTCCCCGCAGTGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTGAGGGACTTTTGa                                    >  2:3300108/1‑150 (MQ=255)

                                                                                                              aGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCCTCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGGTGGGCAGTTCAGGGACTTTTGaa                                   >  2:7067911/1‑150 (MQ=255)

                                                                                                              aGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGaa                                   <  1:6505718/150‑1 (MQ=255)

                                                                                                              aGCGCCGCGCGCTTCGCCCTCCGTTTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGaa                                   <  2:6745120/150‑1 (MQ=255)

                                                                                                               tctccgcgcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCATCGCTCTCTCGCTGTCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGaaa                                  <  2:6179643/145‑1 (MQ=255)

                                                                                                               gctcggcgcgcATCGCCCGTCGTGCCGTAGCCGTCGGCATCGCTCTCTCGCTGGCGCTCGAGTCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGaaa                                  <  2:9213243/145‑1 (MQ=14)

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

                                                                                                               gcgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACCCTCACTCGCTGGCGCGCGAGTCGAACCGTGGTCGTAGCGTCTCCTCCATACCCGCATGGGCAGAATAGTTGATTGCGGATTGGGTGGGGAGGTGAGGGTCGTTTggag                                  >  2:3895963/1‑147 (MQ=2)

                                                                                                                   tgagcgcTTCGCCCTTCGGGTCGTTGCCGGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATTGGCAGAATATTTAATTGCGGATGGGTTGGGAAGTTCAGGGCCTGTTGAAAGTGa                              >  2:811127/4‑150 (MQ=17)

                                                                                                                     cgcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAg                                                                               <  1:6359043/99‑1 (MQ=37)

                                                                                                                     cgcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAg                                                                               >  2:6359429/1‑99 (MQ=18)

                                                                                                                          ttCGCCCTTCGGGTCTTTGCCTTAGGGAACGGTCTCTCGCTGGCGGTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGAtggtggt                       <  1:7376798/150‑1 (MQ=17)

                                                                                                                          ttCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAAtt                                                                    >  2:9720675/1‑105 (MQ=17)

                                                                                                                          ttCGCCCTTCGGGTCGTTGCCGGCGGCAACGCTCTCTCGCTGGCGCACGAGTCGAACCGTGGTCGAAGCTTATCATCCTTCCCCGCAGGGGCAGAATATTTAATAGCGGAATCCTTGGGACGTTCAGGGACGGTTGAAAGtggtggtggg                       >  2:356749/1‑149 (MQ=14)

                                                                                                                              cccTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTgggg                   >  2:11004983/1‑150 (MQ=255)

                                                                                                                                      gTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTg                                                                                                                <  1:6704268/49‑1 (MQ=38)

                                                                                                                                      agcGTTGCCGCAGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCCTGGGCCGAATATTTAATTGCGGATTCGTTGGGCCGTTCAGGGACGGGGGACAGGGATGGTGGTGGGGGACGGAtg           >  2:854612/3‑149 (MQ=14)

                                                                                                                                      agcGTTGCCGCAGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATAGCGGAGTCGTTTGTAAGTTCAGGGACTTTTGAAATGGATGGTGGTGGGTGAAGGAtt           >  2:1678502/3‑150 (MQ=14)

                                                                                                                                         gTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGa        >  1:5879769/1‑150 (MQ=255)

                                                                                                                                          ttGCCTGTGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGaa       <  1:1464743/150‑1 (MQ=255)

                                                                                                                                          ttGCCTGCGGCAACGCTCTCTCGCTGGCGCTCTAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTg                                                                   <  1:3473152/90‑1 (MQ=255)

                                                                                                                                          ttGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCGTGGTCGAAGCTTCACAACCATCCCCCCCGGGGCAGAAAATTTAATTGCGGCATCGTTGGGAAGGTCAGGGACGTTTGAAAGTGATGGGGGTGGGGGACGGCTTCGaa       >  1:10174600/1‑150 (MQ=17)

                                                                                                                                          ttGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGaa                                   <  1:10810068/122‑1 (MQ=255)

                                                                                                                                            gCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAAcc     <  1:10443025/150‑1 (MQ=255)

                                                                                                                                              cTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCtt   >  2:10469338/1‑150 (MQ=255)

CCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTC  >  NC_000913/1286979‑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: `‑`