BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,287,122	T→G	5.9%	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	5.9%	1047.5 / 21.2	426	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 (187/190); new base G (32/3); total (232/194)
Fisher's exact test for biased strand distribution p-value = 5.84e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 440 total aligned reads displayed.

GCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACG‑‑CTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTC  >  NC_000913/1286973‑1287271

gCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCTCCCTTCGGGTCGTTGCCt                                                                                                                                                         >  2:6366277/1‑150 (MQ=255)

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

 ccGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTAGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCTCTTAGCGCTTCGCACTACGGGTCGTGGGCTg                                                                                                                                                        >  1:10079065/1‑150 (MQ=14)

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

    aaCCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGAGTCGCTTGAGAGTTCAGGGACTTTTGAAGGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCgg                                                                                                                                                     >  2:1607912/1‑150 (MQ=255)

          aggcGAAGCTTCTCAGCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCGGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACG‑‑c                                                                                                                                             <  1:3033079/147‑1 (MQ=255)

          aGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCTGGTCGTTGCCTGCGGCAACG‑‑c                                                                                                                                             <  2:2365067/150‑1 (MQ=255)

          aGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTAGCTTGTGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCATCGGGTCGTTGCCTGCGGCaaacc                                                                                                                                               >  1:2950496/1‑147 (MQ=255)

            tCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCATCGCCCTTCGGGTCGTTTCCTGCGGCAAAG‑‑cgc                                                                                                                                           >  2:5030775/1‑148 (MQ=255)

             cGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACG‑‑ctct                                                                                                                                          >  1:7672430/1‑150 (MQ=255)

              gAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACG‑‑ctctc                                                                                                                                         >  1:5099299/1‑150 (MQ=255)

                aGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACG‑‑Ctctctc                                                                                                                                       >  2:8180808/1‑150 (MQ=255)

                          ccTTCCCCGCACTGGCAGAATATTTGATTGCGGATTCGCTTGAGATTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACG‑‑CTCTCTCGCTGGCGCTc                                                                                                                             <  1:176035/150‑1 (MQ=255)

                           cTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACG‑‑CGCACTAGCGGGCGCTCg                                                                                                                            >  2:4088599/1‑150 (MQ=255)

                            ttCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACG‑‑CTCTCTCGCTGGCGCTCGa                                                                                                                           >  1:830780/1‑150 (MQ=255)

                            ttCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACG‑‑CTCTCTCGCTGGCGCTCGa                                                                                                                           <  2:1250687/150‑1 (MQ=255)

                            ttCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCAGCGGCAACG‑‑CACTCACGCTGGCGCTCGa                                                                                                                           >  1:5669749/1‑150 (MQ=255)

                            ttCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCAGCGGCAACC‑‑CTCACTAGCTCgcgcacca                                                                                                                           >  1:10463525/1‑146 (MQ=14)

                            ttCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGACGTTGCCTGGGGCAACG‑‑ATCTCTCGCTTTCGCTCGa                                                                                                                           >  2:1365005/1‑150 (MQ=14)

                                ccGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACG‑‑CTCTCTCGCTGGCGCTCGAGTCg                                                                                                                       >  1:10010213/1‑150 (MQ=255)

                                ccGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACG‑‑CTCTCTCGCTGGCGCGCGAGTCg                                                                                                                       >  2:3663028/1‑150 (MQ=255)

                                        gCAGAATAATTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTTCCTGCGGCAACG‑‑CTCTCTCGCTGGCGCTCGAGTCGAACCTTgg                                                                                                               >  2:1816947/1‑150 (MQ=255)

                                         cAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTTCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACG‑‑CTCTCTCGCTTGCGCTCGACTCGAACCTTGGt                                                                                                              >  2:4430844/1‑150 (MQ=255)

                                         cAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACG‑‑CTCTCTCGCTGGCGCTCGAGTCGAACCTTGGt                                                                                                              >  2:7306049/1‑150 (MQ=255)

                                         cAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCAGCGGCAACC‑‑CTCACACGCTGGCGCTCCAGTCGAACCATGGt                                                                                                              >  2:7684427/1‑150 (MQ=14)

                                          aGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACACAGCGCTGCGCGCATCGCCCCTCGGGTCCTGGCCGCCCCCAACC‑‑CACTCTCGCTCGCGCACGACGCGCACCTTGGTc                                                                                                             >  2:7193818/1‑150 (MQ=14)

                                            aaTATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACG‑‑CTCTCTCGCTGGCGCTCGAGTCGAACCGTGGTCGa                                                                                                           >  1:10614755/1‑150 (MQ=255)

                                                 ttGATTGCGGATTCGCTTGAGATTTCATTGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCATCG‑‑CTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCtt                                                                                                      <  2:2063800/150‑1 (MQ=255)

                                                 ttGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGGGGCAAAG‑‑GTATCTCGCTGGCGCTCGAGTGGAAGCTTGGTCGAATCtt                                                                                                      >  1:8176752/1‑150 (MQ=14)

                                                   gATTGCGGATTCGCTTGAGAGTTCAGGTACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACG‑‑CTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTtct                                                                                                    <  1:5752736/150‑1 (MQ=255)

                                                       gCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACG‑‑CTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATc                                                                                                <  1:9143908/150‑1 (MQ=255)

                                                        cGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACG‑‑CTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATcc                                                                                               <  2:4410633/150‑1 (MQ=255)

                                                          gATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCATGCGGCAACG‑‑CTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTAACATCCtt                                                                                             >  2:3470089/1‑150 (MQ=255)

                                                           aTTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACC‑‑CACACTCGCGGGCGCTCGAGGAGAACCAAGGTCGAAACATCTCaaaactc                                                                                            >  1:4022573/1‑144 (MQ=14)

                                                                 ttGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACG‑‑CTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCa                                                                                      >  2:10837552/1‑150 (MQ=255)

                                                                  tGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACG‑‑CTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCAt                                                                                     <  1:8433189/150‑1 (MQ=255)

                                                                      agTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACG‑‑CTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGc                                                                                 <  2:395352/150‑1 (MQ=255)

                                                                       gTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGTGTCGTTGCCTGCGGCAACG‑‑CTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCAACCTTACCCCCATGGGCa                                                                                >  1:2090991/1‑150 (MQ=255)

                                                                       gTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCATCTCCCTTCGGGACGTTGCCTGCGGCAACG‑‑CTCTCTCGCTGGCGCTCGACTCGAACCTTGGTCGAAGCTTCTCAACCTTCCCCGCAGGGGcc                                                                                >  2:5870020/1‑149 (MQ=14)

                                                                          cAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCTCGCTGGCGCACGAGTCGAACCTTGGTCGTAGCTTCTCATCCTTCCCCGCAAGGGCAg                                                                               >  2:540289/1‑150 (MQ=14)

                                                                           cgggTCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACG‑‑CTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAat                                                                            <  1:8697802/149‑1 (MQ=25)

                                                                           aGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACG‑‑CTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAat                                                                            >  2:9793998/1‑150 (MQ=40)

                                                                            gggACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGAGGCAACG‑‑CTCTCTAGCGGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCAACCATACCCGCATGGGCAGaaga                                                                           >  2:6627123/1‑148 (MQ=14)

                                                                            gggACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCATCGCCCTTCGGGTCGTTGCCGGCGGCAACC‑‑CTCTCACGCAGGCGCTCGAGTCCAACCGTCGTCCAAGCTTAACAAACAACCCCGCAAGGGCAGaaaa                                                                           >  2:6568776/1‑148 (MQ=11)

                                                                             ggACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGGCTGCGGCAACG‑‑CTCTATGGCTGGCGCTGGAGTGGTACCTTGGTGGTAGCTTATAATACTTGACCGCATGGGGAGAatat                                                                          >  1:7104331/1‑150 (MQ=14)

                                                                             ggACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACC‑‑CTCTATCGCTGGCGCGCGAGTCGAACCATGGTCGAAGCGTCACAAACATCACCGCAAGGGCAGAaaat                                                                          >  2:1982081/1‑150 (MQ=14)

                                                                             ggACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCAGCGGCAACC‑‑CACTCACGCGGGCGCGCCGGGCGGAACAGGGGCGAAGCATCTAAACCGTACCCGCAAGGGCAGAatat                                                                          >  1:78535/1‑150 (MQ=2)

                                                                              gACTTTTCAACGTGATGGTGGTGGGGGAAGGATTACTCCGCGCGGCGCGCTTCGCCCTCCGGCTCGTTGCCTGCGGCAACG‑‑CTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATAtt                                                                         <  2:2221762/150‑1 (MQ=11)

                                                                               tcTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACG‑‑CTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATAttt                                                                        <  2:10675622/149‑1 (MQ=37)

                                                                               aCTTTTGAACGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACG‑‑CTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATAttt                                                                        <  2:757657/150‑1 (MQ=38)

                                                                                 ttttGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTTCCTGCGGCAACG‑‑CTCTCTCCCTGGCGCGCGAGTCGAACCTTGGTCGAATCTTCTAATCCATACCCGCATGGGCAGAAAATTTaa                                                                      >  2:334174/1‑150 (MQ=11)

                                                                                 ttttGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACG‑‑CTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTaa                                                                      <  1:9793998/150‑1 (MQ=38)

                                                                                     gAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACG‑‑CTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGc                                                                  <  1:1865218/150‑1 (MQ=39)

                                                                                      gaaGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACG‑‑CTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCg                                                                 <  1:860014/149‑1 (MQ=18)

                                                                                      gaaGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTCGCCTGCGGCAACG‑‑CTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAata                                                                           >  2:10363161/2‑140 (MQ=9)

                                                                                      gaaGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTCGCCTGCGGCAACG‑‑CTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAata                                                                           <  1:10363161/139‑1 (MQ=9)

                                                                                           gATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACG‑‑CTCTCTCGCTGGCGCTCGAGTCGAACCATGGTCGAAGCTTCTCATCCATACCCGCATGGGCAGAATATTTAAATGCGGATTc                                                            >  2:4585735/1‑150 (MQ=255)

                                                                                               gtggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACG‑‑CTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTg                                                        <  1:8994401/150‑1 (MQ=38)

                                                                                               gtggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACG‑‑CTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTg                                                        <  1:4897427/150‑1 (MQ=38)

                                                                                                 ggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACG‑‑CTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTggg                                                      >  2:2697854/1‑150 (MQ=39)

                                                                                                  gtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACG‑‑CTCTCTCGCTGGCGCTCGAGTCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGAt                                                              >  2:8475628/1‑141 (MQ=37)

                                                                                                    gggggAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACG‑‑CTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAg                                                   <  2:6061549/150‑1 (MQ=255)

                                                                                                     ggggAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACG‑‑CTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCACATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGCAGt                                                  >  2:5160055/1‑150 (MQ=255)

                                                                                                     ggggAAGGATTACTCAGCGCTGCGCACTTCGCCCTTCGGGTCGTTGCCTGCGGCAACG‑‑CTCTCTCGCTGGCGCTCGAGTCGAACCTTTGGCGAAGCTTCTCATCCTTACCCCCCAGGGCAGAATATTTAATTTCGGATTAGGTGGGGAGt                                                  >  2:5584075/1‑150 (MQ=11)

                                                                                                     ggggAAGGATTACTCAGCGCTGCGCACTTCGCCCTTCGGGTCGTTGCCTGCGGCAACG‑‑CTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTGGTTGGGAAGt                                                  >  2:9561541/1‑150 (MQ=255)

                                                                                                         aaGGATTACTGAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACG‑‑CTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAg                                              <  1:2855671/150‑1 (MQ=255)

                                                                                                          acgATAACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACG‑‑CTCTCTCGCTGGCGCTCGAGTCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAgg                                             <  1:2638430/148‑1 (MQ=255)

                                                                                                          aGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACG‑‑CTCTCTCGCTGTCGCTCGAGTCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGTAGTTCAgg                                             <  1:110144/150‑1 (MQ=255)

                                                                                                            gATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACG‑‑CTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATc                                                                                                <  1:3687720/97‑1 (MQ=33)

                                                                                                            gATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACG‑‑CTCTCTCGCTGGCGCGCGAGTCGAACCGTGGTCGAAGCTTCACAAACTTCCCCGCATGGGCAGCAAATTTAATTGCGGATTCGTTGGGACGTTCAGGGa                                           >  2:3018616/1‑150 (MQ=14)

                                                                                                             aTTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACG‑‑CTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATcc                                                                                               <  1:4290495/97‑1 (MQ=33)

                                                                                                                 cTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACG‑‑CTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACtttt                                      <  1:2000350/150‑1 (MQ=255)

                                                                                                                 cTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACG‑‑CTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGGTGGGAAGTTCAGGGACtttt                                      >  2:6889303/1‑150 (MQ=255)

                                                                                                                  tCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACG‑‑CTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTg                                     >  2:7913585/1‑150 (MQ=255)

                                                                                                                  tCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACG‑‑CTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTg                                     <  1:9304628/150‑1 (MQ=255)

                                                                                                                  tCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACG‑‑CTCTCTCGCTGACGCTCGAGGCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCCTGGGCAGAATATTTAATTGCGGATTCGGTGGTAAGTTAAGGGACTTTTg                                     >  2:10923353/1‑150 (MQ=11)

                                                                                                                     gcgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACG‑‑CTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGaaa                                  <  1:4585735/150‑1 (MQ=255)

                                                                                                                     gcgcTGCGCGCGTCGCCCTTCGGGTCGTTGCCTGCGGCAACG‑‑CTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGaaa                                  <  1:5167811/150‑1 (MQ=255)

                                                                                                                       gcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACG‑‑CTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGtt                                                         <  1:6688800/125‑1 (MQ=38)

                                                                                                                       gcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACG‑‑CTCTCTCGCTGGCGCTCGAGCCGAACCTTGGTCTAAGCTTATCACCCTTCCCCGCATGGGCATAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGt                                <  1:3938891/150‑1 (MQ=14)

                                                                                                                         tGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACG‑‑CTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGa                                                                              <  1:2404072/102‑1 (MQ=35)

                                                                                                                              gcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACG‑‑CTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGAtggtg                         <  1:9945116/150‑1 (MQ=255)

                                                                                                                                   gCCCTTCGGGTCTTTGCCTGCGGCAACG‑‑CTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTggg                    <  2:8855793/150‑1 (MQ=255)

                                                                                                                                   gCCCTTCGGGTCGTTGCCTGCGGCAACG‑‑CTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTggg                    >  2:7245676/1‑150 (MQ=255)

                                                                                                                                    cccTTCGTGTCGGTGACTGCGTCATCT‑‑CTATCTCGCTGGAGCTCGTGTCGCTGCTATGTAGATGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTgggg                   <  1:6003612/150‑1 (MQ=14)

                                                                                                                                    cccTTCGGGTCGTTGCCTGCGGCAACG‑‑CTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTgggg                   <  2:5636941/150‑1 (MQ=255)

                                                                                                                                    cacTTCGGGTCGTTGCCTGCGGCAACG‑‑CTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTgggg                   >  2:1862767/3‑150 (MQ=255)

                                                                                                                                         cGGGTCGTTGCCTGCGGCAACG‑‑CTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAgg              <  1:1441098/150‑1 (MQ=255)

                                                                                                                                          agagcGTTGCCGCAGGCAACG‑‑CTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGa             >  2:8634921/5‑150 (MQ=255)

                                                                                                                                              cGTTGCCGCAGGCAACG‑‑CTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGtt                                                         <  2:6251546/102‑1 (MQ=255)

                                                                                                                                              cGTTGCCGCAGGCAACG‑‑CTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGtt                                                         >  1:6251546/1‑102 (MQ=255)

                                                                                                                                                ttGCCTGCGGCAACG‑‑CTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGaa       >  2:9007364/1‑150 (MQ=255)

                                                                                                                                                ttGCCTGCGGCAACG‑‑CTCTCTCGCTGGCGCTCGAGTCGAACCATGGTCGAAGCGTCTCAAACATACCCGCAAGGGCCCAAAATTTAATTGAGGAGTCGTTGGGAAGGTCAGGGACGTTTTAAATGGATGGGGGTGGGGGGAGGATTCGaa       >  1:3351888/1‑150 (MQ=17)

                                                                                                                                                ttGCCTGCGGCAACG‑‑CTCTCTCGCAGGCGCTCGAGTCGAACCTTGGTCGAAGCATCTCATCCTTCCCCGCAGGGGCCGGCGATTTAATTGCGGAGTCGGTGGGTAGTGCAGGGGCTTTTGACAGTGAGGGTGGTGGGGGCAGGATTCgga       >  1:10857263/1‑148 (MQ=14)

                                                                                                                                                 tGCCTGCGGCAACG‑‑CTCTCTCGCTGTCGCTCGAGTCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTTATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAAc      <  1:2247844/150‑1 (MQ=255)

                                                                                                                                                    cTGCGGCAACG‑‑CTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCtt   >  2:2565336/1‑150 (MQ=255)

                                                                                                                                                     tGCGGCAACG‑‑CTCTCTCGCTGTCGCTCGAGTCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTTATTGCGGATTCGATTGTTAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTc  <  1:8532439/150‑1 (MQ=11)

                                                                                                                                                     tGCGGCAACG‑‑CTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGaa                                                                             <  2:7829933/75‑1 (MQ=39)

                                                                                                                                                     tGCGGCAACG‑‑CTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTc  <  2:10188842/150‑1 (MQ=255)

                                                                                                                                                     tGCGGCAACG‑‑CTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCa                                               <  1:3186247/105‑1 (MQ=255)

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