BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,234,506	G→A	6.9%	intergenic (‑51/+233)	higB ← / ← rlmG	mRNA interferase toxin of the HigB‑HigA toxin‑antitoxin system/23S rRNA m(2)G1835 methyltransferase, SAM‑dependent

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,234,506	0	G	A	6.9%	650.7 / 26.1	205	intergenic (‑51/+233)	higB/rlmG	mRNA interferase toxin of the HigB‑HigA toxin‑antitoxin system/23S rRNA m(2)G1835 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand): ref base G (88/101); new base A (7/7); total (96/109)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.22e-01

Only 100 of 206 total aligned reads displayed.

TGCCCCCAACATTGAAAACATAATGCTTATCCAGATATTTGAAGTTATCCAGAGATGGGAATACTGCTTTTAATGACTCAGGTTTTTTGAAATATCCCTTAGCAATCGTGTTCCCCAGAGCCACCAACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGAACTGGGCAGACATTGCCGGATGCGACGCTTGCCGCGTCTTATCCGGCCTGGAGTAACGTGCTGAACGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGGTAACTCCCAAT  >  NC_000913/3234263‑3234737

           ttGAAAACATAATGCTTATCCAGATATTTGAAGTTATCCAGAGATGGGAATACTGCTTTTAATGACTCAGGTTTTTTGAAATATCCCTTAGCAATCTTGTTCCACAGGGCCACCAAATCCGTTTTATGTTGCGTGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTTTCTGTGTGACGGAAGAATAAAc                                                                                                                                                                                               <  2:217644/275‑1 (MQ=255)

                              ctaGATATTTGAAGTTATCCAGAGATGGGAATACTGCTTTTAATGACTCAGGTTTTTTGAAATATCCCTTAGCAATCGTGTTCCCCAGAGCCACCAACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACaa                                                                                                                                                                            <  4:169740/273‑1 (MQ=255)

                                       tGAAGTTATCCAGAGATGGGAATACTGCTTTTAATGACTCAGGTTTTTTGAAATATCCCTTAGCAATCGTGTTCCCCAGAGCCACCAACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑Tg                                                                                                                                                                                                                                     >  4:299551/1‑209 (MQ=255)

                                           gTTATCCAGAGATGGGAATACTGCTTTTAATGACTCAGGTTTTTTGAAATATCCCTTAGCAATCGTGTTCCCCAGAGCCACCAACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCaa                                                                                                                                                                                                                                 <  2:299551/209‑1 (MQ=255)

                                              aTCCAGAGATGGGAATACTGCTTTTAATGACTCAGGTTTTTTGAAATATCCCTTAGCAATCGTGTTCCCCAGAGCCACCAACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGTCAATATGATAACATCCGTTGCCAGAT‑TTGCAACGGATGCATTTATCTGTTTGACGGAAGAATAACCGTAGAAAAACGCTAAACAAGAGAATGTAAACGATg                                                                                                                                                            >  2:249195/1‑275 (MQ=255)

                                                                     ttAATGACTCAGGTTTTTTGAAATATCCCTTAGCAATCGTGTTCCCCAGAGCCACCAACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTAt                                                                                                                                     >  2:29248/1‑275 (MQ=255)

                                                                       aaTGACTCAGGTTTTTTGAAATATCCCTTAGCAATCGTGTTCCCCAGAGCCACCAACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTAGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACTATGGTGACCGGGATTTTTTGTTTTATgg                                                                                                                                   >  2:9516/1‑275 (MQ=255)

                                                                                               cccTTAGCAATCGTGTTCCCCAGAGCCACCAACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAg                                                                                                                                                                                                                                          >  5:1354653/1‑149 (MQ=255)

                                                                                                ccTTAGCAATCGTGTTCCCCAGAGCCACCAACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGa                                                                                                                                                                                                                                         >  6:1454804/1‑149 (MQ=255)

                                                                                                ccTTAGCAATCGTGTTCCCCAGAGCCACCAACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGa                                                                                                                                                                                                                                         >  8:657690/1‑149 (MQ=255)

                                                                                                 cTTAGCAATCGTGTTCCCCAGAGCCACCAACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAt                                                                                                                                                                                                                                        <  6:791952/149‑1 (MQ=255)

                                                                                                 cTTAGCAATCGTGTTCCCCAGAGCCACCAACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAt                                                                                                                                                                                                                                        <  8:1067822/149‑1 (MQ=255)

                                                                                                  ttAGCAATCGTGTTCCCCAGAGCCACCAACACCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAt‑t                                                                                                                                                                                                                                      <  5:275299/149‑1 (MQ=255)

                                                                                                        aTCGTGTTCCCCAGAGCCACCAACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGGGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGAACTGGGCAGACATTGcc                                                                                                   >  1:14419/1‑274 (MQ=255)

                                                                                                           gtgtTCCCCAGAGCCACCAACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGa                                                                                                                                                                                                                             <  7:841353/149‑1 (MQ=255)

                                                                                                           gtgtTCCCCAGAGCCACCAACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGa                                                                                                                                                                                                                             <  7:21894/149‑1 (MQ=255)

                                                                                                           gtgtTCCCCAGAGCCACCAACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGa                                                                                                                                                                                                                             >  6:1444704/1‑149 (MQ=255)

                                                                                                              ttCCCCAGAGCCACCAACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGAACTGGGCAGACATTGCCGGATg                                                                                              >  2:330290/1‑273 (MQ=255)

                                                                                                               tCCCCAGAGCCACCAACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAATC‑TGGCAACGGATgtt                                                                                                                                                                                                                         <  5:447333/149‑3 (MQ=255)

                                                                                                               tCCCCAGAGCCACCAACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAATC‑TGGCAACGGATgtt                                                                                                                                                                                                                         <  7:179525/149‑3 (MQ=255)

                                                                                                                ccccAGAGCCACCAACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAAt                                                                                                                                                                                                                        >  5:1410905/1‑149 (MQ=255)

                                                                                                                  ccAGAGCCACCAACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAAttt                                                                                                                                                                                                                      <  6:454313/149‑1 (MQ=255)

                                                                                                                  ccAGAGCCACCAACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAAttt                                                                                                                                                                                                                      <  8:295176/149‑1 (MQ=255)

                                                                                                                      agCCACCAACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCt                                                                                                                                                                                                                  <  6:718381/149‑1 (MQ=255)

                                                                                                                        ccaccaACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCATTATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCtgt                                                                                                                                                                                                                <  8:707263/149‑1 (MQ=255)

                                                                                                                        ccaccaACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑Tg                                                                                                                                                                                                                                     <  8:1398471/128‑1 (MQ=255)

                                                                                                                          accaACTCCGTTTTATATTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCtgtgt                                                                                                                                                                                                              >  8:1388570/1‑149 (MQ=255)

                                                                                                                           ccaACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGAAAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAAAACTGGGCAGACATTGCCGGATGCGAAGTTTGCcgcct                                                                               >  4:315338/1‑273 (MQ=255)

                                                                                                                             aaCTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGAc                                                                                                                                                                                                           <  7:677533/149‑1 (MQ=255)

                                                                                                                             aaCTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGAc                                                                                                                                                                                                           >  8:876485/1‑149 (MQ=255)

                                                                                                                             aaCTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGAc                                                                                                                                                                                                           >  7:935653/1‑149 (MQ=255)

                                                                                                                                tCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGga                                                                                                                                                                                                        >  6:814640/1‑149 (MQ=255)

                                                                                                                                tCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGga                                                                                                                                                                                                        >  6:105691/1‑149 (MQ=255)

                                                                                                                                tCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAATC‑TGGCAACgg                                                                                                                                                                                                                              >  5:719511/1‑127 (MQ=255)

                                                                                                                                 ccGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGgaa                                                                                                                                                                                                       <  5:656332/149‑1 (MQ=255)

                                                                                                                                 ccGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGgaa                                                                                                                                                                                                       <  6:990067/149‑1 (MQ=255)

                                                                                                                                  cGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAGATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGgaag                                                                                                                                                                                                      >  6:974350/1‑149 (MQ=255)

                                                                                                                                   gTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAATC‑TGGCAACGGAt                                                                                                                                                                                                                            >  6:845306/1‑126 (MQ=255)

                                                                                                                                    ttttATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGgaagaa                                                                                                                                                                                                    <  8:192725/149‑1 (MQ=255)

                                                                                                                                    ttttATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGgaagaa                                                                                                                                                                                                    <  6:666341/149‑1 (MQ=255)

                                                                                                                                    ttttATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGgaagaa                                                                                                                                                                                                    >  5:964824/1‑149 (MQ=255)

                                                                                                                                        aTGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATaaa                                                                                                                                                                                                <  7:554628/149‑1 (MQ=255)

                                                                                                                                        aTGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATaaa                                                                                                                                                                                                <  5:1454804/149‑1 (MQ=255)

                                                                                                                                        aTGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATaaa                                                                                                                                                                                                <  5:1444704/149‑1 (MQ=255)

                                                                                                                                         tGTTGCGGGTATTTTTCCGCAGTATCTTATAATGCTTTTTGAGTTATCAGGTGCATTCTTTATCACGTCCGTTTACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAAc                                                                                                                                                                                               <  7:1147722/149‑1 (MQ=255)

                                                                                                                                          gTTGCGGGTATTTTTCCGCAGCCTCTTCCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATACTATCCGTTGCCCGAT‑TGGCAACGGATGAATTTATCTGTGTTACGGAAGAATTAACGTCGAAAAACGCTAAACAATAGAATGTAGACGATGGTGACAGGGATCTTGTGTTTTATGGAAGCGCTATTCTCGAACTGGGCAGACATTGCCGGAATCTACGCATTCCGCGTTTTAACAGGCCTggc                                                                >  4:287073/1‑274 (MQ=255)

                                                                                                                                             gCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAg                                                                                                                                                                                           <  7:1066045/149‑1 (MQ=255)

                                                                                                                                                 gTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGaaaa                                                                                                                                                                                       >  7:769003/1‑149 (MQ=255)

                                                                                                                                                  tATTTTTCCGCAGCATTTTTCAATTCTTTTTGAGTTATCAGGGTCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTTACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGAACTGGGCAGACATTGCCGGATGCGACGCTTGCCGCGTCTTATCCGGCCTGGAGTAACGTg                                                        <  4:226175/275‑1 (MQ=255)

                                                                                                                                                  tATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAAt                                                                                                                                                                                                   >  5:1218946/1‑136 (MQ=255)

                                                                                                                                                   aTTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTTGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAAc                                                                                                                                                                                     <  6:284724/149‑1 (MQ=255)

                                                                                                                                                   aTTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAAc                                                                                                                                                                                     <  8:1454282/149‑1 (MQ=255)

                                                                                                                                                   aTTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAAc                                                                                                                                                                                     <  6:297258/149‑1 (MQ=255)

                                                                                                                                                    tttttCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACg                                                                                                                                                                                    >  8:1349881/1‑149 (MQ=255)

                                                                                                                                                              gcaTCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAaga                                                                                                                                                                          <  5:79707/149‑1 (MQ=255)

                                                                                                                                                              gcaTCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAaga                                                                                                                                                                          <  8:820657/149‑1 (MQ=255)

                                                                                                                                                               caTCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAagag                                                                                                                                                                         >  8:552624/1‑149 (MQ=255)

                                                                                                                                                               caTCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGAACTGGGCAGACATTGCCGGATGCGACGCTTGCCGCGTCt                                                                             >  2:139777/1‑241 (MQ=255)

                                                                                                                                                                 tCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAAc                                                                                                                                                                                     <  7:1390358/135‑1 (MQ=255)

                                                                                                                                                                   tttCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGAACTGGGCAGACATTGCCGGATGCGACGCTTGCCGCGTCTTATc                                                                         <  4:139777/241‑1 (MQ=255)

                                                                                                                                                                         tGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAaga                                                                                                                                                                          <  7:1075244/138‑1 (MQ=255)

                                                                                                                                                                            tttttGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAAt                                                                                                                                                                                                                        <  6:1286349/89‑1 (MQ=255)

                                                                                                                                                                            tttttGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATg                                                                                                                                                            <  7:125607/149‑1 (MQ=255)

                                                                                                                                                                                 gAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGAc                                                                                                                                                       >  8:685217/1‑149 (MQ=255)

                                                                                                                                                                                      aTCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAaga                                                                                                                                                                          >  7:187273/1‑125 (MQ=255)

                                                                                                                                                                                       tCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGAt                                                                                                                                                 >  6:263906/1‑149 (MQ=255)

                                                                                                                                                                                           gTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCTGTTGCCGGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGAttttt                                                                                                                                             >  6:864452/1‑149 (MQ=255)

                                                                                                                                                                                            tGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGAtttttt                                                                                                                                            <  6:1114720/149‑1 (MQ=255)

                                                                                                                                                                                              cATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGt                                                                                                                                          >  6:1174875/1‑149 (MQ=255)

                                                                                                                                                                                               aTTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGtt                                                                                                                                         <  5:597030/149‑1 (MQ=255)

                                                                                                                                                                                                ttcttcATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGttt                                                                                                                                        >  5:634202/1‑149 (MQ=255)

                                                                                                                                                                                                 tcttcATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGAACTGGGCAGACATTGCCGGATGCGACGCTTGCCGCGTCTTATCTGGCCTGGAGTCACGTCCTGAACGTAGGCCGGATAAGGCGTTTCCGCCGCACCCGGCAGGTaa          >  2:36310/1‑274 (MQ=255)

                                                                                                                                                                                                  cttcATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTa                                                                                                                                      <  7:122847/149‑1 (MQ=255)

                                                                                                                                                                                                   ttcATCACGTCCGTTGACAAATTGGCAATATGATAACATCTGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGAtttttt                                                                                                                                            >  8:1337056/1‑142 (MQ=255)

                                                                                                                                                                                                   ttcATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTAt                                                                                                                                     <  5:814640/149‑1 (MQ=255)

                                                                                                                                                                                                   ttcATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTAt                                                                                                                                     <  8:557733/149‑1 (MQ=255)

                                                                                                                                                                                                       tcaCGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGaa                                                                                                                                 <  7:326129/149‑1 (MQ=255)

                                                                                                                                                                                                         aCGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAaga                                                                                                                                                                          <  8:1028442/106‑1 (MQ=255)

                                                                                                                                                                                                         aCGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAaga                                                                                                                                                                          >  7:1028442/1‑106 (MQ=255)

                                                                                                                                                                                                         aCGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAgc                                                                                                                               <  8:1262269/149‑1 (MQ=255)

                                                                                                                                                                                                         aCGTCCGTTGACAAATTGGCAATATGATAACAGCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGAACTGGGCAGACATTGCCGGATGCGACGCTTGCCGCGTCTTATCCGGCCTGGAGTAACGTGCTGAACGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGGTAACTCCCAAt  <  4:29248/274‑1 (MQ=255)

                                                                                                                                                                                                                 tGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTcc                                                                                                                       <  5:263906/149‑1 (MQ=255)

                                                                                                                                                                                                                    cAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGa                                                                                                                    <  6:403783/149‑1 (MQ=255)

                                                                                                                                                                                                                    cAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGa                                                                                                                    <  5:1289625/149‑1 (MQ=255)

                                                                                                                                                                                                                     aaaTTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGaa                                                                                                                   <  6:1314408/149‑1 (MQ=255)

                                                                                                                                                                                                                                 tGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGAACTGGGCAGACAt                                                                                                       >  6:548169/1‑149 (MQ=255)

                                                                                                                                                                                                                                  gATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGAACTGGGCAGACAtt                                                                                                      >  8:1008057/1‑149 (MQ=255)

                                                                                                                                                                                                                                   aTAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGAACTGGGCAGACATTg                                                                                                     >  8:118557/1‑149 (MQ=255)

                                                                                                                                                                                                                                    tAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGCCGGAAGAATAAACGTAGAAAATCGCTAAACAAGAGCATGTAGACGATGGTGACAGGGATTTTTTGTCTTATGGAAGCGCGATTCCAGAACTGGGCAGACATTGc                                                                                                    <  8:371823/149‑1 (MQ=255)

                                                                                                                                                                                                                                    tAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACTGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGAACTGGGCAGACATTGc                                                                                                    >  6:1219989/1‑149 (MQ=255)

                                                                                                                                                                                                                                    tAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGAACTGGGCAGACATTGc                                                                                                    <  7:1219174/149‑1 (MQ=255)

                                                                                                                                                                                                                                     aaCATCCGTTGCCAGAT‑TGGCAACGGATGTATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGAACTGGGCAGACATTGcc                                                                                                   >  8:1229952/1‑149 (MQ=255)

                                                                                                                                                                                                                                        aTCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGAACTGGGCAGACATTGCCGGa                                                                                                <  7:876485/149‑1 (MQ=255)

                                                                                                                                                                                                                                        aTCCGTTGCCAATC‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGt                                                                                                                                          >  6:578953/1‑107 (MQ=255)

                                                                                                                                                                                                                                        aTCCGTTGCCA‑ATCTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGt                                                                                                                                          <  5:578953/107‑1 (MQ=255)

                                                                                                                                                                                                                                         tCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGAACTGGGCAGACATTGCCGGAt                                                                                               <  7:1331408/149‑1 (MQ=255)

                                                                                                                                                                                                                                         tCCGTTGCCAATC‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGAACTgg                                                                                                               >  5:1083923/1‑133 (MQ=255)

                                                                                                                                                                                                                                          ccGTTGCCAATC‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTAt                                                                                                                                     <  7:460926/110‑1 (MQ=255)

                                                                                                                                                                                                                                          ccGTTGCCA‑ATCTGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTAt                                                                                                                                     >  8:460926/1‑110 (MQ=255)

                                                                                                                                                                                                                                                ccAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGAACTGGGCAGACATTGCCGGATGCGACGc                                                                                        <  6:619143/149‑1 (MQ=255)

TGCCCCCAACATTGAAAACATAATGCTTATCCAGATATTTGAAGTTATCCAGAGATGGGAATACTGCTTTTAATGACTCAGGTTTTTTGAAATATCCCTTAGCAATCGTGTTCCCCAGAGCCACCAACTCCGTTTTATGTTGCGGGTATTTTTCCGCAGCATCTTTCAATGCTTTTTGAGTTATCAGGTGCATTCTTCATCACGTCCGTTGACAAATTGGCAATATGATAACATCCGTTGCCAGAT‑TGGCAACGGATGAATTTATCTGTGTGACGGAAGAATAAACGTAGAAAAACGCTAAACAAGAGAATGTAGACGATGGTGACAGGGATTTTTTGTTTTATGGAAGCGCGATTCCAGAACTGGGCAGACATTGCCGGATGCGACGCTTGCCGCGTCTTATCCGGCCTGGAGTAACGTGCTGAACGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGGTAACTCCCAAT  >  NC_000913/3234263‑3234737

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