Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,004,132	0	G	A	30.9%	57.9 / 16.1	43	A16V (GCT→GTT)	rarD	putative transporter RarD
Reads supporting (aligned to +/- strand):  ref base G (9/20);  new base A (13/0);  total (23/20)
Fisher's exact test for biased strand distribution p-value = 2.25e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTG‑AAGTACGCTGGCGCTAT‑ACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTG‑CATCCATCGGGAATACTCATTTTTAATTGGTAA‑CAGCAG‑TTT‑ACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGG‑CAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGGCGCAGATAATCA  >  NC_000913/4003890‑4004380                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                              gATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTG‑AAGTACGCTGGCGCTAT‑ACCCCAAATAAAATAAGCGGCAaga                                                                                                                                                                                                                                                       <  1:441532/251‑1 (MQ=255) gATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTG‑AAGTACGCTGGCGCTAT‑ACCCCAAATAAAATAAGCGGCAag                                                                                                                                                                                                                                                        <  2:413105/250‑1 (MQ=255)      ttCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGGAAATCAACTTT‑AAGTAAGCTGGCGCTAT‑ACCCCAAAAAAAATAAACGGCAAGAGCg                                                                                                                                                                                                                                                    >  1:157610/1‑249 (MQ=255)        cACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAAGGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCTGGCACGTAGTAAATCAACCTG‑AAGTACGCTGGCGCTAT‑ACCCCAAATAAAATAAGCGGCAAGAGCGAGTa                                                                                                                                                                                                                                                <  2:142165/251‑1 (MQ=255)             cccAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCCCGGAGTAAATCAAATTG‑AAGTAAGCTGGCGCTAT‑ACCCCACATAAAAAAAGCGGCAAAAGCGCGTGATACg                                                                                                                                                                                                                                           >  1:428437/1‑251 (MQ=255)               cAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGGAGGAAAACAACCTG‑AAGGACGCTGGCGCGAT‑ACCCCAAATAAAATAAACGGCAAGAGCGGGTaaaaacac                                                                                                                                                                                                                                         >  1:154815/1‑245 (MQ=255)                    aaaCAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAAACAACTTG‑AAGTACGCTGGCGCTAT‑ACCCCAAATAAAATAAACGGCAAGAGCGAggaaaaccccccgcc                                                                                                                                                                                                                                    >  1:569645/1‑237 (MQ=255)                               aGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTG‑AAGTACGCTGGCGCTAT‑ACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGtttt                                                                                                                                                                                                                         <  1:553745/251‑1 (MQ=255)                                     caccaATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTG‑AAGTACGCTGGCGCTAT‑ACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTG‑catcc                                                                                                                                                                                                                  <  1:386864/251‑1 (MQ=255)                                          aTCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTG‑AAGTACGCTGGCGCTAT‑ACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTG‑CATCCATCgg                                                                                                                                                                                                             <  1:452509/251‑1 (MQ=255)                                            cAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTG‑AAGTACGCTGGCGCTAT‑ACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTG‑CATCCATCggg                                                                                                                                                                                                            <  2:517918/250‑1 (MQ=255)                                            cAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGGAGGAAATTAAACTG‑AAGTACGCTGGCGCTAT‑AACCCAAATAAAAAAAACGGCCAGAGAGAGTAATACCCCCCCCCCCGTTTGGTTTT‑CATCCATAgggg                                                                                                                                                                                                           >  1:350874/1‑250 (MQ=255)                                                       agaCTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGGAGTAAATCAACTTG‑AAGAACTCTGGCGCGTT‑ACCCCCAAAAAAAAAAACGGCCAGAGGGGGGAAAACGCCCCGCCGCGGGTGTTTTggaacccacggggataatcatttt                                                                                                                                                                                                 >  1:82844/1‑228 (MQ=255)                                                             ccAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGGAGTAAATCAACTTG‑AAGTACGCTGGCGCTAT‑ACCCCAAAAAAAAAAAACGGCAAGAGCGGGTAAAAAGCCCTGCCCCGTTTGGTTTG‑CATCCATCGGGGATACCCATTTTTAAttt                                                                                                                                                                                          >  1:321234/1‑250 (MQ=255)                                                               aGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGGAGTAAAACAACTTG‑TAGGACGCTGGCGCCAT‑ACCCCAAAAAAAAAAAAAGGCCAGAAAGAGTAATACGCCCCGCCCCGGTTGTTTTG‑CATCCAACggggagattccttttttatttgt                                                                                                                                                                                        >  2:94415/1‑231 (MQ=255)                                                                           ttCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGGAGTAAAACAACTTG‑AAGGGCGCTGGGGCTAT‑ACCCCCAATAAAATAAACGGCAAGAGCGGGGAAAACGCCCTGCCCCGTTTGTTTTT‑CCTCCCACGGGGATACTCATTTTTTATTTTTAA‑cagcaggttt                                                                                                                                                                           >  1:70071/1‑247 (MQ=255)                                                                                                  ttttAAATAGGCCCACTGGCGGCAAATGCTCATCAGCCCCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTG‑AAGTACGCTGGCGCTAT‑ACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTG‑CATCCATCGGGAATACTCATTTTTAATTGGTAA‑CAGCAG‑TTT‑ACCTGCTTTTATGTCTTCAACCt                                                                                                                                                   <  2:75909/250‑1 (MQ=255)                                                                                                          aGGACCACTGGCGGCAAATGCTCATCAGCCCCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGCTTTCATCGGCGGGCACGGAATAAATAAACTTT‑AAGTACGCTGGCGCTAT‑ACCCCCAAAAAAAAAAACGGGAAGAGGGAGGAAAACGCCCTGCCCCGCTTGTTTTT‑TATCCCTCGGGAATACTCCTTTTTTATTGGTAA‑AAACAG‑TTT‑TTCTGCtttttttttttttacataccctata                                                                                                                                           >  2:197356/1‑229 (MQ=255)                                                                                                                   ggcggcAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGGAGTAAATCAACTTT‑AAGTACGCTGGCGCTAT‑ACCCCCAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTTTTTTG‑TATCCATCGGGGATACTCATTTTTAATTGGGAA‑AAACAG‑TTT‑ACCTGCTTTTTTTTCTTCAACCTACCCAATAagtggcgggg                                                                                                                                 >  1:531479/1‑249 (MQ=255)                                                                                                                   ggcggcAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCAAGTACTAAATCAAACTA‑AAGAACGCTGGGGCTCT‑ACCCCCCAAAAAATAAACGGCAAGGGCGGGGAAAAAGCCCTCCCCCGCTTTTTTTT‑CATCCCTAGGGAATACACATTTTTTATTAGGAA‑AAGCAGGTTT‑AC‑TGCTTTTATGTTTTCacaccacactattgggggggggg                                                                                                                                 >  2:336665/1‑229 (MQ=255)                                                                                                                            tGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGAAAGAAAGAAAACAAATTG‑AAGAACGCTGGCGCCAT‑ACCCCCAATAAAAAAAAAGGCAAGAGCGGGGAAAACGCCCCGCCGCGTTTGTTTTG‑CATCCATCCGGGATAATCATTTTTTATTGGTAA‑CAGCAG‑TTT‑TCCTGCCGTTTTTTTTTCCACCTACCAtataaggtgcgggggcaccggc                                                                                                                         >  2:336940/1‑232 (MQ=255)                                                                                                                                catcaGCACCACCATAAAGAATAACGACCAGATCACGCGATGCGTCAGTATTTCATCGGCGGGCTCGTAGGAAATCAACTTG‑AAGTCCTCTGGCGCTAT‑TCCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTG‑CATCCATTGGGAATACTCCTTTTTAATTGGTAA‑CAGCAG‑TTT‑ACCTGCTTTTATGTCTTCCACTTACCATATAGGTGGCGGTGGCACTTG‑CAata                                                                                                                    <  2:523559/250‑1 (MQ=255)                                                                                                                                       accaccATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTG‑AAGTACGCTGGCGCTAT‑ACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTG‑CATCCATCGGGAATACTCATTTTTAATTGGTAA‑CAGCAG‑TTT‑ACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGG‑CAATATAAAGCTg                                                                                                            <  1:526970/251‑1 (MQ=255)                                                                                                                                            caTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGCAGTAAATCAACATG‑AAGTAAGCTGGCGCTAT‑ACCCCCAATAAAATAAACGGCAAAAGAGGGGAAAACGCCCTGCCGCGTTTGTTTTG‑CATACATCGGGGATACTCATTTTTTATTGGTAA‑CAACAG‑TTT‑ACCCGCTTTTATGTCTTCAACCTAACCTATAGGTGGCGGGGGCACTGG‑CCATATAAAACCGttctt                                                                                                       >  2:470484/1‑251 (MQ=255)                                                                                                                                                  gAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGGCATAAATCAAATTGGAAG‑ACGCCGGCGCTATTAACCCAAA‑AAAAAAAGCAGCAAGAGCGGGGAATACGCCCTCCCGCGGTTGTTTTG‑CATCCCTCGGGAATACTCATTTTTAATTTTGAA‑AAGCAG‑TTT‑ACCCTCTTTTTTGTCTTCAAACTACCATATAAGGGGGGGGGGGAATGG‑CAATATAAAGACGGTCTTCAtttt                                                                                                 >  2:509820/1‑249 (MQ=255)                                                                                                                                                         gACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCCCGGAGTAAATTAAATTG‑AAGAAAGCTGGGGCTAT‑ACCCCCAATAAAAAAAGAGGCAAGAGCGAGTAATAAGCCCTCCCGCGTTTGTTTTT‑GATCCATCCGGAATACTCATTTTTAATTGGGAA‑CATGAG‑ATT‑ACCTGCTTTTATGTTTTCAACCTACCCTATATGAGGCGGGGGCACCGG‑CAATATAAAGCTGTTCTTCTTTTTGTAATTc                                                                                          >  1:363186/1‑251 (MQ=255)                                                                                                                                                           ccAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTG‑AAGTACGCTGGCGCTAT‑ACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTG‑CATCCATCGGGAATACTCATTTTTAATTGGTAA‑CAGCAG‑TTT‑ACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGG‑CAATATAAAGCTGTTCTTCATTGTGTAATTCaa                                                                                        <  2:435483/251‑1 (MQ=255)                                                                                                                                                            cAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGCAGGAAAACAACTTC‑AAGTACGCCGGCGCCAT‑ACCCCCAATAAAATAAGCCGGAAAAGCGGGGAATAAGCCCCGCCGCGTTTGTTTTT‑TATCCATCCGGGATACTCATTTTTAATTGGTAAACA‑CAG‑GTT‑ACCTGCTTTTTTGTCTTCAAACTACCCTATAAGTGGCGGGGGCACTGGGCAA‑ATAAAGCTGttctttctttttgaattccac                                                                                       >  1:437153/1‑235 (MQ=255)                                                                                                                                                               aTCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGAAAAACAACTTG‑AAGTAAGCCGGCGCTAT‑ACCCCCAATAAAAAAAAAGGCAAGAGCGAGGAATACGCCCTCCCCCGGTTGTTTTTTCTTCC‑TCGGGGAGACTCATTTTTTATTGGGAA‑CAGCAG‑TTTTACC‑GCTTTTTTTTCTTCAAACCACCATATAAGGGGCGGGGGGCCTGGGAAATATAAAG‑TGttcttccttttgttatttaaagcg                                                                                    >  2:233403/1‑233 (MQ=255)                                                                                                                                                                    gcgATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTG‑AAGTACGCTGGCGCTAT‑ACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTG‑CATCCATCGGGAATACTCATTTTTAATTGGTAA‑CAGCAG‑TTT‑ACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGG‑CAATATAAAGCTGttctt                                                                                                       <  2:403016/227‑1 (MQ=255)                                                                                                                                                                          cGTCAGGGTGTCATCAGCGGGCACGTAGTAAATCAACTTG‑AGGTACGCGGGCGCTAG‑ACCGCAAATAAGATAAGCGGCAATAGCGAGTAATACGCCCTGCCGCGTTTGTTTTG‑CATCCATCGGGAATACTCATTTTTAATTGGGAA‑CAGCAG‑TTT‑ACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGG‑CAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGggcgacggcga                                                                         <  2:80025/251‑1 (MQ=255)                                                                                                                                                                             cAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTG‑AAGTACGCTGGCGCTAT‑ACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTG‑CATCCATCGGGAATACTCATTTTTAATTGGTAA‑CAGCAG‑TTT‑ACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGG‑CAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACt                                                                       <  1:475006/250‑1 (MQ=255)                                                                                                                                                                                gATTTCATCGGCGGGCACGTAGTAAATCAACTTG‑AAGTACGCTGGCGCTAT‑ACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTG‑CATCCATCGGGAATACTCATTTTTAATTGGTAA‑CAGCAG‑TTT‑ACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGG‑CAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTAt                                                                   <  1:366140/251‑1 (MQ=255)                                                                                                                                                                                gATTTCATCGGCGGGCACGTAGTAAATAAACCTG‑CAGTAAGCCGGCGCTAT‑TACCCAAATAAAATAACCGGCCAGAGCGGGGAAAACCCCCTCCCCCGTTTGTTTTT‑CATCCATCGGGAATACTCATTTTTTATTGGTAA‑CAGCAG‑TTT‑ACCTGCTTTTATGTCTTTAACCTACCATATAGGGGGCGGGGGCCCTGG‑CAATATCAAGCTGTTCTTTATTGTGTGATTCAACCCCGGCGGCCGCGGCTTTAt                                                                   >  1:493208/1‑251 (MQ=255)                                                                                                                                                                                                        aaTCAACTTG‑AAGTACGCTGGCGCTAT‑ACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTG‑CCTCCATCGGGAATACTCATTTTTAATTGGTAA‑CAGCAG‑TTT‑ACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGG‑CAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCTGGCGACGGCGACTTTATTGCCTGCACGCACCAGGctact                                             >  2:350693/1‑249 (MQ=255)                                                                                                                                                                                                                tG‑AAGTACGCTGGCGCTAT‑ACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTG‑CATCCATCGGGAATACTCATTTTTAATTGGTAA‑CAGCAG‑TTT‑ACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGG‑CAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAg                                     <  2:437153/249‑1 (MQ=255)                                                                                                                                                                                                                tG‑AAGTACGCTGGCGCTAT‑ACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCCCGTTTGTTTTG‑CATCCATCGGGAATACTCATATTTAAATGGTAA‑CAGCAG‑TTT‑ACCTGCTTTTATGTCTTCAACCTCACATATAGGTGGCGGTGGCACTGG‑CAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGAGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGt                                    <  2:291203/250‑1 (MQ=255)                                                                                                                                                                                                                        cGCTGGCGCTAT‑ACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTG‑CATCCATCGGGAATACTCATTTTTAATTGGTAA‑CAGCAG‑TTT‑ACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGG‑CAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCt                            <  1:484692/251‑1 (MQ=255)                                                                                                                                                                                                                         gCTGGCGCTAT‑ACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTG‑CATCCATCGGGAATACTCATTTTTAATTGGTAA‑CAGCAG‑TTT‑ACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGG‑CAATATAAAGCTGTTCTTCATTGTGTAATTCAAc                                                                                       >  2:566273/1‑191 (MQ=255)                                                                                                                                                                                                                         gCTGGCGCTAT‑ACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTG‑CATCCATCGGGAATACTCATTTTTAATTGGTAA‑CAGCAG‑TTT‑ACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGG‑CAATATAAAGCTGTTCTTCATTGTGTAATTCAAc                                                                                       <  1:566273/191‑1 (MQ=255)                                                                                                                                                                                                                                          aaatacaataAGCGGCACGAGCGAGTAATACGCCCTGCCGCGTTTGCTTTG‑CATCCATCGGGAATACTCATTTTTACTTGGTAA‑CAGCAG‑TTT‑ACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTTGCACTGG‑CCATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTCTTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGgcgc           >  2:440188/1‑251 (MQ=255)                                                                                                                                                                                                                                              aaaataaGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTG‑CATCCATCGGGAATACTCATTTTTAATTGGTAA‑CAGCAG‑TTT‑ACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGG‑CAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGGCGCAGAt       <  1:458204/251‑1 (MQ=255)                                                                                                                                                                                                                                                    aGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTG‑CATCCATCGGGAATACTCATTTTTAATTGGTAA‑CAGCAG‑TTT‑ACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGG‑CAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGGCGCAGATAATCa  >  2:66479/1‑250 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                              GATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTG‑AAGTACGCTGGCGCTAT‑ACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTG‑CATCCATCGGGAATACTCATTTTTAATTGGTAA‑CAGCAG‑TTT‑ACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGG‑CAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGGCGCAGATAATCA  >  NC_000913/4003890‑4004380

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