BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	954,374	C→G	22.8%	P31P (CCG→CCC)	focA ←	formate channel

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	954,374	0	C	G	22.8%	124.4 / 12.1	71	P31P (CCG→CCC)	focA	formate channel
Reads supporting (aligned to +/- strand): ref base C (19/32); new base G (2/14); total (23/48)
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.61e-01

AACTGACCCCAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAG‑AGAGAAGCAAATGCCGCCAACCAGTTTT‑GCC‑ATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGA‑TATAGAAGACGAATGCGATTGAGATGAAAACAC‑CGGCGGTAATCGCCAGAT‑AG‑AAAGTCTTAAGCGGATG‑TTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGCAAATACGGAGT  >  NC_000913/954138‑954595

aaCTGACCCCAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACATCACAAAGAATCAGCCCCAG‑AGAGAAGCAAATGCCGCCAACCAGTTTT‑GCC‑ATGCCGAAGGTCATTGTGCGTGTGCCA‑GTGGTTGCTGTGA‑TATAGAAGACGAATGCGATTGAGATGAAAACAC‑CGGCGGTAATCGCCAGAT‑GG‑AAAGTCTTAAGCgg                                                                                                                                                                                                                              <  2:383322/239‑1 (MQ=255)

  cTGACCCCAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCACCAG‑AGAGAAGCAAATGCCGCCAACCAGTTTT‑GCT‑ATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGA‑TATAGAAGACGAATGCGATTGAGATGAAACCGC‑AGGCGGTAATCGCCAGAT‑AG‑AAAGTCTTAAGCGGAt                                                                                                                                                                                                                            <  2:606206/239‑1 (MQ=255)

      ccccAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAG‑AGAGAAGCAAATGCCGCCAACCAGTTTT‑GCC‑ATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGA‑TATAGAAGACGAATGCGATTGAGATGAAAACAC‑CGGCGGTAATCGCCAGAT‑AG‑AAAGTCTTAAGCGGATG‑ttt                                                                                                                                                                                                                       <  2:981502/239‑1 (MQ=255)

        ccAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAG‑AGAGAAGCAAATGCCGCCAACCAGTTTT‑GCC‑ATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGA‑TATAGAAGACGAATGCGATTGAGATGAAAACAC‑CGGCGGTAATCGCCAGAT‑AG‑AAAGTCTTAAGCGGATG‑TTTCg                                                                                                                                                                                                                     <  1:1070199/239‑1 (MQ=255)

        ccAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAG‑AGAGAAGCAAATGCCGCCAACCAGTTTT‑GCC‑ATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGA‑TATAGAAGACGAATGCGATTGAGATGAAAACAC‑CGGCGGTAATCGCCAGAT‑AG‑AAAGTCTTAAGCGGATG‑TTTCg                                                                                                                                                                                                                     <  2:233335/239‑1 (MQ=255)

        ccAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAG‑AGAGAAGCAAATGCCGCCAACCAGTTTT‑GCC‑ATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGA‑TATAGAAGACGAATGCGATTGAGATGAAAACAC‑CGGCGGTAATCGCCAGAT‑AG‑AAAGTCTTAAGCGGATG‑TTTCg                                                                                                                                                                                                                     <  2:106477/239‑1 (MQ=255)

        ccAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAG‑AGAGAAGCAAATGCCGCCAACCAGTTTT‑GCC‑ATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGA‑TATAGAAGACGAATGCGATTGAGATGAAAACAC‑CGGCGGTAATCGCCAGAT‑AG‑AAAGTCTTAAGCGGATG‑TTTCg                                                                                                                                                                                                                     <  1:960056/239‑1 (MQ=255)

        ccAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAG‑AGAGAAGCAAATGCCGCCAACCAGTTTT‑GCC‑ATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGA‑TATAGAAGACGAATGCGATTGAGATGAAAACAC‑CGGCGGTAATCGCCAGAT‑AG‑AAAGTCTTAAGCGGATG‑TTTCg                                                                                                                                                                                                                     <  1:883888/239‑1 (MQ=255)

        ccAGGTGATGCGCACCGTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGAACGGCTGCGCAGACAACACAAAGAATCAGCCCCAG‑AGAGAAGCAAATGCCGCCATCCAGTTTT‑GCC‑ATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGA‑TATAGAAGACGAATGCGATTGAGATGAAAACAC‑CGGCGGTAATCGCCAGAT‑AG‑AAAGTCTTAAGCGGATG‑TTTCg                                                                                                                                                                                                                     <  1:594935/239‑1 (MQ=255)

        ccAGGAGATGCGCCCACTCGCCTTAGCAAGAACAAGCAACACAGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAG‑AGAGAAGCAAATGCCGCCAACCAGTTTT‑GCC‑ATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGA‑TATAGAAGACGAATGCGATTGAGATGAAAACAC‑CGGCGGTAATCGCCAGAT‑AG‑AAAGTCTTAAGCGGAGG‑TTTCg                                                                                                                                                                                                                     <  2:1025068/239‑1 (MQ=255)

           ggTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAG‑AGAGAAGCAAATGCCGCCAACCAGTTTT‑GCC‑ATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGA‑TATAGAAGACGAATGCGATTGAGATGAAAACAC‑CGGCGGTAATCGCCAGAT‑AG‑AAAGTCTTAAGCGGATG‑TTTCGtt                                                                                                                                                                                                                   >  1:1111108/1‑238 (MQ=255)

            gTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAG‑AGAGAAGCAAATGCCGCCAACCAGTTTT‑GCC‑ATGCCGAAGGGCATTGCGCCTGTGCCA‑GTGGTTGCTGTGA‑TATAGAAGACGAATGCGATTGAGATGAAAACAC‑CGGCGGTTATCGCCCGAT‑AG‑AAAGGCTTAACCGGATG‑TTTCCtt                                                                                                                                                                                                                   >  1:899742/1‑237 (MQ=255)

               aTGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAG‑AGAGAAGCAAATGCCGCCAACCAGTTTT‑GCC‑ATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGA‑TATAGAAGACGAATGCGATTGAGATGAAAACAC‑CGGCGGTAATCGCCAGGT‑AG‑AAAGGCTTAAGCGGATG‑TTTCCGTGCttt                                                                                                                                                                                                              >  2:782293/1‑239 (MQ=255)

                 gcgcCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAG‑AGAGAAGCAAATGCCGCCAACCAGTTTT‑GCC‑ATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGA‑TATAGAAGACGAATGCGATTGAGATGAAAACAC‑CGGCGGTAATCGCCAGAT‑AG‑AAAGACTTAAGCGGATG‑TTTCGTTGCTTtat                                                                                                                                                                                                            >  1:539205/1‑239 (MQ=255)

                 gcgcCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAG‑AGAGAAGCAAATGCCGCCAACCAGTTTT‑GCC‑ATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGA‑TATAGAAGACGAATGCGATTGAGATGAAAACAC‑CGGCGGTAATAGCCAGAT‑AG‑AAAGGCTTAAGAGGGTG‑GTTCGTTgtgtttt                                                                                                                                                                                                            >  2:819249/1‑233 (MQ=255)

                 gcgcCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAG‑AGAGAAGCAAATGCCGCCAACCAGTTTT‑GCC‑ATGCCGAAGGGCATTGCGCCTGTGCCA‑GTGGTTGCTGTGA‑TATAGAAGACGAATGCGATTGAGATGAAAACAC‑CGGCGGGAATCACCACGT‑AA‑AAAGGCTTAAGAGGATGATTTCGtgacttat                                                                                                                                                                                                             >  2:1117218/1‑232 (MQ=255)

                    cccACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAG‑AGAGAAGCAAATGCCGCCAACCAGTTTT‑GCC‑ATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGA‑TATAGAAGACGAATGCGATTGAGATGAAAACAC‑CGGCGGTAATCGCCAGAT‑AG‑AAAGTCTTAAGCGGATG‑TTTCGTTGCTTTATAg                                                                                                                                                                                                          >  1:801427/1‑238 (MQ=255)

                                   aaaaaaaaacAACACGGAGAAAAAAAAGAGACCGGCCCCGCAAACAAAACAAAGAACCAGCCCCAG‑AGAGAACCAAATCCCGCAAAACAGTTTG‑GAC‑AGGCAGAAGGGAATTGGGCCTGGGCCG‑GTGGTGGCTGATA‑TATAAAAAACGAATGCGATTGGGATAAAACCCC‑CGGCGGTAATCCCCAGAG‑AA‑AAATTCTTAAGGGGGTG‑TTTGGTTGCTTTATAGACACCCGCCTCTTCgg                                                                                                                                                                                          <  2:661349/230‑1 (MQ=255)

                                     caacaaTCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAG‑AGAGAAGCAAATGCCGCCAACCAGTTTT‑GCC‑ATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGA‑TATAGAAGACGAATGCGATTGAGATGAAAACAC‑CGGCGGTAATTGCCCGAT‑AG‑AAAGTCTTAAGGGGAAG‑TTTCGTTGCTTTTTAGACACCCGCCGCTTTGGcc                                                                                                                                                                                        >  2:436184/1‑239 (MQ=255)

                                     caacaaTCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAG‑AGAGAAGCAAATGCCGCCAACCAGTTTT‑GCC‑ATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGA‑TATAGAAGACGAATGCGATTGAGATGAAAACAC‑CGGCGGTAATCGCCAGAT‑AG‑AAAGTCTTAAGCGGATG‑TTTCGTTGCTTTATAGACACCCGCCTCTTCGGc                                                                                                                                                                                         >  1:508630/1‑238 (MQ=255)

                                     caacaaTCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAG‑AGAGAAGCAAATGCCGCCAACCAGTTTT‑GCC‑ATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGA‑TATAGAAGACGAATGCGATTGAGATGAAAACAC‑CGGCGGTAATCGCCAGAT‑AG‑AAAATCTTAAGCGGATG‑TTTCGTTGCTTTATAGACACCCGCCTCTTCggg                                                                                                                                                                                         >  1:730962/1‑237 (MQ=255)

                                                      gaagaaAAGAGACCGCCCCCGCAAACAACAAAAAGAATCAGCCCCAG‑AGAGAGGAAATTGCCG‑CAACCCGTTTTGGCC‑AGGCGGAGGGGTATTTTGCCTG‑GCCAGGTGGTTGCTGTAT‑TAAAAAAGACAAATGCTTTTGGAATAAAACCCC‑GGGCGTTAACCCCCAGAA‑AA‑AAAGTCTTAAGCGGTTG‑TTTCTTTGCTTTAAAGACCCCCGCCTCTTCGGCCACTTTGGCCATTgcag                                                                                                                                                                        <  2:449266/238‑1 (MQ=255)

                                                                gaTCGGCTCCGCAGACAACACAAAGAATCAGCCCCAG‑AGAGAAGCAAATGCCGCCAACCAGTTTT‑GCC‑ATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGA‑TATAGAAGACGAATGCGATTGAGATGAAAACAC‑CGGCGGTAATCGCCAGAT‑AG‑AAAGTCTTAAGCGGATG‑TTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCCTTGCCGCAGGAAGTaaa                                                                                                                                                             >  1:859825/1‑239 (MQ=255)

                                                                gaTCGGCTCCGCAGACAACACAAAGAATCAGCCCCAG‑AGAGAAGCAAATGCCGCCAACCAGTTTT‑GCC‑ATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGA‑TATAGAAGACGAATGCGATTGAGATGAAAACAC‑CGGCGGTAATCGCCAGAT‑AG‑AAAATATTAAGGGGATG‑TTTCGTTGCTTTTTAGACCCCCGCCCCTTCGGCCCCTTTTGGCATTGCCGCCGGAAGTaaa                                                                                                                                                             >  2:328116/1‑239 (MQ=255)

                                                                gaTCGGCCCCGCAAACAACACAAAAAACCAGCCCCAG‑AGAGAAGCAAATGCCGCCAACCAGTTTG‑GCC‑ATGCCGAGGGGAATTGTGCCTGGGCCA‑GTGGTTGCTGTGA‑TATAGAAGACGAATGCGATGGAGATAAAAACCC‑CGGCGGTAATCGCCAGAT‑AA‑AAAGTCTTAAGCGGATG‑TTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTaa                                                                                                                                                              <  1:237237/238‑1 (MQ=255)

                                                                 aTCGGCTCCGCAGACAACACAAAGAATCAGCCCCAG‑AGAGAAGCAAATGCCGCCAACCAGTTTT‑GCC‑ATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGA‑TATAGAAGACGAATGCGATTGAGATGAAAACAC‑CGGCGGTAATCGCCAGAT‑AG‑AAAGTCTTAAGCGGATG‑TTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTaaaa                                                                                                                                                            >  1:7927/1‑239 (MQ=255)

                                                                 aTCGGCCCCGCAGAAAACACAAAGAATCAGCCCAAG‑AGAGAAGCAACTGCCGCCAACGATTTTT‑GCC‑ATGCCGAAGGGTGTTGTGCCTGTGCCA‑GTGGTTGCTGTAA‑TATAGAAGCCGAAGGAGATTGAGATAAAAACCC‑CGGCGGTAAGCGCAAGAT‑AG‑AAAGTCTTAGGGGGATG‑TTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTaaaa                                                                                                                                                            <  2:7927/239‑1 (MQ=255)

                                                                    ggCTCCGCAGACAACACAAAGAATCAGCCCCAG‑AGAGAAGCAAATGCCGCCAACCAGTTTT‑GCC‑ATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGA‑TATAGAAGACGAATGCGATTGAGATGAAAACAC‑CGGCGGTAATCGCCAGAT‑AG‑AAAGTCTTAAGCGGATG‑TTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGa                                                                                                                                                          >  1:1028848/1‑238 (MQ=255)

                                                                                aacaCAAAGAATCAGCCCCAG‑AGAGAAGCAAATGCCGCCAACCAGTTTT‑GCC‑ATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGA‑TATAGAAGACGAATGCGATTGAGATGAAAACAC‑CGGCGGTAATCGCCAGAT‑AG‑AAAGTCTTAAGCGGATG‑TTTCGTTGCTTTATAGAAAACCGCCTCTTCCGCCACTTTGGCCATTGCAGCAGGAAATAAAAGGTCAAAAGGGGTgc                                                                                                                                             >  1:1093780/1‑238 (MQ=255)

                                                                                aaaaaaaaaaaTCA‑CCCCAGAAGAGAAGAAAATGCC‑CGAAAGAGTTTT‑TCCGACGCAGAAGGGATTTTTGCCTTTGCCA‑GTGTTGGTTGTAA‑TAAAAAAGAGAATTGCGATTGAAAAGAAACCCC‑CGGCGGTAATCGCCAAAT‑AG‑AAATTCTTAAGCGGATG‑TTTCGTTGTTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGt                                                                                                                                             <  2:1028848/234‑1 (MQ=255)

                                                                                      aaaaaccAGCCCCAG‑AGAGAAGCAAAGGCCGCCAACAATTTTT‑GCA‑ATGCCGAAGGGATTTGTGCCTGGGCCA‑GTGGTTGCTGTGA‑TATAGAAGACGAATGCGATTGAGATGAAACCCC‑GGGCGGTAACCGCCAGAT‑AA‑AAAGTCTAAAGGGGATG‑TTTCGTTGCTTTATAGACCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCtt                                                                                                                                       <  1:955607/233‑1 (MQ=255)

                                                                                          aTCAGCCCCAG‑AGAGAAGCAAATGCCGCCAACCAGTTTT‑GCC‑ATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGA‑TATAGAAGACGAATGCGATTGAGATGAAAACAC‑CGGCGGTAATCGCCAGAT‑AG‑AAAGTCTTAAGCGGATG‑TTTCGTTGCTTTATAGACACCCGCCTCCTCGGCCACTTTGGc                                                                                                                                                                                >  1:799152/1‑194 (MQ=255)

                                                                                            cAGCCCCAG‑AGAGAAGCAAATGCCGCCAACCAGTTTT‑GCC‑ATGCCGAAGGGCATTGTGACTGTGCCA‑GTGGTTGCTGTGA‑TATAGAAGACGAATGCGATTGAGATGAAAACAC‑CGGCGGTTAATGCCAGAT‑AG‑AAAGGCTTAAACGGGTG‑TTTCGTTGCTTTATAGAAACCCGCCCCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGGTGGTAGCTTTCacac                                                                                                                                 >  1:559591/1‑239 (MQ=255)

                                                                                             aGCCCCAG‑AGAGAAGCAAATGCCGCCAACCAGTTTT‑GCC‑ATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGA‑TATAGAAGACGAATGCGATTGAGATGAAAACAC‑CGGCGGTAATCGCCAGAT‑AG‑AAAGGCTTAAGCGGAGG‑ATTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACCTTGGCCCTTGCAGCCGGAAGTAAAAGATCAAAACGGGTGTTAGCTTTCACACt                                                                                                                                >  1:1165095/1‑239 (MQ=255)

                                                                                              accgcgcgaaggaAACAAGTGCGCCCACCCAGTTTT‑ACC‑ATGCGGAAGGGCTTTGTGCCTGTACAG‑GGTGTTGCTGTAA‑TATAAAAAGATAAAGCGGTGAAGAATAAACCCC‑GGGCGGTAACCGCCAAAT‑AG‑AAAGTCTGAAGGGGATT‑TTTCTTTGTTTTATAAAAACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGTAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACt                                                                                                                                <  2:767822/228‑1 (MQ=255)

                                                                                                 ccAG‑AGAGAAGCAAAGGCCCCCACCCAGTTTT‑GCA‑ATGCCGAAGGGCTTTGTGCCTGGGCCA‑GTGGTTGCTGTGA‑TATAGAAGACGAATGCGATTGAGATGAAACCCC‑CGGCGGTAATCGCCAGAT‑AG‑AAATTCTTAAGCGGATG‑TTTGGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAAct                                                                                                                            <  1:29956/239‑1 (MQ=255)

                                                                                                 acAG‑AGAAAAGCAAATGCCGCCAACCAGTTTA‑CCC‑ATGCGGAGGGGCATTGTGCCTGTGCCG‑GGGGTGGCTGTAA‑TATAGAAGAAGAAGGCGTTTGAGAAAAAACCCC‑GGGCGTTAACCGCCAGAT‑AG‑AAAGTCTTAAGCGTATG‑TTTTGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAAct                                                                                                                            <  2:1027389/238‑1 (MQ=255)

                                                                                                       gagaGGCAAGCGCCGCCAACTAGTTGG‑ACC‑ATGCGAAGGGGCGTTGGACTGGTGCGA‑GTGGTTGCTGTGA‑TATAGAAGACGAATGTGAATGAGATGAAACAGC‑GGGCGGTAACCGCCAAAT‑AT‑AAAGTCTTAAGGGGATG‑TTTCGTTGCTTTAAAGACAACCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAActctct                                                                                                                        <  2:274928/238‑1 (MQ=255)

                                                                                                       gagaAGCAAATGCCGCAAACCAGTTTT‑GCC‑ATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGA‑TATAGAAGACGAATGCGATTGAGATAAAACCAC‑CGGCGGTAATCGCCAGAT‑AG‑AAAGTCTTAAGCGGATT‑TTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAActctct                                                                                                                        <  1:808132/238‑1 (MQ=255)

                                                                                                          aaGCAAATGCCGCCAACCAGTTTT‑GCC‑ATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGA‑TATAGAAGACGAATGCGATTGAGATGAAAACAC‑CGGCGGTAATCGCCAGAT‑AG‑AAAGTCTTAAGCGGATG‑TTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAAATCTCTCttt                                                                                                                    >  1:66596/1‑239 (MQ=255)

                                                                                                                             ttttt‑GCC‑ATGCCGAAGGGCATTGTGCCTGGGCCA‑GGGGTTGCTGTGA‑TAAAAAAGACAAATGCGATTGAGATGAAAACCC‑CGGCGGTAATCGCCAGAT‑AG‑AAAGTCTTAAGCGGATG‑TTTCTTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGAt                                                                                                 <  1:4916/238‑1 (MQ=255)

                                                                                                                              gttgcgCC‑ATGCGGAAGTGCATTGCGCATGTGCCA‑TAGGTTGCT‑TGATTATAAAAGCCGAATGCTATTGAGATAAAACCAC‑CGGCGTTAACCGCCAGAT‑AG‑AAATTCTAAAGCGGATT‑TTTCTTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGAt                                                                                                 <  2:376142/234‑1 (MQ=255)

                                                                                                                               tgt‑GCT‑ATACCGAAGGGGATGGTGCGTGGCCCA‑GGGGTTGCTGTAA‑TAAAGAAGCAGATTTCTATTGAGAGAAAAACAC‑CGGCGTAACTCGCGAGAG‑AA‑AAGTTCTCAAGGGGATT‑TTTCGTTGCTTTATAAACCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACgaga                                                                                                  <  2:539205/234‑1 (MQ=255)

                                                                                                                                 tggtgaaTGCCGAAGGTGATGGGGCCGGTGCCG‑GTGGTGGCGTTAA‑AAAAGAAGACGAAGGGGATTGTAATGAAAACAC‑CGGCGGTAATCGCCAGAT‑AA‑AAAGTCTTAGGGGGAGG‑TTTCGTTGCTTTAAAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATa                                                                                                <  2:730962/232‑1 (MQ=255)

                                                                                                                                     c‑ATGCCGAAGGGTTTTGTGCCTGTGCAA‑GTGGTTGCTGTGA‑TATAGAAGACAAATGCGATTGAGATGAAAACCC‑CGGCGGTAATCGCCAGAT‑AG‑AAAGTCTTAAGCGGATT‑TTTCGTTGTTTTAAAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACa                                                                                          <  1:328116/239‑1 (MQ=255)

                                                                                                                                      cctgacggcgggCATTGGGCCTGTGCAG‑ATGGTTGGTGTGT‑AAAGGGAGACAAATGAGATTGAGATGAAAACGC‑CGGATTCAATCTCCGGAA‑AA‑AAGTTATGACGCGGTTT‑TTTCGTTCCTTTGAATACACCCGCCTTTTCGGCCACTTTGGGCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCCCTCTTTATTAAATCGGCGGCGAGATGCTAACa                                                                                          <  2:725102/230‑1 (MQ=255)

                                                                                                                                      catcacGAGGGGCTTTGTGCTGTTGCAG‑GTGGTTCTTTTGA‑TATAGAAAACAAAGGCGATTGAGATAAAACCCC‑CGGCGTTATTCGCCAGAT‑AA‑AATGCTTTACGCGGATT‑TTTCTTTTTTTTAGAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACa                                                                                          <  2:1093780/234‑1 (MQ=255)

                                                                                                                                           ggaggggTATGGTGCGTGTGCAA‑GGGTTTGCTGTGA‑TATAGAAGACGAAGGCGATTGAGATAAAAACCC‑CGGCGGTAACCGCCAGAT‑AG‑AAATTCTTAAGCGGTTG‑TTTCGTTTCTTTATAGACCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCa                                                                                     <  1:998147/235‑1 (MQ=255)

                                                                                                                                              agggagtggtgcctgggccA‑GGGGTTGCTGTAA‑TATAAAAGAGGAAGGCGATTGAGATGAAAACAC‑CGGCGGTAATCGCCAGAT‑AG‑AAAGTCTAAAGCGGATG‑TTTCGTTTCTTTATAGACCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATta                                                                                  <  1:402393/231‑1 (MQ=255)

                                                                                                                                                    gtgtgcgtgtgcaA‑GTGGTGTCTGTAA‑TTGAGAAGGCAAATGCGAAACAGTTGAAAGCAC‑CGGGTGAAACCGACGAAT‑AG‑AAGGTGTTAATGGCAGT‑GTTGGTTGCTTTAGTGACTCCCGCCGCTGCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTAAGCTTTCACACTAACTCTCGCTTTATTAAGGCGGAGACGAGATACTAACAAAGCATGATAGAt                                                                             <  2:904878/237‑1 (MQ=255)

                                                                                                                                                          ttggacta‑GTGTTGGCGGTGA‑TATAAAAGCCGAATGGTATTGAAAGGAAAACAC‑CGGCGGTAACCCCCAAAT‑AG‑AAAGTCTTAAGGGGATG‑TTTAGTTGTTTTATACACAACCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGaaa                                                                       <  1:836115/231‑1 (MQ=255)

                                                                                                                                                           tgggacA‑GTGGTTGTGGTGA‑TAAAAAAGACGAACGCGATAGAAATAAAACCCC‑GGGCGGTAACCCCCAGAA‑AG‑AAAGTCTAAAGCGGTTG‑TTTCGTGGTTTTAAAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAAt                                                                      <  2:508630/233‑1 (MQ=255)

                                                                                                                                                                 atacgagaagggtgA‑TATAGAAGGAGAATGGGATTGAGATGAAAACAC‑CGGCGGTAACGGCCAGTA‑AG‑AACGGTTAAAGCGGATG‑TTTCGTTTCTATAAAAAAGCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGatat                                                                <  2:910235/228‑1 (MQ=255)

                                                                                                                                                                        tGCTGTGA‑TATAGAAGACGAATGCGATTGAGATGAAAACAC‑CGGCGGTAATCGCCAGAT‑AG‑AAAGTCTTAAGCGGATG‑TTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCat                                                         >  1:355290/1‑239 (MQ=255)

                                                                                                                                                                               a‑taaAGAAGACGAATGCGATGGAGATAAAAACCC‑CGGCGGTAATCGCCAGAT‑AG‑AAAGTCGTAAGCGGTTG‑TTTTGTTAGTTTAAAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGc                                                  <  1:425231/239‑1 (MQ=255)

                                                                                                                                                                                ataaaaaagaCGAATGCGATGGAAATAAAACCCC‑CGGCGTAAAACGCAAGAT‑AG‑AAGGCCTAAGGAGGGTT‑TTTCGTTGTTTAAAAGACCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGc                                                  <  1:24071/233‑1 (MQ=255)

                                                                                                                                                                                 tatagacacagAATGCGATTAAGATAAAAACCC‑CGGC‑TTAAACCCCATAA‑GGAAAAGTTTTAGGGGTTTG‑TTTTGTGGTTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGcc                                                 <  2:66596/229‑1 (MQ=255)

                                                                                                                                                                                 tataGAAGACGAATGCGATTGAGATGAAAACAC‑CGGCGGTAACCGCCAGAT‑AG‑AAAGTCTTAAGCGATTG‑TTTCTTTGCTTTATAGCCACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGcc                                                 <  1:1130907/239‑1 (MQ=255)

                                                                                                                                                                                 tataGAAGACGAAGGCGATTGAGATGAAAACAC‑CGGCGGAAATCGCCAGAT‑AG‑AAAGTTTAAAGGGGATT‑TTTCGTTGTTTAAGAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGcc                                                 <  1:27582/239‑1 (MQ=255)

                                                                                                                                                                                 tataGAAAGAGAATGGAATTGAGATAAAAACGC‑CGGCGTAAATGCCCGAAT‑AA‑AAAGTTTGAAGCGTATT‑TTTCGTTGCTTTAGAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACGCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGcc                                                 <  2:1165095/239‑1 (MQ=255)

                                                                                                                                                                                  ataaagaaagAATGCGATTAAGATAAAAACCC‑CGGGGGTAACCGCC‑GATAAG‑AAGGTCTTAGGCGGTTT‑TTTGTTTGCTTTATACACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCt                                                <  2:355290/230‑1 (MQ=255)

                                                                                                                                                                                   taGAAGAGGAATTGGATGGAGAT‑AAAACACGCGGGGTTAACCGCAAGAT‑AG‑AAAGTTTCAAGGTGTGT‑TTTCGTTCTTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGcc                                                 <  2:488192/237‑1 (MQ=255)

                                                                                                                                                                                        gaCGAATGCGATTGAGATGAAAACAC‑CGGCGGAAATCGCCAGAT‑AG‑AAAGTCTTAAGCGGTTG‑TTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCCCTTTTGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCCGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAACCATTATAGATGAGAAACTGATATAAATCATAtc                                                      >  1:910235/1‑227 (MQ=255)

                                                                                                                                                                                          cGAATGCGATTGAGATGAAAACAC‑CGGCGGTAATCGCCAGAT‑AG‑AAAGTCTTAAGCGGATG‑TTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTTGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCACCTGGCttttt                                        >  1:436002/1‑236 (MQ=255)

                                                                                                                                                                                          cGAATGCGATTGAGATGAAAACAC‑CGGCGGTAATCGCCAGAT‑AG‑AAAGTCTTAAGCGGATG‑TTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATtata                                                                                >  1:531565/1‑199 (MQ=255)

                                                                                                                                                                                          cGAAGGCAATTGAAATGAAACCCC‑CGGCGTAAACCGCCAGAT‑AG‑AAAGTCTAAAGGGGATG‑TTTCGTTGTTTTATAGCCACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTtata                                        <  1:688388/239‑1 (MQ=255)

                                                                                                                                                                                           taaTGCGATTAATATGAAAACCC‑CGGCGATAATCGCCGGAT‑AC‑AGAGCTTGAAGTGTATG‑TTTTGTTGCTTTATAGACACCCGCCGCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAg                                       <  2:285664/238‑1 (MQ=255)

                                                                                                                                                                                                                   cggcggTAATCGCCAGAT‑AG‑AAGGTCTTAAGCGGATG‑TTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTa                <  1:436184/239‑1 (MQ=255)

                                                                                                                                                                                                                           accGCCAGAT‑AA‑AAATTCTAAAGCGGATG‑TTTCTTTGCTTTATAGACCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGCAAATa        <  1:782293/237‑1 (MQ=255)

                                                                                                                                                                                                                           aacgcaaaaT‑AG‑AAAGTCGTAAAGGGATG‑TTACGTTTCTTAAAAGACCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGCAAATa        <  1:76609/231‑1 (MQ=255)

                                                                                                                                                                                                                                    tgaaaaattgtttAAGCTGATC‑TTTCGTTGCTTTATAGACCACCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGCAAATACGGAGt  <  2:436002/227‑1 (MQ=255)

AACTGACCCCAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAG‑AGAGAAGCAAATGCCGCCAACCAGTTTT‑GCC‑ATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGA‑TATAGAAGACGAATGCGATTGAGATGAAAACAC‑CGGCGGTAATCGCCAGAT‑AG‑AAAGTCTTAAGCGGATG‑TTTCGTTGCTTTATAGACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGCAAATACGGAGT  >  NC_000913/954138‑954595

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