BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_007779	954,815	T→A	17.2%	Y25F (TAT→TTT)	focA ←	formate transporter FocA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007779	954,815	0	T	A	17.2%	154.3 / 8.0	73	Y25F (TAT→TTT)	focA	formate transporter FocA
Reads supporting (aligned to +/- strand): ref base T (44/16); new base A (2/11); total (46/27)
Fisher's exact test for biased strand distribution p-value = 1.66e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTAT‑AGA‑CA‑CCCGCCTCTTCGGCCACTTTGGCCATTGCAG‑CAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGCAAATACGGAGTAAATATTTGATTATCCAAATAAAAATA  >  NC_007779/954577‑955044

gcgcCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGCTTtat                                                                                                                                                                                                                                           >  1:1258229/1‑239 (MQ=255)

   cccACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTAT‑Ag                                                                                                                                                                                                                                        >  1:79757/1‑238 (MQ=255)

        tCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATTCCCAGAAAGAAAAGCTTTATCGGGTGTTTCGTTGCTTTTA‑TGA‑CA‑cc                                                                                                                                                                                                                                 >  2:1140798/1‑238 (MQ=255)

                    caacaaTCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCCGAAAGAAAGACCTAAGGGGGAGTTTCGTTGGTTTTT‑AGA‑CA‑CCCGCCCCTTCGGc                                                                                                                                                                                                                     >  2:449443/1‑238 (MQ=255)

                    caacaaTCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGATGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAAAACCGGCGGGAAACGCCAGAAAGAAAGACTTTAGCAGGGGTTTCGTTGCTGTTA‑ATA‑CA‑CCCACCTCTTctgac                                                                                                                                                                                                                    >  2:1487278/1‑235 (MQ=255)

                      acaaTCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGGAAACCCCAACTAGAAAATCTTTAAGGGAATTTTCGTTTCCTTTT‑AAA‑CA‑ACCGCCCCCTCGGccc                                                                                                                                                                                                                   >  2:758934/1‑237 (MQ=255)

                                  gTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGAAAGGCTTAAGCGGATGTTTTGTTGCTTTAT‑AGACCC‑CCCGCCCCTTCGCCCA‑TTTTGCCAttt                                                                                                                                                                                                       >  2:850665/1‑237 (MQ=255)

                                  gTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATCCCCAGATAGAAAGTCTTTAGCGGGGGTTTCGGTGCTTTAT‑AGG‑CA‑CCCTCCTCTTCTCCCACTTTCGCCACTgc                                                                                                                                                                                                      >  1:1434722/1‑239 (MQ=255)

                                   tGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTAT‑AGA‑CA‑CCCGCCTCTTCGGCCACTTTGGCCATTgca                                                                                                                                                                                                     >  1:35946/1‑239 (MQ=255)

                                    ggAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTAT‑AGA‑CA‑CCCGCCTCTTCGGCCACTTTGGCCATTtca                                                                                                                                                                                                     >  1:1683215/1‑238 (MQ=255)

                                        gTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTAT‑AGA‑CA‑CCCGCCTCTTCGGCCACTTTGGCCATTGCAG‑CAgg                                                                                                                                                                                               >  1:831730/1‑239 (MQ=255)

                                           aaGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTAT‑AGA‑CA‑CCCGCCTCCTCGGCCACTTTGGCCATTGCAG‑CAGGAAg                                                                                                                                                                                            >  1:1739222/1‑239 (MQ=255)

                                             gagaTCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTAT‑AGA‑CA‑CCCGCCTCTCCGGCCCCTTTGGCCATTGCAG‑CACGAAgaa                                                                                                                                                                                          >  2:442336/1‑237 (MQ=255)

                                               gaTCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATCGGCCGATAAAAAGTATTAAGAGGAGGTTTCGTTGCTTTAT‑TAG‑CA‑CCCGCCTACTCCGGCACCTTGGCCATTTCTG‑CATGAAGTaaa                                                                                                                                                                                        >  2:1129878/1‑239 (MQ=255)

                                               gaTCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTTT‑AGA‑CA‑CCCGCCTCTTCGGCCACTTTTGCCATTGCAG‑CAGGAACTaaa                                                                                                                                                                                        >  1:1441814/1‑239 (MQ=255)

                                               caTCGGCCCCGAAGAAAACAAAAGGAATCAGCCCCAGGAAGAACAAGAGGCCCCCAACCATTTTTGCAAGGCGGAGGGGATTTGGGCCTGTGCGAGGGGTTGTTGGAAAATAAAAGACGAATGGGATTGAGATAAAACCCCCGGCGTAAATCGCCAGAAAGAAATTTTAACGGGGATGTTTCGTTGTTTTAT‑AGC‑CC‑CCCGCCTCTTCGGCCACTTTGGCCATTGCAG‑CAGGAAGTaaa                                                                                                                                                                                        <  2:1314592/238‑1 (MQ=255)

                                                  cGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGGAATCACCAGATAGAAAGACTTAAGGGGGTGTTTTGTTGCTTTATAAGA‑AA‑CCCGCC‑CTTCCGCCCCTTTTGCCATTCCAGCCGGGAAGAaaaaa                                                                                                                                                                                      >  2:1009148/1‑238 (MQ=255)

                                                     cTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTAT‑AGA‑CA‑CCCGCCTCTTCGGCCACTTTGGCCATTGCAG‑CAGGAAGTAAAAGATCa                                                                                                                                                                                  >  1:1799827/1‑239 (MQ=255)

                                                          cAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTAT‑AGA‑CA‑CCCGCCTCTTCGGCCACTTTGGCCATTGCAG‑CAGGAAGTAAAAGATCAAAAgg                                                                                                                                                                             >  1:1174695/1‑239 (MQ=255)

                                                                acacAAAGAACCAGCCCAAGAGAGAAGAAAATGCCGCCAACCAGTTTTGCCATGCGGAAGGGCTTTGTGCCTGGGCCAGTGGTTGCTGTAATATAGAAGACGAATGCGTTTGAGATAAAACCCCGGGCGGTAATCGCAAGATAGAAAGTCTTAAGCGGATGTTTCTTTGCTTTAT‑AGA‑CC‑CCCGCCTCTTCGGCCACTTTGGCCATTGCAG‑CAGGAAGTAAAAGATCAAAAGGGTTGTc                                                                                                                                                                       <  1:850665/239‑1 (MQ=255)

                                                                  acaAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGACAACACCGGCGGTAATAGCCAGATAGAAAAGCTTAAGCGGGTGTTTTGTTTGTTTAT‑TGA‑CC‑CCCGGCTCTTTGGCCAACTTGGGCCTTTCCG‑TCGGAAAGAAAAGATCAAAAaggggcgta                                                                                                                                                                      >  2:874086/1‑232 (MQ=255)

                                                                      aGAACCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCAGGCCGAGGGGCTTTGTGCCTGTGCCAGTGTTTGTTGTAATATAGAAGACGAATGCTATTGAGATAAAAACCCCGGCGGTAATCGCCAGATAGAAATTCTTAAGCGGATTTTTCGTTGCTTTAT‑AGA‑CA‑CCCGCCTCTTCGGCCACTTTGGCCATTGCAG‑CAGGAAGTAAAAGATCAAAAGGGTTGTCAGCttt                                                                                                                                                                 <  1:1129878/239‑1 (MQ=255)

                                                                         aTCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAAATAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTTT‑AGA‑AA‑CCCCCCTCTTCCGCCACTTTGGCCATTTCAG‑CAAGAAGTAAAAAACACAAAGGGTTTGCAGCTTTtac                                                                                                                                                              >  2:217126/1‑239 (MQ=255)

                                                                           cAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTAT‑AGA‑CA‑CCCCCCCCTTCGGCCCCTTTGGCCATTGCAG‑CAGGAAGTAAAAAACCAAAAGGGTTGTCGGCTTTcac                                                                                                                                                              >  1:1470140/1‑237 (MQ=255)

                                                                            aGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTAT‑AGA‑CA‑CCCGCCTCTTCGGCCACTTTGGCCATTGCAG‑CAGGAAGTAAAAGATCAAAAg                                                                                                                                                                              >  1:1009900/1‑220 (MQ=255)

                                                                            aGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTAT‑AGA‑CA‑CCCGCCTCTTCGGCCACTTTGGCCATTGCAG‑CAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCAAACt                                                                                                                                                           >  1:1635802/1‑239 (MQ=255)

                                                                              ccccAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTAT‑AGA‑CA‑CCCGCCTCTTCGGCCACTTTGGCCATTGCAG‑CAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTa                                                                                                                                                          >  1:748896/1‑238 (MQ=255)

                                                                                caagagagAAGAAAATGCCGCCAACCAGTTTGGCCATGCCGAGGGGCATGGTGCCTGTGCCAGTGGTTGCTGTAATAAAAAAGACAAAGGCGATTGAGATAAAAACCCCGGCGGTAATCGCCAGAAAGAAAGTCTAAAGCGGTGGTTTCTTTGTTTTAT‑AGA‑CA‑CCCGCCTCTTCGGCCACTTTGGCCATTGCAG‑CAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAAct                                                                                                                                                       <  1:675304/237‑1 (MQ=255)

                                                                                  agagagAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGGTAGAAAGTCTTAAGCGGGAGTTTCGTTGGTTtat‑a                                                                                                                                                                                                                                         >  2:1470855/1‑158 (MQ=255)

                                                                                        aaGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTAT‑AGA‑CA‑CCCGCCTCTTCGGCCACTTTGGCCATTGCAG‑CAGGAAGTAAAAGATCAAAAGGGGTGTCAGCTTTCACACTAACTCTCTCttt                                                                                                                                               >  1:1804381/1‑239 (MQ=255)

                                                                                            aaaTGCCCCAAACAAGTTTGGCCAGGCGGAAGGGCAGGGGGCCTGGGCCAGTGTTTGTTGGAATAAAGAAGACGAAGGCGAGTGAGAGGAAACCCCGGGCGGAAATCGCCAGATAGAAAGTTTAAGGGGGATGTTTCGTTGCTTTAA‑AGA‑CA‑CCCGCCTCTTCGGCCACTTTGGCCATTGCAG‑CAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTttat                                                                                                                                             <  1:138002/237‑1 (MQ=255)

                                                                                                                gCCATCCCGAAGGGCTTTGTGCCTGGGCGGGTGGTTGTTGAGAAGAAGAAGACGATTGCGATGGAAATAAAAACCCCGGCGATAACCGCAAGAAAGAAATTTTTAAGCGGATTTTTCGTTGCTTTAA‑AGA‑CG‑CCCGCCTCTTCGGCCACTTTGGCCATTGCAG‑CAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTaa                                                                                                                       <  2:386759/239‑1 (MQ=255)

                                                                                                                      ccGAGGGGTTTTGTGCCTGTGCCGGTGGTTGCTGTGATATAGAAGACGAATGCGTTGGAGATAAAAACCCCGGCGGTAATCGCCAGATAGAAAGTTTAAAGCGGTTTTTTCTTTGTTTTAA‑AGA‑CC‑CCCGCCTCTTCGGCCACTTTGGCCATTGCAG‑CAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAg                                                                                                                  <  1:995992/238‑1 (MQ=255)

                                                                                                                       cGAAGGGCATTGTGCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACCCCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGGTTGTTTCGTTGTTTTAT‑AGA‑CA‑CCCGCCTCTTCGGCCACTTTGGCCATTGCAG‑CAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCa                                                                                                                <  1:449443/239‑1 (MQ=255)

                                                                                                                                ttGTGCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTAT‑AGA‑CA‑CCCGCCTCTTCGGCCACTTTGGCCATTGCAG‑CAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATg                                                                                                       >  1:580042/1‑239 (MQ=255)

                                                                                                                                ttGTGCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTAT‑AGA‑CA‑CCCGCCTCTTCGGCCACTTTGGCCATTGCAG‑CAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATg                                                                                                       >  1:1758788/1‑239 (MQ=255)

                                                                                                                                 tgtgccaggcccGGTGGTGGCTGTATTAAAGAAGAGGAAGGGGATTGGAATGAAACCCCCGGCGGAAATCGCCAGATAGAAATTCTTAAGCGGATGTTTCGTTGTTTTAT‑AGA‑CA‑CCCGCCTCTTCGGCCACTTTGGCCATTGCAG‑CAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATga                                                                                                      <  2:1758788/239‑1 (MQ=255)

                                                                                                                                  gtgcctgtgccAGTGGTTGCTGTAATATAGAAGACGAATGCGATTGAGATAAAACCCCCGGCGGTAATCGCCAGAAAGAAAGTCTTAAGCGGTTGTTTCGTTGCTTTAA‑AGA‑CC‑CCCGCCTCTTCGGCCACTTTGGCCATTGCAG‑CAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATga                                                                                                      <  1:1067823/238‑1 (MQ=255)

                                                                                                                                      ctgtgccGGTGTTTGCTGTAATATAGAAGACAAATGCGATTGAGATAAAAACCCCGGCGGTAACCGCCAGAAAGAAAGTTTAAAGGGGTTGTTTCGTTGTTTTAA‑AGC‑CC‑CCCGCCTCTTCGGCCACTTTGGCCATTGCAG‑CAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGaaa                                                                                                  <  1:330731/238‑1 (MQ=255)

                                                                                                                                      ctgtgccAGTGTTTGCTGTGATATAGAAGACGAATGCGATTGAAATGAAACCCCCGGCGGTAAGCGCCAGAGAGAAAGTCTTAAGCGGTTGTTTCGTTGTTTTAT‑AGA‑CA‑CCCGCCTCTTCGGCCACTTTGGCCATTGCAG‑CAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAAt                                                                                                 <  1:947156/239‑1 (MQ=255)

                                                                                                                                      ctgggccAGTGGTTGTTGTGATATAGAAGACGAATGGGATTGAGATAAAAACCCCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGTTTTAA‑AGA‑CA‑CCCGCCTCTTCGGCCACTTTGGCCATTGCAG‑CAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGaaa                                                                                                  <  1:1115787/238‑1 (MQ=255)

                                                                                                                                        gtgccAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTAT‑AGA‑CA‑CCCGCCTCTTCGGCCACTTTGGCCATTGCAG‑CAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTg                                                                                               >  1:537149/1‑239 (MQ=255)

                                                                                                                                         tgccAGTGGTTGCTGTAATAAAGAAGACGAATGGGATTGAGATGAAAACCCCGGCGGTAATCGCCAGATAGAAAGTCTAAAGGGGTTGTTTCGTTGTTTTAT‑AGC‑CC‑CCCGCCTCTTCGGCCACTTTGGCCATTGCAG‑CAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTcaca                                                                                                                                                             <  1:1369613/176‑1 (MQ=255)

                                                                                                                                            cgggggTTGCTGTGATAAAAAAGACGAATGCGATTGAGATGAAAACCCCGGCGGTAATCGCAAGATAGAAAGTCTTAAGCGTATGTTTCGTTGCTTTAT‑AGA‑CC‑CCCGCCTCTTCGGCCACTTTGGCCATTGCAG‑CAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGata                                                                                            <  1:10805/234‑1 (MQ=255)

                                                                                                                                             aGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGGGTGTTTCGTTGCTTTAT‑AGA‑CA‑CCCCCCTCTTCGGCCACTTTGGCCATTGCAG‑CAGGAAGGAAAAGATCAAAAGGGGTGTCAGCTTTCACACTAACTCTCTCCTTATTAAGTTGGCGACGAGATACTAACAAAGAAATATAGATGAGAAAATGata                                                                                            >  2:1067823/1‑237 (MQ=255)

                                                                                                                                                   tGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGAAAGGCTTAAGCGGGTGGTTTGTTGCTTTAT‑AGA‑CA‑CCCGCCTCTTCGGCCACTTTGGCCATTGGAG‑CAGGAAATAAAAAATCAAAAGGGGTGTCAGCTTTCCACCTAActctct                                                                                                                                                   >  2:1369613/1‑176 (MQ=255)

                                                                                                                                                   tGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATAGCCAGATAGAAAGGCTTAAGCGGAAGTTTCGGTGCTTTAT‑AGA‑CA‑CCCGCCTCTTCGGCCACTTTGGCCATTTCAG‑CAGGAAGTAAAAGATCCAAAGGGGTGTCAGCCTTCCCACTAACTCTCTCTTTATTAAGTCGGCGAAGAGATACTATCAAAGCATTATAAATGAGAAAATGATATAAAACat                                                                                    >  2:10805/1‑239 (MQ=255)

                                                                                                                                                      tgtgATATAGAAGACGAATGCGATTGAGATGAAAACCCCGGCGGTAACCGCCAGATAGAAAGTCTTAAGCGGTTTTTTCTTTGTTTTAT‑AGA‑CC‑CCCGCCTCTTCGGCCACTTTGGCCATTGCAG‑CAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCatat                                                                                  <  1:442336/238‑1 (MQ=255)

                                                                                                                                                        tgATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTAT‑AGA‑CA‑CCCGCCTCTTCGGCCCCTTTGGCCATTGCAG‑CAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATtata                                                                                                           >  1:1107141/1‑211 (MQ=255)

                                                                                                                                                         gATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGAAAGTCTTTAGCGGGTGTTTCGTTGCTTTAT‑AGA‑CA‑CCCGCCTCTTCGGCCACTTTGGCCATTGCAG‑CAGGAAGTAAAAGATCAAAAGGGGTGTCAGCTTTCACACTAACTCTCTCTTTTTTAAATCGGCGACGAGAAACTAACAAAGAATTATAGATTATAAAATGATATAGATCATAtcta                                                                               >  2:98764/1‑237 (MQ=255)

                                                                                                                                                           tataGAAGACGAATGGGATTGAAATAAAACCCCCGGCGGTAATCGCCAGAAAGAAATTCTAAGGCGGTTGTTTCGTTGTTTAAA‑AGA‑CC‑CCCGCCTCTTCGGCCACTTTGGCCATTGCAG‑CAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGc                                                                             <  1:1690818/238‑1 (MQ=255)

                                                                                                                                                           taaaaaagaagAATGCGATGGAGATAAAAACCCCGGCGGAAACCCCAAAATAAAAAGTTTAAAGGGGATTTTTCTTTGTTTTAA‑AGA‑CA‑CCCGCCTCTTCGGCCACTTTGGCCATTGCAG‑CAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGcc                                                                            <  1:758934/234‑1 (MQ=255)

                                                                                                                                                           aataGAAGACGAATGCGATTGAGATAAAAACCCCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGGTTGTTTCGTTGTTTTAT‑AGA‑CA‑CCCGCCTCTTCGGCCACTTTGGCCATTGCAG‑CAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGcc                                                                            <  1:1140798/238‑1 (MQ=255)

                                                                                                                                                             aaaaagaCAAATGCGTTTGAGATGAAACCCCCGGCGGTAATCGCCAGAAAGAAAGTCTTAAGGGGTTTTTTGTTTGTTTTAA‑AGA‑CC‑CCCGCCTCTTCGGCCACTTTGGCCATTGCAG‑CAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTg                                                                          <  1:874086/236‑1 (MQ=255)

                                                                                                                                                               gaagaCGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTAT‑AGA‑CA‑CCCGCCTCTTCGGCCACTTTGGCCATTGCAG‑CAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCa                                                                                     >  1:681106/1‑226 (MQ=255)

                                                                                                                                                                aagaCGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTAT‑AGA‑CA‑CCCGCCTCTTCGGCCACTTTGGCCATTGCAG‑CAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCt                                                                       >  1:162622/1‑239 (MQ=255)

                                                                                                                                                                  gaCGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGGGTGTTTCGTTGCTTTAT‑AGA‑CA‑CCCGCCTCTTCGGCCACTTTTGCCATTTCAG‑CAGGGAGTAAAAGATCAAAAGGGTTGGCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACCAACAAAGCCTTATAGATTGGAAATTGATATAGAACATATCTTGCCTGGGttt                                                                     >  2:685473/1‑238 (MQ=255)

                                                                                                                                                                    cGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGACAGAAAGTCTTAACCCGATGCTTCGTTGCCTTAT‑AGA‑CA‑CACGCCTCTTCCGCCACTTTGGTCATTGCAG‑CAAGAAGCAAAAGATCCAAAGGGTTGCCAGCTTTCACCATAAATTTCTCTTTATTAACTAGGCGACGAGATACTAAAAAAACCTTATAGATGAGAAATTGAAATATATCATATCTCCCCTCGCttcaa                                                                   >  1:1277073/1‑236 (MQ=255)

                                                                                                                                                                     gAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGAAAGACTTAAGCGGATGGTTTGTTGCTTTAT‑AGA‑CA‑CCCGCCTCTTCGGCCACTTTGGCCCTTGCAG‑CAGGAAGTAAAAGATCCAAAAGGTTGTCCGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAg                                                                  >  1:1602427/1‑239 (MQ=255)

                                                                                                                                                                       aTGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTAT‑AGA‑CA‑CCCGCCTCTTCGGCCACTTTGGCCATTGCAG‑CAGGAAGTAAAAGa                                                                                                                                                                                     >  1:1139941/1‑122 (MQ=255)

                                                                                                                                                                       aTGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGAAAGGCTTAAGCGGATGTTTCGTTGCTTTAT‑AGA‑CA‑CCCGCCTCTTCGGCCACTTTGGCCATTGCAG‑CAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAACATtata                                                                                                           >  1:749805/1‑196 (MQ=255)

                                                                                                                                                                            aTTGAGATGAAAACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTAT‑AGA‑CA‑CCCGCCTCTTCGGCCACTTTGGCCATTGCAG‑CAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAg                                                                                                                                         >  1:58393/1‑161 (MQ=255)

                                                                                                                                                                                       aaCACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTAT‑AGA‑CA‑CCCGCCTCTTCGGCCACTTTGGCCATTGCAG‑CAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGata                                                                                            >  1:440919/1‑195 (MQ=255)

                                                                                                                                                                                       aaCACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTAT‑AGA‑CA‑CCCGCCTCTTCGGCCACTTTGGCCATTGCAG‑CAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGt                                                >  1:1807338/1‑239 (MQ=255)

                                                                                                                                                                                              ggggTAAACGCCAAGTAAAAAATTTATAGGGGATGGTTTGTTTGTTTAA‑AAA‑CACCCCGCCTCTTCGGCCACTTTGGCCATTGCAG‑CAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTa                                           <  2:748896/236‑1 (MQ=255)

                                                                                                                                                                                               cggTAATCGCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTAT‑AGA‑CA‑CCCGCCTCTTCGGCCACTTTGGCCATTGCAG‑CAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATg                                         >  1:1626868/1‑238 (MQ=255)

                                                                                                                                                                                                    accGCAAGATAAAAATTTTAAGGCGGTTTTTTCTTTGTTTTAT‑AGA‑CA‑CCCGCCTCTTCGGCCACTTTGGCCATTGCAG‑CAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGCAAATa                                   <  1:685473/237‑1 (MQ=255)

                                                                                                                                                                                                      cGCAAGATAAAAAGTCTTAAGGGGTTGTTTTGTTGTTTTTA‑AGC‑CC‑CCCGCCTCTTCGGCCACTTTGGCCATTGCAG‑CAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGCAAATACg                                 <  1:217126/239‑1 (MQ=255)

                                                                                                                                                                                                          ggATAGAATTTATAAAGGGGATGTTTTTTTGCTTTAT‑AGA‑CC‑CCCGCCTCTTCGGCCACTTTGGCCATTGCAG‑CAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGCAAATACGGAGt                             <  2:537149/238‑1 (MQ=255)

                                                                                                                                                                                                          aaaTAGAAAGTTTTAAGGGGTTTTTTTTTTGTTTTAT‑AGA‑CC‑CCCGCCTCTTCGGCCACTTTGGCCATTGCAG‑CAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGCAAATACGGAGt                             <  1:1487278/237‑1 (MQ=255)

                                                                                                                                                                                                             gaaaaaGTATTAGGTGGTTGTTCTTTTGTTTTAT‑AGA‑CA‑CCCGCCTCTTCGGCCACTTTGGCCATTGCAG‑CAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACGAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGCAAATACGGAGt                             <  1:1382791/233‑1 (MQ=255)

                                                                                                                                                                                                                    tCTAAAGGGGTTGTTTCTTTGTTTAAT‑AGA‑CC‑ACCGCCTCTTCGGCCACTTTGGCCATTGCAG‑CAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGCAAATACGGAGTAAATATTTGa                   <  1:1009148/239‑1 (MQ=255)

                                                                                                                                                                                                                              aTGCTTCGTTGCTTTAT‑AGA‑CA‑CCCGCCTCTTCGGCCACTTTGGCCATTGCAG‑CAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGCAAATACGGAGTAAATATTTGATTATCCaaat         >  1:421571/1‑239 (MQ=255)

                                                                                                                                                                                                                                     gTTGCTTTAT‑AGA‑CA‑CCCGCCTCTTCGGCCACTTTGGCCATTGCAG‑CAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGCAAATACGGAGTAAATATTTGATTATCCaaataaaaata  >  2:214331/1‑239 (MQ=255)

                                                                                                                                                                                                                                        gCTTTAT‑AGA‑CA‑CCCGCCTCTTCGGCCACTTTGGCCATTGCAG‑CAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGc                                                                >  1:660788/1‑174 (MQ=255)

GCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAGATGAAAACACCGGCGGTAATCGCCAGATAGAAAGTCTTAAGCGGATGTTTCGTTGCTTTAT‑AGA‑CA‑CCCGCCTCTTCGGCCACTTTGGCCATTGCAG‑CAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGCAAATACGGAGTAAATATTTGATTATCCAAATAAAAATA  >  NC_007779/954577‑955044

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