Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,336,897	0	T	C	37.8%	35.3 / 13.4	39	G45G (GGA→GGG)	adiC	arginine:agmatine antiporter
Reads supporting (aligned to +/- strand):  ref base T (15/8);  new base C (0/14);  total (15/24)
Fisher's exact test for biased strand distribution p-value = 7.00e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GATGGTTAATACCAATGGATCTTTCAGAATCGGGAAGAAGTAACTTAAATATCCTACGCCAATGACCACCATGGCGATATTGCCGATCCAGCAGGCCAGCCAGTAGAGGACGTTGGT‑TTG‑ATAACCGAGAAACGGG‑CCAAA‑G‑CAGCGGCGGGCGTAAGCGTAAGAACCACCA‑GGACTTGGGTCGAGGAACGACATTTTGGCGTATACCAT‑CGAGAGCCCCAGTGCACCGATAATCGTCACC‑AACC‑ATC‑CATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGTTAT  >  NC_000913/4336651‑4337139                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                    gATGGTTAATACCAATGGATCTTTCAGAATCGGGAAGAAGTAACTTAAATATCCTACGCCAATGACCACCATGGCGATATTGCCGATCCAGCAGGCCAGCCAGTAGAGGACGTTGGT‑TTG‑ATAACCGAGAAACGGG‑CCAAA‑G‑CAGCGGCGGGCGTAAGCGTAAGAACCACCA‑GGACTTGGGTCGAGGAACGACATTTTGGCGTATACCAT‑CGAGAGCCCCAGTGCACCGATAATCGTCACC‑AAcc‑atc‑cat                                                                                                                                                                                                                                                <  1:71349/251‑1 (MQ=255)     gTTAATACCAATGGATCTTTCAGAATCGGGAAGAAGTAACTTAAATATCCTACGCCAATGACCACCATGGCGATATTGCCGATCCAGCAGGCCAGCCAGTAGAGGACGTTGGT‑TTG‑ATAACCGAGAAACGGG‑CCAAA‑G‑CAGCGGCGGGCGTAAGCGTAAGAACCACCA‑GGACTTGGGTCGAGGAACGACATTTTGGCGTATACCAT‑CGAGAGCCCCAGTGCACCGATAATCGTCACC‑AACC‑ATC‑CATAAAt                                                                                                                                                                                                                                            <  2:378453/251‑1 (MQ=255)     gTTAATACCAATGGATCTTTCAGAATCGGGAAGAAGTAACTTAAATATCCTACGCCAATGACCACCATGGCGATATTGCCGATCCAGCAGGCCAGCCAGTAGAGGACGTTGGT‑TTG‑ATAACCGAGAAACGGG‑CCAAA‑G‑CAGCGGCGGGCGTAAGCGTAAGAACCACCA‑GGACTTGGGTCGAGGAACGACATTTTGGCGTATACCAT‑CGAGAGCCCCAGTGCACCGATAATCGTCACC‑AACC‑ATC‑CATAAAt                                                                                                                                                                                                                                            <  2:368900/251‑1 (MQ=255)     gTTAATACCAATGGATCTTTCAGAATCGGGAAGAAGTAACTTAAATATCCTACGCCAATGACCACCATGGCGATATTGCCGATCCAGCAGGCCAGCCAGTAGAGGACGTTGGT‑TTG‑ATAACCGAGAAACGGG‑CCAAA‑G‑CAGCGGCGGGCGTAAGCGTAAGAACCACCA‑GGACTTGGGTCGAGGAACGACATTTTGGCGTATACCAT‑CGAGAGCCCCAGTGCACCGATAATCGTCACC‑AACC‑ATC‑CATAAAt                                                                                                                                                                                                                                            >  1:220413/1‑251 (MQ=255)                                                  tatCCTACGCCAATGACCACCATGGCGATATTGCCGATCCAGCAGGCCAGCCAGTAGAGGACGTTGGT‑TTG‑ATAACCGAGAAACGGG‑CCAAA‑G‑CAGCGGCGGGCGTAAGCGTAAGAACCACCA‑GGACTTGGGTCGAGGAACGACATTTTGGCGTATACCAT‑CGAGAGCCCCAGTGCACCGATAATCGTCACC‑AACC‑ATC‑CATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGGAACAGAAAAac                                                                                                                                                                                               >  1:171848/1‑251 (MQ=255)                                                       cccgcaaatgcccacgagggggATATGGCCGAACAAGCAGGCCAGCAAGTAGAGGACGTTGGT‑TTG‑AAACCCGA‑AAACGGG‑GCAAA‑GGCGGGGGCGGGCGAAAGGGAAAGAACCACCA‑GGACTGGGGCGGAGGACCGACATTGTGGCGTATCCCAG‑CGAGAGCCCCAGGGCACCGATACGCGCCCCC‑AACC‑ATC‑CAAAAAGGACAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTg                                                                                                                                                                                          <  2:86913/230‑1 (MQ=255)                                                           ccAATGCCCCCCATGTCGATATGGCCGATCCAGCAGACCAGCCGGTAGAGGACGTTGTT‑TTA‑AAACCCGAAAAAGGGG‑CCAAA‑G‑CAGCGGCGGGCGTAAGCGTAAGACCCACCG‑GGACTTGGGTGGAGGAACGACTTTTTGGCGTATCCCAC‑CGAGACCCCCAGTGCCCAGATAATCTTCCCC‑ACCC‑ACC‑CAAAAATGACAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGAccc                                                                                                                                                                                      <  1:197456/251‑1 (MQ=255)                                                             cggatctaccaGGGCGATATTGCGAATCCAGCGGGCCAGCAGGTAGGGGACGTTGTT‑TGG‑AAAACCGAGAACGGGG‑CAAAA‑G‑CAGCGGGGGGGGAAGGGGTAAGAACCCCCA‑GGACTGGGGTGGAGAAACGACTTTTGGGCGTATACCAT‑GAAGACCCCCAGTGCACCGATAATCGCCCCC‑ACCC‑ACC‑CAAAAAGGGCAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCa                                                                                                                                                                                    <  1:2007/244‑1 (MQ=255)                                                             aaTGACCACCATGGCGATATTGCCGATCCAGCAGGCCAGCCAGTAGAGGACGTTGGT‑TTG‑ATAACCGAGAAACGGG‑CCAAA‑G‑CAGCGGCGGGCGTAAGCGTAAGAACCACCA‑GGACTTGGGTCGAGGAACGACATTTTGGCGTATACCAT‑CGAGAGCCCCAGTGCACCGATAATCGTCACC‑AACC‑ATC‑CATAACTGGCAATCCCGCCAGTAGCGGCCAGGTTTCGAGGTAAACGCAAAACACCTGAacccc                                                                                                                                                                                    >  2:168141/1‑250 (MQ=255)                                                              atgccccccagggCGATATGGCCGACCCAGCAGCCCAGCCAGTAGAGGACGTTGGT‑TTG‑ATACCCGAAAAACGGG‑CCAAA‑G‑CAGCGGCGGGCGTAAGCGAAAGAACCCCCA‑GGACTTGGGTCGAGGAACGACATTTTGGCGTATACCAT‑CGAGAGCCCCGGTCCACCGATAACCGTCCCC‑ACCC‑ACC‑CAAAAATGGCAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCat                                                                                                                                                                                   <  1:289230/240‑1 (MQ=255)                                                                    accaTGGCGATATTGCCGATCCAGCAGGCCAGCCAGTAGAGGACGTTGGT‑TTG‑ATAACCGAGAAACGGG‑CCAAA‑G‑CAGCGGCGGGCGTAAGCGTAAGAACCACCA‑GGACTTGGGTCGAGGCACGACATTTTGGCGTATACCAT‑CGAGAGCCCCAGTGCACCGATAATCGTCACC‑AACC‑ATC‑CATAAATGGCAATCCCGCCAGAAGAGGCCAGGTTTGCAGGTAAAAGAAAAACACCTAACCCCa                                                                                                                                                                                    >  2:2007/1‑244 (MQ=255)                                                                           ggttagtGCCGACCAAGCGGGCCAGCCAGTAGGGGTCGTTGTT‑TAG‑ATAACCGAGAACGGGG‑CAAAA‑G‑CAGGGGCGGGCGTAAGCGTAAGAACCCCCA‑GGTCTGGGGGCGAGGACCGACATTTGGGTGTATCCCAT‑CGAGAGCCCCAGTGCCCCGATAATCGCCACC‑AACC‑ATC‑AATAAAGGGAAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGac                                                                                                                                                                      <  2:208522/245‑1 (MQ=255)                                                                                 tcccaacccAGTAGGCCAGCAATTAGAGGGCGT‑GGTGTGG‑ATAACAGAAAAAGGGG‑CCAAA‑G‑CAGGGGGGGGCTAAAGCGTAAAACCCACCA‑GGACTGGGGTCGAGAAACGACATTTTGGCGTAAACAAG‑CGAGGGCCCCAGCGCCCCAAAAATCGTCACC‑AACC‑ATC‑AAAAAATGGAAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCAt                                                                                                                                                                 <  1:144076/243‑1 (MQ=255)                                                                                 gcccgccccaccggcccagcaagTAGAGGGCGTTGGT‑TGG‑AAAACCGA‑AAACGGGGCAAAA‑G‑CGGCGGGGGGCGTAGGCGAAAGACCCACCA‑GGACTTGGGCAGAGGAACGACTTTTTGGCGTATCCCAT‑CGAGCGCCCCAGTGCACCGATAATCG‑CCCC‑AACCCACC‑CATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATc                                                                                                                                                                <  2:220413/236‑1 (MQ=255)                                                                                        ccAGCAGGCCAGCCAGTAGAGGACGTTGGT‑TTG‑ATAACCGAGAAACGGG‑CCAAA‑G‑CAGCGGCGGGCGTAAGCGTAAGAACCACCA‑GGACTTGGGTCGAGGAACGACATTTTGGCGAATACCAC‑CGAGAGCCCCAGTGCACCGATAATCGTCACC‑AACC‑ATC‑CATAAATTGCAATCCCGCCAGTAGAAGCCAGGGTTGCAGGGAACAGAAAAAAAACTGACCCCCTAATATTCACCGACACCactcgagtg                                                                                                                                                          >  2:356421/1‑242 (MQ=255)                                                                                         cagcagCCCGGCAAGTAGAGGA‑GGTGGT‑TTGGAAACCCGAAAAACGGC‑CAAAA‑G‑CAGCGGCGGGCGTAAGCGAAAGAAC‑ACCAAGGCCTTGGTTCGAGGACCGACTTTTTGGCGTATCCCAT‑CGAGAGCCCCAGTGCCCCAATAATCGTCCCC‑AACC‑ATC‑CAAAAATGACAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGAc                                                                                                                                                        <  1:283033/251‑1 (MQ=255)                                                                                                      agTAGAGGACGTTGGT‑TTG‑ATAACCGAGAAACGGG‑ACAAA‑G‑CAGCGGCGGGCGTAAGCGTAAGAACCACCA‑GGACTTGGGTCGAGGAACGACATTTTGGCGTATACCAT‑CGAGCGCCCCAGAGCACCGATAATCGTCAAC‑AACC‑ATC‑CATACAAGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGAAACAGAAAAAAAACTGACCCCATAATATTCCCCGGCACCCTACGGGGGACGGGGTTTAAGccc                                                                                                                                           >  2:333373/1‑251 (MQ=255)                                                                                                       ggagaggcggTTGGT‑TTG‑AAACCCGAAAACCGGG‑CAAAA‑G‑CAGGGGGGGGCGTAAGCGAAAGACCCCCCA‑GGACTGGGGTCGAGGAACGATATTTGGGCGTATACCAT‑CGAGACCCCCAGTGCACCGAAAATCGTCCCC‑ACCC‑ACC‑CAAAAAGGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCa                                                                                                                                          <  1:250484/242‑1 (MQ=255)                                                                                                         ggaggggagtttgtttagaaaaCCGAGAA‑CGGG‑CCAAAAG‑CAGAGGCGGGGAAAAGCGAAAGACCCCCCA‑GGACTGGGGTGGAGAAACGACTTTTTGCCGTATACAGA‑CGAGCGCCCCTGTCCACGAATAAGCCACCCA‑AACC‑AAC‑CAAAAGTGGCAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACGCCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCa                                                                                                                                          <  2:243711/231‑1 (MQ=255)                                                                                                         agagGACGTTGGT‑TTG‑ATAACCGAGAAACGGG‑CCAAA‑G‑CAGCGGCGGGCGTAAGCGTAAGAACCACCA‑GGACTTGGGTCGAGGAACGACATTTTGGCGTATACCAT‑CGAGAGCCCCAGTGCACCGATAATCGTCACC‑AACC‑ATC‑CATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACt                                                                                                                                        >  1:137779/1‑251 (MQ=255)                                                                                                                   ggT‑TTG‑ATAACCGAGAAACGGG‑CCAAA‑G‑CAGCGGCGGGCGTAAGCGTAAGAACCACCA‑GGACTTGGGGCGAGGAACGACATTTTGGCGTATACCAT‑CGCGAGCCCCAGTGAACCGATAATCGTCACC‑AACC‑ATC‑CATAAATGGCAATCCCGACAGTAGAGGCAAGGGTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGAAACCATCAGGGTGACGGGGATTAAGCCCACGTTGTGaggat                                                                                                                              >  2:149476/1‑248 (MQ=255)                                                                                                                                      ggggg‑CCAAA‑G‑CAGGGGGGGGCGAAAGCGAAAGACCCCCCG‑GGATTTGGGTCGAGAAACGCAATTTTGGCTTATCCCAG‑CGAGCGCCCCAGTCCACCGAAAATCGTCCCC‑ACCC‑ACC‑CATAAATGGCAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATa                                                                                                             <  2:245924/249‑1 (MQ=255)                                                                                                                                        gaggggcaatacgagggcgggggataagcgtaagACCCACCA‑GGGCTGGGGGGGAGGAACGCTATTTTGGTTTCT‑CCATACGAGAGCCCTGGTGCCCCGAAAATCGCCACC‑AACC‑‑TCTAAAAAATGGCAACCCCGCCAGGAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACAt                                                                                                              <  2:150194/227‑1 (MQ=255)                                                                                                                                          ggcaaaa‑G‑AAGCGGCGGGCGTAAGCTTAAGACCCACCA‑GGACTTGGGTCGAGAAACGCCTTTTTGGCTTATCCCAT‑GGAGACCCCCGGTGCCCCGATAATCGTCACC‑ACCC‑ACC‑CATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAAtt                                                                                                          <  1:37857/247‑1 (MQ=255)                                                                                                                                                a‑g‑cagcGGCGGGCGTAAGCGTAAGAACCACCA‑GGACTTGGGTCGAGGAACGACATTTTGGCGTATACCAT‑CGAGAGCCCCAGTGCACCGATAATCGTCACC‑AACC‑ATC‑CATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTc                                                                                                    >  1:171329/1‑251 (MQ=255)                                                                                                                                                             ggataagcgtaagACACACCA‑GCACTTGGTTAAAGGACCGACATTTGGGTATATCCCAT‑CGCGGGCCCCAGTGCACCTAAAATCATCACA‑ACCC‑ACC‑CATAAATGACAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCAGCCGAAGACATAATTAAACTCCTGCGAAGGCg                                                                                         <  1:51844/248‑1 (MQ=255)                                                                                                                                                                          aaCCACCA‑GGACTTGGGTCGAGGAACGACATTTTGGCGTATACCAT‑CGAGAGCCCCAGTGCACCGATAATCGTCACC‑AACC‑ATC‑CATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACa                                                                             >  2:372976/1‑250 (MQ=255)                                                                                                                                                                                           gTCGAGGAACGACATTTTGGCGTATACCAT‑CGAGAGCCCCAGTGCACCGATAATCGTCACC‑AACC‑ATC‑CATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGtcatatc                                                            >  1:16227/1‑251 (MQ=255)                                                                                                                                                                                             ggAGGAACGACATTTTGGCGTAACCCAG‑CGAGCGCCCCAGTGCACCAAAAACCG‑CACCCAACC‑ACC‑CAAAAAAGGCAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGtcatatcat                                                          <  2:28712/250‑1 (MQ=255)                                                                                                                                                                                               aGGAACGACATTTTGGCGTATACCAT‑CGAGAGCCCCAGTGCACCGATAATCGTCACC‑AACC‑ATC‑CATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATAc                                                        >  1:229530/1‑251 (MQ=255)                                                                                                                                                                                                          tttggctaaaaCCAA‑CGAAAGCCCCAGTCCCCCAATAACCGCCACA‑AACC‑ACC‑CAAAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTc                                              <  2:248181/240‑1 (MQ=255)                                                                                                                                                                                                           ttgggtgtatccccccgaggcccccAGTGCCCCGATAATCGTACCC‑ACCC‑ACA‑CAAAAAGGGAAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTc                                              <  2:166784/229‑1 (MQ=255)                                                                                                                                                                                                           ttGGCGTATACCAT‑CGAGAGCCCCAGTGCACCGATAATCGTCACC‑AACC‑ATC‑CATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGAGACAATAAATTACGTCATATCATACGCCTGCATTCg                                             >  2:124413/1‑250 (MQ=255)                                                                                                                                                                                                            gggCGTATACCAG‑CGAGACCCCCAGTGCCCAGATACGCGTCCCC‑ACCC‑CAC‑AATAAATGGAAAACCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGa                                            <  1:166621/249‑1 (MQ=255)                                                                                                                                                                                                                     ccAT‑CGAGAGCCCCAGTGCCCCGATAAGCGTCACC‑AACC‑ACC‑CAAAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAatt                                  <  1:138267/251‑1 (MQ=255)                                                                                                                                                                                                                                       gCACCGATAATCGTCACC‑AACC‑ATC‑CATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTAt                 >  2:269707/1‑251 (MQ=255)                                                                                                                                                                                                                                            gATAATCGACACC‑AACC‑ACC‑CAAAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTg            <  1:125770/251‑1 (MQ=255)                                                                                                                                                                                                                                               aaTCGCCACC‑AACC‑ACC‑CATAAAGGCCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATg         <  1:86506/251‑1 (MQ=255)                                                                                                                                                                                                                                                   gTCACC‑AACC‑ATC‑CATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGt     >  2:263534/1‑251 (MQ=255)                                                                                                                                                                                                                                                      ccc‑aacc‑aTC‑AATAAATGCCACACCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACGCCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAAGTATTATCGACCGGGTTATCGCTGATGCCGttat  <  1:193940/250‑1 (MQ=255)                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                    GATGGTTAATACCAATGGATCTTTCAGAATCGGGAAGAAGTAACTTAAATATCCTACGCCAATGACCACCATGGCGATATTGCCGATCCAGCAGGCCAGCCAGTAGAGGACGTTGGT‑TTG‑ATAACCGAGAAACGGG‑CCAAA‑G‑CAGCGGCGGGCGTAAGCGTAAGAACCACCA‑GGACTTGGGTCGAGGAACGACATTTTGGCGTATACCAT‑CGAGAGCCCCAGTGCACCGATAATCGTCACC‑AACC‑ATC‑CATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGTTAT  >  NC_000913/4336651‑4337139

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