Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,411,135	0	A	.	67.7%	26.7 / 15.3	32	coding (583/699 nt)	yjfJ	putative transcriptional regulator effector protein
Reads supporting (aligned to +/- strand):  ref base A (6/4);  new base . (0/21);  total (6/26)
Fisher's exact test for biased strand distribution p-value = 2.85e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.31e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGG‑AAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCG‑TTT‑CGACAGCGGCAG‑AATCCTT‑AAAACGCCTGC‑AAACGCGTCAGGCCG‑AACGTCAGGCTCGCCTGGA‑TGCTGCCGCACAGTTGGAG‑AAAG‑TCGCAGACGGTCGCG‑ACCTTGACG‑AAAAGCT‑GGC‑GGAAGCCGG‑AATTGGC‑GGTAG‑CAAT‑AAAAGTAGCG‑‑CCCAGGATGTATTAGC‑AAGACTG‑CAACGCC‑AACAGGGCG‑AGTAATTTTTTTGCCAGCCTCGCT‑AAAAGGCT‑GGCAACTATTTTAAGG‑ATAAAATATGTCTGGTTTTTTCCAGCGTCTGTTTGGCAAGGATAATAAGCCAGCTATCGCTCGTGGTCCGCTGGGACTTCATCTCAATAGTGGCTTTACGCTCGACACG  >  NC_000913/4410859‑4411403                                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                               ctcATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGGAAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCG‑TTT‑CGACAGCGGCAG‑AATCCTTAAAAACGCCTGCAAAACGCGTCAGGCCG‑AACGTCAGGCTCGCCTGGA‑TGCTGCCGCACAGTTGGAG‑aagtcg                                                                                                                                                                                                                                                                                           <  1:312054/285‑5 (MQ=255)                  gCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGG‑‑AAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCG‑TTT‑CGACAGCGGCAG‑AATCCTT‑AAAACGCCTGCAAAACGCGTCAGGCCG‑AACGTCAGGCTCGCCTGGA‑TGCTGCCGCACAGTTGGAG‑‑AAG‑TCGCAGACGGTCGCG‑ACCTTGACG‑AAAAGCT‑GGC‑GGAAGCCgg                                                                                                                                                                                                                                             <  1:908514/305‑1 (MQ=255)                                           tcCCGTGATGGCGTGG‑AAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTT‑CAGCG‑TTT‑CGACAGCGGCAG‑AATCCTT‑AAAACGCCTGC‑AAACGCGTCAGGCCG‑AACGTCAGGCTCGCCTGGACTGCTGCCGCACAGTTGGAG‑aaa                                                                                                                                                                                                                                                                                              >  1:312058/1‑237 (MQ=255)                                             ccGTGATGGCGTGG‑AAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCG‑TTT‑CGACAGCGGCAG‑AATCCTT‑AAAACGCCTGC‑AAACGCGTCAGGCCG‑AACGTCAGGCTCGCCTGGA‑TGCTGCCGCACAGTTGGAG‑‑AAG‑TCGCAGACGGTCGCG‑ACCTTGACG‑aaa                                                                                                                                                                                                                                                               <  1:908519/262‑1 (MQ=255)                                              cGTGATGGCGTGG‑AAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCG‑TTT‑CGACAGCGGCAG‑AATCCTT‑AAAACGCCTGC‑AAACGCGTCAGGCCG‑AACGTCAGGCTCGCCTGGA‑TGCTGCCGCACAGTTGGAG‑‑AAG‑TCGCAg                                                                                                                                                                                                                                                                                      <  1:312057/240‑1 (MQ=255)                                                   tGGCGTGG‑AAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGC‑AATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCG‑TTT‑CGACAGCGGCAG‑AATCCTT‑AAAACGCCTGC‑AAACGCGTCAGGCCG‑AACGTCAGGCTCGCCTGGA‑TGCTGCCGCACAGTTGGAG‑‑AAG‑TCGCAGACGGTCGCG‑ACCTTGACG‑AAAAGCT‑GGC‑GGAAGCCGG‑AATTGGC‑GGTAG‑CAAT‑AAAAGTAGCG‑‑ccccag                                                                                                                                                                                                       <  1:908521/303‑4 (MQ=255)                                                                              cagcgcagcgTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCG‑TTT‑CGACAGCGGCAG‑AATCCTT‑AAAACGCCTGC‑AAACGCGTCAGGCCG‑AACGTCAGGCTCGCCTGGA‑TGCTGCCGCACAGTTGGAG‑‑AAG‑TCGCAGACGGTCGCG‑ACCTTGACG‑AAAAGCT‑GGC‑g                                                                                                                                                                                                                                                     <  1:908524/238‑1 (MQ=255)                                                                                gcgcagcgTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCG‑TTT‑CGACAGCGGCAG‑AATCCTTAAAAACGCCTGC‑AAACGCGTCAGGCCG‑AACGTCAGGCTCGCCTGGA‑TGCTGCCGCACAGTTGGAG‑‑AAG‑TCGCAGACGGTCGCG‑ACCTTGACG‑AAAAGCT‑GGC‑GGAAGCCGG‑AATTGGC‑GGTAG‑CAAT‑aa                                                                                                                                                                                                                       <  1:908525/263‑1 (MQ=255)                                                                                                             aaTGGAAGTCGTT‑AAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCG‑TTT‑CGACAGCGGCAG‑AATCCTT‑AAAACGCCTGC‑AAACGCGTCAGGCCG‑AACGTCAGGCTCGCCTGGA‑TGCTGCCGCACAGTTGGAG‑‑AAG‑TCGCAGACGGTCGCG‑ACCTTGACG‑AAAAGCT‑GGC‑GGAAGCCGG‑AATTGGC‑GGTAG‑CAAT‑AAAAGTAGCGCCCCCAGGATGTATTAGC‑AAg                                                                                                                                                                                         <  1:908527/261‑1 (MQ=255)                                                                                                                    gTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCG‑TTT‑CGACAGCGGCAG‑AATCCTT‑AAAACGCCTGCAAAACGCGTCAGGCCG‑AACGTCAGGCTCGCCTGGA‑TGCTGCCGCACAGTTGGAGAAAAG‑TCGCAGACGGTCGCG‑ACCTTGACG‑aaa                                                                                                                                                                                                                                                               <  1:908529/195‑1 (MQ=255)                                                                                                                    gTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTT‑GCGCATCTTCCAGCG‑TTT‑CGACAGCGGCAGAAATCCTT‑AAAACGCCTGC‑AAACGCGTCAGGCCG‑AACGTCAGGCTCGCCTGGA‑TGCTGCCGCACAGTTGGAG‑‑AAG‑TCGCAGACGGTCGCG‑ACCTTGACG‑AAAAGCT‑GGC‑GGAAGCCGG‑AATTGGC‑GGTAG‑CAAT‑AAAAGTAGCGC‑CCCAGGATGTATTAGC‑AAGACTG‑CAACGCC‑AACAGGGCG‑AGTAATTTTTTTGCCAGCCTCGCT‑AAAAGGCT‑GGCAACTATTTTAAGG‑ATAAAATATGTCTGGTTTTTTCCAGCGTCTGTTTGGCAAGGATAATAAGCCAGCTATCGCt                                                  <  1:908530/383‑1 (MQ=255)                                                                                                                                        cATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCG‑TTT‑CGACAGCGGCAG‑AATCCTT‑AAAACGCCTGC‑AAACGCGTCAGGCCG‑AACGTCAGGCTCGCCTGGA‑TGCTGCCGCACAGTTGGAG‑‑AAG‑TCGCAGACGGTCGCG‑ACCTTGACG‑AAAAGCT‑GGC‑GGAAGCCGG‑a                                                                                                                                                                                                                                           <  1:908533/189‑1 (MQ=255)                                                                                                                                              gCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCG‑TTT‑CGACAGCGGCAG‑AATCCTT‑AAAACGCCTGC‑AAACGCGTCAGGCCG‑AACGTCAGGCTCGCCTGGA‑‑GCTGCCGCACAGTTGGAG‑‑AAG‑TCGCAGACGGTCGCG‑ACCTTGACg                                                                                                                                                                                                                                                                   <  1:312063/162‑1 (MQ=255)                                                                                                                                                    acaacaGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCG‑TTT‑CGACAGCGGCAGAAATCCTT‑AAAACGCCTGCAAAACGCGTCAGGCCG‑AACGTCAGGCTCGCCTGGA‑TGCTGCCGCACAGTTGGAG‑‑AAG‑TCGCAGACGGTCGCG‑ACCTTGACG‑AAAAGCT‑GGC‑GGAAGCCGG‑AATTGGCAGGTAG‑CAAT‑AAAAGTAGCG‑‑CCCAGGATg                                                                                                                                                                                                    <  1:312064/214‑1 (MQ=255)                                                                                                                                                               tAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTT‑CGACAGCGGCAG‑AATCCTT‑AAAACGCCTGC‑AAACGCGTCAGGCCG‑AACGTCAGGCTCGCCTGGA‑TGCTGCCGCACAGTTGGAG‑AAAGATCGCAGAC‑GTCGCG‑ACCTTGACG‑AAAAGCT‑GGC‑GGAAGCCGG‑AATTGGC‑GGTAG‑CAAT‑a                                                                                                                                                                                                                        >  1:908535/1‑184 (MQ=255)                                                                                                                                                                        ctACCGTTGGCGCATCTTCCAGCG‑TTT‑CGACAGC‑GCAGAAATCCTT‑AAAACGCCTGC‑AAACGCGTCAGGCCG‑AACGTCAGGCTCGCCTGGA‑TGCTGCCGCACAGTTGGAG‑‑AAG‑TCGCAGA‑GGTCGCG‑ACCTTGACG‑AAAAGCT‑‑GC‑GGAAGCCGG‑AATTGGC‑GGTAG‑CAAT‑AAAAGTAGCG‑‑CCCAGGATGTATTAGC‑AAGACTG‑CAACGCC‑AACAGGGCG‑AGTAATTTTTTTGCCAGCCTCGCT‑AAAAGGCT‑GGCAACTATTTTAAGG‑ATaaaa                                                                                                         <  1:312065/273‑1 (MQ=255)                                                                                                                                                                           ccGTTGGCGCA‑CTTCCAGCG‑TTTACGACAGCGGCAG‑AATCCTTAAAAACGCCTGCAAAACGCGTCAGGCCG‑AACGTCAGGCTCGCCTGGA‑TGCTGCCGCACAGTTGGAG‑aagtcg                                                                                                                                                                                                                                                                                           <  1:1220391/115‑5 (MQ=255)                                                                                                                                                                                        ttCCAGCGTTTT‑CGACAGCGGCAG‑AATCCTT‑AAAACGCCTGC‑AAACGCGTCAGGCCGAAACGTCAGGCTCGCCTGGA‑TGCTGCCGCACAGTTGGAG‑AAAGATCGCAGAC‑GTCGCG‑ACCTTGACG‑AAAAGCT‑GGC‑GGAAGCCGG‑AATTGGC‑GGTAG‑CAAT‑AAAAGTAGCG‑‑CCCAGGATGTATTAGC‑AAGACTG‑CAACGCC‑aa                                                                                                                                                                          >  1:312068/1‑201 (MQ=255)                                                                                                                                                                                        ttCCAGCGTTTT‑CGACAGCGGCAG‑AATCCTT‑AAAACGCCTGC‑AAACGCGTCAGGCCG‑AACGTCAGGCTCGCCTGGA‑TGCTGCCGCACAGTTGGAG‑AAAGATCGCAGAC‑GTCGCG‑ACCTTGACG‑AAAAGCT‑GGC‑GGAAGCCGG‑AATTGGC‑GGTAG‑CAAT‑AAAAGTAGCG‑‑CCCAGGATGTATTAGC‑AAGACTG‑CAACGCC‑AACAGGGCG‑AGTAA‑TTTTTTG‑CAGCCTCGCT‑AAAAGGCT‑GGCAACTATTTTAAGGAATAAAATATGTCTGGTTTTTTCCAGCGTCt                                                                                 >  1:312066/1‑284 (MQ=255)                                                                                                                                                                                        ttCCAGCG‑TTT‑CGACAGCGGCAG‑AATCCTT‑AAAACGCCTGC‑AAACGCGTCAGGCCG‑AACGTCAGGCTCGCCTGGA‑TGCTGCCGCACAGTTGGAG‑AAAGATCGCAGAC‑GTCGCG‑ACCTTGACG‑AAAAGCT‑GGC‑GGAAGCCGG‑AATTGGC‑GGTAG‑CAATAAAAAGTAGCG‑‑CCCAGGATGTATTAGC‑AAGACTG‑CAACGCC‑AACAGGGCG‑AGTAA‑TTTTTTG‑CAGCCTCGCT‑AAAAGGCT‑GGCAACTA‑TTTAAGG‑At                                                                                                             >  1:312067/1‑254 (MQ=255)                                                                                                                                                                                                        cAGCGGCAG‑AATCCTT‑AAAACGCCTGCAAAACGCGTCAGGCCG‑AACGTCAGGCTCGCCTGGA‑TGCTGCCGCACAGTTGGAG‑‑AAG‑TCGCAGACGGTCGCG‑ACCTTGACG‑AAAAGCT‑GGC‑GGAAGCCGG‑AATTGGC‑GGTAG‑CAAT‑AAAAGTAGCG‑‑CCCAGGATGTATTAGC‑AAGACTG‑CAACGCC‑AACAGGGCG‑AGTAATTTTTTTGCCAGCCTCGCT‑AAAAGGCT‑GGCAACTATTTTAAGG‑ATAAAATATGTCTGGTTTTTTCCAGCGTCTGTTTGGCAAGGat                                                                    <  1:908536/283‑1 (MQ=255)                                                                                                                                                                                                          gCGGCAG‑AATCCTTAAAAACGCCTGCAAAACGCGTCAGGCCG‑AACGTCAGGCTCGCCTGGA‑TGCTGCCGCACAGTTGGAG‑‑AAG‑TCGCAGACGGTCGCG‑ACCTTGACG‑AAAAGCT‑GGC‑GGAAGCCGG‑AATTGGC‑GGTAG‑CAAT‑AAAAGTAGCG‑‑CCCAGGATGTATTAGC‑AAGACTG‑CAACGCC‑AACAGGGCG‑AGTAATTTTTTTGCCAGCCTCGCT‑AAAAGGCTGGGCAACTATTTTa                                                                                                                   <  1:312069/237‑1 (MQ=255)                                                                                                                                                                                                              cAG‑AATCCTT‑AAAACGCCTGC‑AAACGCGTCAGGCCG‑AACGTCAGGCTCGCCTGGA‑TGCTGCCGCACAGTTGGAG‑‑AAG‑TCGCAGACGGTCGCG‑ACCTTGACG‑AAAAGCT‑GGC‑GGAAGCCGG‑AATTGGC‑GGTAG‑CAAT‑AAAAGTAGCG‑‑CCCAGGATGTATTAGC‑AAGACTG‑CAACGCC‑AACAGGGCG‑AGTAATTTTTTTGCCAGCCTCGCT‑AAAAGGCT‑GGCAACTATTTTAAGG‑ATaa                                                                                                           <  1:908537/237‑1 (MQ=255)                                                                                                                                                                                                              cAG‑AATCCTT‑AAAACGCCTGC‑AAACGCGTCAGGCCG‑AACGTCAGGCTCG‑CTGGA‑TGCTGCCGCACAGTTGGAG‑‑AAG‑TCGCAGACGGTCGCG‑ACCTTGACG‑AAAAGCT‑GGC‑GGAAGCCGG‑AATTGGC‑GGTAG‑CAAT‑AAAAGTAGCGC‑CCCAGGATGTATTAGC‑AAGACTGTCAACGCC‑AACAGGGCGAAGTAATTTTTTTGCCAGCCTCGCT‑AAAAGGCT‑GGCAACTATTTTAAGG‑ATAAAATATGTCTGGTTTTTTCCAGCGTCTGTTTGGCAAGGATAATAAGCCAGCTATCGc                                                   <  1:908538/295‑1 (MQ=255)                                                                                                                                                                                                                       ttaaaaaCGCCTGCAAAACGCGTCAGGCCG‑AACGTCAGGCTCGCCTGGA‑TGCTGCCGCACAGTTGGAG‑‑AAG‑TCGCAGA‑GGTCGCG‑ACCTTGACG‑AAAAGCT‑GGC‑GGAAGCCGG‑AATTGGC‑GGTAG‑CAAT‑AAAAGTAGCG‑‑CCCAGGATGTATTAGC‑AAGACTG‑CAACGCC‑AACAGGGCG‑AGTAATTTTTTTGCCAGCCTCGCT‑AAAAGGCT‑GGCAACTATTTTAAGG‑ATa                                                                                                            <  1:312071/226‑1 (MQ=255)                                                                                                                                                                                                                                 ctgcaaaaCGCGTCAGGCCG‑AACGTCAGGCTCGCCTGGA‑TGCTGCCGCACAGTTGGAG‑AAAG‑TCGCAGACGGTCGCG‑ACCTTGACG‑AAAAGCT‑GGC‑GGAAGCCGG‑AATTGGC‑gg                                                                                                                                                                                                                                  <  1:312072/109‑1 (MQ=255)                                                                                                                                                                                                                                      aaaCGCGTCAGGCCG‑AACGTCAGGCTCGCCTGGA‑TGCTGCCGCACAGTTGGAG‑AAAG‑TCGCAGACGGTCGCG‑ACCTTGACG‑AAAAGCT‑GGC‑GGAAGCCGG‑AATTGGC‑GGTAG‑CAAT‑AAAAGTAGCGC‑CCCAGGATGTATTAGC‑AAGACTG‑c                                                                                                                                                                                   <  1:908540/151‑1 (MQ=255)                                                                                                                                                                                                                                            gTCAGGCCG‑AACGTCAGGCTCGCCTGGA‑TGCTGCCGCACAGTTGGAG‑AAAG‑TCGCAGA‑GGTCGCG‑ACCTTGACGAAAAAGCT‑GGCGGGAAGCCGG‑AATTGGC‑GGTAGAGAAT‑AAAAGTAGCG‑‑CCCAGGATGTATTAGC‑AAGACTG‑CAACGCCAAACAGGGCG‑AGTAAttt                                                                                                                                                          <  1:961430/170‑1 (MQ=255)                                                                                                                                                                                                                                                 gCCG‑AACGTCAGGCTCGCCTGGA‑TGCTGCCGCACAGTTGGAG‑AAAG‑TCGCAG‑C‑GTCGCG‑ACCTTGACG‑AAAAGCT‑GGC‑GGAAGCCGG‑AATTGGC‑GGTAG‑CAAT‑AAAAGTAGCG‑‑‑CCAGGATGTATTAGC‑AAGACTG‑CAACGCC‑AACAGGGCG‑AGTAA‑TTTTTTGCCAGCCTCGCT‑AAAAGGCT‑GGCAACTAt                                                                                                                       >  1:908547/1‑191 (MQ=255)                                                                                                                                                                                                                                                  ccG‑AACGTCAGGCTCGCCTGGA‑TGCTGCCGCACAGTTGGAG‑‑AAG‑TCGCAGA‑GGTCGCGAACCTTGACG‑AAAAGCTAGGC‑GGAAGCCGGAAATTGGC‑GGTAG‑CAAT‑AAAAGTAGCG‑‑CCCAGGATGTATTAGCAAAGACTG‑CAACGCCAAACAGGGCGAAGTAATTTTTTTGCCAGCCTCGCT‑aa                                                                                                                                       <  1:1154112/183‑1 (MQ=255)                                                                                                                                                                                                                                                       aCGTCAGGCTCGCCTGGA‑TGCTGCCGCACAGTTGGAG‑‑AAG‑TCGCAGACGGTCGCG‑ACCTTGACG‑AAAAGCT‑GGC‑GGAAGCCGG‑AATTGGC‑GGTAG‑CAAT‑AAAAGTAGCGC‑CCCAGGATg                                                                                                                                                                                                    <  1:908546/119‑1 (MQ=255)                                                                                                                                                                                                                                                           cAGGCTCGCCTGGA‑TGCTGCCGCACAGTTGGAG‑‑AAG‑TCGCAGA‑GGTCGCG‑ACCTTGACG‑AAAAGCT‑GGC‑GGAAGCCGG‑AATTGGC‑GGTAG‑CAAT‑AAAAGTAGCG‑‑CCCAGGATGTATTAGC‑AAGACTG‑CAACGCC‑AACAGGGCG‑AGTAATTTTTTTGCCAGCCTCGCTAAAAAGGCT‑‑GCAACTATTTTAAGG‑ATAAAATATGTCTGGTTTTTTCCAGCGTCTGTTTGGCAAGGATAATAAGCCAGCTATCGCTCGTGGTCCGCTGGGACTTCATCTCAATAGTGGCTTTACGCTCGACACg  <  1:1136374/300‑1 (MQ=255)                                                                                                                                                                                                                                                                   ccTGGA‑TGCTGCCGCACAGTTGGAA‑GAAG‑TCGCAGACGGTCGCG‑A‑CTTGACG‑AAAAGCT‑‑GC‑GGAAG‑CGG‑AATTGGC‑GGTAG‑CAAT‑AAAAGTAGCG‑‑CCCAGGATGTATTAGC‑aa                                                                                                                                                                                          <  1:1256457/113‑1 (MQ=255)                                                                                                                                                                                                                                                                    cTGGA‑TGCTGCCGCACAGTTGGAG‑‑AAG‑TCGCAGACGGTCGCG‑ACCTTGACG‑AAAAGCT‑GGC‑GGAAGCCGG‑AATTGGC‑GGTAG‑CAAT‑AAAAGTAGCG‑‑CCCAGGATGTATTAGC‑AAGACTG‑CAACGCC‑AACAGGGCG‑AGTAATTTTTTTGCCAGCCTCGCT‑AAAAGGCT‑GGCAACTATTTTAAGG‑ATaaaa                                                                                                         <  1:908548/189‑1 (MQ=255)                                                                                                                                                                                                                                                                     tGGA‑TGCTGCCGCACAGTTGGAG‑AAAG‑TCGCAGAC‑GTCGCG‑ACCTTGACG‑AAAAGCT‑GGC‑GGAAGCCGG‑AATTGGC‑GGTAG‑C‑AT‑AAAAGt                                                                                                                                                                                                                   >  1:312075/1‑90 (MQ=255)                                                                                                                                                                                                                                                                                            g‑aAAG‑TCGCAGACGGTCGCG‑ACCTTGACG‑AAAAGCT‑GGC‑GGAAGCCGG‑AATTGGC‑GGTAG‑CAAT‑AAAAGTAGCG‑‑CCCAGGATGTATTAGC‑AAGACTG‑CAACGCC‑AACAGGGCG‑AGTAA‑TTTTTTGCCAGCCTCGCT‑AAAAGGCT‑GGCAACTATTTTAAGG‑ATAAAATATGTCTGGTTTTTTCCAGCGTCTGTTTGGCAAGGATAATAAGCCAGCTATc                                                     >  1:908550/1‑218 (MQ=255)                                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                               CTCATTACCGCTGAAGAGCAAGTGTTGTCGAATCTGGAAGTCTCCCGTGATGGCGTGG‑AAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCG‑TTT‑CGACAGCGGCAG‑AATCCTT‑AAAACGCCTGC‑AAACGCGTCAGGCCG‑AACGTCAGGCTCGCCTGGA‑TGCTGCCGCACAGTTGGAG‑AAAG‑TCGCAGACGGTCGCG‑ACCTTGACG‑AAAAGCT‑GGC‑GGAAGCCGG‑AATTGGC‑GGTAG‑CAAT‑AAAAGTAGCG‑‑CCCAGGATGTATTAGC‑AAGACTG‑CAACGCC‑AACAGGGCG‑AGTAATTTTTTTGCCAGCCTCGCT‑AAAAGGCT‑GGCAACTATTTTAAGG‑ATAAAATATGTCTGGTTTTTTCCAGCGTCTGTTTGGCAAGGATAATAAGCCAGCTATCGCTCGTGGTCCGCTGGGACTTCATCTCAATAGTGGCTTTACGCTCGACACG  >  NC_000913/4410859‑4411403

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