BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,388,285	A→G	5.5%	intergenic (+24/+20)	tyrR → / ← tpx	aromatic amino acid biosynthesis and transport regulon transcriptional regulator; autorepressor; ATPase; phosphatase/lipid hydroperoxide peroxidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,388,285	0	A	G	5.5%	729.4 / 21.3	222	intergenic (+24/+20)	tyrR/tpx	aromatic amino acid biosynthesis and transport regulon transcriptional regulator; autorepressor; ATPase; phosphatase/lipid hydroperoxide peroxidase
Reads supporting (aligned to +/- strand): ref base A (113/96); new base G (3/9); total (116/106)
Fisher's exact test for biased strand distribution p-value = 7.27e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.33e-02

Only 100 of 238 total aligned reads displayed.

TGGACGGTTATGAGCTGCGTCCACAGGATATTTTGTTGCCGGATTATGACGCCGCAACGGTAGCCGTGGGCGAAGATGCGATGGAAGGTTCGCTGGACGAAATCACCAGCCGTTTTGAACGCTCGGTATTAACCCAGCTTTATCGCAATTATCCCAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAA‑CACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACGGCACGCGCTGCCAGACCTTTCAGTGGGCCATCAGCAATTGCCACACCGTAAGCTTGCAGAAATTCAGCGTTACGGAAAGTGGAGAGGGTGATAACGTTGTTCAGACCTTCTGCGCCGCAGAAACGAGACTGGGCGAACGGCAGATCGGCA  >  NC_000913/1388008‑1388560

   aCGGTTATGAGCTGCGTCCACAGGATATTTTGTTGCCGGATTATGACGCCGCAACGGTAGCCGTGGGCGAAGATGCGATGGAAGGTTCGCTGGACGAAATCACCAGCCGTTTTGAACGCTCGGTATTAACCCAGCTTTATCGCAATTATCCCAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCaa                                                                                                                                                                                                                                                                                      <  2:292982/275‑1 (MQ=255)

                 cGTCCACAGGATATTTTGTTGCCGGATTATGATGCCGCAACGGTAGCCGTGGGCGAAGATGCGAACGAAGGTTCTATGGACGAAATCACCAAACGAGTTTAACGCTCGGTAATAACCCAACTTTTTCGCCATTATCCCAGCACGCGCAACCTTGCAAGACGTATCGGCGTTTCACATACCGCGATTGCCGACAAGTTCCGGGAAGACGGTCGGAGGCAGAAGGAGAACGGAGGGGAGGCGGGAAAGTGCATGGTTGTGCAa‑caaccgcgagaatg                                                                                                                                                                                                                                                                       >  4:38832/1‑263 (MQ=255)

                       caGGGTATTTGGGTGCCGGATGAGTAGGGCTAAAGGGTAGGCGTGGGGGAAGATGCGATGGAAGGGTCGCTTGACGAAAGCACCAGCCGGTTTGCACGCTAGGTATTAACCCAGCTATATCGAAAGTACCCCAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATAACGCGATTGCCAATAAGTTGCGAGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGAGCAA‑CACCATCAGGCATATTAAAt                                                                                                                                                                                                                                                                 <  4:55858/275‑1 (MQ=255)

                                           ttATGACGCCGCAACGGTAGCCGTGGGCGAAGATGCGATGGAAGGTTCGCTGGACGAAATCACCAGCCGTTTTGAACGCTCGGTATTAACCCAGCTTTATCGCAATTATCCCAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAA‑CACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCag                                                                                                                                                                                                                                             <  2:502165/275‑1 (MQ=255)

                                                 cgccgcAACGGTAGCCGTGGGCGAAGATGCGATGGAAGGTTCGCTGGACGAAATCACCAGCCGTTTTGAACGCTCGGTATTAACCCAGCTTTATCGCAATTATCCCAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGCAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAA‑CACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTc                                                                                                                                                                                                                                    >  1:230498/1‑278 (MQ=255)

                                                    cgcAACGGTAGCCGTGGGCGAAGATGCGATGGAAGGTTCGCTGGACGAAATCACCAGCCGTTTTGAACGCTCGGTATTAACCCAGCTTTATCGCAATTATCCCAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTTCAA‑CACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTaa                                                                                                                                                                                                                                >  1:209698/1‑279 (MQ=255)

                                                                         aGATGCGATGGAAGGTTCGCTGGACGAAATCACCAGCCGTTTTGAACGCTCGGTATTAACCCAGCTTTATCGCAATTATCCCAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAA‑CACCATCAGGCATAtt                                                                                                                                                                                                                                                                     <  2:624779/221‑1 (MQ=255)

                                                                              cGATGGAAGGTTCGCTGGACGAAATCACCAGCCGTTTTGAACGCTCGGTATTAACCCAGCTTTATCGCAATTATCCCAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGCATATGGTTTGAGTCAGAAGAAGAACGAAGAGGAAGCGCGAATATGACTGATGGTGCAA‑CACCAACACGAATATTAAATTATGCTTTCAGTAACGCCAGAGCTGCTTCGTATTCCGGCTCGGTGGTGAttttct                                                                                                                                                                                                          >  4:550513/1‑272 (MQ=255)

                                                                                  ggaggtttgccTGGAAGAAATCGCCGGCGGTTTTAACCGTGAGGTATTAACCCAGCTTTATCGCAATTATCCCAGCAGACNAAAACTGGCAAAACGTTGCGGCGTTTCACATACCGGGTTTGCCAATAAGTTGGGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAA‑CACCATCAGGCATATTAAAGTATGCTTTCAGTACAGCCAGAGCTGCGTCGTAATCCGGCTCGGTGGGGATTTCATCCACCAg                                                                                                                                                                                                   <  3:230498/268‑1 (MQ=255)

                                                                                         tCGCTGGAAGCAATCACAAGCCGGTTTTAACGCTCGGTATTAACCCAGATTTATCGAAATTATCCCAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAA‑CACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAgctgg                                                                                                                                                                                               <  4:130874/275‑1 (MQ=255)

                                                                                              ggACGAAATCACCAGCCGTTTTGAACGCTCGGTATTAACCCAGCTTTATCGCAATTATCCCAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAA‑CACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGaa                                                                                                                                                                                          >  2:606082/1‑275 (MQ=255)

                                                                                                                ttttTAACGCTCGTTATTACCCCAGCTTTAGCGCAATTATCCCAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGACAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAA‑CACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCa                                                                                                                                                                    <  1:27733/279‑1 (MQ=255)

                                                                                                                  ttGAACGCTCGGTATTAACCCAGCTTTATCGCAATTATCACAGCACGCGCAAAATGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCAGAGGAAGAAGAAGAACGAAGAGTAAGCGCGAATATGCGTGATGGTGAAACCACCA‑CAGGCAGATCAAATTATGTTTTCAGTACAGCCAAGACTGCCTCGTACACGCGATCTATGGTGATTGTATCCACAAGCTGACTGGAAATCAACTTGTCATTTTCGt                                                                                                                                                                      >  4:143795/1‑275 (MQ=255)

                                                                                                                    gAAAGCTCGGTATTACCCCAGCTTTATCGCAATTATCCCCGCACGCGCAAACTGGCAAAACGTACATGCGTTGCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAA‑CACCATCAGGCATATTAAATTAGGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCa                                                                                                                                                                    <  4:606078/275‑1 (MQ=255)

                                                                                                                                  aaCCCAGCTTTATCGCAATTATCCCAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCaa                                                                                                                                                                                                                                                                                      >  5:820969/1‑148 (MQ=255)

                                                                                                                                     ccAGCTTTATCGCAATTATCCCAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGACAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAA‑CAcc                                                                                                                                                                                                                                                                                 >  6:1271755/1‑149 (MQ=255)

                                                                                                                                        gCTTTATCGCAATTATCCCAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAA‑CACCATc                                                                                                                                                                                                                                                                              <  6:832076/149‑1 (MQ=255)

                                                                                                                                        gCTTTATCGCAATTATCCCAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAA‑CACCATCAGGCATATTAAAtt                                                                                                                                                                                                                                                                >  1:115617/1‑163 (MQ=255)

                                                                                                                                        gCTTTATCGCAATTATCCCAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAA‑CACCATCAGGCATATTAAAtt                                                                                                                                                                                                                                                                <  3:115617/163‑1 (MQ=255)

                                                                                                                                        gCTTTATCGCAATTATCCCAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAA‑CACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACGGCAcgcg                                                                                                                                                 >  2:551596/1‑274 (MQ=255)

                                                                                                                                        gCTTTATCGCAATTACCCCAGCACGCGCAACCTGGCAGAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAA‑CACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCt                                                                                                                                                                                                                         <  1:148650/202‑1 (MQ=255)

                                                                                                                                          tttATCGCAATTATCCCAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGTTG‑CACCATCAg                                                                                                                                                                                                                                                                            <  6:1007696/149‑1 (MQ=255)

                                                                                                                                          tttATCGCAATTATCCCAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGTTG‑CACCATCAg                                                                                                                                                                                                                                                                            <  5:1175756/149‑1 (MQ=255)

                                                                                                                                          tttATCGCAATTATCCCAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGTTG‑CACCATCAg                                                                                                                                                                                                                                                                            <  6:716950/149‑1 (MQ=255)

                                                                                                                                          tttATCGCAATTATCCCAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATACGCCTGATGGTGTTG‑CACCATCAg                                                                                                                                                                                                                                                                            <  8:1218291/149‑1 (MQ=255)

                                                                                                                                            tATCGCAATTATCCCAGCACGCGCAAACTGTCAAAACGTCTCGGCAGTTCACATACCGCGATTGCCAATAAGTGGTGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAA‑CACCATTAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACGACGGCACGCACTGcc                                                                                                                                            <  4:84342/275‑1 (MQ=255)

                                                                                                                                            tATCGCAATTATCCCAGCACGCGCAAACTGGCAAAACGTCTCTGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAA‑CACCAt                                                                                                                                                                                                                                                                               <  2:665979/144‑1 (MQ=255)

                                                                                                                                               cGCAATTATCCCAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGTTG‑CACCATCAg                                                                                                                                                                                                                                                                            <  5:1378664/144‑1 (MQ=255)

                                                                                                                                               cGCAATTATCCCAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAA‑CACCATCAGGCata                                                                                                                                                                                                                                                                       >  5:1041119/1‑149 (MQ=255)

                                                                                                                                                gCAATTATCCCAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAA‑CACCATCAGGCatat                                                                                                                                                                                                                                                                      >  5:34697/1‑149 (MQ=255)

                                                                                                                                                  aaTTATCCCAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAA‑CACCATCAGGCATATTa                                                                                                                                                                                                                                                                    >  8:50227/1‑149 (MQ=255)

                                                                                                                                                    ttATCCCAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAA‑CACCATCAGGCATATTaaa                                                                                                                                                                                                                                                                  <  8:272492/149‑1 (MQ=255)

                                                                                                                                                       tCCCAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAA‑CACCATCAGGCATATTAAATTa                                                                                                                                                                                                                                                               <  8:1503691/149‑1 (MQ=255)

                                                                                                                                                         ccAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAA‑CACCATCAGGCATATTAAATTATg                                                                                                                                                                                                                                                             >  8:12811/1‑149 (MQ=255)

                                                                                                                                                            gCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAA‑CACCATCAGGCATATTAAATTATGCtt                                                                                                                                                                                                                                                          >  6:524434/1‑149 (MQ=255)

                                                                                                                                                                   cAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAA‑CACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACGGCACGCGCTGCCAGACCTTTCAGTGGGCCATCAGc                                                                                                                     <  2:38832/275‑1 (MQ=255)

                                                                                                                                                                    aaaCTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGATGCGGGAATATGGGCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAA‑CACCATCAGGCATATTAAATTATGCTTTCAGTACa                                                                                                                                                                                                                                                  >  8:313154/1‑149 (MQ=255)

                                                                                                                                                                     accTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAA‑CACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACGGCACGCGCTGCCAGACCTTTCAGTGGGCCATCAGCaa                                                                                                                   <  2:143797/273‑1 (MQ=255)

                                                                                                                                                                      aCTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGTTG‑CACCATCAg                                                                                                                                                                                                                                                                            >  6:666391/1‑121 (MQ=255)

                                                                                                                                                                      aCTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGTTG‑CACCATCAg                                                                                                                                                                                                                                                                            <  5:666391/121‑1 (MQ=255)

                                                                                                                                                                        tGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAA‑CACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCa                                                                                                                                                                                                                                              <  5:577253/149‑1 (MQ=255)

                                                                                                                                                                               aCGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAA‑CACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAgctgc                                                                                                                                                                                                                                       >  7:49641/1‑149 (MQ=255)

                                                                                                                                                                                  tctcGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAA‑CACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTc                                                                                                                                                                                                                                    >  5:142534/1‑149 (MQ=255)

                                                                                                                                                                                   ctcGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAA‑CACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCg                                                                                                                                                                                                                                   <  5:1047606/149‑1 (MQ=255)

                                                                                                                                                                                    tcGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAA‑CACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGt                                                                                                                                                                                                                                  >  8:647999/1‑149 (MQ=255)

                                                                                                                                                                                      ggCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAA‑CACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTaa                                                                                                                                                                                                                                <  8:1158177/149‑1 (MQ=255)

                                                                                                                                                                                         gTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAA‑CACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATcc                                                                                                                                                                                                                             <  5:1534823/149‑1 (MQ=255)

                                                                                                                                                                                         gTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAA‑CACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATcc                                                                                                                                                                                                                             <  6:1027466/149‑1 (MQ=255)

                                                                                                                                                                                          tttCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAA‑CACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATcc                                                                                                                                                                                                                             >  5:501599/1‑148 (MQ=255)

                                                                                                                                                                                          tttCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAA‑CACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCg                                                                                                                                                                                                                            >  6:479201/1‑149 (MQ=255)

                                                                                                                                                                                          tttCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAA‑CACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTCCGTAATCCg                                                                                                                                                                                                                            >  7:1376885/1‑149 (MQ=255)

                                                                                                                                                                                                aTACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAA‑CACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCgg                                                                                                                                                                                                                      <  5:1603232/149‑1 (MQ=255)

                                                                                                                                                                                                 tACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAA‑CACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCggt                                                                                                                                                                                                                     <  7:313154/149‑1 (MQ=255)

                                                                                                                                                                                                 tACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAA‑CACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCggt                                                                                                                                                                                                                     <  5:685397/149‑1 (MQ=255)

                                                                                                                                                                                                      cgATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAA‑CACCATCAg                                                                                                                                                                                                                                                                            <  3:196671/89‑1 (MQ=255)

                                                                                                                                                                                                          agaacagaaGGTGCGGGAGAATGGTCTGAGTCAGAAGAAGAACGAGAGGTAAGCGCAAAAAGGCCTGATGGTGCAA‑CACCATCAGGCAGAGTAAATTATGTTTTCAGTACAGCCAGAGCGGCTTCGTAATCCGGCTCGGGGGTGATTTCATCCACCAGCTGGCTGAAAATAACATTGTCATTTTCGTCAATAACCACAACGGCACGCGCTGCAAGACCTTTCAGTGGGCCATCAGCAATTGCCACACCGTAAGCTTGCAGAAATTCAGCGTTAc                                                                               <  4:351455/267‑1 (MQ=255)

                                                                                                                                                                                                           gCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAA‑CACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCa                                                                                                                                                                                                           <  5:205939/149‑1 (MQ=255)

                                                                                                                                                                                                           gCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAA‑CACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCa                                                                                                                                                                                                           <  6:256173/149‑1 (MQ=255)

                                                                                                                                                                                                             cAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAA‑CACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATc                                                                                                                                                                                                         >  7:139965/1‑149 (MQ=255)

                                                                                                                                                                                                             cAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAA‑CACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATc                                                                                                                                                                                                         >  8:79253/1‑149 (MQ=255)

                                                                                                                                                                                                             cAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAA‑CACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATc                                                                                                                                                                                                         <  5:482669/149‑1 (MQ=255)

                                                                                                                                                                                                              aaTAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAA‑CACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATcc                                                                                                                                                                                                        >  8:1178215/1‑149 (MQ=255)

                                                                                                                                                                                                                   gTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAA‑CACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTc                                                                                                                                                                                                            >  5:588658/1‑140 (MQ=255)

                                                                                                                                                                                                                   gTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAA‑CACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACGGCACGCGCTGCCAGACCTTTCAGTGGGCCATCAGCAATTGCCACACCGTAAGCTTGCAGAAATTCAGCGTTACGGAAAGTGga                                                                     <  2:47843/275‑1 (MQ=255)

                                                                                                                                                                                                                      gCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGCGTAAGCGCGAATATGCCTGATGGTGCAA‑CACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACGGCACGCGCTGCCAGACCTTTCAGTGGGCCATCAGCATTTGCCACACCGTAAGCTTTCAGAATTTTAGCGTTACGGAAAGTGaag                                                                    >  2:351459/1‑273 (MQ=255)

                                                                                                                                                                                                                        gggAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAA‑CACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCgg                                                                                                                                                                                                                      >  7:704995/1‑125 (MQ=255)

                                                                                                                                                                                                                        gggAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAA‑CACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCgg                                                                                                                                                                                                                      <  8:704995/125‑1 (MQ=255)

                                                                                                                                                                                                                                    tGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAA‑CACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACAt                                                                                                                                                                                  >  5:1120339/1‑149 (MQ=255)

                                                                                                                                                                                                                                    tGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAA‑CACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACAt                                                                                                                                                                                  >  8:1431374/1‑149 (MQ=255)

                                                                                                                                                                                                                                     gAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAA‑CACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTTGTGATTTCATCCCCCAGCTGGCTGGAAATCACATTGTCATTTTCGTCAAAAACCACAACAGCACGCGCTGCCCGCCCTTTCAGTGGGCCATCAGCACTTGCCACACCTTAAGCTTACAGACATTCAGCGTTACGGAAAGTGGAGAGGGTGATAACGTTGTTAAg                                                >  1:186133/1‑278 (MQ=255)

                                                                                                                                                                                                                                       gTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAA‑CACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACgg                                                                                                                                                       >  4:199208/1‑173 (MQ=255)

                                                                                                                                                                                                                                        tCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAA‑CACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTc                                                                                                                                                                              <  8:139965/149‑1 (MQ=255)

                                                                                                                                                                                                                                        tCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAA‑CACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTc                                                                                                                                                                              >  7:1608607/1‑149 (MQ=255)

                                                                                                                                                                                                                                          agaggaagaACGAAGAGGAAGCGCGAAGAGGCCTGAGGGTACAA‑CACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACGGCATGCGCTGCCAGACCTTTCAGTGGGCCATCAGCAATTGCCACACCGTATGCTTGCAGAAATTCAGCGTTAc                                                                               <  2:110905/242‑1 (MQ=255)

                                                                                                                                                                                                                                               aagaaCGAAGAGTAAGCGCGAATATGCCTGATGGTGCAA‑CACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCg                                                                                                                                                                       >  6:696198/1‑149 (MQ=255)

                                                                                                                                                                                                                                                  aaCGAAGAGTAAGCGCGAATATGCCTGATGGTGCAA‑CACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCa                                                                                                                                                                    >  6:1548427/1‑149 (MQ=255)

                                                                                                                                                                                                                                                   aCGAAGAGTAAGCGCGAATATGCCTGATGGTGCAA‑CACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTGATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCaa                                                                                                                                                                   <  7:1461958/149‑1 (MQ=255)

                                                                                                                                                                                                                                                     gAAGAGTAAGCGCGAATATGCCTGATGGTGCAA‑CACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATa                                                                                                                                                                 >  5:1146463/1‑149 (MQ=255)

                                                                                                                                                                                                                                                      aaGAGTAAGCGCGAATATGCCTGATGGTGCAA‑CACCATCAGGCATATTAAATTATGTTTTCAGTACAGCCAGAGCTGCTTCGTACTCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATaa                                                                                                                                                                <  6:1663382/149‑1 (MQ=255)

                                                                                                                                                                                                                                                         agTAAGCGCGAATATGCCTGATGGTGCAA‑CACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACGGCACGCGCTGCCAGACCTTTCAGTGGGCCATCAGCAATTGCCACACCGTAAGCTTGCAGAAATTCAGCGTTACGGAAAGTGGAGAGGGTGATAACGTTGTTCAGACCTTCTGCGCCGCAgc                               >  2:184547/1‑274 (MQ=255)

                                                                                                                                                                                                                                                          gTAAGCGCGAATATGCCTGATGGAGCAA‑CACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACGGCACGCGCTGCCAGACCTTTCAGTGGGCCATCAGCAATTGCCACACCGTAAGCTTGCAGAAATTCAGCGTTACGGaa                                                                           <  2:421622/230‑1 (MQ=255)

                                                                                                                                                                                                                                                           tAAGCGCGAATATGCCTGATGGTGCAA‑CACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACcaca                                                                                                                                                           >  6:1589405/1‑149 (MQ=255)

                                                                                                                                                                                                                                                            aaGCGCGAATATGCCTGATGGTGCAA‑CACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACaa                                                                                                                                                          <  7:209376/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                   aaTATGCCTGATGGTGCAA‑CACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACGGCACGCGCTGCCAGACCTTTCAGTGGGCCATCAGCAATTGCCACACCGTAAGCTTGCAGAAATTCAGCGTTACGGAAAGTGGAGAGGGTGATAACGTTGTTCAGACCTTCTGCGCCGCAGAAACGAGACt                      >  4:640279/1‑274 (MQ=255)

                                                                                                                                                                                                                                                                    atatGCCTGATGGTGCAA‑CACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACGGCAcgc                                                                                                                                                  <  6:142534/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                      atGCCTGATGGTGCAA‑CACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACGGCACgcgc                                                                                                                                                <  7:79253/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                      atGCCTGATGGTGCAA‑CACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACGGCACgcgc                                                                                                                                                <  6:978277/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                        gCCTGATGGTGCAA‑CACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACGGCACGCGCTg                                                                                                                                              <  6:1135865/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                         ccTGATGGTGCAA‑CACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATaa                                                                                                                                                                >  6:577259/1‑130 (MQ=255)

                                                                                                                                                                                                                                                                          cTGATGGTGCAA‑CACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACGGCACGCGCTGcc                                                                                                                                            >  5:266279/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                           tGATGGTGCAA‑CACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACGGCACGCGCTGCCa                                                                                                                                           >  7:1540617/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                             aTGGTGCAA‑CACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAgctggc                                                                                                                                                                                              <  5:1060029/96‑1 (MQ=255)

                                                                                                                                                                                                                                                                             aTGGTGCAA‑CACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAgctggc                                                                                                                                                                                              >  6:1060029/1‑96 (MQ=255)

                                                                                                                                                                                                                                                                              tGGTGCAA‑CACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCTTTTTCGTCAATAACCACAACGGCACGCGCTGCCAGAc                                                                                                                                        >  8:614485/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                              tGGTGCAA‑CACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACGGCACGCGCTGCCAGAc                                                                                                                                        <  7:581061/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                                 tGCAA‑CACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACgg                                                                                                                                                       >  2:716010/1‑131 (MQ=255)

                                                                                                                                                                                                                                                                                 tGCAA‑CACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACGGCACGCGCTGCCAGACCtt                                                                                                                                     >  5:721881/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                                  gCAA‑CACCATCAGGCATATTAAATTAGGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAAGAACCACAACGGCACGCGCTGCCAGACCTTTCAGTGGGCCATCAGCAATTGCCACACCGTAAGCTTGCAGAAATTCAGCGTTACGGAAAGTGGAGAGGGTGATAACGTTGTTCAGACCTTCTGCGCCGCAGAAACGAGACTGGGCGAACGGCAGATCGGCa  <  1:72543/279‑1 (MQ=255)

                                                                                                                                                                                                                                                                                    aa‑CACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACGGCACGCGCTGCCAGACCTTTCa                                                                                                                                  >  8:1027436/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                                     a‑cacCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACGGCAcg                                                                                                                                                   <  4:716005/131‑1 (MQ=255)

TGGACGGTTATGAGCTGCGTCCACAGGATATTTTGTTGCCGGATTATGACGCCGCAACGGTAGCCGTGGGCGAAGATGCGATGGAAGGTTCGCTGGACGAAATCACCAGCCGTTTTGAACGCTCGGTATTAACCCAGCTTTATCGCAATTATCCCAGCACGCGCAAACTGGCAAAACGTCTCGGCGTTTCACATACCGCGATTGCCAATAAGTTGCGGGAATATGGTCTGAGTCAGAAGAAGAACGAAGAGTAAGCGCGAATATGCCTGATGGTGCAA‑CACCATCAGGCATATTAAATTATGCTTTCAGTACAGCCAGAGCTGCTTCGTAATCCGGCTCGGTGGTGATTTCATCCACCAGCTGGCTGAAAATCACATTGTCATTTTCGTCAATAACCACAACGGCACGCGCTGCCAGACCTTTCAGTGGGCCATCAGCAATTGCCACACCGTAAGCTTGCAGAAATTCAGCGTTACGGAAAGTGGAGAGGGTGATAACGTTGTTCAGACCTTCTGCGCCGCAGAAACGAGACTGGGCGAACGGCAGATCGGCA  >  NC_000913/1388008‑1388560

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