BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,318,080	A→T	5.5%	intergenic (‑74/+133)	rimP ← / ← metY	ribosome maturation factor for 30S subunits/tRNA‑Met

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,318,080	0	A	T	5.5%	964.8 / 31.6	275	intergenic (‑74/+133)	rimP/metY	ribosome maturation factor for 30S subunits/tRNA‑Met
Reads supporting (aligned to +/- strand): ref base A (133/125); new base T (5/10); total (139/136)
Fisher's exact test for biased strand distribution p-value = 1.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.38e-01

Only 100 of 276 total aligned reads displayed.

CACAGCACTTACCTGGTGGCTCACATCAGCACAATCATCAACATTGATGCCATCTTCACTATCAATATAGATGCGCAGTGTGGATGTGCGACCGCGAATAAATTCGATGCCAACCAGTTCAAAACCCAGGGCCTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGACTCTAAGGGAAAGTGGTTGCGGGGGCCGGATTTGAACCGACGATCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGCCTGAAACGTGGCAAATTCTACTCGTTTTGGGTAAAAAATGCAAAT  >  NC_000913/3317824‑3318335

cacagcacTTACCTGGTGGCTCACATCAGCACAATCATCAACATTGATGCCATCTTCACTATCAATATAGATGCGCAGTGTGGATGTGCGACCGCGAATAAATTCGATGCCAACCAGTTCAAAACCCAGGGCCTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAAc                                                                                                                                                                                                                                               >  2:309903/1‑275 (MQ=255)

                           atacagatcatcTACATTGATGCCATCTTCACTATCAATATAGATGAGCAGTGTGGATGTGCGACCGAGAATAAATTCGATGCCAACCAGTTCAAAACCCATGGCCTCAACTGGCGCAGTAATCATCTCTGTTAAGTTTTGCTCTAATGTGGACAAGCCCACCCACAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAg                                                                                                                                                                                                                <  1:288487/273‑1 (MQ=255)

                                            tGATGCCATCTTCACTATCAATATAGATGCGCAGTGTGGATGTGCGACCGCGAATAAATTCGATGCCAACCAGTTCAAAACCCAGGGCCTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATc                                                                                                                                                                                                                                                               >  2:132882/1‑216 (MQ=255)

                                                                         cgcAGTGTGGATGTGCGACCGCGAATAAATTCGATGCCAACCAGTTCAAAACCCAGGGCCTCAACTGGCGCAGTAACCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCNTCAAGCACAGCTATTCTTTAGCCAGAGAACAAAAACCACAGATAAAACGGA‑GATTTATATAAATGAACCGAAAANAAGAAACAGCGAACTGGAGAAAGTTACAAANNGATTTTTTAACAACACACAAAGAAGGCCaaana                                                                                                                                                                  >  3:100265/1‑277 (MQ=255)

                                                                           cAGTGTGGATGTGCGACCGCGAATAAATTCGATGCCAACCAGTTCAAAACCCAGGGCCTCAACTGGCGCAGTAAACATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTTAGATAACAAAAAACACCGATAAATCGGG‑GCTTTATATAACTGAAACCTATAAACGCAACTGAGGGCTGGAACAATTTCCAGAAGGCATTTTTCAAAACCCAATACGCAAGCCCaa                                                                                                                                                                    >  4:143767/1‑275 (MQ=255)

                                                                              tgtgGATGTGCGACCGCGAATAAATTCGATGCCAACCAGTTCAAAACCCAGGGCCTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCt                                                                                                                                                                                                                    >  4:453193/1‑224 (MQ=255)

                                                                              tgtgGATGTGCGACCGCGAATAAATTCGATGCCAACCAGTTCAAAACCCAGGGCCTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTATGGAGCACTTTCCAGAAGGATTTTTTCAAATCCAACTACGAAGGCCGAAGTc                                                                                                                                                                 >  2:515626/1‑275 (MQ=255)

                                                                                tgGATGTGCGACCGCGAATAAATTCGATGCCAACCAGTTCAAAACCCAGGGCCTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCtt                                                                                                                                                               <  2:577607/275‑1 (MQ=255)

                                                                                        cGACCGCGAATAAATTCGATGCCAACCAGTTCAAAACCCAGGGCCTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGta                                                                                                                                                        >  2:635897/1‑274 (MQ=255)

                                                                                                    aaTTCGATGCCAACCAGTTCAAAACCCAGGGCCTCAACTGGCGCATTAATCATCTCTGTTAATTTTTGCTCTAATATGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGACAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAg                                                                                                                                            <  2:161002/274‑1 (MQ=255)

                                                                                                         gATGCCAACCAGTTCAAAACCCAGGGCCTCAACTTGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGAct                                                                                                                                        <  1:11179/273‑1 (MQ=255)

                                                                                                                     ttCAAAACCCAGGGCCTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATTTATCggg‑gtt                                                                                                                                                                                                                                                        <  6:1407760/149‑3 (MQ=255)

                                                                                                                     ttCAAAACCCAGGGCCTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATTTATCggg‑gtt                                                                                                                                                                                                                                                        <  6:179193/149‑3 (MQ=255)

                                                                                                                     ttCAAAACCCAGGGCCTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATTTATCggg‑gtt                                                                                                                                                                                                                                                        <  5:1234496/149‑3 (MQ=255)

                                                                                                                     ttCAAAACCCAGGGCCTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATTTATCggg‑gtt                                                                                                                                                                                                                                                        <  7:1123379/149‑3 (MQ=255)

                                                                                                                     ttCAAAACCCAGGGCCTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCt                                                                                                                                                                                                                                                        >  6:892827/1‑149 (MQ=255)

                                                                                                                      tCAAAACCCAGGGCCTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATTTATCGGG‑Gttt                                                                                                                                                                                                                                                       <  5:1074913/149‑1 (MQ=255)

                                                                                                                      tCAAAACCCAGGGCCTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCtt                                                                                                                                                                                                                                                       <  5:348899/149‑1 (MQ=255)

                                                                                                                       cAAAACCCAGGGTCTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTATTGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATTTATCGGG‑Gtttt                                                                                                                                                                                                                                                      <  5:686887/149‑1 (MQ=255)

                                                                                                                        aaaaCCCAGGGCCTCAACTGGCGCAGTAATCGTCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAGGGGCCTAATGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTta                                                                                                                                                                                                                                                     >  8:1035459/1‑149 (MQ=255)

                                                                                                                        aaaaCCCAGGGCCTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTta                                                                                                                                                                                                                                                     <  5:265319/149‑1 (MQ=255)

                                                                                                                        aaaaCCCAGGGCCTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTta                                                                                                                                                                                                                                                     >  6:62045/1‑149 (MQ=255)

                                                                                                                          aaCCCAGGGCCTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTtata                                                                                                                                                                                                                                                   <  7:592530/149‑1 (MQ=255)

                                                                                                                                gggCCTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTg                                                                                                                                                                                                                                             <  8:1400761/149‑1 (MQ=255)

                                                                                                                                    cTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCtt                                                                                                                                                                                                                                                       <  5:436389/135‑1 (MQ=255)

                                                                                                                                    cTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGAAcc                                                                                                                                                                                                                                         >  7:28343/1‑149 (MQ=255)

                                                                                                                                     tCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGAAccc                                                                                                                                                                                                                                        <  6:326285/149‑1 (MQ=255)

                                                                                                                                       aaatgggcaaGTAATCATTTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCAAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGACTCTAAGGGAAAGTGGTTGCGGGGGCCGGATTTg                                                                                                        <  4:635892/266‑1 (MQ=255)

                                                                                                                                         cTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGAACCCtata                                                                                                                                                                                                                                    >  6:962943/1‑149 (MQ=255)

                                                                                                                                             cgcAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACAT‑AAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAAc                                                                                                                                                                                                                                               >  4:150704/1‑133 (MQ=255)

                                                                                                                                                    atcatcTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGc                                                                                                                                                                                                                         >  7:709332/1‑149 (MQ=255)

                                                                                                                                                    atcatcTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGc                                                                                                                                                                                                                         <  5:1532186/149‑1 (MQ=255)

                                                                                                                                                       atcTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGt                                                                                                                                                                                                                      >  7:456716/1‑149 (MQ=255)

                                                                                                                                                          tctGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCt                                                                                                                                                                                                                    >  2:359248/1‑148 (MQ=255)

                                                                                                                                                            tGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGACTCTAAGGGAAAGTGGTTGCGGGGGCAGGATTTGAACCGACGATCTTCGGGTTAt                                                                                   >  2:68080/1‑275 (MQ=255)

                                                                                                                                                              ttAATTTTTGCTGTAATGTGAACAAGCCCACCCCCCAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAACCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACGGCGGTCTGGAg                                                                                                                                                                                                                <  4:359244/148‑1 (MQ=255)

                                                                                                                                                                aaTTTTTGCACCAATGTGGACAAGCCCACCCCCAAGGCATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTATATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAg                                                                                                                                            <  4:58382/214‑1 (MQ=255)

                                                                                                                                                                 aTTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACt                                                                                                                                                                                                            >  8:1191175/1‑149 (MQ=255)

                                                                                                                                                                  ttgttGCTCTAATGTGAACACGCCCACCCCCAAAACAGAAAAAAAGGGCCTAAGGCCCAGTGATTCTGTAGTCAGATAACAAAAACCCCCGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGACTCTAAGGGAAAGTGGTTGCGGGGGCCGGATTTGAACCGACGATCTTCGGGTTATGAGCCcgacg                                                                         <  1:100265/276‑1 (MQ=255)

                                                                                                                                                                    tttGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTc                                                                                                                                                                                                         <  8:887975/149‑1 (MQ=255)

                                                                                                                                                                    tttGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTc                                                                                                                                                                                                         <  8:468753/149‑1 (MQ=255)

                                                                                                                                                                    tttGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTc                                                                                                                                                                                                         <  6:26621/149‑1 (MQ=255)

                                                                                                                                                                     ttGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTcc                                                                                                                                                                                                        <  7:514298/149‑1 (MQ=255)

                                                                                                                                                                         tctAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGaa                                                                                                                                                                                                    <  7:1550046/149‑1 (MQ=255)

                                                                                                                                                                         tctAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGaa                                                                                                                                                                                                    >  5:527637/1‑149 (MQ=255)

                                                                                                                                                                           tAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCT‑GAGCACTTTCCAGAAGGa                                                                                                                                                                                                 >  8:599549/1‑149 (MQ=255)

                                                                                                                                                                             aTGTGGACGAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGaa                                                                                                                                                                                                                                           <  5:891367/106‑1 (MQ=255)

                                                                                                                                                                             aTGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTAGATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGAt                                                                                                                                                                                                >  5:1555553/1‑149 (MQ=255)

                                                                                                                                                                              tgtgGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGAtt                                                                                                                                                                                               <  8:1003094/149‑1 (MQ=255)

                                                                                                                                                                                tgGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGAtttt                                                                                                                                                                                             <  7:1304825/149‑1 (MQ=255)

                                                                                                                                                                                  gACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGAtttt                                                                                                                                                                                             >  5:1674336/1‑147 (MQ=255)

                                                                                                                                                                                     aaGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCaa                                                                                                                                                                                        <  7:1604650/149‑1 (MQ=255)

                                                                                                                                                                                       gCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑Gc                                                                                                                                                                                                                                                         >  2:764480/1‑82 (MQ=255)

                                                                                                                                                                                               ccAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTAc                                                                                                                                                                              >  6:1651135/1‑149 (MQ=255)

                                                                                                                                                                                               ccAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAACACTTTCCAGAAGGATTTTTTCAAATCCCACTAc                                                                                                                                                                              >  7:1379641/1‑149 (MQ=255)

                                                                                                                                                                                                   gACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAg                                                                                                                                                                          >  8:865161/1‑149 (MQ=255)

                                                                                                                                                                                                      aTAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGt                                                                                                                                                                                                                      >  6:1116778/1‑102 (MQ=255)

                                                                                                                                                                                                      aTAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGt                                                                                                                                                                                                                      <  5:1116778/102‑1 (MQ=255)

                                                                                                                                                                                                       tAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTta                                                                                                                                                                                                                                                     >  3:189575/1‑70 (MQ=255)

                                                                                                                                                                                                       tAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTta                                                                                                                                                                                                                                                     <  1:189575/70‑1 (MQ=255)

                                                                                                                                                                                                       tAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCg                                                                                                                                                                      >  5:32237/1‑149 (MQ=255)

                                                                                                                                                                                                        aaaaaaaGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGa                                                                                                                                                                     >  8:229728/1‑149 (MQ=255)

                                                                                                                                                                                                        aaaaaaaGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGACTCTAAGGGAAAGTGGTTGCGGGGGCCGGATTTGAAc                                                                                                     >  4:266103/1‑213 (MQ=255)

                                                                                                                                                                                                           aaaaGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGt                                                                                                                                                                  >  6:853574/1‑149 (MQ=255)

                                                                                                                                                                                                           aaaaGGGCCTAAAGCCCAGTTATTCTGTAGCCAGATAACAAAAAACCCCGATAAAACGGG‑GTTTTATATAACTGAACCCTATAACGTCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACACCAAAGGCCGAAGt                                                                                                                                                                  >  7:861795/1‑149 (MQ=255)

                                                                                                                                                                                                               gggCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCTGGCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGcc                                                                                                                                                                       >  6:158504/1‑140 (MQ=255)

                                                                                                                                                                                                                 gCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTcac                                                                                                                                                            >  7:209810/1‑149 (MQ=255)

                                                                                                                                                                                                                  ccTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTcaca                                                                                                                                                           <  8:1343820/149‑1 (MQ=255)

                                                                                                                                                                                                                   cTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAg                                                                                                                                                          >  6:814294/1‑149 (MQ=255)

                                                                                                                                                                                                                    tAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGa                                                                                                                                                                                                                 <  6:730588/93‑1 (MQ=255)

                                                                                                                                                                                                                    tAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGa                                                                                                                                                                                                                 >  5:730588/1‑93 (MQ=255)

                                                                                                                                                                                                                    tAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGt                                                                                                                                                                  >  7:1021061/1‑140 (MQ=255)

                                                                                                                                                                                                                      aaGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAg                                                                                                                                                                                                                <  1:140130/92‑1 (MQ=255)

                                                                                                                                                                                                                      aaGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAg                                                                                                                                                                                                                >  3:140130/1‑92 (MQ=255)

                                                                                                                                                                                                                           cAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAACGGACCCTAAGGGAAAGTGGTTGCGGGGGCCGGATTTGAACCGACGATCTTCGGGTTATGAGCCGGACGAGCTACCAGGCTGCTCCACCCCGCGCCTGAAACGTGGCAAATTCTCCTCGtat                    >  2:3966/1‑273 (MQ=255)

                                                                                                                                                                                                                                aTTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAACTCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGaaaaa                                                                                                                                             <  5:647113/149‑1 (MQ=255)

                                                                                                                                                                                                                                aTTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGaaaaa                                                                                                                                             <  6:427915/149‑1 (MQ=255)

                                                                                                                                                                                                                                 ttCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAg                                                                                                                                            <  8:709332/149‑1 (MQ=255)

                                                                                                                                                                                                                                   cctgaaacacaTAAAAAAAAACCCC‑AGAAANNGGGTGCTGTATATAACTGAAACTTATAACCNCAACTGCGGACTGGCGCACTGTCCAGAATGATTTTTTCAAAGCCCACTACGAGGGCCGATGTCTTCACAGTATANTTGAAAAAGGACTCTAAGGGAAAGTGGTTGCGGGGGCAGGATTTGAACCGACTATCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGCCTGAAACGTGGCAAATTCTACTCGTTTTGGGTaaaa           <  3:118346/266‑1 (MQ=37)

                                                                                                                                                                                                                                   cTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTcc                                                                                                                                                                                                        >  8:1604238/1‑87 (MQ=255)

                                                                                                                                                                                                                                   cTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGa                                                                                                                                          <  5:1201749/149‑1 (MQ=255)

                                                                                                                                                                                                                                     gTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGAct                                                                                                                                        <  6:1352286/149‑1 (MQ=255)

                                                                                                                                                                                                                                        gTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGaa                                                                                                                                                                    >  6:1006719/1‑118 (MQ=255)

                                                                                                                                                                                                                                        gTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGaa                                                                                                                                                                    <  5:1006719/118‑1 (MQ=255)

                                                                                                                                                                                                                                           aGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTCTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGACTCTAAgg                                                                                                                                  >  5:232022/1‑149 (MQ=255)

                                                                                                                                                                                                                                             aTAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCAATACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGACTCTAAGGGAAAGTGGTTGCGGGGGCCGGATTGGAACCGACGATCTTCGGGTTATGATCTCGACGAGCTCCAAGGCTGCTCCACCCAGCCCCTGAAACTGGGCAAATTCCACTCGTTTTGGGTaaaaaaggaaact  >  4:102082/1‑268 (MQ=255)

                                                                                                                                                                                                                                               aaCAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACt                                                                                                                                                                                                                           <  4:760712/56‑1 (MQ=255)

                                                                                                                                                                                                                                               aaCAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGACTCTAAGGGaaa                                                                                                                              <  8:54854/149‑1 (MQ=255)

                                                                                                                                                                                                                                                  aaaaaaCCCCGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGAtttt                                                                                                                                                                                             <  5:585174/83‑1 (MQ=255)

                                                                                                                                                                                                                                                  aaaaaaCCCCGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGACTCTAAGGGAAAGTg                                                                                                                           >  7:960083/1‑149 (MQ=255)

                                                                                                                                                                                                                                                    aaagCCCCGATTTATCGGG‑GTTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGACTCTAAGGGAAAGTGGt                                                                                                                         >  8:931632/1‑149 (MQ=255)

                                                                                                                                                                                                                                                    aaagCCCCGATTTATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGACTCTAAGGGAAAGTGGt                                                                                                                         >  7:684286/1‑149 (MQ=255)

                                                                                                                                                                                                                                                    aaagCCCCGATTTATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGACTCTAAGGGAAAGTGGt                                                                                                                         >  5:1623014/1‑149 (MQ=255)

                                                                                                                                                                                                                                                          ccGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGaaaaa                                                                                                                                             <  8:410443/123‑1 (MQ=255)

                                                                                                                                                                                                                                                          ccGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGACTCTAAGGGAAAGTGGTTGCggg                                                                                                                   >  6:36119/1‑149 (MQ=255)

                                                                                                                                                                                                                                                          ccGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGACTCTAACGGAAAGTGGTTGCggg                                                                                                                   >  6:505436/1‑149 (MQ=255)

                                                                                                                                                                                                                                                           cGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGaa                                                                                                                                                                                                    >  8:685758/1‑67 (MQ=255)

                                                                                                                                                                                                                                                           cGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGACTCTAAGGGAAAGTGGTTGCgggg                                                                                                                  <  5:836682/149‑1 (MQ=255)

                                                                                                                                                                                                                                                            gATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGaaaaa                                                                                                                                             >  5:395441/1‑121 (MQ=255)

                                                                                                                                                                                                                                                                aaTCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGACACTAAGGGAAAGTGGTTGCGGGCGTCgg                                                                                                             >  8:1088617/1‑149 (MQ=255)

CACAGCACTTACCTGGTGGCTCACATCAGCACAATCATCAACATTGATGCCATCTTCACTATCAATATAGATGCGCAGTGTGGATGTGCGACCGCGAATAAATTCGATGCCAACCAGTTCAAAACCCAGGGCCTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGACTCTAAGGGAAAGTGGTTGCGGGGGCCGGATTTGAACCGACGATCTTCGGGTTATGAGCCCGACGAGCTACCAGGCTGCTCCACCCCGCGCCTGAAACGTGGCAAATTCTACTCGTTTTGGGTAAAAAATGCAAAT  >  NC_000913/3317824‑3318335

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