BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,318,079	A→T	5.4%	intergenic (‑73/+134)	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,079	0	A	T	5.4%	979.6 / 32.1	275	intergenic (‑73/+134)	rimP/metY	ribosome maturation factor for 30S subunits/tRNA‑Met
Reads supporting (aligned to +/- strand): ref base A (134/126); 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 = 2.25e-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)

                                                                                                               aaCCAGTTCAAAACCCAGGGCCTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATc                                                                                                                                                                                                                                                               >  6:340596/1‑149 (MQ=255)

                                                                                                                  cAGTTCAAAACCCAGGGCCTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCggg                                                                                                                                                                                                                                                            <  6:220173/149‑1 (MQ=255)

                                                                                                                     ttCAAAACCCAGGGCCTCAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATTTATCggg‑gtt                                                                                                                                                                                                                                                        <  5:1234496/149‑3 (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)

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

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

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

                                                                                                                                      cAACTGGCGCAGTAATCATCTCTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGAACCCt                                                                                                                                                                                                                                       <  6:1380788/149‑1 (MQ=255)

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

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

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

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

                                                                                                                                                        tctcTGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATTTATCGGG‑Gtttt                                                                                                                                                                                                                                                      <  7:893240/116‑1 (MQ=255)

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

                                                                                                                                                          tctGTTAATTTTTGCTCTAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTg                                                                                                                                                                                                                   >  7:1003094/1‑149 (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:468753/149‑1 (MQ=255)

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

                                                                                                                                                                         tctAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGTG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGaa                                                                                                                                                                                                    >  6:836682/1‑149 (MQ=255)

                                                                                                                                                                         tctAATGTGGACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGaa                                                                                                                                                                                                    >  6:792259/1‑149 (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)

                                                                                                                                                                             aTGTGGACGAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGaa                                                                                                                                                                                                                                           >  6:891367/1‑106 (MQ=255)

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

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

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

                                                                                                                                                                                  gACAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGAtttt                                                                                                                                                                                             <  6:1674336/147‑1 (MQ=255)

                                                                                                                                                                                   aCAAGCCCACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTc                                                                                                                                                                                          >  6:1067248/1‑149 (MQ=255)

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

                                                                                                                                                                                         ccACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATcc                                                                                                                                                                                    >  6:1573123/1‑149 (MQ=255)

                                                                                                                                                                                          cACCCCCAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATccc                                                                                                                                                                                   >  7:660112/1‑149 (MQ=255)

                                                                                                                                                                                             ccccAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACt                                                                                                                                                                                >  7:1343820/1‑149 (MQ=255)

                                                                                                                                                                                               ccAAGACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCa                                                                                                                                                                                  >  7:15260/1‑145 (MQ=255)

                                                                                                                                                                                                   gACATAAAAAAAGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAg                                                                                                                                                                          >  8:865161/1‑149 (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‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACg                                                                                                                                                                             >  6:725500/1‑138 (MQ=255)

                                                                                                                                                                                                            aaaGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTc                                                                                                                                                                 >  6:1588275/1‑149 (MQ=255)

                                                                                                                                                                                                              aGGGCCTAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCtt                                                                                                                                                               >  7:54854/1‑149 (MQ=255)

                                                                                                                                                                                                                 gCCTAAAGCCCATTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTcac                                                                                                                                                            >  8:1057800/1‑149 (MQ=255)

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

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

                                                                                                                                                                                                                    tAAAGCCCAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGt                                                                                                                                                                  <  8:1021061/140‑1 (MQ=255)

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

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

                                                                                                                                                                                                                          ctaGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAGTCGGG‑GCTTTAAATAACTGAACCCTAGAACCGCAACTGCGGTCAGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATAttt                                                                                                                                                   >  6:1215210/3‑149 (MQ=255)

                                                                                                                                                                                                                          ccAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATAttt                                                                                                                                                   <  6:496842/149‑1 (MQ=255)

                                                                                                                                                                                                                          ccAGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATAttt                                                                                                                                                   >  7:1194774/1‑149 (MQ=255)

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

                                                                                                                                                                                                                            aGTTATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGa                                                                                                                                                 >  7:479954/1‑149 (MQ=255)

                                                                                                                                                                                                                              ttATTCTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGaaa                                                                                                                                               <  7:254481/149‑1 (MQ=255)

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

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

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

                                                                                                                                                                                                                                   cTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTcc                                                                                                                                                                                                        <  7:1604238/87‑1 (MQ=255)

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

                                                                                                                                                                                                                                   cTGTAGTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGa                                                                                                                                          <  7:1673841/149‑1 (MQ=255)

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

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

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

                                                                                                                                                                                                                                        gTCAGATAACAAAAAACCCCGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTca                                                                                                                                                             >  6:310098/1‑125 (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                                                                                                                                                                                             >  6:585174/1‑83 (MQ=255)

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

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

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

                                                                                                                                                                                                                                                    aaaaCCCCGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGACTCTAAGGGAAAGTGGt                                                                                                                         <  7:963935/149‑1 (MQ=255)

                                                                                                                                                                                                                                                          ccGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGaaaaa                                                                                                                                             >  7:410443/1‑123 (MQ=255)

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

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

                                                                                                                                                                                                                                                           cGATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGaa                                                                                                                                                                                                    <  7:685758/67‑1 (MQ=255)

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

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

                                                                                                                                                                                                                                                            gATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGaaaaa                                                                                                                                             <  6:395441/121‑1 (MQ=255)

                                                                                                                                                                                                                                                            gATAAATCGGG‑GCTTTATATAACTGAACCCTATAACCGCAACTGCGGTCTGGAGCACTTTCCAGAAGGATTTTTTCAAATCCCACTACGAAGGCCGAAGTCTTCACAGTATATTTGAAAAAGGACTCTAAGGGAAAGTGGTTGCgg                                                                                                                    <  6:908099/146‑1 (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` `< 24 ≤` `ATCG/ATCG` `< 36 ≤` `ATCG/ATCG` `< 40 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`