BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,093,153	G→C	14.2%	S95R (AGC→AGG)	rhaM ←	L‑rhamnose mutarotase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,093,153	0	G	C	14.2%	157.7 / 2.1	69	S95R (AGC→AGG)	rhaM	L‑rhamnose mutarotase
Reads supporting (aligned to +/- strand): ref base G (29/29); new base C (9/2); total (38/31)
Fisher's exact test for biased strand distribution p-value = 9.55e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.16e-01

AGTGAAAGGTCAGATTTGCGGAGTAATGCACATAATGGTTATTTAAATAAACCACATGAATCATTAATGGGTTATTCATTATTTTTGTGATTTACTCACAATGTATCA‑GGGAAATATAACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGC‑AGC‑TCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATA‑ACATCGGTCAT‑ATA‑TTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGC‑ATTCCAGCGTTCTTCAGATTCAATCTCTACCATGGCAAACAG‑CAGATTACGCGCTTTGTCGAGATAGATGGCGTAGTTATGCGCACCGTGA‑GATTTCAGCACTGCTTCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTC  >  NC_000913/4092921‑4093389

agTGAAAGGTCAGATTTGCGGAGTAATGCACATAATGGTTATTTAAATAAACCACATGAATCATTAATGGGTTATTCATTATTTTTGTGATTTACTCACAATGTATCA‑GGGAAATATAACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGGAAAAAAATCCTTTC‑GGC‑CCGCGACta                                                                                                                                                                                                                                              >  2:633356/1‑238 (MQ=255)

         tCAGATTTGCTGAGTAATGCACTTAATGGTTAGTTAAATAAACCACATGAATCATTAATGGGTTATTCATTATTTTTGTGATTTACTCACAATGTATCA‑GGGAAATATAACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGC‑AGC‑TCGCTACTCACCGGGc                                                                                                                                                                                                                                       <  2:51016/237‑1 (MQ=255)

          cAGATTTGCGGAGTAATGCACATAATGGTTATTTAAATAAACCACATGAATCATTAATGGGTTATTCATTATTTTTGTGATTTACTCACAATGTATCA‑GGGAAATATAACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGGAAAAAACTTCTTTC‑AAC‑TCGCTACTCCCCCGGgtgt                                                                                                                                                                                                                                    >  2:1047193/1‑239 (MQ=255)

           aGATTTGCGGAGTAATGCACATAATGGTTATTTAAATAAACCACATGAATCATTAATGGGTTATTCATTATTTTTGTGATTTACTCACAATGTATCA‑GGGAAATATAACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAAAAACTCTTGC‑AGC‑TCGCCACTCCCCCGGCTggt                                                                                                                                                                                                                                   >  2:749247/1‑237 (MQ=255)

              tttgggggatAATGGACATAATGGTTATTTAAATAAACCACATGAATCATTAAGTGGTTATTTATTATTTTGTTGATTTACTCACAATGTATCA‑GGGAAATATAACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGC‑AGC‑TCGCTACTCACCGGGCTGTTATc                                                                                                                                                                                                                                <  1:186277/230‑1 (MQ=255)

                                 aaTGGTTATTTAAATAAACCACATTAATCATTAATGGGTTATTCATTATTTTTGTGATTTACTCACAATGTATCA‑GGGAAATATAACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGC‑AGC‑TCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATA‑ACa                                                                                                                                                                                                            <  2:1232572/239‑1 (MQ=255)

                                 aaTGGTTATTTAAATAAACCACATGAATCATTAATGGGTTATTCATTATTTTTGTGATTTACTCACAATGTATCA‑GCGAAATATAACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAGATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTTCGTCCGGGTAAACACGTATTCC‑AGC‑GCGCGACTACCCCGGCTGTTATCCCGGGTCCCGGGCATA‑Ac                                                                                                                                                                                                             >  2:776966/1‑238 (MQ=255)

                                  tttgTTTTTTAATTTACCCACATTAATCATTAATTGTTTTTTCATTATTTTTGTGATTTACTCACAATGTATCA‑GGGAAATATAACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGC‑AGC‑TCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATA‑ACAt                                                                                                                                                                                                           <  1:1252135/236‑1 (MQ=255)

                                  aTGGTTATTTAAATAAACCACATGAATCATTAATGGGTTATTCATTAGTTTTGTGATTTACCCACAATGTATCA‑GGGAAATATAACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGC‑AGC‑TCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATA‑ACAt                                                                                                                                                                                                           <  2:574771/239‑1 (MQ=255)

                                    ggTTATTTAAATAAACCACATGAATCATTAATGGGTTATTCATTATTTTTGTGATTTACTCACAATGTATCA‑GGGAAATATAACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGGAAAACACTTCTTTC‑CAC‑TCCCTACTCACCCGGCTGTTATTCCGGGTAGCGGGGATa‑aaaacg                                                                                                                                                                                                         >  2:1183169/1‑234 (MQ=255)

                                        aTTTAAATAAACCACATGAATCATTAATGGGTTATTCATTATTTTTGTGATTTACTCACAATGTATCA‑GGGAAATATAACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAAACATTCTTGC‑AGC‑TCGCTACTCAACGGGCCGGTCTCCCGGGTCACGTGCATA‑AAATcaggcc                                                                                                                                                                                                     >  2:1144297/1‑234 (MQ=255)

                                         tttAAATTATCCATTGGAATTATTAGTGTGTTTTTCATTCTTTTTGTGATTTACTAATAATGTAACA‑GGGAAATAGAACTTACCGTCTGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCTTCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACGTCTTTC‑AGG‑TCGCTACTCACCGGTCTGTTATCCGAGTTATCGGGCATA‑ACATCGGTCat                                                                                                                                                                                                    <  1:40273/239‑1 (MQ=255)

                                         tttAAATAAACCCAAGTAATCATTAATGGGTTATTTATTATTTTTTGGATTTACTTACAAAGTTTCA‑TGGAAATATAACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCCCATTGGGCAATTACGGCAGGTAAAACACTTCTTGC‑AGC‑TCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATA‑ACATCGGTCat                                                                                                                                                                                                    <  1:1047193/239‑1 (MQ=255)

                                         tttAAATAAACCACATGAATCATTAATGGGTTATTCATTATTTTTGTGATTTACTCACAATGTATCA‑GGGAAATATAACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTTAAAAACTTCTTTC‑AGC‑CCGCCACTCACCCGGCTGGTTTCCGGGTTAGCGGGCCTa‑aaaatgggcca                                                                                                                                                                                                    >  2:734921/1‑229 (MQ=255)

                                         tttAAATAAACCACATGAATCATTAATGGGTTATTCATTATTTTTGTGATTTACTCACAATGTATCA‑GGGAAATATAACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTAC‑AGC‑TCCCTACTCAACCGGCTGTTTTCCGGGGTAGCGGGCCTC‑AAATAGGTcct                                                                                                                                                                                                    >  2:76247/1‑237 (MQ=255)

                                         tttAAATAAACCACATGAATCATTAATGGGTTATTCATTATTTTTGTGATTTACTCACAATGTATCA‑GGGAAATATAACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAAAACTTCTTGG‑AGC‑TCGCTACTCACCGGGCTGTTTTCCCGGGTAAGGGGGATA‑AAATCGGGCat                                                                                                                                                                                                    >  2:854905/1‑239 (MQ=255)

                                           taaataaaCCACATGAATCATTAATGGGTTATTCATTATTTTTGTGATTTACTCACAATGTATCA‑GGGAAATATAACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGACAGGTAAAAAACATCTTAC‑AGC‑ACGTTACTACCCGGCTTGTTATCCAGGTAAggaggggaaaaaacgggacaaa                                                                                                                                                                                                   >  1:574771/1‑218 (MQ=255)

                                                       aTGAATCATTAATGGGTTATTCATTATTTTTGTGATTTACTCACAATGTATCA‑GGGAAATATAACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGGAAAAAACCTCTTGC‑AGC‑CCGCCACTTACCCGGCTGTTTTCCCGGTTTCGGGGCATA‑AAATCCGGCCA‑ATA‑TTTCCCCCAAc                                                                                                                                                                                    >  2:618752/1‑239 (MQ=255)

                                                              aTTAATGGGTTATTCATTATTTTTGTGATTTACTCACAATGTATCA‑GGGAAATATAACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGC‑AGC‑TCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATA‑ACATCGGTCAT‑ATA‑TTTCCACCAACGTTGGCa                                                                                                                                                                             <  2:125437/239‑1 (MQ=255)

                                                              aTTAATGGGTTATTCATTATTTTTGTGATTTACTCACAATGTATCA‑GGGAAATATAACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGGAAAAAACTTCTTGC‑AAA‑TCGCTACCCCCCCGGCTGTTTTCCGGGTTTACCGGCAAA‑AACTCGGGCAA‑TAA‑TTTTCACCCACCCTGGCa                                                                                                                                                                             >  2:486171/1‑239 (MQ=255)

                                                                     ggTTATTCATTATTTTTGTGATTTACTCACAATGTATCA‑GGGAAATATAACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGTCAGGGAAAACACTTTTTTC‑CCG‑TCCCCATCTACCCGGCTGTTATCCTGGTTAGCGCGCATACACA‑CCGTCCT‑ATT‑TTTTCACCACCGTTGGCaaaactca                                                                                                                                                                      >  2:651974/1‑234 (MQ=255)

                                                                                 tttttGTGATTTACTCACAATGTATCA‑GGGAAATATAACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAAAAATTCTTGC‑AGA‑TAGCTAATCACCGGGCTGGTTTACGGGTTAACGGGCATG‑AAATACGGCAT‑AAA‑TTTTCACCCACCTTgggacaaaaacagttctgcaacg                                                                                                                                                          >  1:23192/1‑218 (MQ=255)

                                                                                         aTTTACTCATAATGT‑TCATGGGAATAATAAATTACCGGGAGATGTAAATTGTTTAATTTTTCAAATTGAAATTTAAAATTTTGTGCCGGAGGCATCTCTGGCACATTGGGCACTGTCGGCAGGTAAAACACTTCTTGC‑AGC‑TCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATA‑ACATCGGTCAT‑ATA‑TTTCCGACAACGTTGGCAAACATCAGTGCTGGCTACCGC‑ATTCCa                                                                                                                                                 <  1:76247/239‑1 (MQ=255)

                                                                                            tataaaagacagaTCA‑GGGAAATATAACTTACCGGGAGATGCAATTGGTTTAAATTTTCAAATCGAAATTTAAAATCTTGTGCCGGAGTCATCTCTGGCCAATTGGGCAATTACGGCAGGTAAAACTCTTCTTGC‑AGC‑TCGCTACTCACCGGGCTGGTATCCGGGTTAGCGGGCATA‑ACATCGGTCAT‑ATA‑TTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGC‑ATTCCAGCg                                                                                                                                              <  1:820152/227‑1 (MQ=255)

                                                                                                caaaaTGTATCA‑GGGAAATATAACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGC‑AGC‑TCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATA‑ACATCGGTCAT‑ATA‑TTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGC‑ATTCCAGCGttct                                                                                                                                          <  2:23192/236‑1 (MQ=255)

                                                                                                 acaATGTATCA‑GGGAAATATAACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGC‑AAC‑TCGCTACTCACCGGGCTGTTATCCGGGCTAGCGGGCATT‑ACATCGGTCATAATA‑TTT‑CACCAACGTTGGCAAAAATCAGGTATGGCACCCGC‑ATTCCACCGttct                                                                                                                                          >  2:683183/1‑238 (MQ=255)

                                                                                                   aaTGTATCA‑GGGAAATATAACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAACAAACTCTTTC‑AGT‑GTCCCAGTAACCGGGCTATTCAACGAGGTAGCGGGAATA‑ACATCAGGCAT‑AAATTTTACCACAACG‑TGGCCAACACAAATACAGGCAACCGA‑ACTCAAGCGttcgtac                                                                                                                                       >  1:1006009/1‑235 (MQ=255)

                                                                                                      gTATCA‑GGGAAATATAACTTACCGGGAAATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGC‑AGC‑TCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATA‑ACATCGGTCAT‑ATA‑TTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGC‑ATTCCAGCGTTCTTCAGAt                                                                                                                                    <  2:547735/239‑1 (MQ=255)

                                                                                                         tCA‑GGGAAATATAACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTGAAAAACTTCTTGC‑CGC‑TCCCTACTCTCCTGGCTGTTTTCCCGGTTTGCCGGCCTG‑ACATCGGTCCT‑ATA‑TTTTCACCAACGTTGGCAACCATCAACGCTGGCAAACGC‑ATTCCAGCGTTTTTTCGATTCa                                                                                                                                 >  2:78555/1‑239 (MQ=255)

                                                                                                             gggAAATATAACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAAAACTTCTTGC‑AGCACCGCTA‑TCACCGGGCTTTTATTCGGGTTTGAGCGGAGA‑ACATCGGTCAT‑ATA‑TTTCCCCCCACGTTGGGCAACATAAAGATTGGTAGCCGC‑AATCCAGCggtatttagtatagaaa                                                                                                                              >  1:271894/1‑223 (MQ=255)

                                                                                                              ggAAATATAACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTTC‑AGC‑TCCCTACTCACCCGGCTGTTATCCGGGTTAGCGGGCCTA‑ACATCGGGCAC‑ATA‑TTTCCACCCACGTTGGCCAACATCAGTGCTGGCAACCGC‑CTTCCAGCGttctttcgagttcatc                                                                                                                              >  2:728933/1‑227 (MQ=255)

                                                                                                                    ataACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGC‑AGC‑TCGCTACTCACCGGGCTGTTATCCGGGTTTGCGGGCATA‑ACATCGGTCAT‑TTA‑TTTCCACCAACGTTGGCAAACATCCATGCTGGGCACCCC‑CTTCCCACGTTCTTCACATTCTGTctctttt                                                                                                                        >  2:820152/1‑235 (MQ=255)

                                                                                                                        cTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGGTAAACACTTCCTGC‑AGC‑TCCCTACTTCCCGCGCTGTTATCCGGGTTAGAGGGCATA‑AAAACGGGCAT‑ATA‑TTTTCACCCACGTTTGGAAACATCAGTACTGGCCACCGCAATTCCA‑CGGTCTTCAGTATCAAAATCTACCAAggg                                                                                                                   >  1:621168/1‑238 (MQ=255)

                                                                                                                               ggAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGC‑AGC‑TCGCTACTCACCGGGCTGTTATTCGGGGTAGCGGGCCTA‑AAATCGGTCCT‑ATA‑TTTCCACCAACGTTTGCAAACATCAGTGCTGGCAAACGC‑ATTCCAGCCTTCTTCAGATTCAATCTCTACCATTGGAAAcaac                                                                                                            >  2:717785/1‑237 (MQ=255)

                                                                                                                                      tAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGC‑AGC‑TCGCTACTCCCCCGGCTGTTATCCGGGTTAGCGGGCATA‑ACATCGGTCAT‑ATA‑TTTCCACCAACGTTTGCAAACATCCGTGCTGGCAAccc                                                                                                                                                         >  2:171069/1‑187 (MQ=255)

                                                                                                                                           tGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGC‑AGC‑TCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATA‑ACATCGGTCAT‑ATA‑TTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGC‑ATTCCAGCGTTCTTCAGATTCAATCTCTACCATGGCAAACAG‑CAGATTACGCGCt                                                                                               <  2:923418/239‑1 (MQ=255)

                                                                                                                                             tttAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGC‑AGC‑TCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATA‑ACAt                                                                                                                                                                                                           >  2:955625/1‑133 (MQ=255)

                                                                                                                                             tttAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGC‑AGC‑TCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATA‑ACAt                                                                                                                                                                                                           <  1:955625/133‑1 (MQ=255)

                                                                                                                                             tttAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACATCTTTC‑AGA‑CCCCTACTCACCAGGCTGTTTTGCGGGTTAGCGGGCTTA‑ACATCGGTGAT‑ATC‑TTTCCACCACCGTTGGCGAACATCAGTGCTGGCAACCGC‑ATTCCAACGGTCTTCAAATTTAAACTCAACCCCGGCaaaaacacaatatcccgcttt                                                                                              >  1:441147/1‑221 (MQ=255)

                                                                                                                                              ttAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGC‑AGC‑TCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATA‑ACATCGGTCAT‑ATA‑TTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGC‑ATTCCAGCGTTCTTCAGATTCAATCTCTACCATGGCAAACAG‑CAGATTACGCGCTTTg                                                                                            <  2:291074/239‑1 (MQ=255)

                                                                                                                                              ttAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGC‑AGC‑TCGCTACTCACCCGGCTGTTATCCGGGGTAGCGGGCATA‑ACATCGGTCAT‑ATA‑TTTCCACCAACGTTTGCAAACATCAGTGCTGGCAAACGC‑ATTCCAGCGTTCTTCCGATTCCATTTCTACCATGGCAAACAG‑CAGATTACACGCTTTg                                                                                            >  2:307659/1‑239 (MQ=255)

                                                                                                                                              ttAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGGAAAAAAACTCTTTCAACC‑TCCTCATTCCCCGGCATGTT‑TCCCGGGTAGCGGGGCTA‑ACATAGGGCAT‑ATA‑TTTTCAACAACCTTGGGAACCATCCGGGCTGGAAAACGC‑ATTCCAGCGTTCTTTAGATTAAAATCCTACCAAGGCCAAAAGACACATAACGCGtttg                                                                                             >  1:923418/1‑238 (MQ=255)

                                                                                                                                                aaTTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGC‑AGC‑TCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATA‑ACATCGGTCAT‑ATA‑TTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGC‑ATTCCAGCGTTCTTCAGATTCAATCTCTACCATGGCAAACAG‑CAGATTACGCGCTTTGTc                                                                                          <  2:441147/239‑1 (MQ=255)

                                                                                                                                                      gcAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACCCTTCTTGC‑AGC‑TCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATA‑ACATCGGTCAT‑ATA‑TTTCCACCAACGTTGGAAAACATCAGTGCTGGCAACCGC‑ATTCCAGCGTTCTTCAGATTCAATCTCTACCATGGCAAACAG‑CAGATTACGCGCTTTGTCGagata                                                                                    <  2:45994/238‑1 (MQ=255)

                                                                                                                                                                           tgtgCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGC‑AGC‑TCGCTACTCACCGGGGTGTTATCCGGGTTAGCGGGCATA‑ACATCGGTCAT‑ATA‑TTTCCACCAACGTTGGCAAACATCAGTGCTGGAAACCGC‑ATTCCAGCGTTCTTCAGATTCAATTTCTACCCTGGCAAACAC‑CAGATTACGCGCTTTGTCGGGTTAACTGGGGTACTTATTCGCAcc                                                               >  2:1231344/1‑239 (MQ=255)

                                                                                                                                                                           tgtgCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGGAAAACACTTCTTGC‑AGC‑TCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATA‑ACATCGGGCAT‑ATA‑TTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGC‑ATTCCAGCGGTCTTCAGATTCAATCTCTACCCTGGCAAACAG‑CAGATTACGCGCTTTGTCGAGATAGATGGCGTAGTTATGCGCAcc                                                               >  2:245102/1‑239 (MQ=255)

                                                                                                                                                                           tgtgCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGGAAAACACTTCTTGC‑AGC‑TCGCTACTCACCGGGCTGTTATCCGGGGTAACGGGCATA‑ACATCGGTCCT‑ATA‑TTTCCACCCACGTTGGCAAACATCAGTGCTGGCAACCGC‑ATTCCAGCGTTCTTCAGATTCAATCTCTACCATGGCAAACAG‑CCGATTACGCGCTTTGTTGGGGTAGATTGCGTGGTTATTCGCAcc                                                               >  2:591624/1‑239 (MQ=255)

                                                                                                                                                                                         ttttGGCTCTTTGGGCAATTACGGGAGGTAAACCTCTTTTTGG‑AGC‑CCGCTACGCGCCGGGCTGTTATCCGGGTTAGCCGGCATA‑ACATCGGTCAT‑ATA‑TTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGC‑ATTCCAGCGTTCTTCAGATTCAATCTCTACCATGGCAAACAG‑CAGATTACGCGCTTTGTCGAGATAGATGGCGTAGTTATGCGCACCGTGA‑GATTTCAGCa                                                <  1:617448/236‑1 (MQ=255)

                                                                                                                                                                                         ctctTGGACCTTGGGGAATGACGGCAGGTAAAACCCTTCCTGC‑AGC‑TCCCTTCTCACCGGGCTGTTATCCGGGTTAGCGGGCTTA‑ACATCGGTCAT‑ATA‑TGTCCACCAACGTTGGCAACCATCAGTGCTGGCAACCGC‑ATTCCAGCGTTCTTCAGATTCAATCTCTACCATGGCAAACAG‑CAGATTACGCGCTTTGTCGAGATAGATGGCGTAGTTATGCGCACCGTGA‑GATTTCAGCa                                                <  1:854905/239‑1 (MQ=255)

                                                                                                                                                                                         ctctGGCACATTGGGCAATTACGGCAGGTAAAACCCTTCTTGC‑AGC‑TCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATA‑ACATCGGTCAT‑ATA‑TTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGC‑ATTCCAGCGTTCTTCAGATTCAATCTCTACCATGGCAAACAG‑CAGATTACGCGCTTTGTCGAGATAGATGGCGTAGTTATGCGCACCGTGA‑GATTTCAGCa                                                <  2:271894/239‑1 (MQ=255)

                                                                                                                                                                                           ctgtgcccttGGGCTATTCCGGCCGGTAAAAACCTTCTTTC‑AAC‑TCCATACCCTCCGGGCTGTTATCCGGGTTAGCGGGCATA‑ACATCGGGCAT‑ATA‑TTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGC‑ATTCCAGCGTTCTTCAGATTCAATCTCTACCATGGCAAACAG‑CAGATTACGCGCTTTGTCGAGATAGATGGCGTAGTTATGCGCACCGTGA‑GATTTCAGCACt                                              <  1:734921/231‑1 (MQ=255)

                                                                                                                                                                                                    tGGGCAATTACGGCAGGGAAAACACTTCTTGC‑AGC‑CCGCTACCCACCCGGCTGTTATCCGGGGTAGCGGGCATA‑ACATTGGTCAT‑ATA‑TTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGC‑ATTCCAGCGTTCTTCAGATTCAATCTCCACCATGGCAAACAG‑CAGATTTAGAGCTTTGTCGAAGTTATTGGCGGAGTTATCCGCACCGCGG‑GATTTCAGCCCTGCTTCCaag                                     >  2:313091/1‑237 (MQ=255)

                                                                                                                                                                                                           gataagagacagaaaCACTTCTTGC‑AAC‑TCGCTACTCACCGGGCTTTTATCCGGGTTAGCGGGCATA‑ACATCGGTCAT‑ATA‑TTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGC‑ATTCCAGCGTTCTTCAGATTCAATCTCTACCATGGCAAACAG‑CAGATTACGCGCTTTGTCGAGATAGATGGCGTAGTTATGCGCACCGTGA‑GATTTCAGCACTGCTTCCAGTTCTGGc                               <  2:1061561/226‑1 (MQ=255)

                                                                                                                                                                                                              cgggaagTAAAACCCTTCTTGC‑AGC‑TCGCTTCGCACCGGGCTGTTATCCGGGTTAGCGGGCATA‑ACATCGGTCAT‑ATA‑TTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGC‑ATTCCAGCGTTCTTCAGATTCAATCTCTACCATGGCAAACAG‑CAGATTACGCGCTTTGTCGAGATAGATGGCGTAGTTATGCGCACCGTGA‑GATTTCAGCACTGCTTCCAGTTCTGGCCAGa                           <  2:284205/233‑1 (MQ=255)

                                                                                                                                                                                                              cGGCAGGTAAAACACTTCTTGC‑AGC‑TCGGTACTCCCCGGGCTGTTAACCGGGTTAGCGGGCATA‑ACATCGGTCAT‑ATA‑TTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGC‑ATTCCAGCGTTCTTCAGATTCAATCTCTACCATGGCAAACAG‑CAGATTACGCGCTTTGTCGAGATAGATGGCGTAGTTATGCGCACCGTGA‑GATTTCAGCACTGCTTCCAGTTCTGGCCAGa                           <  2:593531/239‑1 (MQ=255)

                                                                                                                                                                                                              cGGCAGGGAAAACTCTTTCTGC‑AGC‑TCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATA‑ACATCGGTCAT‑ATA‑TTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGC‑ATTCCAGCGTTCTTCAGATTCAATCTCTACCATGGCAAACAG‑CAGATTACGCGCTTTGTCGAGATAGATGGCGTAGTTATGCGCACCGTGA‑GATTTCAGCACTGCTTCCAGTTCTGGCCAGa                           <  2:1006009/239‑1 (MQ=255)

                                                                                                                                                                                                                     tAAAACCCTTCTTGC‑AGC‑TCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATA‑ACATCGGTCAT‑ATA‑TTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGC‑ATTCCAGCGTTCTTCAGATTCAATCTCTACCATGGCAAACAG‑CAGATTACGCGCTTTGTCGAGATAGATGGCGTAGTTATGCGCACCGTGA‑GATTTCAGCACTGCTTCCAGTTCTGGCCAGATGGGa                      <  2:1231224/237‑1 (MQ=255)

                                                                                                                                                                                                                            cttcttGC‑AGC‑TCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATA‑ACATCGGTCAT‑ATA‑TTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGC‑ATTCCAGCGTTCTTCAGATTCAATCTCTACCATGGCAAACAG‑CAGATTACGCGCTTTGTCGAGATAGATGGCGTAGTTATGCGCACCGTGA‑GATTTCAGCACTGCTTCCAGTTCTGGCCAGATGGGATTATGCCGa             <  2:976261/239‑1 (MQ=255)

                                                                                                                                                                                                                              tcttGC‑AGC‑TCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATA‑ACATCGGTCAT‑ATA‑TTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGC‑ATTCCAGCGTTCTTCAGATTCAATCTCTACCATGGCAAACAG‑CAGATTACGCGCTTTGTCGAGATAGATGGCGTAGTTATGCGCACCGTGA‑GATTTCAGCACTGCTTCCAGTTCTGGCCAGATGGGATTATGCCGACg           <  2:770158/239‑1 (MQ=255)

                                                                                                                                                                                                                              tcttGC‑AGC‑TCGCTACTCACCGGGCTGTTATCCGGGATAGCGGGCATA‑ACATCGGTCAT‑ATA‑TTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACACC‑ATTCCAGCGTTCTTCAGATTCAAACTCTACCATGGCAAAAAG‑CAGCTTACGCGCTTTGACGACATAGATGGCCTCAGTACGCGC‑CCCTGACGATTTCCGCACCCTCTCCAGTTCTGACCAGACGGGAATATGcgaag            >  1:605278/1‑234 (MQ=255)

                                                                                                                                                                                                                               cttGC‑AGC‑TCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATA‑ACATCGGTCAT‑ATA‑TTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGC‑ATTCCAGCGTTCTTCAGATTCAATCTCTACCATGGCAAACAG‑CAGATTACGCGCTTTGTCGAGATAGATGGCGTAGTTATGCGCACCGTGA‑GATTTCAGCACTGCTTCCAGTTCTGGCCAGATGGGATTATGCCGAc            >  2:1174527/1‑237 (MQ=255)

                                                                                                                                                                                                                               cttGC‑AGC‑TCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATA‑ACATCGGTCAT‑ATA‑‑TTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGC‑ATTCCAGCCTTCTTCAGATTAAATCTCTACCATGGCAAACAG‑CAGATTACCCGCTTTGTCGAGATAGTTGGCTTAGTTCTGCGCACCGTGT‑GATTTCAGACCTGCTTCCCGTTCTGCCACGATggactttgtccacacat         >  1:10717/1‑224 (MQ=255)

                                                                                                                                                                                                                                  gC‑AGC‑TCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATA‑ACATCGGTCAT‑ATA‑TTTCCACCAACGTTGGCAAACCTCAGTGCTGGCAACCGC‑ATTCCAGCGTTCTTCAGATTCAATCTCTACCATGGCAAACAG‑CAGCTTACGCCCTTTTTCGAGATAGATGGCTTCGTACGTCGCCACGTCA‑GATTTCAGCACTGCTTCCATTACTGCCCAGATGGGaattgttctcttaata       >  1:1061561/1‑224 (MQ=255)

                                                                                                                                                                                                                                  gC‑AGC‑TCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATA‑ACATCGGTCAT‑ATA‑TTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGC‑ATTCCAGCGTTCTTCAGATTCAATCTCTACCATGGCAAACAG‑CAGATTACGCGCTTTGTCGAGATAGATGGCGTAGTTATGCGCACCGTGA‑GATTTCAGCACTGCTTCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGa       <  2:420225/239‑1 (MQ=255)

                                                                                                                                                                                                                                      gC‑TCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATA‑ACATCGGTCAT‑ATA‑TTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGC‑ATTCCAGCGTTCTTCAGATTCAATCTCTACCATGGCAAACAG‑CAGATTACGCGCTTTGTCGAGATAGATGGCGTAGTTATGCGCACCGTGA‑GATTTCAGCACTGCTTCCAGTTCTGGc                               >  2:1064027/1‑212 (MQ=255)

                                                                                                                                                                                                                                         tCGCTACTCCCCGGGCTGTTATCCGGGTTAGCGGGCATA‑ACATCGGTCAT‑ATA‑TTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGC‑ATTCCAGCGTTCTTCAGATTCAATCTCTACCATGGCAAACAG‑CAGATTACGCGCTTTGTCGAGATAGATGGCGTAGTTATGCGCACCGTGA‑GATTTCAGCACTGCTTCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATa     <  2:282568/236‑1 (MQ=255)

                                                                                                                                                                                                                                         tCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATA‑ACATCGGTCAT‑ATA‑TTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGC‑ATTCCAGCGTTCTTCAGATTCAATCTCTACCATGGCAAACAG‑CAGATTACGCGCTTTGTCGAGATAGATGGCGTAGTTATgcg                                                                   <  1:775413/175‑1 (MQ=255)

                                                                                                                                                                                                                                         tCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATA‑ACATCGGTCAT‑ATA‑TTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGC‑ATTCCAGCGTTCTTCAGATTCAATCTCTACCATGGCAAACAG‑CAGATTACGCGCTTTGTCGAGATAGATGGCGTAGTTATgcg                                                                   >  2:775413/1‑175 (MQ=255)

                                                                                                                                                                                                                                         tCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATA‑ACATCGGTCAT‑ATA‑TTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGC‑ATTCCAGCGTTCTTCAGATTCAATCTCTACCATGGCAAACAG‑CAGATTACGCGCTTTGTCGAGATACATGGCGTATTTATCCCCCCCTGGC‑GATTTCACCCATGCTTCCCGTTCTGCCCCAATGGGATTTTGCCGACCCTgcgacct  >  1:1126132/1‑233 (MQ=255)

AGTGAAAGGTCAGATTTGCGGAGTAATGCACATAATGGTTATTTAAATAAACCACATGAATCATTAATGGGTTATTCATTATTTTTGTGATTTACTCACAATGTATCA‑GGGAAATATAACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGC‑AGC‑TCGCTACTCACCGGGCTGTTATCCGGGTTAGCGGGCATA‑ACATCGGTCAT‑ATA‑TTTCCACCAACGTTGGCAAACATCAGTGCTGGCAACCGC‑ATTCCAGCGTTCTTCAGATTCAATCTCTACCATGGCAAACAG‑CAGATTACGCGCTTTGTCGAGATAGATGGCGTAGTTATGCGCACCGTGA‑GATTTCAGCACTGCTTCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTC  >  NC_000913/4092921‑4093389

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