BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,287,122	T→G	11.4%	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	11.4%	539.5 / 20.0	254	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 (108/102); new base G (31/2); total (149/105)
Fisher's exact test for biased strand distribution p-value = 9.31e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 258 total aligned reads displayed.

GCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTC  >  NC_000913/1286973‑1287271

gCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGATTGAGAGTTCCGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCt                                                                                                                                                       <  2:3714658/150‑1 (MQ=16)

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

 ccGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCGTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTg                                                                                                                                                      >  2:3479654/1‑150 (MQ=255)

  cGAACCTTAGTCGAAGCTTCTCTTCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGc                                                                                                                                                     <  2:3338204/150‑1 (MQ=255)

  cGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTTCCTGc                                                                                                                                                     >  1:4372693/1‑150 (MQ=255)

    aaCCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCgg                                                                                                                                                   >  2:1121199/1‑150 (MQ=255)

       cTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCaa                                                                                                                                                >  2:5701670/1‑150 (MQ=255)

       cTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCaa                                                                                                                                                >  1:2693288/1‑150 (MQ=255)

            tCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGGGGGGGAAGGATTACTCCGCGCGGCGCGCTGCGCCCTTCGGGTCGTGGCCGGCGGCACCGctc                                                                                                                                           >  2:4708348/1‑150 (MQ=14)

                    tctcATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTgg                                                                                                                                                   <  2:6263699/134‑1 (MQ=14)

                    tctcATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTg                                                                                                                                   >  1:2130192/1‑150 (MQ=255)

                    tctcATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGACGTTGCCTGCGGCAACGCTCTCTCGcag                                                                                                                                   >  2:2877243/1‑148 (MQ=255)

                      tcATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCGACGCTCTCTCGCTGgc                                                                                                                                 >  1:289065/1‑150 (MQ=255)

                       cttCCTTCCGCGCATGGGCAGATTATTTGATTGCGGTTTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTGAGCGCTGCGCGGTTCGCCCATCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGgcg                                                                                                                                <  1:2580102/148‑1 (MQ=14)

                         gcattCCCCGCCTGGTCTGAATATTTGATTGCGGATTCGATTGAGATTTCAGGGTCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCt                                                                                                                              <  1:2778581/147‑1 (MQ=11)

                          ccTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTc                                                                                                                             >  2:1838937/1‑150 (MQ=255)

                            ttCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGAGCTGCGCGCTTCGCCCTTCGGGTCGTTTCCTGCGGCAACGCTCTCTCGCTGGCGCGCGa                                                                                                                           >  1:4829283/1‑150 (MQ=255)

                                  gcttGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGaa                                                                                                                     <  2:752927/147‑1 (MQ=255)

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

                                     tGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCt                                                                                                                  <  2:4195534/150‑1 (MQ=255)

                                        gCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTgg                                                                                                               <  2:6342638/150‑1 (MQ=255)

                                        gCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCACACGCAGGCGCTCGAGTCGAACCTTgg                                                                                                               >  1:4074032/1‑150 (MQ=255)

                                        gCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCGCTCTCGCTGGCGCTCGAGTCGAACCATgg                                                                                                               >  2:387625/1‑150 (MQ=255)

                                        gCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACCCACTCTCGCAGGCGCTCGAGTCGAACCTTgg                                                                                                               >  1:6714421/1‑150 (MQ=255)

                                         cAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGt                                                                                                              <  2:4414338/150‑1 (MQ=255)

                                         cAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCGGCGCGCTTCGCCCGACGGGTCCATGCCGGCCGCACCGCCCTCACGCTCGCGCCCCCGTCGAAccggggg                                                                                                              >  2:8158786/1‑145 (MQ=14)

                                                    aTTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTtctc                                                                                                   <  1:1838937/150‑1 (MQ=255)

                                                     ttGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCa                                                                                                  >  2:5253347/1‑150 (MQ=255)

                                                       gCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATc                                                                                                <  2:6016691/150‑1 (MQ=255)

                                                       gCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATc                                                                                                <  1:3022781/150‑1 (MQ=255)

                                                       gCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACCCTCTCTCGCGGGCGCGCGAGTCGCACCTTGGTCGGAGCGTCTCaac                                                                                                >  2:3595870/1‑148 (MQ=14)

                                                        cGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATcc                                                                                               <  1:6368218/150‑1 (MQ=255)

                                                         ggATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCACGCTGGCGCTCGAGTCGAAACGTGGTCGAATCTTCTCATcat                                                                                              >  2:77703/1‑148 (MQ=255)

                                                         ggATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACGCACACTAGCTGGCGCTCCAGTAGCACCATGGTAGAAGCATAACATcca                                                                                              >  1:3400873/1‑149 (MQ=14)

                                                           aTTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACGCGCACACGCGGGCGCTCGAGTCGGACCTTGGTCGAAGCGTCTCaaacatc                                                                                            >  1:1388598/1‑144 (MQ=14)

                                                               gCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTccccc                                                                                        >  1:5989366/1‑149 (MQ=255)

                                                                 ttGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTTCCTGCGGCAACGCACTCTCGCTGGCGCGAGAGTCGAACCGTGGTCGAAGCGTCTCATCCATCCCCGCa                                                                                      >  2:6608092/1‑150 (MQ=14)

                                                                 ttGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCAACCTTCCCCGCa                                                                                      >  1:4536250/1‑150 (MQ=255)

                                                                 ttGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCACGCGGGCGCTCGACTCGAAACATGGTCGAAGCGTCTCAACCATACCCCCa                                                                                      >  1:3817424/1‑150 (MQ=14)

                                                                 ttGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCACTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCGTCACAACCTTACCCCCa                                                                                      >  2:6478014/1‑150 (MQ=14)

                                                                 ttGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCACTCACGCGGGCGCTCGAGTCGAACCTTGGTAGAAGCATCTCAAACTTCCCCGCa                                                                                      >  1:6158415/1‑150 (MQ=14)

                                                                 ttGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACCCTCTCTCGCAGTCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCAACCTTACCCGCa                                                                                      >  2:3117087/1‑150 (MQ=14)

                                                                       gTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTTCCTGCGGCAACGCTCTCTAGCTGGCGCTCGAGTCGAACCGTTGTAGAAGCTTCTAAAACATCCCCGCATGGGcc                                                                                >  1:2521295/1‑149 (MQ=14)

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

                                                                        ttCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCACACTCGCTGGCGCACGAGTCGAACCTTGGTCGAATCGTGTCATACATCACCCCAGGGGccg                                                                               >  1:7914588/1‑148 (MQ=14)

                                                                           aGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTTGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCAACCTTCCCCGAATggtgacaaa                                                                            >  2:6817577/1‑143 (MQ=14)

                                                                            gggACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAata                                                                           >  2:544398/1‑150 (MQ=39)

                                                                            gggACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCGCTCTCGCTGGCGCTCGAGTCGAACCTTGGTAGAAGCTTCTCATCCTTCCCCGCATGGGCAGAata                                                                           >  2:2034763/1‑150 (MQ=255)

                                                                              gtcTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATAtt                                                                         <  1:5529726/148‑1 (MQ=37)

                                                                              gACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATAtt                                                                         <  1:1832425/150‑1 (MQ=38)

                                                                                  tttGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCAGGCGCTCGAGTCGAACCTTGGTCGAAGCATCTCAACCATCCCCGCATGGGCAGCATATATAAt                                                                     >  2:4892433/1‑150 (MQ=11)

                                                                                  tttGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCAGCGGCAACGCACTCACGCGGGCGCTCGAGGCGAAACTTGGTCGAAGCGTCTAATACGTACCCCCCAGGGCGGAATATTTAAt                                                                     >  1:7644638/1‑150 (MQ=11)

                                                                                      gaaGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCa                                                                                <  1:8232595/134‑1 (MQ=14)

                                                                                      aaaGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTTCCTGGGGCCACGCTCTCTCGCTTGCGCGCGAGTCGAAACTTTGTCGAAGCGTCTCAACCATCCCCCCCGGGGGCGAATATTTaaaagcg                                                                 >  1:1012812/1‑145 (MQ=11)

                                                                                          tGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACCCTCTCTCGCTGGCGCTCGAGTCGAACCAAGGTCGAAGCTTCTCAACCTTCCCCCCATCGGCAGAAAATTAAATATCGGaat                                                             >  2:4360089/1‑148 (MQ=11)

                                                                                          tGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACCCTCTCACGCTTGCCCTCGAGTGGAAACATGGGCGAAGCATCACAACCAACCCCCCCAGGGCCGAATAGTTATGTGCGGGtt                                                             >  1:5254915/1‑150 (MQ=2)

                                                                                          tGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACGCGCTCTCGCAGGCGCGCGAGTCGAACCTTGGTCGAAGCGTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGttt                                                             >  2:2446472/1‑150 (MQ=11)

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

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

                                                                                                  gtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCACGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCAAGGGCAGAATATTTAATTGCGGATTCGTTGGGa                                                     >  2:8251423/1‑150 (MQ=255)

                                                                                                        gAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCa                                               <  1:568906/150‑1 (MQ=255)

                                                                                                        gAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCa                                               >  2:5426708/1‑150 (MQ=255)

                                                                                                        gAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTGCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCa                                               >  2:3010554/1‑150 (MQ=255)

                                                                                                        gAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCGCTCGCTGGCGCTCGAGTCGAACCGTTGTCGAAGCTTCACAACCATCCCCCCATGGGCAGAATATTTAATTGCGGCTTAGGTGGGAAGttaa                                               >  2:7871388/1‑148 (MQ=14)

                                                                                                         aaGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGTGTCGTTGCCTGCTGCAACGCTCTCTCGCTGTCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAg                                              <  1:6618951/150‑1 (MQ=255)

                                                                                                          aGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAgg                                             <  1:7438681/150‑1 (MQ=255)

                                                                                                            gATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATc                                                                                                >  2:7582693/1‑97 (MQ=33)

                                                                                                            gATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCAACCTTCCCCCCATGGGCACAAAATTTAAAATCGGAGACGTTGGGAAAATAAGGGa                                           >  2:7207229/1‑150 (MQ=11)

                                                                                                             aTTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGGTGCGTGCGGCATCTCTCTCTCGCTGGCGCTCGAGTATAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGAc                                          <  1:3637397/150‑1 (MQ=17)

                                                                                                               tACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCa                                                                                      <  1:7326808/104‑1 (MQ=33)

                                                                                                               tACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCa                                                                                      >  2:7326808/1‑104 (MQ=33)

                                                                                                                aCTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACttt                                       <  1:5152941/150‑1 (MQ=255)

                                                                                                                aCTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACttt                                       <  1:8125435/150‑1 (MQ=255)

                                                                                                                 cTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGTAGCTTCTCATCCTTACCCGCCTGGGCCGAATATTTAATTGCGGATTGGGTGGGAAGTTCAGGGGCtttt                                      >  2:1786798/1‑150 (MQ=14)

                                                                                                                  tCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAAtt                                                                    <  1:6509011/119‑1 (MQ=39)

                                                                                                                  tCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTg                                     >  2:2805063/1‑150 (MQ=255)

                                                                                                                         tGCGCGCTTCGCCCTTTGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGa                              <  2:1591945/150‑1 (MQ=255)

                                                                                                                         tGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTc                                                                                                                        <  1:8303219/60‑1 (MQ=35)

                                                                                                                         tGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTc                                                                                                                        <  1:5298619/60‑1 (MQ=35)

                                                                                                                         tGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGa                              <  2:413366/150‑1 (MQ=255)

                                                                                                                         tGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGa                              <  1:4349658/150‑1 (MQ=255)

                                                                                                                         tGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGa                              <  1:6739482/150‑1 (MQ=255)

                                                                                                                           cgcgcgTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGGCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGAtg                            <  1:2697092/150‑1 (MQ=255)

                                                                                                                                ttCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTg                                                                                                                <  1:1694787/61‑1 (MQ=35)

                                                                                                                                 tCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGAtggtggtg                      <  2:5284146/150‑1 (MQ=255)

                                                                                                                                   gCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCGTTGTCGAAGCTTCTCAAACATCCCCCCAAGGGCAGAATATTTAATTGCGGATTAGTTGGGGAGTTAAGGGACTTTTTTAAGTGATGGTGGTggg                    >  2:3791494/1‑150 (MQ=14)

                                                                                                                                      cTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGaa                                                                                                                     >  2:8368576/1‑50 (MQ=38)

                                                                                                                                      cTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGaa                                                                                                                     <  1:8368576/50‑1 (MQ=38)

                                                                                                                                      cTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATgg                                                                                   <  1:7038889/84‑1 (MQ=38)

                                                                                                                                      cTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATgg                                                                                   >  2:7038889/1‑84 (MQ=38)

                                                                                                                                       ttCGGGTAGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCGTCCCCGCATGGGCATAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGaa                <  1:7550582/150‑1 (MQ=255)

                                                                                                                                        tCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAg               <  1:4667481/150‑1 (MQ=255)

                                                                                                                                            gTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGTCAGAATATTTAATTGCGGATTAGTTTGTAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGAtt           <  1:3291688/150‑1 (MQ=255)

                                                                                                                                            agcGTTGCCGCAGGCAACGCTATCTCGCTGGCGCTCGAGTCGAACCGTGGTGGAAGCATATCAAACTTACCCGCAAGGGCACAATATTTAATTTCGGGTTCGGTGGGAAATTAATGGACTTTTGTAAGGGATGTTGGTGGGGGAAGGAtt           >  1:5771160/3‑150 (MQ=255)

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

                                                                                                                                               gTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCtt                                                                                                                 <  1:7676856/45‑1 (MQ=38)

                                                                                                                                               gTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCtt                                                                                                                 >  2:7676856/1‑45 (MQ=38)

                                                                                                                                                 tGCCTGCGGCATCGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAAc      <  1:6857277/150‑1 (MQ=255)

                                                                                                                                                   ccTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCa                                                                                                  <  1:3218570/56‑1 (MQ=35)

                                                                                                                                                   ccTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCt    >  2:5021529/1‑150 (MQ=255)

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