BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,070,413	T→G	6.6%	intergenic (‑12/+102)	yihR ← / ← yihS	putative sulphoquinovose mutarotase/sulphoquinovose isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,070,413	0	T	G	6.6%	747.3 / 25.7	220	intergenic (‑12/+102)	yihR/yihS	putative sulphoquinovose mutarotase/sulphoquinovose isomerase
Reads supporting (aligned to +/- strand): ref base T (111/93); new base G (6/10); total (117/103)
Fisher's exact test for biased strand distribution p-value = 2.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 224 total aligned reads displayed.

AACTGATATTCCTGCCCCTGATAGCGGTAACAGCCGTTGGCGATCCGGTTCGGCCACGGAATAAGCACTTTACCCAGATGCGCCAGCGGCATCTCTTCGGGTTTGTGCGGGATGACTAAATGACAACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCC‑AAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAATGCAGCAGGTGATAAATATCCTGTTTGCCGTCCCAGACTTTGGTGGTGACCTTATTGTCCGCATCCAGCTCCTGCCACCAGGAACCA  >  NC_000913/4070139‑4070686

             tccccTGTAAGCGGGAACAGCCGTTGGGGATCCGGTTGGGCCACGGAATAAGCACTTTACCCAGATGCGCCAGCGGCATCTCTTCGGGTTTGTGCGGGATGACTAAATGACAACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCAAAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGCCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCa                                                                                                                                                                                                                                                                       <  4:496530/274‑1 (MQ=255)

                                            ccGGTTCGGCCCCGGAATAAGCACGTTACCCAGAGTCGCCAGCGGCATCTCTTCGGGGTTGCGCGGGAGGACGAAATGACAACCCTGAAAGGTTAATTCAGCCAGACCAGCACAGACGGTCACAATGGTGGCGTGATAATCCCCCGCGGCAAGTGAAACAGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCC‑AAAAACTGCACCGGACGGTCCCCTCg                                                                                                                                                                                                                                       <  4:543559/275‑1 (MQ=255)

                                             cGGTTCGGCCACGGAATAAGCACTTTACCCAGATGCGCCAGCGGCATCTCTTCGGGTTTGTGCGGGATGACTAAATGACAACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGt                                                                                                                                                                                                                                                                                     <  2:84806/230‑1 (MQ=255)

                                                               aGCACTTTACCCAGATGCGCCAGCGGCATCTCTTCGGGTTTGTGCGGGATGACTAAATGACAACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCC‑a                                                                                                                                                                                                                                                                >  2:165186/1‑231 (MQ=255)

                                                                                  ccAGCGGCATCTCTTCGGGTTTGTGCGGGATGACTAAATGACAACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCC‑AAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGgaga                                                                                                                                                                                                                         <  2:499899/251‑1 (MQ=255)

                                                                                                                               ccTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTg                                                                                                                                                                                                                                                                                    >  6:420390/1‑149 (MQ=255)

                                                                                                                                 tGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAg                                                                                                                                                                                                                                                                                  >  8:1068082/1‑149 (MQ=255)

                                                                                                                                 tGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAg                                                                                                                                                                                                                                                                                  >  8:1355720/1‑149 (MQ=255)

                                                                                                                                       gTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGGTAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑At                                                                                                                                                                                                                                                                           <  7:1404626/149‑1 (MQ=255)

                                                                                                                                       gTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑At                                                                                                                                                                                                                                                                           <  5:1162400/149‑1 (MQ=255)

                                                                                                                                        gaaaTTCAGCCAGACCAGCAGTGACGGACACAGTGGTGGGGTGATAATCCCCGGCGGCAAGGGAACCGGTCTGTCAACAGCAATGAAGAGTGGTGATTTGCAGACAAGGCATCATAATTAACGACATACAGCCTCCCGTGAGGCTG‑AGGTCAAACCC‑AAGAACGGCACCGGACGGTCCCCTCGCCCCTTTGGTGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCa                                                                                                                                            <  4:283565/272‑1 (MQ=255)

                                                                                                                                         tAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCC‑AAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCCCCCTATTTTCAGAACAGTTATTTCGCATTACAATCCAGCAGACCCGCTGCaa                                                                                                                                           >  2:557548/1‑274 (MQ=255)

                                                                                                                                         tAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACAT‑CAGCCTCACGGGAGGCTgtatgt                                                                                                                                                                                                                                                                         <  7:762131/149‑6 (MQ=255)

                                                                                                                                          aaTTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTc                                                                                                                                                                                                                                                                        <  5:849699/149‑1 (MQ=255)

                                                                                                                                           aTTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCa                                                                                                                                                                                                                                                                       >  6:812617/1‑149 (MQ=255)

                                                                                                                                            ttCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCaa                                                                                                                                                                                                                                                                      >  7:151348/1‑149 (MQ=255)

                                                                                                                                             tCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACAT‑CAGCCTCACGGGAGGCTGTATGTcgt                                                                                                                                                                                                                                                                      <  7:865834/148‑3 (MQ=255)

                                                                                                                                                gCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAAccc                                                                                                                                                                                                                                                                  <  7:1356608/149‑1 (MQ=255)

                                                                                                                                                     aCCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCC‑AAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCCACCGCTGGAGCAAt                                                                                                                              >  2:299636/1‑275 (MQ=255)

                                                                                                                                                       cAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCTGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCC‑AAAAACt                                                                                                                                                                                                                                                          >  5:499518/1‑149 (MQ=255)

                                                                                                                                                          cACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCC‑AAAAACTGCa                                                                                                                                                                                                                                                       <  8:430670/149‑1 (MQ=255)

                                                                                                                                                          cACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCC‑AAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGga                                                                                                                                                                                                                           <  3:207332/177‑1 (MQ=255)

                                                                                                                                                                   tCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCC‑AAAAACTGCACCGGACGGt                                                                                                                                                                                                                                              >  6:1261843/1‑149 (MQ=255)

                                                                                                                                                                   tCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCC‑AAAAACTGCACCGGACGGt                                                                                                                                                                                                                                              >  5:19197/1‑149 (MQ=255)

                                                                                                                                                                       aaTGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCC‑AAAAACTGCACCGGACGGTcccc                                                                                                                                                                                                                                          >  6:1402020/1‑149 (MQ=255)

                                                                                                                                                                        aTGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCC‑AAAAACTGCACCGGACGGTCCCCt                                                                                                                                                                                                                                         >  8:387165/1‑149 (MQ=255)

                                                                                                                                                                           gTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCC‑AAAAACTGCACCGGACGGTCCCCTCGc                                                                                                                                                                                                                                      <  6:966979/149‑1 (MQ=255)

                                                                                                                                                                              gcgTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGa                                                                                                                                                                                                                                                                                   <  8:862895/103‑1 (MQ=255)

                                                                                                                                                                               cgTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCC‑AAAAACTGCACCGGACGGTCCCCTCGCCCCt                                                                                                                                                                                                                                  >  7:426620/1‑149 (MQ=255)

                                                                                                                                                                                gTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCC‑AAAAACTGCACCGGACGGTCCCCTCGCCCCtt                                                                                                                                                                                                                                 >  5:1377408/1‑149 (MQ=255)

                                                                                                                                                                                 tGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCC‑AAAAACTGCACCGGACGGTCCCCTCGCCCCttt                                                                                                                                                                                                                                >  6:1050379/1‑149 (MQ=255)

                                                                                                                                                                                     aaTCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCC‑AAAAACTGCACCGGACGGTCCCCTCGCCCCTTTgggg                                                                                                                                                                                                                            >  5:1181384/1‑149 (MQ=255)

                                                                                                                                                                                       tCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCC‑AAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGga                                                                                                                                                                                                                           >  8:526237/1‑148 (MQ=255)

                                                                                                                                                                                        cccccGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCC‑AAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGgaga                                                                                                                                                                                                                         <  8:154071/149‑1 (MQ=255)

                                                                                                                                                                                         ccccGCGGGAAGTGAAGCGGTCTGTGCCCTGCACTGCATATTGGTGGTTTGCATACAAGGCACGTTAAGTAACGACATACAGCCTCCCGGGAGGCTG‑ATGTCAAACCC‑AAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGagag                                                                                                                                                                                                                        >  6:24788/1‑149 (MQ=255)

                                                                                                                                                                                         ccccGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCC‑AAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGagag                                                                                                                                                                                                                        >  8:1260580/1‑149 (MQ=255)

                                                                                                                                                                                          cccGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCC‑AAAAACTGCACCGGACGGTCCCCTCTCCCCAAAGGGGCGAgg                                                                                                                                                                                                                       <  7:524990/149‑1 (MQ=255)

                                                                                                                                                                                          cccGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCC‑AAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAgg                                                                                                                                                                                                                       <  5:427042/149‑1 (MQ=255)

                                                                                                                                                                                          cccGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCC‑AAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAgg                                                                                                                                                                                                                       >  7:1462289/1‑149 (MQ=255)

                                                                                                                                                                                          cccGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCC‑AAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAgg                                                                                                                                                                                                                       >  5:1465397/1‑149 (MQ=255)

                                                                                                                                                                                          cccGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCC‑AAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAgg                                                                                                                                                                                                                       >  8:1145133/1‑149 (MQ=255)

                                                                                                                                                                                           ccGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCC‑AAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAggg                                                                                                                                                                                                                      <  7:671699/149‑1 (MQ=255)

                                                                                                                                                                                               ggCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCC‑AAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAg                                                                                                                                                                                                                  <  5:990195/149‑1 (MQ=255)

                                                                                                                                                                                                gCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCC‑AAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAgg                                                                                                                                                                                                                 <  5:812617/149‑1 (MQ=255)

                                                                                                                                                                                                gCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCC‑AAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAgg                                                                                                                                                                                                                 <  5:1402020/149‑1 (MQ=255)

                                                                                                                                                                                                  aaGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCC‑AAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGga                                                                                                                                                                                                               <  6:851543/149‑1 (MQ=255)

                                                                                                                                                                                                  aaGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCC‑AAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGga                                                                                                                                                                                                               >  7:1616708/1‑149 (MQ=255)

                                                                                                                                                                                                      gAAACTGTCTGTCCACTGCAATGCATATTGGTGATTAGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCC‑AAAAACTGCACCGGACGGTCCCCTCGGCCCTTTGGGGAGAGGGTTAGGGAGAgg                                                                                                                                                                                                           <  8:282166/149‑1 (MQ=255)

                                                                                                                                                                                                      gAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCC‑AAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAgg                                                                                                                                                                                                           <  7:1260580/149‑1 (MQ=255)

                                                                                                                                                                                                         gcGGTCTGCCCGCGCCAATGCAAATGGGGGATTTGCATACAAGGCACCTAAATTAACGCCATACAGCCTCCCGGGAGGCTG‑ATGTAAAACCC‑AAAAACTGCACCGGACGGCCCCCTCGCCCCTTTGGGGAGAGGGTGAGGGAGAGGGGa                                                                                                                                                                                                        <  8:1549156/148‑1 (MQ=255)

                                                                                                                                                                                                             tCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCC‑AAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAAc                                                                                                                                                                                                    >  5:309831/1‑149 (MQ=255)

                                                                                                                                                                                                               tGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCC‑AAAAACTGCACCGGACGGTccc                                                                                                                                                                                                                                           >  2:44353/1‑108 (MQ=255)

                                                                                                                                                                                                                 tCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCC‑AAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAggg                                                                                                                                                                                                          >  7:1570552/1‑139 (MQ=255)

                                                                                                                                                                                                                 tCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCC‑AAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAggg                                                                                                                                                                                                          <  8:1570552/139‑1 (MQ=255)

                                                                                                                                                                                                                    aCTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCC‑AAAAACTGCACCGGACGGTCCCCTCGCCCCTTTggg                                                                                                                                                                                                                             >  8:1602388/1‑117 (MQ=255)

                                                                                                                                                                                                                      tGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCC‑AAAAACTGCACCGGACGGTCCCCTCGCCCCtt                                                                                                                                                                                                                                 >  1:154139/1‑111 (MQ=255)

                                                                                                                                                                                                                      tGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCC‑AAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCt                                                                                                                                                                                           >  7:28758/1‑149 (MQ=255)

                                                                                                                                                                                                                        cAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCC‑AAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTAt                                                                                                                                                                                         >  7:201534/1‑149 (MQ=255)

                                                                                                                                                                                                                         aaTGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCC‑AAAAACTGCACCGGACGGTCCCCTCGCCCTTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTAtt                                                                                                                                                                                        <  8:780058/149‑1 (MQ=255)

                                                                                                                                                                                                                           tGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCC‑AAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCaa                                                                                                                                           >  2:694420/1‑192 (MQ=255)

                                                                                                                                                                                                                           tGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCC‑AAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAATGCAGCAGGTTATAAATATCCTGTTTTCCGTCCCaaa                                                    >  1:282235/1‑277 (MQ=255)

                                                                                                                                                                                                                             cATATTGTTGATTTGCATACAAGGCCCATTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCC‑AAAACCTGCACCGGACGGCCCCCTCGCCCCTTTGGGGGGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCa                                                                                                                                                                                    <  7:1599614/149‑1 (MQ=39)

                                                                                                                                                                                                                             cATATTGGGGATTTGCAAACAAGGCACCTAAATTAACGACACACAGCCTCCCGTGAGGCTG‑ATGTCAAACCC‑AAAAACTGCACCGGACGGCCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCa                                                                                                                                                                                    <  8:1345926/149‑1 (MQ=39)

                                                                                                                                                                                                                               tatTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCC‑AAAAACTGCACCGGACGGTCCCCTCGcc                                                                                                                                                                                                                                     <  2:308968/98‑1 (MQ=255)

                                                                                                                                                                                                                                atTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCC‑AAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAATGCAGCAGGTGATAAATATCCTGTTTGCCGTCCCCGCCTTtg                                               >  1:208026/1‑279 (MQ=255)

                                                                                                                                                                                                                                 ttGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCC‑AAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATc                                                                                                                                                                                >  7:75216/1‑149 (MQ=255)

                                                                                                                                                                                                                                 tctgtgtTTTGCATACAAGCCACCTTAATTAACGACATACCCCCCCCCGGGAGGCTG‑ATGTCAAACCC‑AAAAACTGCCCCGGACGGCCCCCCCCCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATc                                                                                                                                                                                <  8:743702/146‑1 (MQ=255)

                                                                                                                                                                                                                                    gTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCC‑AAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGt                                                                                                                                                                             >  6:793449/1‑149 (MQ=255)

                                                                                                                                                                                                                                     tGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCC‑AAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGG‑TAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTa                                                                                                                                                                           >  6:58711/1‑149 (MQ=255)

                                                                                                                                                                                                                                      gATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTAAAACCC‑AAAAACGGCACCGGACGGCCCCCTCGCCCTTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTa                                                                                                                                                                           <  6:173913/149‑1 (MQ=255)

                                                                                                                                                                                                                                       aTTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCC‑AAAAACTGCACCGGACGGTccc                                                                                                                                                                                                                                           >  6:142270/1‑84 (MQ=255)

                                                                                                                                                                                                                                           gCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCC‑AAAAACTGCACCGGACGGTCCCCTCGCCCCtt                                                                                                                                                                                                                                 >  7:1687184/1‑90 (MQ=255)

                                                                                                                                                                                                                                           gCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCC‑AAAAACTGCACCGGACGGTCCCCTCGCCCCtt                                                                                                                                                                                                                                 <  8:1687184/90‑1 (MQ=255)

                                                                                                                                                                                                                                           gCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCC‑AAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCg                                                                                                                                                                      >  6:1236300/1‑149 (MQ=255)

                                                                                                                                                                                                                                               aCAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCC‑AAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCAtt                                                                                                                                                                  >  5:1102438/1‑149 (MQ=255)

                                                                                                                                                                                                                                                cAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCC‑AAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTa                                                                                                                                                                 >  8:646457/1‑149 (MQ=255)

                                                                                                                                                                                                                                                   ggCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCC‑AAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCCCCAGGCAATGCAGCAGGTGATAAATAACCTGTTTGCCGTCCCAGACTTTGGTGGTGACCTTAtt                                 >  4:384019/1‑274 (MQ=255)

                                                                                                                                                                                                                                                   ggCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCC‑AAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAATGCAGCAGGTGATAAATATCCTGTTTGCCGTCCCAGACTTTGGTGGTGACCTTAt                                  >  1:11879/1‑273 (MQ=255)

                                                                                                                                                                                                                                                       ccTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATcc                                                                                                                                                           >  5:694804/1‑149 (MQ=255)

                                                                                                                                                                                                                                                       ccTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCC‑AAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCca                                                                                                                                                          >  7:85663/1‑149 (MQ=255)

                                                                                                                                                                                                                                                        cTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCC‑AAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCcag                                                                                                                                                         >  5:1634505/1‑149 (MQ=255)

                                                                                                                                                                                                                                                        cTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCC‑AAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCcag                                                                                                                                                         >  7:34031/1‑149 (MQ=255)

                                                                                                                                                                                                                                                        cTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCC‑AAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCcag                                                                                                                                                         >  6:1326232/1‑149 (MQ=255)

                                                                                                                                                                                                                                                        cTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCC‑AAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCcag                                                                                                                                                         >  7:553222/1‑149 (MQ=255)

                                                                                                                                                                                                                                                         ttaattaaCGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCC‑AAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGCTATTTCGCATTAATATCcagc                                                                                                                                                        >  7:1560599/1‑149 (MQ=255)

                                                                                                                                                                                                                                                           aattaaCGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCC‑AAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCcagcag                                                                                                                                                      >  5:1299457/1‑149 (MQ=255)

                                                                                                                                                                                                                                                            attaaCGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCC‑AAAAACGGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGa                                                                                                                                                     <  8:593428/149‑1 (MQ=255)

                                                                                                                                                                                                                                                             ttaaCGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCC‑AAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGAc                                                                                                                                                    >  6:730652/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                aCGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCC‑AAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCg                                                                                                                                                 >  6:113287/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                aCGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCC‑AAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCg                                                                                                                                                 >  5:576399/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                 cGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCC‑AAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGc                                                                                                                                                >  6:1530164/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                  gACATCCCGCCTCCCGTGAGGCTG‑ATGTCAAACCC‑AAAACCTGCACCGGACGGTCCCCTCGCCCTTTTGGGGAGAGGGTTGGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCt                                                                                                                                               <  6:1634505/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                     aTACAGCCTCCCGTGAGGCTG‑ATGTCAAACCC‑AAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCa                                                                                                                                            >  5:530781/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                      tACAGCCTCCCGTGAGGCTG‑ATGTCAAACCC‑AAAAACTGCACCGGACGGTc                                                                                                                                                                                                                                             <  5:1367300/51‑1 (MQ=255)

                                                                                                                                                                                                                                                                      atcAGCCTCACGGGAGGCTGTATGTCAAACCC‑AAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCa                                                                                                                                            >  7:1430324/3‑149 (MQ=39)

                                                                                                                                                                                                                                                                      atcAGCCTCACGGGAGGCTG‑ATGTCAAACCC‑AAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCaa                                                                                                                                           >  6:541496/3‑149 (MQ=255)

                                                                                                                                                                                                                                                                          gCCTCCCGTGAGGCTG‑ATGTCAAACCC‑AAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGc                                                                                                                                       >  8:852052/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                           ccTCCCGTGAGGCTG‑ATGTCAAACCC‑AAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCt                                                                                                                                      <  8:28758/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                            cTCCCGTGAGGCTG‑ATGTCAAACCC‑AAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGCTCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTg                                                                                                                                     >  6:1447337/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                                  tGAGGCTG‑ATGTCAAACCC‑AAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGcc                                                                                                                                >  5:1008682/1‑148 (MQ=255)

AACTGATATTCCTGCCCCTGATAGCGGTAACAGCCGTTGGCGATCCGGTTCGGCCACGGAATAAGCACTTTACCCAGATGCGCCAGCGGCATCTCTTCGGGTTTGTGCGGGATGACTAAATGACAACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCC‑AAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAATGCAGCAGGTGATAAATATCCTGTTTGCCGTCCCAGACTTTGGTGGTGACCTTATTGTCCGCATCCAGCTCCTGCCACCAGGAACCA  >  NC_000913/4070139‑4070686

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