BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,489,485	T→G	25.2%	T234P (ACC→CCC)	yfdV ←	putative transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,489,485	0	T	G	25.2%	91.3 / 9.6	61	T234P (ACC→CCC)	yfdV	putative transporter
Reads supporting (aligned to +/- strand): ref base T (20/24); new base G (15/2); total (35/26)
Fisher's exact test for biased strand distribution p-value = 3.31e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCGACTGACATAAATCCACAAGGGAGCCGTGACGACAAAGCCCAGTACGCTCACAGCCAATGACGCTGTA‑CCAGTACGGGTGTAGACATTAAACC‑GGCTGGCAATAATGATTCCGGAGAATG‑CCGGCGGTAATGCGCCTGCCAGTACCATCATCTGCAGA‑TGTTCGCTGTTCAAATGACATGCCATACCAACGAGAAGCA‑GTGCCAGTGGCATCAGAATCAGCTTCAGGAAGGTGTTATAAGCAATTTCGGCACTGAACTCGAATTTATGTGCAGCCAGAGTCAACCCCGCAGCGAATACCGC‑TA‑CCCCTGAGTTAGCCTTCGCAATCAGATTAAAGGTTGGGTCCCATGCTGCCGGAATTTTTACCCCAACCAACACCAGGATCGTTGCCAGAACAGGTGCCCATACTACTGGCTCCTTTGCCGCAGAAATTAATGCGCTCAGATTACTATTCTTCTTACCATCC  >  NC_000913/2489251‑2489717

