BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,778,123	G→T	6.8%	intergenic (+341/+23)	ypjF → / ← ypjA	CP4‑57 prophage; toxin of the YpjF‑YfjZ toxin‑antitoxin system/adhesin‑like autotransporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,778,123	0	G	T	6.8%	786.1 / 30.5	229	intergenic (+341/+23)	ypjF/ypjA	CP4‑57 prophage; toxin of the YpjF‑YfjZ toxin‑antitoxin system/adhesin‑like autotransporter
Reads supporting (aligned to +/- strand): ref base G (103/109); new base T (6/10); total (110/119)
Fisher's exact test for biased strand distribution p-value = 4.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 229 total aligned reads displayed.

AAAAAACCTGACATTTTCTTTTAGGACCAACAATGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCAGCCCACCAATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGAATGAGGCTTTAATACCTGCCTGGTATATTCCACGTCTGCCCGACAAATAATTGACGAAATTACCGTCACTATTCACTTTCACCCGGTTATCATCGACAAATTCTTTGCGCACAGCCGCCTTCAGCCACGGCTC  >  NC_000913/2777868‑2778379

aaaaaaCCTGACATTTTCTTTTAGGACCAACAATGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCAGCCCACCAATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGATCCTGGTAACATGAAAGCCCACAAATGTGGCATTTTTCGTTgg                                                                                                                                                                                                                                               >  2:329188/1‑274 (MQ=255)

                                aTGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCAGCCCACCAATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACcg                                                                                                                                                                                                               <  2:90989/275‑1 (MQ=255)

                                                          aaGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCAGCCCACCAATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACCCCAGCCACCGGCTTCCCCGGGGGTTCCCAACCgga                                                                                                                                                                                      >  4:317522/1‑273 (MQ=255)

                                                                       atTTAACAACAAGTTACACACCCAATTCAGACTCCGCCAGCCCACCAATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACCCCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATacc                                                                                                                                                                        <  2:269833/275‑1 (MQ=255)

                                                                                                    gACTCCGCCAGCCCACCAATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTAACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTActg                                                                                                                                           >  4:33932/1‑275 (MQ=255)

                                                                                                    gACTCCGCCAGCCCACCAATCATGATTGGACGGTGTAAAGACAACAACAACAAAGACAGGAAGTTAGAAGTCGAAGCAAAGCACCGACCAGACGGTGGGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTActg                                                                                                                                           <  2:159386/275‑1 (MQ=255)

                                                                                                                ccaccaATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGAtgtg                                                                                                                                                                                                                                                             >  6:148175/1‑149 (MQ=255)

                                                                                                                 caccaATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTgg                                                                                                                                                                                                                                                            <  7:368733/149‑1 (MQ=255)

                                                                                                                  accaATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGgggg                                                                                                                                                                                                                                                           >  6:1392492/1‑149 (MQ=255)

                                                                                                                   ccaATCATGATTGGACGGTGTAGGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCTGCAAAACACCGACCAGACGGTGGGGAGACATAAAAGGATACGCAAAGGAACCGCGGCTCCTGGTAACATGAAAGCCCACAGATGCGGGc                                                                                                                                                                                                                                                          >  8:1228008/1‑149 (MQ=255)

                                                                                                                     aaTCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACATCTGTGGGCtt                                                                                                                                                                                                                                                        <  6:1248747/149‑1 (MQ=255)

                                                                                                                      aTCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCttt                                                                                                                                                                                                                                                       <  5:1461150/149‑1 (MQ=255)

                                                                                                                       tCAAGAGTGGACGGTGGAAGGACAAAACCAACAAAAACAGGAAGTTAGAAGGCTCAGAAAAACACCGAGCAGACGGGGAGGGGACATAAAAGGATAGTCAAAGGAGCCGCGGGTCCTGGTAACATGAAAGCCCACAGATGGGGGCGTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGAATGAGGCTTTAATAc                                                                                                                        <  4:160114/275‑1 (MQ=25)

                                                                                                                          tGATTGGATGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCg                                                                                                                                                                                                                                                   >  5:931808/1‑149 (MQ=255)

                                                                                                                           gATTGGAGGGTGTAAGGGAAACACCAACCAAAAGAGTAAGTTAGAGAACTAAACAAAACACCGGCCAGAGGGTGAGGAGACATAAAAGGTTACGAAAAGGAGCCGCGGCGCCTGGTAAAATGAACGCCCACCGATGTGGGCGTTTTCGTTGATGGTCAGACCGACTAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCAGCTTAACGTACTGCTGAATGAGGCTGTAATAcctgc                                                                                                                    <  4:91570/275‑1 (MQ=21)

                                                                                                                           gATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGt                                                                                                                                                                                                                                                  >  5:1035398/1‑149 (MQ=255)

                                                                                                                           gATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGt                                                                                                                                                                                                                                                  <  5:545020/149‑1 (MQ=255)

                                                                                                                           gATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGt                                                                                                                                                                                                                                                  <  7:570865/149‑1 (MQ=255)

                                                                                                                                gACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATg                                                                                                                                                                                                                                             >  7:196975/1‑149 (MQ=255)

                                                                                                                                         aGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACg                                                                                                                                                                                                                                    >  7:880294/1‑149 (MQ=255)

                                                                                                                                           gACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGAc                                                                                                                                                                                                                                  <  7:236665/149‑1 (MQ=255)

                                                                                                                                            acaacaCCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGAcc                                                                                                                                                                                                                                 <  6:535958/149‑1 (MQ=255)

                                                                                                                                            acaacaCCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGAcc                                                                                                                                                                                                                                 >  8:194213/1‑149 (MQ=255)

                                                                                                                                              aacaCCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAg                                                                                                                                                                                                                               <  7:628185/149‑1 (MQ=255)

                                                                                                                                              aacaCCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAg                                                                                                                                                                                                                               <  6:375511/149‑1 (MQ=255)

                                                                                                                                               acacCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGt                                                                                                                                                                                                                              >  6:262557/1‑149 (MQ=255)

                                                                                                                                               acacCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCg                                                                                                                                                                                        >  2:205891/1‑187 (MQ=255)

                                                                                                                                                cacCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGtt                                                                                                                                                                                                                             >  5:439682/1‑149 (MQ=255)

                                                                                                                                                cacCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGtt                                                                                                                                                                                                                             <  6:198425/149‑1 (MQ=255)

                                                                                                                                                 acCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACATCTGTGGGCtt                                                                                                                                                                                                                                                        <  6:609388/121‑1 (MQ=255)

                                                                                                                                                   caacaaAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTgg                                                                                                                                                                                                                                                            <  7:1442005/115‑1 (MQ=255)

                                                                                                                                                   caacaaAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGt                                                                                                                                                                                                                                                  <  7:858033/125‑1 (MQ=255)

                                                                                                                                                   caacaaAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGt                                                                                                                                                                                                                                                  >  8:858033/1‑125 (MQ=255)

                                                                                                                                                   caacaaAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTcac                                                                                                                                                                                                                          >  7:1088419/1‑149 (MQ=255)

                                                                                                                                                   caacaaAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGGGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTcac                                                                                                                                                                                                                          >  8:1122004/1‑149 (MQ=255)

                                                                                                                                                    aacaaAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGa                                                                                                                                                                                                                                               <  5:961411/127‑1 (MQ=255)

                                                                                                                                                    aacaaAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGa                                                                                                                                                                                                                                               >  6:961411/1‑127 (MQ=255)

                                                                                                                                                    aacaaAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTcaca                                                                                                                                                                                                                         <  8:1082480/149‑1 (MQ=255)

                                                                                                                                                    aacaaAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTcaca                                                                                                                                                                                                                         <  8:575931/149‑1 (MQ=255)

                                                                                                                                                    aacaaAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTcaca                                                                                                                                                                                                                         >  5:492910/1‑149 (MQ=255)

                                                                                                                                                     acaaAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCacac                                                                                                                                                                                                                        <  5:1392492/149‑1 (MQ=255)

                                                                                                                                                     acaaAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCacac                                                                                                                                                                                                                        >  8:282772/1‑149 (MQ=255)

                                                                                                                                                      caaAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACAcc                                                                                                                                                                                                                       <  7:687788/149‑1 (MQ=255)

                                                                                                                                                       aaaaaCAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACTCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCa                                                                                                                                                                                                                      >  7:1116648/1‑149 (MQ=255)

                                                                                                                                                        aaaaCAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTgg                                                                                                                                                                                                                                                            >  5:930574/1‑110 (MQ=255)

                                                                                                                                                        aaaaCAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTgg                                                                                                                                                                                                                                                            <  6:930574/110‑1 (MQ=255)

                                                                                                                                                         aaaCAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGc                                                                                                                                                                                                                    >  8:496296/1‑149 (MQ=255)

                                                                                                                                                            cAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACATCTGTGGGCtt                                                                                                                                                                                                                                                        <  7:96073/110‑1 (MQ=255)

                                                                                                                                                               gAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGAATGAGGCTTTAATACCTGCCTGGTATATTCCACGTCTGCCCGACAAATAATTg                                                                                  >  2:43288/1‑273 (MQ=255)

                                                                                                                                                                 aGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGt                                                                                                                                                                                                            <  8:26887/149‑1 (MQ=255)

                                                                                                                                                                   ttAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTc                                                                                                                                                                                                          >  6:371222/1‑149 (MQ=255)

                                                                                                                                                                      gAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCAc                                                                                                                                                                                                       <  8:13224/149‑1 (MQ=255)

                                                                                                                                                                        aGTCTCAGCAAAACACCGACCAGACGGTGAGGAGCCATAAAAGCATACCCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACgg                                                                                                                                                                                                     <  8:1370450/149‑1 (MQ=255)

                                                                                                                                                                        aGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACgg                                                                                                                                                                                                     <  5:381517/149‑1 (MQ=255)

                                                                                                                                                                         gTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACggg                                                                                                                                                                                                    >  8:1285999/1‑149 (MQ=255)

                                                                                                                                                                              aGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTc                                                                                                                                                                                               >  5:670613/1‑149 (MQ=255)

                                                                                                                                                                              aGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTc                                                                                                                                                                                               <  5:513493/149‑1 (MQ=255)

                                                                                                                                                                               gCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTcc                                                                                                                                                                                              <  6:1448665/149‑1 (MQ=255)

                                                                                                                                                                                cAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCca                                                                                                                                                                                             <  7:1391805/149‑1 (MQ=255)

                                                                                                                                                                                 aaaaCATCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCcac                                                                                                                                                                                            >  6:1045719/1‑149 (MQ=255)

                                                                                                                                                                                      acCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCgg                                                                                                                                                                                       >  6:512545/1‑149 (MQ=255)

                                                                                                                                                                                       ccGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGc                                                                                                                                                                                      >  5:259680/1‑149 (MQ=255)

                                                                                                                                                                                         gACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCAc                                                                                                                                                                                    >  7:817956/1‑149 (MQ=255)

                                                                                                                                                                                           ccAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTa                                                                                                                                                                                  <  6:100104/149‑1 (MQ=255)

                                                                                                                                                                                                   tGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCtata                                                                                                                                                                          <  5:34316/149‑1 (MQ=255)

                                                                                                                                                                                                    gaggagACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATac                                                                                                                                                                         <  8:374868/149‑1 (MQ=255)

                                                                                                                                                                                                       gagaCATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATacccc                                                                                                                                                                      >  6:25284/1‑149 (MQ=255)

                                                                                                                                                                                                          aCATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATaccccacc                                                                                                                                                                   >  6:1295787/1‑149 (MQ=255)

                                                                                                                                                                                                            aTAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATaccccacccc                                                                                                                                                                 >  6:86075/1‑149 (MQ=255)

                                                                                                                                                                                                              aaaaGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGtt                                                                                                                                                                                                           <  8:1448328/105‑1 (MQ=255)

                                                                                                                                                                                                               aaaGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACggg                                                                                                                                                                                                    <  8:155414/111‑1 (MQ=255)

                                                                                                                                                                                                               aaaGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACggg                                                                                                                                                                                                    >  7:155414/1‑111 (MQ=255)

                                                                                                                                                                                                               aaaGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAg                                                                                                                                                              >  7:1096226/1‑149 (MQ=255)

                                                                                                                                                                                                                aaGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGtt                                                                                                                                                                                                           >  8:1423629/1‑103 (MQ=255)

                                                                                                                                                                                                                aaGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGa                                                                                                                                                             <  7:399093/149‑1 (MQ=255)

                                                                                                                                                                                                                aaGGATACGCAAAGCAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGa                                                                                                                                                             >  5:1111755/1‑149 (MQ=255)

                                                                                                                                                                                                                     tACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGccc                                                                                                                                                        <  8:1096226/149‑1 (MQ=255)

                                                                                                                                                                                                                        gCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGtt                                                                                                                                                                                                                                                 >  5:596191/1‑57 (MQ=255)

                                                                                                                                                                                                                        gCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCAc                                                                                                                                                                                                       >  7:1227837/1‑99 (MQ=255)

                                                                                                                                                                                                                        gCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCt                                                                                                                                                     >  8:318718/1‑149 (MQ=255)

                                                                                                                                                                                                                          aaaGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTa                                                                                                                                                   >  6:75846/1‑149 (MQ=255)

                                                                                                                                                                                                                          aaaGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTa                                                                                                                                                   >  6:521373/1‑149 (MQ=255)

                                                                                                                                                                                                                           aaGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTaa                                                                                                                                                  >  6:240212/1‑149 (MQ=255)

                                                                                                                                                                                                                            aGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAAc                                                                                                                                                 <  7:501327/149‑1 (MQ=255)

                                                                                                                                                                                                                            aGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAAc                                                                                                                                                 <  8:846339/149‑1 (MQ=255)

                                                                                                                                                                                                                            aGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAAc                                                                                                                                                 <  5:86075/149‑1 (MQ=255)

                                                                                                                                                                                                                            aGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAAc                                                                                                                                                 <  5:737538/149‑1 (MQ=255)

                                                                                                                                                                                                                             ggAGCCGCGGCACCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACg                                                                                                                                                >  7:1261217/1‑149 (MQ=255)

                                                                                                                                                                                                                              gAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATg                                                                                                                                                           <  6:448691/137‑1 (MQ=255)

                                                                                                                                                                                                                                gCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTAc                                                                                                                                             <  6:582837/149‑1 (MQ=255)

                                                                                                                                                                                                                                    cgGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTActgct                                                                                                                                         <  5:1295787/149‑1 (MQ=255)

                                                                                                                                                                                                                                     ggCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGAATGAGGCTTTAATACCTGCCTGGTATATTCCACGTCTGCCCGACAAATAATTGACGAAATTACCGTCACTATTCACTTTCACCCGGTTATCATCGACAAATTCTTTGCGCACAGCCGCCTTCAGc          <  4:43288/275‑1 (MQ=255)

                                                                                                                                                                                                                                         ccagatAACATGAAAGCCCACATAGTTGGGATTATTCGTTGATGCGCAGAAGGACCAGCTCACACCAGCTAGTGCGTTCCACGGGGAGTCCCCACCGGCACGATGGCTATACCCCACCCCAAGATGCACGCTTAACGCACTCCTCAATg                                                                                                                                    >  8:1226125/6‑149 (MQ=255)

                                                                                                                                                                                                                                             gTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCCGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGAATGAGGc                                                                                                                                >  8:714194/1‑149 (MQ=255)

                                                                                                                                                                                                                                              tAACATAAAAGCCCTCTGATGTGGGCTTTTTCGTTGAGGGTAAGAACGACCAGTTCAAACCGTCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGAATGAGGCTTTAATACCTGCCTGGTATATTCCACGTCTGCCCGACAAATAATTGACGAAATTACCGTCACTATTCACTTTCACCCGGTTATCATCGACAAATTCTTTGCGCACAGCCGCCTTCAGCCACGGCTc  <  2:126342/274‑1 (MQ=255)

                                                                                                                                                                                                                                                    gAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGAATGAGGCTTTAATa                                                                                                                         >  7:72314/1‑149 (MQ=255)

                                                                                                                                                                                                                                                       aGCCCACATCTGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATgg                                                                                                                                                                               <  7:909065/92‑1 (MQ=255)

                                                                                                                                                                                                                                                       aGCCCACATCTGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGAATGAGGCTTTAATAcct                                                                                                                      >  8:50358/1‑149 (MQ=255)

                                                                                                                                                                                                                                                          ccACAGATGTGGGATTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGAATGAGGCTTTAATAcctgcc                                                                                                                   <  8:1242017/149‑1 (MQ=255)

                                                                                                                                                                                                                                                              aGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGAATGAGGCTTTAATACCTGCCTGGt                                                                                                               <  7:282772/149‑1 (MQ=255)

AAAAAACCTGACATTTTCTTTTAGGACCAACAATGGGACCAAAATGAAAATTGAACTGAAGATTACAATCAATTTAACAACAAGTTACACAACCAATTCAGACTCCGCCAGCCCACCAATCATGATTGGACGGTGTAAGGACAACACCAACAAAAACAGGAAGTTAGAAGTCTCAGCAAAACACCGACCAGACGGTGAGGAGACATAAAAGGATACGCAAAGGAGCCGCGGCTCCTGGTAACATGAAAGCCCACAGATGTGGGCTTTTTCGTTGATGGTCAGAACGACCAGTTCACACCAGCTACCGCGTTCCACGGGGATTCCACACCGGCACTATGGCTATACCCCACCCCAAGATGCCCGCTTAACGTACTGCTGAATGAGGCTTTAATACCTGCCTGGTATATTCCACGTCTGCCCGACAAATAATTGACGAAATTACCGTCACTATTCACTTTCACCCGGTTATCATCGACAAATTCTTTGCGCACAGCCGCCTTCAGCCACGGCTC  >  NC_000913/2777868‑2778379

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