BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	445,570	T→C	21.8%	D366G (GAC→GGC)	yajR ←	putative transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	445,570	0	T	C	21.8%	107.6 / 5.6	57	D366G (GAC→GGC)	yajR	putative transporter
Reads supporting (aligned to +/- strand): ref base T (25/17); new base C (4/9); total (29/28)
Fisher's exact test for biased strand distribution p-value = 1.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGCGATTCGTCACTTTGCTGTCGATTTTCACATAAGCTGAATGTTCTTCTTCTGCAATCAACACTTCTTTGATGCCTTCAGTTTCTAGCAAACGCACTTTTAACG‑CCTCGTTTGCGGCAATGTTCGCCGGGATTT‑C‑AATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAG‑CCACAC‑TGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCT‑GA‑CCG‑T‑CAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAA‑GAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGATAAGTGAGGGCAGGAGGGCTTCCATCAAATTAAACGCCACAAAGAAAAGCTGCACGCCGACCACCAGTTGCCAGAACTGCGTTTGCGCGTTCCACAACACAATTTCCGCAACCACGATCA  >  NC_000913/445333‑445807

aaGCGATTCGTCACTTTGCTGTCGATTTTCACATAAGCTGAATGTTCTTCTTCTGCAATCAACACTTCTTTGATGCCTTCAGTTTCTAGCAAACGCACTTTTAACG‑CCTCGTTTGCGGCAATGTTCGCCGGGATTT‑C‑AATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGGACTGGCGAGCGTGAG‑CCACACACACGGCCGCCATT‑CGCCAGCGAGGAGAACCCCCT‑GG‑CCG‑T‑c                                                                                                                                                                                                                                               >  2:646280/1‑239 (MQ=255)

     aTTCGTCACTTTGCTGTCGCTTTTCACATAAGCTGAATGTTCTTCTTCTGCAATCAACACTTCTTTGATGCCTTCAGTTTCTAGCCAACGCACTTTTAACG‑CCTCGTTTGCGGCAATGTTCGCCGGGATTT‑C‑AATGCGCAATCTGCTGACATACGGCGGTTCTTTCATGGGACCGGGGTCTTTTTA‑CCCCAC‑TGCGGGCCGAATTTCGCCCGCGGGAAAAACCCCCT‑GC‑CCG‑C‑CCAACa                                                                                                                                                                                                                                          >  1:1127004/1‑239 (MQ=255)

     aTTCGTCACTTTGCTGTCGATTTTCACATAAGCTGAATGTTCTTCTTCTGCAATCAACACTTCTTTGATGCCTTCAGTTTCTAGCAAACGCACTTTTAACG‑CCTCGTTTGCGGCAATGTTCGCCGGGATTT‑C‑AATGCGCAAACTGCTGACATACGGCGGTTCTTTCATAGGACTTGAGACGGTCAG‑CCACAC‑TGCGGCCAGCAATGCGCCAGCGAGAAATAACCCCC‑GA‑CCG‑C‑CAAACa                                                                                                                                                                                                                                          >  1:1038170/1‑239 (MQ=255)

      ttCGTCACTTTGCTGTCGATTTTCACATAAGCTGAATGTTCTTCTTCTGCAATCAACCCTGCTTTGATGCCTTCAGTTTCTAGCAAACGCACTTTTAACG‑CCGCGTTTGCGGCAATGTTCGCCGGGATTT‑C‑AATGCGCAAACTGCTGACATACGGCGGGTCTTTCATGGTACTGGCGACTGTCAG‑CCACAC‑TGCGGCCAGCATTGCGCCGGCGAGAAATACCCCCT‑GA‑CCG‑T‑CAAACa                                                                                                                                                                                                                                          <  2:795816/238‑1 (MQ=255)

      ttCGTCACTTTGCTGTCGATTTTCACATAAGCTGAATGTTCTTCTTCTGCAATCAACACTTCTTTGATGCCTTCAGTTTCTAGCAAACGCACTTTTAACG‑CCTCGTTTGCGGCAATGTTCGCCGGGATTT‑C‑AATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAG‑CCACAC‑TGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCT‑GA‑CCG‑T‑CAAACAt                                                                                                                                                                                                                                         <  1:1041951/239‑1 (MQ=255)

                tGCTGTCGATTTTCACATAAGCTGAATGTTCTTCTTCTGCAATCAACACTTCTTTGATGCCTTCAGTTTCTAGCAAACGCACTTTTAACG‑CCTCGTTTGCGGCAATGTTCGCCGGGATTT‑C‑AATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAG‑CCACAC‑TGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCT‑GA‑CCG‑T‑CAAACATGCCGTTAATc                                                                                                                                                                                                                               <  1:1127040/239‑1 (MQ=255)

                        aTTTTCACATAAGCTGAATGTTCTTCTTCTGCAATCAACACTTCTTTGATGCCTTCAGTTTCTAGCAAACGCACTTTTAACG‑CCTCGTTTGCGGCAATGTTCGCCGGGATTT‑C‑AATGCGCAAACTGCTGACATACGGCGGTTCTTTCCTGGTACTTGCGGCCGTCAG‑CCCCAC‑CCCGGCCAGCCTTTCGCCAGCGGGAAATAACCCCC‑GA‑CCG‑T‑CCAAAAACCCGTTTATCCAgcccc                                                                                                                                                                                                                        >  1:105874/1‑236 (MQ=255)

                                              cttcttctGCAATCAACACTTCTTTGATGCCTTCAGTTTCTAGCAAACGCACTTTTAACG‑CCTCGTTTGCGGCAATGTTCGCCGGGATTT‑C‑AATGCGCAAACTGCTGACATACGGCGGTTCTTTCCTGGTACTGGCGACTGTCAG‑CCACAC‑TTCGGCCCGCCTTTCGCCAGCGAGAAATACCCCCC‑GA‑CCC‑C‑CCACACTTCCGCTTATCCACCCGCCCCGCGaaacgccaaacgcccc                                                                                                                                                                                                  >  2:426537/1‑236 (MQ=255)

                                               ttcttctGCAATCAACACTTCTTTGATGCCTTCAGTTTCTAGCAAACGCACTTTTAACG‑CCTCGTTTGCGGCAATGTTCGCCGGGATTT‑C‑AATGCGCAAACTGCTGACATACGGCGGTTCTTTCTTGGTTCTTGCGACGGTCAG‑CCCCAC‑TTCGGCCAGCCTTTCCCCAGCGAGAAATACCCCCC‑GA‑CCG‑T‑CCAACATCCCGTTAATTCACACCCCCCGCGCACCGCCAATcgccccc                                                                                                                                                                                                 >  2:240386/1‑235 (MQ=255)

                                                 cttctGCAATCAACACTTCTTTGATGCCTTCAGTTTCTAGCAAACGCACTTTTAACG‑CCTCGTTTGCGGCAATGTTCGCCGGGATTT‑C‑AATGCGCAAACTGCTGACATACGGCGTTTCTTTCATGGTACTGGCGACTGTCAG‑CCACAC‑TGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCT‑GA‑CCG‑T‑CAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATcgccacgcca                                                                                                                                                                                              <  1:952231/239‑1 (MQ=255)

                                                 cttccGCAATAAACACTCTTTGGAGCCCTTCAGTTTCTGGCAAACACCCTTTAACCG‑CCCTGTTTGCGGCAATGTCCCCCGGGATTT‑C‑AATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTCCTGGCGCCTGCCGC‑CCACAC‑TGCGGCAAGCATTGCCCCAGCGAGAAATACCCCCT‑GA‑CCG‑T‑CAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATcgccacgcca                                                                                                                                                                                              <  2:860855/239‑1 (MQ=255)

                                                       cAATCAACACTTCTTTGATGCCTTCAGTTTCTAGCAAACGCACTTTTAACG‑CCTCGTTTGCGGCAATGTTCGCCGGGATTT‑C‑AATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGGGACTGTCAG‑CCACCC‑TGCGGCCAGCATTTCGCCAGCGAGAAAAACCCCCC‑GA‑CCG‑T‑CAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCaa‑aaaaa                                                                                                                                                                                       >  1:466255/1‑238 (MQ=255)

                                                           cAACCCTCCTTTGATGCCTTCATTTTCTAGCAACCGCACTTTTAACG‑CCTCGTTTGCGGCAATGTTCGCCGGGATTT‑C‑AATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAG‑CCACAC‑TGCGGCCAGCATTGCGCCAGCGAGAAATCCCCCCT‑GA‑CCG‑T‑CAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAA‑GAAActggc                                                                                                                                                                                   <  1:459165/239‑1 (MQ=255)

                                                                   cttTGATGCCTTCAGTTTCTAGCAAACGCACTTTTAACG‑CCTCGTTTGCGGCAATGTTCGCCGGGATTT‑C‑AATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAG‑CCACAC‑TGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCT‑GA‑CCG‑T‑CAAACa                                                                                                                                                                                                                                          <  1:1072176/177‑1 (MQ=255)

                                                                    tttGATGCCTTCAGTTTCTAGCAAACGCACTTTTAACG‑CCTCGTTTGCGGCAATGTTCGCCGGGATTT‑C‑AATGCGCAAACTGCTGACATACGGCGGTTTTTTTATGGGACTTGCGGCTGTCCA‑GCACAC‑AGCCGCCCGCATTTGGCCAGCGAGTAATACCCCCC‑GA‑CCGCC‑CAAACA‑GCCGTTAATCCAGGCGCCCCGCGGAACGCCCAACGCCAAGACAACAAAACGGGCgggtgggg                                                                                                                                                                           >  2:921691/1‑237 (MQ=255)

                                                                            tttCTGTTTCTAACAAACCCACTTTTACCG‑CCTGTTTGGCGCAAATGTT‑GGCGGGATTT‑CAAATGCGCAACCGGCGGACACCGGGCGTTTATTGCATGGAACGGGCGACTGTCGC‑CCACAC‑GGCGGACAGCATGCCGCCAGCGAAAAATCCCCCCT‑GA‑CCG‑T‑AAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAA‑GAAACTGGCTGGTGGAGTAAACAccc                                                                                                                                                                  <  2:212411/238‑1 (MQ=255)

                                                                             ttCAGTTTTAAGAAAACGCACTTTAACCG‑CCTCGTTTGCGGCAATGTTCGCCGGGATTT‑C‑AATGCGCAAACTGCTGACATACGGCGGTTCTTCCATGGTACGGGCGACTGTCAG‑CCACAC‑GGCGGCCAGCATTGCGCCAGCGAAAAATCCCCCCG‑GC‑CCG‑T‑AAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAA‑GAAACTGGCTGGTGGAGTAAACAcc                                                                                                                                                                   <  1:1092657/237‑1 (MQ=255)

                                                                             ttCAGTTTCTAGCAAACGCACTTTTAACG‑CCTCGTTTGCGGCAATGTTCGCCGGGATTT‑C‑AATGCGCAAACTGCTGACATACGGCGGTTCTTTCCTGGGTCTTGCGGCTGTCAG‑CCACCC‑CTCCGCCCGCCTTGCGCCCGCGAGAAATAACCCCC‑CA‑CCC‑T‑CCAACATGCCGTTCATTCCACCGCCCCCCGCACCCCCAACCGCCACCCCAC‑ACAACTAGCCGGGGGGGGAAACAccc                                                                                                                                                                  >  2:751384/1‑238 (MQ=255)

                                                                               cAGTTTCCAGCAAAACCCCGTTTGACGTCGTCGTTGGCGGAATTGTCCGCGGGGTTTC‑A‑AATGCGAAAACGGCGGACATCGGGGGTTTCTTTCATGGTACTGGCGCCGTTCGG‑CCACAC‑GGCGGCCAGCATGGCGCCAGCGAAAAATCCCCCCG‑GC‑CCG‑T‑CAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAA‑GAAACTGGCTGGTGGAGTAAACACCCAt                                                                                                                                                                <  2:53357/239‑1 (MQ=255)

                                                                                    ccTAGAAAACCCCCTTATACCG‑CCCCGTTGGCGGAAATGTACGCCGGGATTT‑A‑AAGGCGCAACCTGCTGACATACGGCGTTTCTTTCATGGTCCTGCCAACTGTCGC‑CCACAC‑TGCGGCCAGGTTTCCGCCAGCGAGAAATCCCCCCG‑GA‑CCG‑T‑AAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAA‑GAAACTGGCTGGTGGAGTAAACACCCATCGCCGt                                                                                                                                                          <  2:603659/238‑1 (MQ=255)

                                                                                       aGCAAACGCACTTTTAACG‑CCTCGTTTGCGGCAATGTTCGCCGGGATTT‑C‑AATGCGCAAACTGCTGACATACGGCGGTTTTTTTATGGTACTGGCGACTGTCAG‑CCACAC‑TGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCT‑GA‑CCG‑T‑CCAACATGACGTTAATCCAGCCGCCCCGCGAAACGCCAATCGCCCCGCCCA‑GAAACTGGCAGTGGGGGTAAACAccccaccccccgc                                                                                                                                                        >  2:1092657/1‑228 (MQ=255)

                                                                                                                tttGCGACAATGTCCCCCGGGATTC‑A‑AAGGCGCAAACTCCGGACTAACGGGGTTTTTTTCATGGT‑CTTGCGGCCTTCCGCCCACAC‑GGCGGCCAGCATGGCGCCAGCGAAAAATCCCCCCT‑GC‑CCG‑C‑CAACCATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAA‑GAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCttt                                                                                                                               <  2:375456/239‑1 (MQ=255)

                                                                                                                                cccggATTTAC‑AATGCGCAACCTGGTGACATTGGGCGGTTCTTTGATGGTACGGGCGACTGCCAG‑CCCCAC‑TGCGGCCAGCATGGCGCCAGCGAGAAATCCCCCGT‑GA‑CCG‑T‑CAAACATGCGGTTAATCCAGCCGCCCGGCGAACCGCCAATCGCCACGCCAA‑GAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGATAAGTGAggg                                                                                                                <  2:975064/236‑1 (MQ=255)

                                                                                                                                    gATTT‑C‑AATCCGCAAACTGCTGACATACGGCGGTTCTTTCATGTAACTGGCGACTGTCAC‑CCACCC‑GGCGGCCAGCATTGCGCCAGCAAAAAATCCCCCCT‑GA‑CCG‑T‑AAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAA‑GAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGATAAGTGAGGGCagg                                                                                                            <  1:940225/238‑1 (MQ=255)

                                                                                                                                    gATTC‑C‑AATGCGAAACCTGCTGACATCCGGGGGTTCTTCCAGGGTACTGGCGACTGTCAG‑CCACCC‑TGCGGCCAGCATTGCGCCAGCGAAAAATCCCCCCG‑GA‑CCG‑T‑CAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAA‑GAAActggctgg                                                                                                                                                                                <  1:748413/170‑1 (MQ=255)

                                                                                                                                          c‑AATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGCCAG‑CCACAC‑TGCGGCCAGCATTGCGCCAGGGAAAAAT‑CCCCCC‑GA‑CCC‑TCCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAA‑GAAACTGGCTGGTGGAGTAAACACCCATCGCCg                                                                                                                                                           <  1:1065014/186‑1 (MQ=255)

                                                                                                                                                      tggctgACACGCGGGGGTTTTTTCATGGTACGGGAGCCTGTCAG‑CCACAC‑GCCGGCAATCATTGCGCCAGCGAAAAATCCCCCCT‑GC‑CCT‑C‑CAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAA‑GAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGATAAGTGAGGGCAGGAGGGCTTCCATCAAAtt                                                                                           <  2:466255/237‑1 (MQ=255)

                                                                                                                                                          tgACATACGGGGGTTCTTCAATGGTCCTGGCGATGGTCCA‑CAAAAC‑TGCGCCCAG‑ATTACCCCCGGCAGAAAATCCCCCTGGA‑CGG‑A‑CACACATGCCGTTAACCAAGCCGCCAAGGAAGCCGCCAAGCGCCACGCCAA‑GAAACTGGTTGGTGGATTAAACACCCAACGCCGTACCATTGTAACCTGCTGGCGACTCTTTACTGATAAGTGAGGGCAGGAGGGCTTCCATCAAATTaaa                                                                                        <  2:490289/238‑1 (MQ=255)

                                                                                                                                                               tACGGCGTTTCTTTCATGGTACTGGCGACTGTCAG‑CCACCC‑GGGGGCCAGCATGGCGCCAGCGAAAAATCCCCCCT‑GC‑CCG‑C‑AAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAA‑GAAACTGGCTGATGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGATAAGTGAGGGCAGGAGGGCTTCCATCaaa                                                                                             <  1:462597/228‑1 (MQ=255)

                                                                                                                                                                      cttcttcCAGGGTGCGGGCGACGGCCAG‑CCAAAC‑TGCGGGCCGCATTGCGCCAGTGAAAAACCCCCCCT‑GC‑CCG‑C‑AAAACATGCCGTTAACCCAGCCGCCCAGCGACCCGCCAATCGCCACGCCAA‑GAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGATAAGTGAGGGCAGGAGGGCTTCCATCAAATTAAACGCCACAAAGaa                                                                            <  2:566488/237‑1 (MQ=255)

                                                                                                                                                                                 gTACTGGCGACTGTCAG‑CCACAC‑TGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCT‑GA‑CCG‑T‑CAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAA‑GAAACTGGCTGGTGGAGTAAacac                                                                                                                                                                    >  1:128238/1‑139 (MQ=255)

                                                                                                                                                                                 gTACTGGCGACTGTCAG‑CCACAC‑TGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCT‑GA‑CCG‑T‑CAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAA‑GAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCTACTCTTTACTACTAACTGAGGGCAGGCGGGCTTCCATCAAATTACACGCCACAACAACAAGTCGCACGcc                                                                >  2:294987/1‑239 (MQ=255)

                                                                                                                                                                                    cTGGCGACTGTCAG‑CCACAC‑TGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCT‑GA‑CCG‑T‑CAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAA‑GAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGATAAGTGAGGGCAGGAGGGCTTCCATCAAATTAAACGCCACAAAGAAAAGCTGCACGCCGa                                                              >  1:191104/1‑238 (MQ=255)

                                                                                                                                                                                     tGGCGACTGTCAG‑CCACAC‑TGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCT‑GA‑CCG‑T‑CAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAA‑GAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGATAAGTGAGGGCAGGAGGGCTTCCATCAAATTAAACGCCACAAAGAAAAGCTGCACGCCGac                                                             >  1:110948/1‑238 (MQ=255)

                                                                                                                                                                                            tGTCAG‑CCACAC‑TGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCT‑GA‑CCG‑T‑CAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAA‑GAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGATAAGTGAGGGCAGGAGGGCTTCCATCAAATTAAACGCCACAAAGAAAAGCTGCACGCCGACCACCAGtt                                                     >  2:886850/1‑239 (MQ=255)

                                                                                                                                                                                            tGTCAG‑CCACAC‑TGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCT‑GA‑CCG‑T‑CAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAA‑GAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGATAAGTGAGGGCAGGAGGGCTTCCATCAAATTAAACGCCACAAAGAAAAGCTGCACGCCGACCACCAGtt                                                     >  1:576874/1‑239 (MQ=255)

                                                                                                                                                                                            tGTCAG‑CCACAC‑TGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCT‑GA‑CCG‑T‑CAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAA‑GAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGATAAGTGAGGGCAGGAGGGCTTCCATCAAATTAAACGCCACAAAGAAAAGCTGCACGCCGACCACCAGtt                                                     >  1:559657/1‑239 (MQ=255)

                                                                                                                                                                                            tGTCAG‑CCACAC‑TGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCT‑GA‑CCG‑T‑CAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAA‑GAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGATAAGTGAGGGCAGGAGGGCTTCCATCAAATTAAACGCCACAAAGAAAAGCTGCACGCCGACCACCAGtt                                                     >  1:548435/1‑239 (MQ=255)

                                                                                                                                                                                            tGTCAG‑CCACAC‑TGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCT‑GA‑CCG‑T‑CAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAA‑GAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGATAAGTGAGGGCAGGAGGGCTTCCATCAAATTAAACGCCACAAAGAAAAGCTGCACGCCGACCACCAGtt                                                     >  1:522379/1‑239 (MQ=255)

                                                                                                                                                                                            tGTCAG‑CCACAC‑TGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCT‑GA‑CCG‑T‑CAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAA‑GAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGATAAGTGAGGGCAGGAGGGCTTCCATCAAATTAAACGCCACAAAGAAAAGCTGCACGCCGACCACCAGt                                                      >  1:233532/1‑238 (MQ=255)

                                                                                                                                                                                            tGTCAG‑CCACAC‑TGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCT‑GA‑CCG‑T‑CAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAA‑GAAACTGGCTGGTGGAGTAAACACCCATCGACGTACCTTTGTAACATACTGGCGACTCTTTACTGATAAGTGAGGGGCGGAGGGCTTCCATCAAATTAAACGCCACAAAGAAAAGCTCCACGCCGaccacccgac                                                     >  2:756886/1‑235 (MQ=255)

                                                                                                                                                                                                tccgccccccggcgcCCAGCATTGCGCCACCGAAAAACACCCCCT‑GC‑CCG‑T‑AAACCATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAA‑GAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGATAAGTGAGGGCAGGAGGGCTTCCATCAAAATAAACGCCACAAAGAAAAGCTGCACGCCGACCACCAGTTGc                                                   <  1:646280/228‑1 (MQ=255)

                                                                                                                                                                                                      cac‑TGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCT‑GA‑CCG‑T‑CAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAA‑GAAACCGGCTGGGGGAGTAAACACCCATCGCCGGACCTTTGTAACCTGCTGGCGACTCTTTACTGATAAGAGAGGGAAGGAGGGGTTCCATCAAATTAAACGCCCCAAAGAAAAGCTGCACGCCGACCACAAGtgacccgc                                               >  2:896660/1‑229 (MQ=255)

                                                                                                                                                                                                          tGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCT‑GA‑CCG‑T‑CACACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAA‑GAAACTGGCTGGTGGAGGAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGAGAAGTGAGGGCAGGAGGCCTTCCATCAAATTAAACGCCACAAAGGAAAGCTGCACGACGCCCACCAGTTGCCTGAACTGCg                                         >  2:645707/1‑239 (MQ=255)

                                                                                                                                                                                                          tGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCT‑GA‑CCG‑T‑CAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAG‑GAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGATAAGTGAGGGCAGGAGGGCTTCCATCAAATTAAACGCCACAAAGAAAAGCTGCACGCCGACCACCAGTTGCCAGAACTGCg                                         >  1:396049/1‑239 (MQ=255)

                                                                                                                                                                                                          tGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCT‑GA‑CCG‑T‑CAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAA‑GAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGATAAGTGAGGGCAGGAGGGCTTCCATCAAATTAAACGCCACAAAGAAAAGCTGCACGCCGACCACCAGTTGCCAGAACTGCg                                         >  1:213192/1‑239 (MQ=255)

                                                                                                                                                                                                          tGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCT‑GA‑CCG‑T‑CAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAA‑GAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGATAAGTGAGGGCAGGAGGGCTTCCATCAAATTAAACGCCACAAAGAAAAGCTGCACGCCGACCACCAGTTGCCAGAACTGCg                                         >  1:416220/1‑239 (MQ=255)

                                                                                                                                                                                                            cGGCCAGATTGGCGCAAGCGAAAAATCCCCCCG‑GA‑CCG‑C‑AAAACATGCCTTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAA‑GAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGATAAGTGAGGGCAGGAGGGCTTCCATCAAATTAAACGCCACAAAGAAAAGCTGCACGCCGACCACCAGTTGCCAGAACTGCGt                                        <  1:751384/238‑1 (MQ=255)

                                                                                                                                                                                                               ccAGCATTGCGCCAGCGAGAAATACCCCCT‑GA‑CCG‑T‑CAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAA‑GAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGAt                                                                                                                         >  1:18677/1‑154 (MQ=255)

                                                                                                                                                                                                                cggaatggcgcCAGC‑AAAAATCCCCCCC‑GAGCCG‑A‑AAAAAATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAA‑GAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGATAAGTGAGGGCAGGAGGGCTTCCATCAAATTAAACGCCACAAAGAAAAGCTGCACGCCGACCACCAGTTGCCAGAACTGCGTTTgcg                                   <  1:1058238/232‑1 (MQ=255)

                                                                                                                                                                                                                  ccactGCGCCCGCGAAACACCCCCCCT‑CA‑CCC‑C‑CACACATCCCTTTAACCCCCCCGCCCACCGAACCTCCAATCGCCACGCCAA‑GAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGATAAGTGAGGGCAGGAGGGCTTCCATCAAATTAAACGCCACAAAGAAAAGCTGCACGCCGACCACCAGTTGCCAGAACTGCGTTTGcgcg                                 <  2:213192/235‑1 (MQ=255)

                                                                                                                                                                                                                   cATGGCGCCAGCGAAAAAACCCCCCT‑GC‑CCG‑C‑CAACCATGCCGTTAATCCAGCCGCCCACCGAACCGCCAATCGCCACGCCAA‑GAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGATAAGTGAGGGCAGGAGGGCTTCCATCAAATTAAACGCCACAAAGAAAAGCTGCACGCCGACCACCAGTTGCCAGAACTGCGTTTgcg                                   <  1:921691/236‑1 (MQ=255)

                                                                                                                                                                                                                    aTTGCGCCGGCGAAAAATCCCCCCT‑GC‑CCG‑T‑AAACCATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAA‑GAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGATAAGTGAGGGCAGGAGGGCTTCCATCAAATTAAACGCCACAAAGAAAAGCTGCACGCCGACCACCAGTTGCCAGAACTGCGTTTGcgcg                                 <  1:426537/237‑1 (MQ=255)

                                                                                                                                                                                                                             gCGAGAAATACCCCCT‑GA‑CCG‑T‑CAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAA‑GAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGATAAGTGAGGGCAGGAGGGCTTCCATCAAATTAAACGCCACAAAGAAAAGCTGCACGCCGACCACCAGTTGCCAGaa                                              >  1:316015/1‑215 (MQ=255)

                                                                                                                                                                                                                                          ccT‑GA‑CCT‑C‑AAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAA‑GAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGATAAGTGAGGGCAGGAGGGCTTCCATCAAATTAAACGCCACAAAGAAAAGCTGCACGCCGACCACCAGTTGCCAGAACTGCGTTTGCGCGTTCCACAACACAATTTCCGCAAcc         <  1:1075673/239‑1 (MQ=255)

                                                                                                                                                                                                                                          ccT‑GA‑CCG‑C‑CAAACAGGCCGTAAATCCAGCCGCCCAGCGACCCGCCAATCGCCACGCCAA‑GAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGATAAGTGAGGGCAGGAGGGCTTCCATCAAATTAAACGCCACAAAGAAAAGCTGCACGCCGACCACCAGTTGCCAGAACTGCGTTTGCGCGTTCCACAACACAATTTCCGCAAcc         <  1:602558/239‑1 (MQ=255)

                                                                                                                                                                                                                                              gA‑CCG‑T‑CAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAA‑GAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGATAAGTGAGGGCAGGAGGGCTTCCATCAAATTAAACGCCACAAAGAAAAGCTGCACGCCGACCACCAGTTGCCAGAACTGCGTTTGCGCGTTCCACAACACAATTTCCGCAACCACg      >  2:508963/1‑239 (MQ=255)

                                                                                                                                                                                                                                                     t‑CAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAA‑GAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGATAAGTGCGGGCAGGAGGGCTTCCATCAAATTAAACGCCACAAAGAAAAGCTGCACGCCGACCACCAGTTGCCAGAACTGCGTTTGCGCGTTCCACAACACAATTTCCGCAACCACGATCa  >  1:125490/1‑238 (MQ=255)

                                                                                                                                                                                                                                                     t‑CAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAA‑GAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGATAAGTGAGGGCAGGAGGGCTTCCATCAAATTAAACGCCACAAAGAAAAGCTGCACGCCGACCACCAGTTGCCAGAACTGCGTTTGCGCGTTCCACAACACAATTTCCGCAACCACGATCa  >  2:63255/1‑238 (MQ=255)

                                                                                                                                                                                                                                                     t‑CAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAA‑GAAACTGGCTGGTGGAGGAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGATAAGTGAGGGCAGGAGGGCTTCCATCAAATTAAACGCCACAAAGAAAAGCTGCACGCCGACCACCAGTTGCCAGAACTGCGTTTGCGCGTTCCACAACCCAATTTCCGCAACCACGATCa  >  1:908900/1‑238 (MQ=255)

AAGCGATTCGTCACTTTGCTGTCGATTTTCACATAAGCTGAATGTTCTTCTTCTGCAATCAACACTTCTTTGATGCCTTCAGTTTCTAGCAAACGCACTTTTAACG‑CCTCGTTTGCGGCAATGTTCGCCGGGATTT‑C‑AATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAG‑CCACAC‑TGCGGCCAGCATTGCGCCAGCGAGAAATACCCCCT‑GA‑CCG‑T‑CAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAA‑GAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGATAAGTGAGGGCAGGAGGGCTTCCATCAAATTAAACGCCACAAAGAAAAGCTGCACGCCGACCACCAGTTGCCAGAACTGCGTTTGCGCGTTCCACAACACAATTTCCGCAACCACGATCA  >  NC_000913/445333‑445807

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