gTCGACTGACATAAATCCACAAGGGAGCCGTGACGACAAAGCCCAGTACGCTCACAGCCAATGACGCTGTA‑CCAGTACGGGTGTAGACATTAAACC‑GGCTGGCAATAATGATTCCGGAGAATG‑CCGGCGGTAATGGGCCGGCCAGTACCATCCTCTGCAGA‑AGTTCGCTGTTGAAATGAACAACCATACCAACGAAAAACAGGACCCAGGGCCATCA‑AAACAGGTTAAGGAAGGGGTta                                                                                                                                                                                                                                        >  1:559515/1‑239 (MQ=255)

   caccgACATAAATCCCAAAGGGACCCGTGACGAAAAAGCCCAGTACCCTCACAGCCAAGGACGCTGTA‑CCAGTACGGGTGTAGACATTAAACC‑GGCTGGCAATAATGATTCCGGAGAATG‑CCGGCGGTAATGCGCCTGCCAGTACCATCATCTGCAGA‑TGTTCGCTGTTCAAATGACATGCCATACCAACGAGAAGCA‑GTGCCAGTGGCATCAGAATCAGCTTCAGGAAGGTGTTATaa                                                                                                                                                                                                                                     <  2:655447/235‑1 (MQ=255)

             aaTCCACAAGGGAGCCGTGACGACAAAGCCCAGTACGCTCACAGCCAATGACGCTGTA‑CCAGTACGGGTGTAGACATTAAACC‑GGCTGGCAATAATGATTCCGGAGAATG‑CCGGCGGTAATGCGCCTGCCACGACCATCATCTTCAGA‑GGTTCGCGGTTCAAATGAAATGCCATAACAACGGGAGACA‑CTGCCAGTGGGATACAAAACAGCTCCAGGAAGGGGGTATaaaaaaatttgg                                                                                                                                                                                                                           >  1:661187/1‑229 (MQ=255)

                  acaAGGGAGCCGTGACGACAAAGCCCAGTACGCTCACAGCCAATGACGCTGTA‑CCAGTACGGGTGTAGACATTAAACC‑GGCTGGCAATAATGATTCCGGAGAATG‑CCGGCGGTAATGCGGCTGCCCGTACCATCATCTGCAGA‑TGTTCGCTGTTCAAATGACATGCCCTACCAACGAGCAGCA‑TTGCCAGTGGCCTCAGAATCAGCTTCAGGAAGGGGTTATAAGCAAATTTGGCACTg                                                                                                                                                                                                                      >  2:806171/1‑239 (MQ=255)

                    aaGGGAGCCGTGACGACAAAGCCCAGTACGCTCACAGCCAATGACGCTGTA‑CCAGTACGGGTGTAGACATTAAACC‑GGCTGGCAATAATGATTCCGGAGAATG‑CCGGCGGTAATGCGCCCTCCCGCACCACCCTCTTCAAATTGTCCCGGGTACAAATG‑CATGCCCTACCGACAGGAAACA‑AGCTCCGCGGGATCAGAAAAAGCCTCAGGAAAGGGTTATAAGAAAATTGGGCACTGaa                                                                                                                                                                                                                    >  1:669864/1‑239 (MQ=255)

                     agggggcctgtaccaaaaaGCCCAGTCCGCTCACAGCCAATGACGCGTTA‑CCAGTACGGTTGTAGAAATTAAACC‑GGCGGGCAATAAGGATTCCGGAGAATG‑CCGGCGGTAATGCGCCTGCCAGTACCATCATCTGCAGA‑TGTTCGCTGTTCAAATGACATGCCATACTAACGAGAAGCA‑GTGCCAGTGGCATCAGAATCAGCTTCAGGAAGGTGTTATAAGCAATTTCGGCACTGAAc                                                                                                                                                                                                                   <  2:730468/223‑1 (MQ=255)

                              tGACGACAAAGCCCAGTACGCTCACAGCCAATGACGCTGTA‑CCAGTACGGGTGTAGACATTAAACC‑GGCTGGCAATAATGATTCCGGAGAATG‑CCGGCGGTAATGCGCCTGCCAGTACCATCATCTGCAGA‑TGTTCGCTGGTCAAATGACATGCCATACCAACGAGAAGCA‑GTGCCAGGGGCATCAGAATCAGCTTCAGGAAGGGGTTTTAAGCAAATTCGGCACTGAACTCGAATTTa                                                                                                                                                                                                          >  2:71038/1‑239 (MQ=255)

                                  ccaaagccccgcctctcctcagcccAATGAC‑CTGTGCCCAGTCGGGGTGGAAATATTAACCC‑GGCTGGCAATTATGATCGCGGA‑AATGCCCGGCTGAAATGCGCCTGCCAGTACCATCATCTGCAGA‑TGTTCGCTGTTCAAATGACATGCCATACCAACGAGAAGCA‑GTGCCAGTGGCATCAGAATCAGCTTCAGGAATGTGTTATAAGCAATTTCGGCACTGAACTCGAATTTAtgtg                                                                                                                                                                                                      <  1:15228/216‑1 (MQ=255)

                                    aaaaGCCCAGTACCCTCCCAGCCAATGACGCTGTA‑CCAGTACGGGTGTAGACATTAAACC‑GGCTGGCAATAATGATTCCGGAGAATG‑CCGGCGGTAATGCGCCTGCCAGTACCATCATCTGCAGA‑TGTTCGCTGTTCAAATGACATGCCATACCAACGAGAAGCA‑GTGCCAGTGGCATCAGAATCAGCTTCAGGAAGGTGTTATAAGCAATTTCGGCACTGAACTCGAATTTAtgtg                                                                                                                                                                                                      <  2:347417/236‑1 (MQ=255)

                                                     caccccagatgccggtttccgcgTAGGGGTGTAGCCCTTAACCC‑GGCTTCACAAAGTGATTCCGGAGAATG‑CCGGGGGTAAGCCCCCTGCCAGTACCAGCATCTGCAGA‑TGTTCTCTGTTCAAATGACATGCCATACCAACGAGAAGCA‑GTGCCAGTGGCATCAGAATCAGCTTCAGGAAGGTGTTATAAGCAATTTCGGCACTGAACTCGAATTTATGTGCAACCAGAGTCAACCCCg                                                                                                                                                                                    <  1:470738/217‑1 (MQ=255)

                                                                 gCTGTA‑CCAGTACGGGTGTAGACATTAAACC‑GGCTGGCAATAATGATTCCGGAGAATG‑CCGGCGGTAATGCGCCTGCCAGTACCATCATCTGCAGA‑TGTTCGCTGTTCAAATGACATGCCATACCCACGAGAAGCC‑GTGCCAGTGGCATCAGAATCAGCTTCAGGAAGGTGTTATAAGAAATTTCGGCACTGAAATCGAATTTATGTGCAGCCAGAGGCAACCCCCCCACGAATAccc                                                                                                                                                                        >  2:766636/1‑237 (MQ=255)

                                                                     tA‑CCAGTACGGGTGTAGACATTAAACC‑GGCTGGCAATAATGATTCCGGAGAATG‑CCGGCGGTAATGCGCCTGCCAGTACCATCATCTGCAGA‑TGTTCGCTGTTCAAATGACATGCCATACAAACGAGAAGCA‑GTGCCAGTGGCATCATAATCAGCTTCAGGAAGGGGTTTTAAGCAATTTCGGCACTGAAATCTAATTTTTGTGCAGCCCAAGTCAACCCCGCAAGCAATACCCC‑Taa                                                                                                                                                                   >  2:397396/1‑237 (MQ=255)

                                                                       cccgtgaCGGGTGTAAAACTTAAACC‑GGCTGGCAATATAGATTCCGGAGAATG‑CCGGCGGTAATGCCCCTGCCAGTACCATCATCTGCAGA‑TGTTCGCTGGTCAAATGACATGCCATACCAACGAGAAGCA‑GTGCCAGTGGCATCAGAATCAGCTTCAGGAAGGTGTTATAAGCAATTTCGGCACTGAACTCGAATTTATGTGCAGCCAGAGTCAACCCCGCAGCGAATACCGC‑TA‑ccc                                                                                                                                                                <  1:89475/233‑1 (MQ=255)

                                                                        ccatgaggggTGTGAACATTAAACC‑GGCTGGCAATAATGATTCCGGAGAATG‑CCGGCGGTAATGCGCCTGCCAGTACCATCATCTGCAGA‑TGTTCGCTGTTCAAATGACATGCCATACCAACGAGAAGCA‑GTGCCAGTGGCATCAGAATCAGCTTCAGGAAGGTGTTATAAGCAATTTCGGCACTGAACTCGAATTTATGTGCAGCCAGAGTCAACCCCGCAGCGAATACCGC‑TA‑ccc                                                                                                                                                                <  2:829438/231‑1 (MQ=255)

                                                                           gTACGGGTGTAGACATTAAACC‑GGCTGGCAATAATGATTCCGGAGAATG‑CCGGCGGTAATGCGCCTGCCAGTACCATTATCTGCAGA‑TGTTCGCTGTTCAAATGACATGCCATACCAACGAGCAGCA‑GTGCCAGTGGCATCAGAATCAGATTTAGGAAGGGGGTATAAGCAAATTCGGAACTGAACTCGAATTTATGTGCATCCAGAGTCACACCCCCACCGAAAACCCC‑TC‑CCCCCGa                                                                                                                                                            >  2:861601/1‑239 (MQ=255)

                                                                            tACGGGTGTAGACATTAAACC‑GGCTGGCAATAATGATTCCGGAGAATG‑CCGGCGGTAATGCGCCTGCCAGTACCATCATCTGCAGA‑TGTTCGCTGTTCAAATGACATGCCATACCAACGAGAAGCA‑GTGCCAGTGGCATCAGAATCAGCTTCAGGAAGGGGGTATAAGCAATTTCGGCACTGAACTCGAATTTATGTGCAGCCCGAGTCAACCCCCCAGCGAATACCGC‑TA‑cccccggg                                                                                                                                                           >  2:687718/1‑235 (MQ=255)

                                                                            tACGGGTGTAGACATTAAACC‑GGCTGGCAATAATGATTCCGGAGAATG‑CCGGCGGTAATGCGCCTGCCAGTACCATCATCTGCAGA‑TGTTCGCTGTTCAAATGACATGCCATACCAACGAGAAGCA‑GTGCCAGTGGCATCAAAATCAGCTTCAGGAAGGGGTTATAAGCAATTTCGGCCCTGAACTCGAATTTATGTGCAGCCCGAGTCAACCCCCCCGCGAAAACCCC‑TA‑CCCCCGAg                                                                                                                                                           >  2:8457/1‑239 (MQ=255)

                                                                              ccgggttaagcaatAAACCGGGCTGGCAAAAAGGATTCCGGAGAATG‑CCGGCGGTAATGCGCCTGCCAGTACCATCATCTGCAGA‑TGTTCGCTGTTCAAATGACATGCCATACCAACGAGAAGCA‑GTGCCAGTGGCATCAGAATCAGCTTCAGGAAGGTGTTATAAGCAATTTCGGCACTGAACTCGAATTTATGTGCAGCCAGAGTCAACCCCGCAGCGAATACCGC‑TA‑CCCCTGAg                                                                                                                                                           <  2:456808/225‑1 (MQ=255)

                                                                                ggTGTAGACATTAAACC‑GGCTGGCAATAATGATTCCGGAGAATG‑CCGGCGGTAATGCGCCTGCCAGTACCATTATCTGCAGA‑TGTTCGCTGTTCAAATTACCTGCCATACCAACGAGAAGAA‑GTGCCAGTGGCATCAAAATAAGATTCAGGAAGGGGTTATAAACAATTTCGGCACTTCAGTCCAATTTATGTTAAGACAGAGTCACACCCCCCGCGCAAACCGC‑CA‑acccccagtaaa                                                                                                                                                       >  1:529742/1‑231 (MQ=255)

                                                                                ggTGTAGACATTAAACC‑GGCTGGCAATAATGATTCCGGAGAATG‑CCGGCGGTAATGCGCCTGCCAGTACCATCATCTGCAGA‑TGTTCGCTGTTCAAATGACATGCCATACCAACGAGAAGCC‑GTGCCAGTGGCATCAGAATCAGCTTCAGGGAGGTGTTATAAGCAATTTCGGCACTGAACTCGAATTTATGTGCAGCCAGAGTCAACCCCGCAGCCAATACCGC‑TA‑CCCCTTAGTTa                                                                                                                                                        >  2:346650/1‑238 (MQ=255)

                                                                                ggTGTAGACATTAAACC‑GGCTGGCAATAATGATTCCGGAGAATG‑CCGGCGGTAATGCGCCTGCCAGTACCATCATCTGCAGA‑TGTTCGCTGTTCAAATGACATGCCATACCAACGAGAAGCA‑GTGCCAGTGGCATCAGAATCAGCTTCAGGAAGGGGTTATAAGCAATTTCGGCACTGAACTCGAATTTATTGTCAGCCAGGAGCAAACCCGCAGCCGATAACGC‑CACCCCCTGAGtt                                                                                                                                                         >  2:548298/1‑238 (MQ=255)

                                                                                ggTGTAGACATTAAACC‑GGCTGGCAATAATGATTCCGGAGAATG‑CCGGCGGTAATGCGCCTGCCAGTACCATCATCTGCAGA‑TGTTCGCTGTTCAAATGACATGCCATACCAACGAGAAGCA‑GTGCCAGTGGCATCAGAATCAGCTTCAGGAAGGGGTTATAAGCAATTTCGGACCTGAACTCGAATTTATTTGCAGCCAAGGTCAACCCCCCCACGAAAACCCC‑AA‑CCCCTGAGttg                                                                                                                                                        >  2:677037/1‑237 (MQ=255)

                                                                                ggTGTAGACATTAAACC‑GGCTGGCAATAATGATTCCGGAGAATG‑CCGGCGGTAATGCGCCTGCCAGTACCATCATCTGCAGA‑TGTTCGCTGTTCAAATGACATGCCATACCAACGAGAAGCA‑ATGCCAGTGGCATCAGAATCAGCGTCAGGAAGGTGTTATAAGCAATTTCGGCACTGAACTCGAATTTATTGGCAGCCCGAGTGAACCCCCCAACCAAAAACCC‑AA‑CCCCTCAGTTa                                                                                                                                                        >  2:581681/1‑238 (MQ=255)

                                                                                ggTGTAGACATTAAACC‑GGCTGGCAATAATGATTCCGGAGAATG‑CCGGCGGGAATGCGCCTGCCAGTACAATCATCTGCAGA‑TGGTCGCTGTTTAAATGACATGCCATACCAACGAGCAGGA‑GTGCCAGGGGCATCAGAATCAGCTTCAGGAAGGTGTTATAAGCAATTTTGGCACTGAACTCGAATTTATATGCAGCCAGAGACAACACCCCAGCGAATACCCT‑CA‑CCCCTGGGTTa                                                                                                                                                        >  2:60321/1‑238 (MQ=255)

                                                                                ggTGTAGACATTAAACC‑GGCTGGCAATAATGATTACGGAGAATG‑CCCGCGTAAATTCGGCTGCCAGTACCACCATCTGTCGA‑TTTCCTCTGTTCAAATGACATGCACTACACAAGCGCAGAT‑GGGGCGCGTGCATACGAATAAGCTCCAGGAAGGTGGTATATGCGAGTTCGGCACTGGACCGGAATTTATGTTCAGCCAGAGTCAACCCCGCAGCgcaacacgcccccccagtgctaa                                                                                                                                                         >  1:868029/1‑217 (MQ=255)

                                                                                      gACATTAAACC‑GGCTGGCAATAATGATTCCGGAGAATG‑CCGGCGGTAATTCGCCTGCCCGTACCATCATCTGCAGA‑TGTTCGCCGTTCAAATGACCTTCCATACCAACGAGAAGCA‑GTGCCAGTGGCATACGAAACAGCTTCCGGAAAGGGTTATAAGCAATTTTGCCACCGAACCCGAATTTTTGCTCAGCCAAAGTCAACCCCCACACGAACACCGCATA‑ccccgcgttaacactc                                                                                                                                                   >  1:539187/1‑226 (MQ=255)

                                                                                         aTTAAACC‑GGCTGGCAATAATGATTCCGGAGAATG‑CCGGCGGTAATGCGCCTGCCAGTACCATCATCTGCAGA‑TGTTCGCTGTTCAAATGACATGCCATACCAACGAGAAGCA‑GTGCCAGAGGCATCAGAATCAGCTTCAGGAAGGGGTTATAAGCAATTTTGGCACTGAACTCGAATTTATTTGACGCCAAAGTCAAACCCCCACCGAATACCCC‑TA‑CCCCTCAGTTAGCCTTTGCaa                                                                                                                                              >  2:136302/1‑239 (MQ=255)

                                                                                          ttAAACC‑GGCTGGCAATAATGATTCCGGACAATG‑CCGGCCGTAAGGCGCCTGCCAGTACCATCATCGGCAGA‑TGTTCGCTGTTCAAATGACATGCCATACCAACGAGAAGCA‑GTGCCAGTGGCATCAGAATCAGCTTCAGGAAGGTGTTATAAGCAATTTCGGCACTGAACTCGAATTTATGTGCAGCCAGAGTCAACCCCGCAGCGAATACCGC‑TA‑CCCCTGAGTTAGCCTTCGCAAt                                                                                                                                             <  2:518393/239‑1 (MQ=255)

                                                                                           tAAACC‑GGCTGGCAATAATGATTCCGGAGAATG‑CCGGCGGGAATGCGCCTGCCCGTACCATCATCTGCAGA‑TGTTCGCTGTTCAAATGACATGCCATACCAACGAGAAGCA‑GTGCCAGTGGCATCAGAATCAGCTTCAGGAAGGTGGTTTAAGCAATTTCGGCCCTGAAATCGAATTTTTGTGCAGCCACAGTCAACCCCGCAGCGAATACCGC‑TC‑CCCCCGAGTTAGCCCTTGCacac                                                                                                                                            >  1:227588/1‑236 (MQ=255)

                                                                                                                      gagaATG‑CCGGCGGTAATGCGCCTGCCCGTACCATCATCTGCAGA‑TGTTCGCTGTTCAAATGACATGCCATACCAACGAGAAGCC‑GTGCCAGTGGCATCAGAATCAGCTTCAGGGAGGTGTTTTAAGCAATTTCGGCCCTGAACTCGCATTTATTGTCCGCCAGAGGCAAACCCCCCCCGGATACCGC‑AA‑CCCCTGAGTTAGCCTTTCCAATCCGCTTAAAGAGTGGGGccccccgtg                                                                                                                   >  1:720704/1‑232 (MQ=255)

                                                                                                                           tG‑CCGGCGGTAATGCGCCTGCCAGTACCATCATCTGCAGA‑TGTTCGCTGTTCAAATGACATGCCATACCAACGAGAAGCA‑GAGCCAGTGGCATCAGAATCAGCTTCAGGAAGGGGTTATAAGCAATTTCGGCACTGAACTCGAATTTATGTGCAGCCAGAGTCAACCCCGCAGCGAATACCGC‑TA‑CCCCTGAGTTAGCCTTCGCAATCAGATTAAAGGATGGGTCCCATGCCGCCGGaa                                                                                                             >  2:680767/1‑239 (MQ=255)

                                                                                                                           tG‑CCGGCGGTAATGCGCCTGCCAGGACCATCATCTGCAGA‑TGTTCGCTGTTCAAATGACATGCCATACCAACGAGAAGCA‑GTGCCAGTGGCATCAGAATCAGCTTCAGGAAGGTGTTATAAGCAATTTCGGGACTGAACTCGAATTTATGTGCAGCCCAAGGCAAACCCCCAGCAGATACCCC‑TA‑CCCCTGAGTTAGCCTTCGCAATCAGATTAAAAGGTGGGTGCCACGCTGCCGGaa                                                                                                             >  1:67856/1‑239 (MQ=255)

                                                                                                                                    gttgtGCCCCTGCCGGTCCCATCATCTGAGGA‑TTGTCGCTGTTAAAATGACAGGCCATACCAACGAGAACCA‑GTGCCAGTGGCATCAGAATCAGCTTCAGGAAGGTGTTATAAGCAATTTCGGCACTGAACTCGAATTTATGTGCAGCCAGAGTCAACCCCGCAGCGAATACCGC‑TA‑CCCCTGAGTTAGCCTTCGCAATCAGACTAAAGGTTGGGTCCCATGCTGCCGGAATTTTTAcc                                                                                                     <  1:766636/235‑1 (MQ=255)

                                                                                                                                     tatggcgcCTGCCAGTGCCATCATCTGCAGA‑TGTTCGCTGTTCAAATGACATGCCATACCAACGAGAAGCA‑GTGCCAGTGGCATCAGAATCAGCTTCAGGAAGGTGTTATAAGCAATTTCGGCACTGAACTCGAATTTATGTGCAGCCAGAGTCAACCCCGCAGCGAATACCGC‑TA‑CCCCTGAGTTAGCCTTCGCAATCAGATTAAAGGTTGGGTCCCATGCTGCCGGAATTTTTAcc                                                                                                     <  2:837222/234‑1 (MQ=255)

                                                                                                                                      tattccccTCCCAGTTCCATCATGTGTAGA‑TGTTCGCTGTTCAAATGACATGCAATACCAACTAGAAGCA‑GTGCCAGTGGCAGCAGAAGCCGCTTCAGGAAGGTGTTATAAGCCAGTCCGGCCCTGAACTCGAATTTATGTGCAGCCAGAGTCAACTCCGCAGCGAATACCGC‑TA‑CCCCTGAGTTAGCCTTCGCAATCAGATTAAAGGTTGGGTCCCATGCTGCCGGAATTTTTAcccc                                                                                                   <  1:426663/233‑1 (MQ=255)

                                                                                                                                           gcCTGCCAGTACCATCATCTGCAGA‑TGTTCGCTGTTCAAATGACATGCCATACCAACGAGAAGCA‑GTGCCAGTGGCATCAGAATCAGCTTCAGGAAGGTGTTATAAGCAATTTCGGCACTGAACTCGAATTTATGTGCAGCCAGAGTCAACCCCGCAGCGAATACCGC‑TA‑CCCCTGAGTTAGCCTTCGCAATCAGATTAAAGGTTGGGTCCCATGCTGCCGGAATTTTTACCccaacc                                                                                               <  2:41251/238‑1 (MQ=255)

                                                                                                                                            ccTGCCAGTACCATCATCTGCAGA‑TGTTCGCTGGTCAAATGACATGCCATACCAACGAGAAGCA‑GTGCCAGGGGCATCAGAATCAGCTTCAGGAAGGTGTTATAAGCAATTTCGGCACTGAACTCGAATTTATGTGCAGCCAGAGTCAAACCCGCCGCGAATACCGC‑TA‑CCCCTGAGTTAGCCTTCGCAAACAACTTAAAAGTTGGGGCCCCAGCGGCCGGAAATTTTACCccaaacaa                                                                                             >  1:107633/1‑239 (MQ=255)

                                                                                                                                                     gcccacatgggCAGA‑TGTTGGCTGTGAAAATGACATGCCATGCCAACGAGAAGCA‑GTGGCAGTGGCATCAGAATCAGCTTCAGGAAGGTGTTATAAGCAATTTCGGCACTGAACTCGAATTTATGTGCAGCCAGAGTCAACCCCGCAGCGAATACCGC‑TA‑CCCCTGAGTTAGCCTTCGCAATCAGATTAAAGGTTGGGTCCCATGCTGCCGGAATTTTTACCCCAACCAACACCAGGa                                                                                     <  1:346650/228‑1 (MQ=255)

                                                                                                                                                      ccacattctgaggA‑TGTTCGCTTTTCAAATGACCTGCCATACCAGCGAGAAGCA‑GTGCCAGTGGCATCCGAATCAGCTTCAGGAAGGTGTTATAAGCAATTTCGGCAATGAGCTCGAAGTTATGTGCAGCCAGAGTCAACCCCGCAGCGAATACCGC‑TA‑CCCCTGAGTTAGCCTTCGCAATCAGATTAAAGGTTGGGTCCCATGCTGCCGGAATTTTTACCCCAACCAACACCAGGATc                                                                                   <  1:282642/227‑1 (MQ=255)

                                                                                                                                                      ccATCAACTGCAGA‑GGTTCGCGGTTCAAATGCCATGCCATTCCGACGGGAAGCA‑GTGCCAGTGGCATCAGAATCAGCTTCAGGAAGGTGTTAAAAGCCATTTCGGCACTGAACTCGAATTTATGTTCAGCCAGAGGCAACCCCGCAGCGAATACCGC‑TA‑CCCCTGAGTTAGCCTTCGCAATCAGATTAAAGGTTGGGTCCCATGCTGCCGGAATTTTTACCCCAACCAACACCAGGATc                                                                                   <  1:373768/239‑1 (MQ=255)

                                                                                                                                                           cccttcAGG‑TGTTCCTTGTTAAAATGACCTGCCATACCAACGAGAAGCA‑GTGCCAGTGGCATCAGAATCAGCTTCAGGAAGGTGTTATAAGCAATTTCGGCACTGAACTCGAATTTATGTGCAGCCAGAGTCAACCCCGCAGCGAATACCGC‑TA‑CCCCTGAGTTAGCCTTCGCAATCAGATTAAAGGTTGGGTCCCATGCTGCCGGAATTTTTACCCCAACCAACACCAGGATCGTTGc                                                                              <  1:687718/234‑1 (MQ=255)

                                                                                                                                                                         cGCTGTTCAAATGACATGCCATACCAACGAGAAGCA‑GTGCCAGTGGCATCAGAATCAGCTTCAGGAAGGTGTTATAAGCAATTTCGGCACTGAACTCGAATTTATGTGAAGCCAGAGTCAACCCCGCAGCGAATACCGC‑TA‑CCCCTGAGTTACCCTTCGCAATCATATTAAAGCGTGGGCCCCATGCCGCCGCAATTTTTACCCAAACCAACACCCGGATCGttcaaaaaaaccgggc                                                                  >  1:692906/1‑224 (MQ=255)

                                                                                                                                                                                        aTGCCATACCAACGAGAAGCA‑GTGCCAGTGGCATCAGAATCAGCTTCAGGAAGGTGTTATAAGCAGTTTAGGCACGGAACTCTAATTTATGTGCTACCAGCGTCAACCCCGCACCGAATACCGC‑TA‑CCCCTCAGTTAGCCTTCGCATCCAGCTTAAAGTCTGGGCCCCACCCTGCCGGAATTTTTACCCCAACAAACACCATGACCTTTACACACACAGCGGCCCACACTACTgcac                                                    >  1:290017/1‑234 (MQ=255)

                                                                                                                                                                                        aTGCCATACCAACGAGAAGCA‑GTGCCAGTGGCATCAGAATCAGCTTCAGGAAGGTGTTATAAGCAATTTCGGCACTGAACTCGAATTTATGTGCAGCCAGAGTCAACCCCGCAGCGAATACCGC‑TA‑CCCCTGAGTTAGACTTCGCAATCAGATTAAAGGTTGGGTCCCATGCTGCCGGAATTTTTACCCCAACCAACACCAGGATCGTTGCCAGAAAAGGGGCCCaaaaaaaaggcc                                                    >  1:774599/1‑226 (MQ=255)

                                                                                                                                                                                               aCCAACGAGAAGCA‑GTGCCAGTGGCATAAGAATCAGCTTCAGGAAGGTGTTATAAGCAATTTCGGCACTGAACTCGAATTTATGTGCAGCCAGAGTCAACCCCGCAGCGAATACCGC‑TA‑CCCCTGAGTTAGCCTTCGCAATCATATTAAAGGTTGGGTCCCATGCGTCCGGACTTTTTTCCCCAACCAACACCAAGAACGTTTCCAGAAAAGGGGccccacctaccgggccccttaccg                                           >  1:656192/1‑218 (MQ=255)

                                                                                                                                                                                                    cGAGAAGCA‑GTGCCAGTGGCATCAGAATCAGCTTCAGGAAGGTGTTATAAGCAATTTCGGCACTGAACTCGAATTTATGTGCAGCCAGAGTCAACCCCGCAGCGAATACCGC‑TA‑CCCCTGAGTTAGCCTTCGCAATCAGATTAAAGGTTGGGTCCCATGCTGCCGGAATTTTTACCCCAACCAACACCAGGATCGTTGCCAGAACAGGTGCCCATACTACTGGCTCCTTTGCCGCAGaa                                      >  2:886659/1‑239 (MQ=255)

                                                                                                                                                                                                     gagaAGCA‑GTGCCAGTGGCATCAGAATCAGCTTCAGGAAGGTGTTATAAGCAATTTCGGCACTGAACTCGAATTTATGTGCAGCCAGAGTCAACCCCGCAGCGAATACCGC‑TA‑CCCCTGAGTTAGCCTTCGCAATCAGATTAAAGGTTGGGTCCCATGCTGCCGGAATTTTTACCCCAACCAAACCCAGAATCGTTGCCAGAACAGTGGCCACAACACCAGGCTCCTTTACCGAAGaaa                                     >  1:362986/1‑239 (MQ=255)

                                                                                                                                                                                                             cgTGCCATGGGCATCAAAATCAGTTTCAGGATGGTGTTATAAGCAATTTCGGCACTGAACTCGAATTTAGGTGCAGCCAGAGTCAACCCTGCAGCGAATACCGC‑TT‑CCCCTGAGTTAACCTTCGCAATCAGATTAAAGGCTGTGTCCCATGCTGCCGGAATTTTTACCCCAACCAACACCAGGATCGGTGCCAGAACAGGTGCCCATACTACTTGATCCTTTGCCGCAGAAATTAATgc                              <  1:215475/238‑1 (MQ=255)

                                                                                                                                                                                                               tGCCAGTGGCATCAGAATCAGCTTCAGGAAGGTGTTATAAGCAATTTCGGCACTGAACTCGAATTTATGTGCAGCCAGAGTCAACCCCGCAGCGAATACCGC‑TA‑CCCCTGAGTTAGCCTTCGCAATCAGATTAAAGGTTGGGTCCCATGCTGCCGGAATTTTTACCCCAACCAACACCAGGATCGTTGCCAGAACAGGTGCCCATACTACTGGCTCCTTTGCCGCAGAAATTAATgc                              >  2:244626/1‑237 (MQ=255)

                                                                                                                                                                                                                gCCAGTGGCCTCAAAATCATCTTCAGGAAGGGGTTATAAGCAATTGCGGCACTGAACTCGAATTTATGTGCAGCCAGAGTCAACCCCGCAGCGAATACCGC‑TA‑CCCCTGAGTTAGCCTTCGCAATCAGATTAAAGGTTGGGTCCCATGCTGCCGGAATTTTTACCCCAACCAACACCAGGATCGTTGCCAGAACAGGTGCCCATACTACTGGCTCCTTTGCCGCAGAAATTAATGCGCt                           <  1:677037/239‑1 (MQ=255)

                                                                                                                                                                                                                gCCAGTGGCATCAGAATCAGCTTCAGGAAGGTGTTATAAGCAATTTCGGCACTGAACTCGAATTTATGTGCAGCCAGAGTCAACCCCGCAGCGAATACCGC‑TA‑CCCCTGAGTTAGCCTTCGCAATCAGATTAAAGGTTGGGTCCCATGCTGCCGGAATTTTTACCCCAACCAACACCAGGATCGTTGCCAGAACAGGTGCCCATACTACTGGCTCCTTTGCCGCAGAAATTAATGCGCt                           <  2:774599/239‑1 (MQ=255)

                                                                                                                                                                                                                gCCAGTGGCATCAGAATCAGCTTCAGGAAGGTGTTATAAGCAATTTCGGCACTGAACTCGAATTTATGTGCAGCCAGAGTCAACCCCGCAGCGAATACCGC‑TA‑CCCCTGAGTTAGCCTTCGCAATCAGATTAAAGGTTGGGTCCCATGCTGCCGGAATTTTTACCCCAACCAACACCAGGATCGTTGCCAGAACAGGTGCCCATACTACTGGCTCCTTTGCCGCAGAAATTAATGCGCt                           <  2:432927/239‑1 (MQ=255)

                                                                                                                                                                                                                  cggtgctatccaaaacAGTTTGAGGAAGGGCTTATAAGCAATTTCGGCACTGAACCCGAATTTATGTGCAGCCATAGTCAACCCCGCAGCGAATGCCGC‑TA‑CCCCTGAGTTAGCCTTCGCAATCAGATTAAAGGTTGGGACCCATGCTGCCGGAATTTTTACCCCAACCAACACCAGGATCGTTGCCATAACAGGTGCCCATACTACTGGCTCCTTTGCCGCAGAAATTAATGCGCTCa                         <  1:861601/224‑1 (MQ=255)

                                                                                                                                                                                                                    gTGGCATCAGAATCAGCTTCAGGAAGGTGTTATAAGCAATTTCGGCACTGAACTCGAATTTATGTGCAGCCAGAGTCAACCCCGCAGCGAATACCGC‑TA‑CCCCTGAGTTAGCCTTCGCAATCAGATTAAAGGTTGGGTCCCATGCTGCCGGAATTTTTACCCCAACCAACACCAGGATCGTTGCCAGAACAGGGGCCCATACTCCTGGCCCCTTTGCCGCCGAAATTAATGCGCTCAGa                       >  2:8048/1‑239 (MQ=255)

                                                                                                                                                                                                                      ggCATCAGAATCAGCTTCAGGAAGGTGTTATAAGCAATTTCGGCACTGAACTCGAATTTATGTGCAGCCAGAGTCAACCCCGCAGCGAATACCGC‑TA‑CCCCTGAGTTAGCCTTCGCAATCAGATTAAAGGTTGGGTCCCATGCTGCCGGAATTTTTACCCCAACCAACACCAGGATCGTTGCCAGAACAGGTGCCCATACTACTGGCTCCTTTGCCGCAGAAATTAATGCGCTCAGAt                      <  2:529742/238‑1 (MQ=255)

                                                                                                                                                                                                                       gCATCAGAATCAGCTTCAGGAAGGTGTTATAAGCAATTTCGGCACTGAACTCGAATTTATGTGCAGCCAGAGTCAACCCCGCAGCGAATACCGC‑TA‑CCCCTGAGTTAGCCTTCGCAATCAGATTAAAGGTTGGGTCCCATGCTGCCGGAATTTTTACCCCAACCAACACCAGGATCGTTGCCAGAACAGGGGCCCATACTACTGGCTCCTATACCGCAAAAATTAATGCGCTCAGttt                     >  2:471429/1‑238 (MQ=255)

                                                                                                                                                                                                                          tCAGAATCAGCTTCAGGAAGGTGTTATAAGCAATTGCGGCACTGAACTCGAATTTATGTGCAGCCAGAGTCAACCCCGCAGCGAATACCGC‑TA‑CCCCTGAGTTAGCCTTCGCAATCAGATTAAAGGTTGGGTCCCATGCTGCCGGAATTTTTACCCCAACCAACACCAGGATCGTTGCCAGAACAGGTGCCCATACTACTGGCTCCTTTGCCGCAGAAATTAATGCGCTCAGATTACTa                 <  2:868029/239‑1 (MQ=255)

                                                                                                                                                                                                                                tCAGCTTCAGGAAGGTGTTATAAGCAATTTCGGCACTGAACTCGAATTTATGTGCAGCCAGAGTCAACCCCGCAGCGAATACCGC‑TA‑CCCCTGAGTTAGCCTTCGCAATCAGATTAAAGGTTGGGTCCCATGCTGCCGGAATTTTTACCCCAACCAACACCAGGATCGTTGCCAGAACAGGTGCCCATACTACTGGCTCCTTTGCCGCAGAAATTAATGCGCTCAGATTACTAttcttc           <  2:91308/239‑1 (MQ=255)

                                                                                                                                                                                                                                   gCTTCAGGAAGGTGTTATAAGCAATTTCGGCACTGAACTCGAATTTATGTGCAGCCAGAGTCAACCCCGCAGCGAATACCGC‑TA‑CCCCTGAGTTAGCCTTCGCAATCAGATTAAAGGTTGGGTCCCATGCTGCCGGAATTTTTACCCCAACCAACACCAGGATCGTTGCCAGAACAGGTGCCCATACTACTGGCTCCTTTGCCGCAGAAATTAATGCGCTCAGATTACTATTcttctt         >  2:546867/1‑238 (MQ=255)

                                                                                                                                                                                                                                   gCTTCAGGAAGGTGTTATAAGCAATTTCGGCACTGAACTCGAATTTATGTGCAGCCAGAGTCAACCCCGCAGCGAATACCGC‑TA‑CCCCTGAGTTAGCCTTCGCAATCAGATTAAAGGTTGGGTCCCATGCTGCCGGAATTTTTACCCCAACCAACACCAGGATCGTTGCCAGAACAGGTGCCCATACTACTGGCTCCTTTGCCGCAGAAATTAATGCGCTCAGATTACTATTcttctt         >  1:89580/1‑238 (MQ=255)

                                                                                                                                                                                                                                         gggaGGTGTTATAAACAATTTCGGCACTGAACTCGAACTTAGGTGCAGCCAGAGTCAACCCCGGAGCGAATACCGC‑TA‑CCCCTGAGTTAGCCGTCGCAATCAGAGTAAAGGTTGGGTCCCATGCTGCCGGAATTTTTACCCCAACCAACACCAGGATCGTTGCCAGAACAGGTGCCCATACTACTGGCTCCTTTGCCGCAGAAATTAATGCGCTCAGATTACTATTCTTCTTACCATcc  <  1:806171/236‑1 (MQ=255)

GTCGACTGACATAAATCCACAAGGGAGCCGTGACGACAAAGCCCAGTACGCTCACAGCCAATGACGCTGTA‑CCAGTACGGGTGTAGACATTAAACC‑GGCTGGCAATAATGATTCCGGAGAATG‑CCGGCGGTAATGCGCCTGCCAGTACCATCATCTGCAGA‑TGTTCGCTGTTCAAATGACATGCCATACCAACGAGAAGCA‑GTGCCAGTGGCATCAGAATCAGCTTCAGGAAGGTGTTATAAGCAATTTCGGCACTGAACTCGAATTTATGTGCAGCCAGAGTCAACCCCGCAGCGAATACCGC‑TA‑CCCCTGAGTTAGCCTTCGCAATCAGATTAAAGGTTGGGTCCCATGCTGCCGGAATTTTTACCCCAACCAACACCAGGATCGTTGCCAGAACAGGTGCCCATACTACTGGCTCCTTTGCCGCAGAAATTAATGCGCTCAGATTACTATTCTTCTTACCATCC  >  NC_000913/2489251‑2489717

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: `‑`