Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,336,901	0	T	A	29.0%	101.6 / 21.2	62	Y44F (TAT→TTT)	adiC	arginine:agmatine antiporter
Reads supporting (aligned to +/- strand):  ref base T (35/9);  new base A (0/18);  total (35/27)
Fisher's exact test for biased strand distribution p-value = 2.53e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TGGTTAATACCAATGGATCTTTCAGAATCGGGAAGAAGTAACTTAAATATCCTACGCCAATGACCACCATGGCGATATTGCCGATCC‑AGCAGGCCAGCCAGTAGAGGA‑CGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACC‑AGGACTTGGGTCGAGGAACGACATT‑TTG‑GCGT‑ATACCATCGAGAGCCCCAGTGCACCGATAATCGTC‑ACCAACCATCCATAAATGG‑CAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGTT  >  NC_000913/4336653‑4337137                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                              tGGTTAATACCAATGGATCTTTCAGAATCGGGAAGAAGTAACTTAAATATCCTACGCCAATGACCACCATGGCGATATTGCCGATCC‑AGCAGGCCAGCCAGTAGAGGA‑CGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACC‑AGGACTTGGGTCGAGGAACGACATT‑TTG‑GCGT‑ATACCATCGAGAGCCCCAGTGCACCGATAATCGTC‑ACCAAACATCCATa                                                                                                                                                                                                                                              >  2:311697/1‑250 (MQ=255)     tAATACCAATGGATCTTTCAGAATCGGGAAGAAGTAACTTAAATATCCTACGCCAATGACCACCATGGCGATATTGCCGATCC‑AGCAGGCCAGCCAGTAGAGGA‑CGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACC‑AGGACTTGGGTCGAGGAACGACATT‑TTG‑GCGT‑ATACCATCGAGAGCCACAGTGCACCGATAATCGTC‑ACCAACCATCCATAAATg                                                                                                                                                                                                                                          <  2:214403/250‑1 (MQ=255)      aaTACCAATGGATCTTTCAGAATCGGGAAGAAGTAACTTAAATATCCTACGCCAATGACAACCATGGCGATATTGCCGATCC‑AGCAGGCCAGCCAGTAGAGGA‑CGATGGTTTGATAACCGAGAAACGGGCCACAGCAGCGGCGGGGGGAAGCGGAAAAAACAAC‑AGAACCTGGGGCGAGGGACGACATT‑TTG‑G‑GTGAAAACCTCGAGAGCCCCAGGGCACCGATAACAGTA‑ACAAACCATCCATAACTGG‑c                                                                                                                                                                                                                                       >  2:377470/1‑251 (MQ=255)         aCCAATGGATCTTTCAGAATCGGGAAGAAGTAACTTAAATATCCTACGCCAATGACCACCATGGCGATATTGCCGATCC‑AGCAGGCCAGCCAGTAGAGGA‑CGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACC‑AGGACTTGGGTCGAGGAACGACATT‑TTG‑GCGT‑ATACCATCGAGAGCCCCAGTGCACCGATAATCGTC‑ACCAACCATCCATAAATGG‑CAAt                                                                                                                                                                                                                                    <  1:441896/251‑1 (MQ=255)         aCCAATGGATCTTTCAGAATCGGGAAGAAGTAACTTAAATATCCTACGCCAATGACCACCATGGCGATATTGCCGATCC‑AGCAGGCCAGCCAGTAGAGGA‑CGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACC‑AGGACTTGGGTCGAGGAACGACATT‑TTG‑GCGT‑ATACCATCGAGAGCCCCAGTGCACCGATAATCGTC‑ACCAACCATCCATAAAAGG‑Caa                                                                                                                                                                                                                                     >  2:40492/1‑250 (MQ=255)              tGGATCTTTCAGAATCGGGAAGAAGTAACTTAAATATCCTACGCCAATGACCACCATGGCGATATTGCCGATCC‑AGCAGGCCAGCCAGTAGAGGA‑CGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACC‑AGGACTTGGGTCGAGGAACGACATT‑TTG‑GCGT‑ATACCATCGAGAGCCCCAGTGCACCGATAATCGTC‑ACCAACCATCCATAAATGG‑CAAtggac                                                                                                                                                                                                                                >  1:172094/1‑246 (MQ=255)                        aGAATCGGGAAGAAGTAACTTAAATATCCTACGCCAATGACCACCATGGCGATATTGCCGATCC‑AGCAGGCCAGCCAGTAGAGGA‑CGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACC‑AGGACTTGGGTCGAGGAACGACATT‑TTG‑GCGT‑ATACCATCGAGAGCCCCAGTGCACCGATAATCGTC‑ACCAACCATCCATAAATGG‑CAATCCCGCCAGTAGAGGc                                                                                                                                                                                                                     >  1:157706/1‑251 (MQ=255)                             gggaaaaaagtacttaaaaTACCCTACGCCAATGACCCCCTGGGCGATAGGGCC‑ATCCCGGCAGGCCAGCGAGGAGAGGA‑CGTTGGTTTGATACCCGAGAAACGGGCCAAAGCAGGGGCGGGCGTAAGCGTAAAACCCCCC‑AGGACTTGGGTCGAGGACCGACATT‑TTG‑GCGT‑ATCCCAGCGAGAGCCCCAGGGCCCCGAAAATCGCC‑ACCACCCATCCAAAAATGG‑CAATCCCGCCAGTAGAGGCCAGGt                                                                                                                                                                                                                <  2:157706/235‑1 (MQ=255)                             cGGGAAGAAGTAACTTAAATATCCTACGCCAATGACCACCATGGCGATATTGCCGATCC‑AGCAGGCCAGCCAGTAGAGGA‑CGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTATGAACCACC‑AGGACTTTGGTCGAGGAACGACATT‑TTG‑GCGA‑ATACCATCAAGAGCCCCAGCGCACCGCTAATCGTC‑ACCAACCAAACATAAATGG‑CAAACCCGCCAGTAGAGGcccggg                                                                                                                                                                                                                >  2:57924/1‑247 (MQ=255)                                   gaagTAACTTAAATATCCTACGCCAATGACCACCATGGCGATATTGCCGATCC‑AGCAGGCCAGCCAGTAGAGGA‑CGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACC‑AGGACTTGGGTCGAGGAACGACATT‑TTG‑GCGT‑ATACCATCGAGAGCCCCAGTGCACCGATAATCGTC‑ACCAACCATCCATAAATGG‑CAATCCCGCCAGTAGAGGCCAGGTTTGc                                                                                                                                                                                                            >  2:596227/1‑249 (MQ=255)                                    aagTAACTTAAATATCCTACGCCAATGACCACCATGGCGATATTGCCGATCC‑AGCAGGCCAGCCAGTAGAGGA‑CGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACC‑AGGACTTGGGTAGAGGAACGACATT‑TTG‑GCGT‑ATACCATCGAGAGCCCCAGTGCACCGAGAATCGTC‑AACAAACATCCATAAATGG‑CAATCCCCCCAGAAGAGGCCAGGTTTGCAgg                                                                                                                                                                                                         >  2:87277/1‑251 (MQ=255)                                      ntAACTTAAATATCCTACGCCAATGACCACCATGGCGATATTGCCGATCC‑AGCAGGCCAGCCAGTAGAGGA‑CGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACC‑AGGACTTGGGTCGAGGAACGACATT‑TTG‑GCGT‑ATACCATCGAGAGCCCCAGTGCACCGCTAATCGTC‑ACCAACCATCCATAAATGG‑CAATCCCGCCAGTAAAGGCCAGGTTTGCAGGt                                                                                                                                                                                                        >  2:277091/2‑250 (MQ=255)                                        aaCTTAAATATCCTCCGCCAATGACCACCATGGCGATATTGCCGATCC‑AGCAGGCCAGCCAGTAGAGGA‑CGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACC‑AGGACTTGGGTCGAGGAACGACATT‑TTG‑GCGT‑ATACCATCGAGAGCCCCAGTGCACCGATAATCGTC‑ACCAACCATCCATAAATGG‑Caa                                                                                                                                                                                                                                     <  2:109102/219‑1 (MQ=255)                                        aaCTTAAATATCCTACGCCAATGACCACCATGGCGATATTGCCGATCC‑AGCAGGCCAGCCAGTAGAGGA‑CGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACC‑AGGACTTGGGTCGAGGAACGACATT‑TTG‑GCGT‑ATACCATCGAGAGCCCCAGTGCACCGATAATCGTC‑ACCAACCATCCATAAATGG‑Caa                                                                                                                                                                                                                                     >  1:109102/1‑219 (MQ=255)                                           ttAAATATCCTACGCCAATGACCACCATGGCGATATTGCCGATCC‑AGCAGGCCAGCCAGTAGAGGA‑CGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACC‑AGGACTTGGGTCGAGGAACGACATT‑TTG‑GCGT‑ATACCATCGAGAGCCCCAGTGCACCGATAATCGTC‑ACCAACCATCCATAAATGG‑CAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAAc                                                                                                                                                                                                     >  1:535003/1‑248 (MQ=255)                                             aaaTATCCTACGCCAATGACCACCATGGCGATATTGCCGATCC‑AGCAGGCCAGCCAGTAGAGGA‑CGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACC‑AGGACTTGGGTCGAGGAACGACATT‑TTG‑GCGT‑ATACCATCGAGAGCCCCAGTGCACCGATAATCGTC‑ACCAACCATCCCTAAAAGG‑CCATCCCGCCAGGAGAGGCCAGGTTTTCAGGTAACaaaa                                                                                                                                                                                                 >  2:248854/1‑247 (MQ=255)                                                     tACGCCAATGACCACCATGGCGATATTGCCGATCC‑AGCAGGCCAGCCAGTAGAGGA‑CGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACC‑AGGACTTGGGTCGAGGAACGACATT‑TTG‑GCGT‑ATACCATCGAGAGCCCCAGTGCACCGATAATCGTC‑ACCAACCATCCATAAATGG‑CAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTg                                                                                                                                                                                        >  1:286761/1‑251 (MQ=255)                                                          aaaTGACCACCAGGGCGATATTGCCGATCC‑AGCAGGCCGGCCAGGAGAGGA‑CGTTGGTTTGATAACCGAAAAAGGGGCAAAAGCAGGGGCGGGCGTAAGGGAAAGACCCACC‑GGGACTGGGGTCGAGGAACGACATT‑TTG‑GCGT‑ATCCCATCGAGAGCCCCAGTGCCCCGATAATCGTC‑ACCAACCATCCAAAAATGG‑CAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGAcccc                                                                                                                                                                                   <  1:250518/250‑1 (MQ=255)                                                               accaccATGGCGATATTGCCGATCC‑AGCAGGCCAGCCAGTAGAGGA‑CGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACC‑AGGACTTGGGTCGAGGAACGACATT‑TTG‑GCGT‑ATACCATCGAGAGCCCCAGTGCACCGATAATCGTC‑ACCAACCATCCATAAATGG‑CAATCCCGCCAGAAGAGGCCAGGTTTGCAGGCAACAGAAAAACACCTGACCCCataat                                                                                                                                                                              >  2:175754/1‑251 (MQ=255)                                                                       ggCGATATTGCCGATCC‑AGCAGGCCAGCCAGTAGAGGA‑CGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACC‑AGGACTTTGGGCGAGGAACGACATT‑TTG‑GCGT‑ATACCATCGAGAGCCCCAGTGCACCGATAATCGTC‑ACCAACCATCCATAAATGG‑CAATCCCGCCCGTAGAGGCCAGGGTTTCAGGGAACATAAAAACACCTGACCCCATAATATTCACCg                                                                                                                                                                      >  2:315892/1‑251 (MQ=255)                                                                                  ccgacccgccagGCCAGCCAG‑AGGGGAACGTTGGTTTAATAACCGAAAAACGGCCAAAAGCGGCGGGGGGCGTAAGCGTAAGACCCACC‑AGGACTGGGGTGGAGGAACGAAATT‑TTG‑GCTT‑ATACCATCAAGGGCCCCAGTGCACGGATAATCGCC‑CCCACCCATCCAAAAATGG‑CAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAg                                                                                                                                                            <  1:363050/242‑1 (MQ=255)                                                                                         agcagcCCAGCCGGTAGAGGT‑CGTGGGTTTTATAACCGAGAACCGGGCTAAAGCAGCGGGGGGGGTAAGCGAAAAACCCCCC‑AGGACTTGGGTCGGGGAACGACATT‑TTG‑GCGA‑ATACCATCGAGACCCCCGGGGCCCCAAAAATCGAC‑CCCACCCACCCAAAAATGG‑CAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACg                                                                                                                                                     <  1:160387/251‑1 (MQ=255)                                                                                           cagGCCAGCCAGTAGAGGA‑CGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACC‑AGGACTTGGGTCGAGGAACGACCTT‑TTG‑GCGT‑ATACCATCGAGAGCCCCCGTGCACCGATAATCGTC‑ACCAACCAGCCATAAATGG‑CAATCCCGCCACTAGAGGCCAGGATTGCCGGTAACAGAAAAACACCTGAAACCATAATATTCACCGGCACCATCCGGGTGACggg                                                                                                                                                   >  2:188534/1‑251 (MQ=255)                                                                                             gcccagccaaGAGAGGA‑TGTTTGTTGGATACCCAAGAAACGGACAAAAGCGGCGGCGGGCGTAAGCTTAAGACCCACA‑AGGCCGTGGTTCGAGAAACGCCTTT‑TTG‑CCTT‑AACCCACCGAGAGCCCCGGTGCCCCAATAATCGTC‑CCCAACCACCAAAAAATGG‑AAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGa                                                                                                                                                 <  2:47087/249‑1 (MQ=255)                                                                                                      gTAGAGGA‑CGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACC‑AGGACTTGGGTCGAGGAACGACATT‑TTG‑GCGT‑ATACCATCGAGAGCCCCAGTGCACCGATAATCGTC‑ACCAACCATCCATAAATGG‑CAATCCCGCCAGTAGAGGCCAGGTTTGc                                                                                                                                                                                                            >  1:560871/1‑183 (MQ=255)                                                                                                       tAGAGGA‑CGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACC‑AGGACTTGGGTCGAGGAACGACATT‑TTG‑GCGT‑ATACCATCGAGAGCCCCAGTGCACCGATAATCGTC‑ACCAACCATCCATAAATGG‑CAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTACCCGACACCATCAGGGTGACGGGGATTAAGCCCAc                                                                                                                                       >  1:534178/1‑251 (MQ=255)                                                                                                        agagGA‑CGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGGGTAAGTGTAAGAACCACC‑AGGACTTGGGGCGAGGAACGACCTT‑TTG‑GCGT‑ATAACATCGAGAGCCCCAGGGCACCGATAATCGCC‑ACCAACCATCCATAAAAGG‑CAATCCCCCCAGTAGGGGCCAGGTTTGCAGGGAACAGAAAAACACCTAATCCCATCATATTCCCCGACACCATCCGGGTGACGGGGGTTAAGCCCAcc                                                                                                                                      >  2:325045/1‑250 (MQ=255)                                                                                                         gagGA‑CGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACC‑AGGACTTGGGTCGAGGAACGACATT‑TTG‑GCGT‑ATACCATCGAGAGCCCCCGTGCACCGATAATTGTC‑AACAACCATCCATAAATGG‑CAATCCCGCCAGTAGAAGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTTCCCGACACCATCAGGGGGACGGGGGTTAAGCCCCCtt                                                                                                                                     >  2:389441/1‑251 (MQ=255)                                                                                                                   ggTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACC‑AGGACTTGGGTCGAGGAACGACATT‑TTG‑GCGT‑ATACCATCGAGAGCCCCAGTGCACCGATAATCGTC‑ACCAACCATCCATAAATGG‑CAATCCCGCCAGTAGAGGCCAGGATTGCAGGAAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGAGACGGGGATTAAGCCCACTTTGTGagcat                                                                                                                            >  1:173616/1‑251 (MQ=255)                                                                                                                     tttGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACC‑AGGACTTGGGTCGAGGAACGACATT‑TTG‑GCGT‑ATACCATCGAGAGCCCCAGTGCACCGATAATCGTC‑ACCAACCATCCATAAATGG‑CAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGagcatc                                                                                                                           >  1:238766/1‑250 (MQ=255)                                                                                                                      ttGACAACCAAGAAGCGGGCCAAAGCGGCGGGGGGCGTAAG‑GTAAAAACCCCCCAGGATGTGGTGCGAGGAACGACTTT‑TTG‑GCGT‑ATACCACCGAGAGCCCCAGTGCACCGAAAATCGAC‑ACCAACCATCCATAAATGG‑CAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGagcatca                                                                                                                          <  2:535003/250‑1 (MQ=255)                                                                                                                       aaacaacggagaAACGGGCAAAAGCGGCGGGGGGCGAAAGCGAAAGAACCCCC‑AGGACTTGGGTCGAGGAACGA‑ATTGTGG‑GCGT‑ATCCAACCGAGAGCCCCAGTGCACCGAAAACCGCC‑ACCACCCACCCAAAAATGG‑AAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGagcatcagc                                                                                                                        <  2:593653/243‑1 (MQ=255)                                                                                                                               ggaaaaCGGCCAAAAGCGGCGGGGGGCTTAAGCGAAAGAACCCCC‑AGGATTTGGGTCGAGAAACGACATT‑TGG‑GCTT‑ATACCATCGAGAGCCCCAGTGCACCGATAATCGTC‑ACCACCCATCAAAAAAGGG‑CAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAg                                                                                                                <  1:326370/248‑1 (MQ=255)                                                                                                                                   accggaacAAAGCAGCGGGGGGCGTAAGCTTAAGAACCACC‑GGCACTGGGGTCGAGGAACGATATT‑TTG‑GCGT‑ATACCCCCGAGAGCCCCAGTGCACCGATAAT‑GTCCACCAACCACCCATAA‑TGGACACTCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACAt                                                                                                            <  1:389441/244‑1 (MQ=255)                                                                                                                                   aaCGGGACACAGCAGCGGCGGGCGTAAGGGTAAGAACCACC‑AGGACTTGGGTCGAGGAACGACATT‑TAT‑GCGA‑ATACCAAAGAGAACCCCAGTGCACCGATAATCGTC‑AACAAACATCCATAAATGG‑CAATCCCGCCAGAAGGGGCCAGGGTTGCGGGTAACAGCAATACACCTGACCCCCTAATATTCCCCGACACCATCAGGGAGACGGGGGTTAAGCCCACTTTGTGggcgtcagcatacgtaaaagt                                                                                                            >  2:68434/1‑241 (MQ=255)                                                                                                                                       ggCCAATGCAGCGGCGGGCGTAAGCGTAAAACCCCCC‑AGGACTGGGGTCGAGAAACGACATT‑TTG‑GCGT‑ATACCCTCGAGCGCCCCAGGGACCCGATAACCTTC‑ACCAACCACCCAAAAAGGG‑CAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAAtt                                                                                                        <  2:238766/251‑1 (MQ=255)                                                                                                                                        ccaaaaGCAGCGGCGGGCGTAAGCGTAAGAACCCCC‑AGGTCTGGGGTCGAGGAACGACATT‑TTG‑GCGT‑ATACCATCGAGAGCCCCAGTGCCCCGATAATCGCC‑CCCAACCATCCAAAAATGG‑CAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTa                                                                                                       <  1:563282/248‑1 (MQ=255)                                                                                                                                           aacgcagcGGCGGGCGTAAGCGTAAGACCCCCC‑AGGTCTTGGGTCGAGGAACGACATT‑TTG‑GCGT‑ATACCATCGAAACCCCCGGTGCACCGATAATCGTC‑ACCAACCATCCAAAAATGG‑CAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAAc                                                                                                    <  1:401727/248‑1 (MQ=255)                                                                                                                                           aaaGCAGCGGCGGGCGTAAGCGTAAGAACCACC‑AGGACCTGGGTCGAGGAACGACATT‑TTG‑GCGT‑ATACCATCGAGAGCCCCAGTGCACCGATAATCGTC‑ACCAACCATCCATAAATGG‑CAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGGGACGGGGATTAAGCCCACTTTGTGAGAATCAGCAACCGAAGACATAATTaaa                                                                                                     >  2:432975/1‑250 (MQ=255)                                                                                                                                              gcagcGGCGGGCGTAAGCTTAAAACCCACC‑AGGACTTGGGTCGAGGAACGACATT‑TTG‑GCGT‑ATACCATCGAGAGCCCCAGTGCACCGATAATCGTC‑CCCACCCACCCAAAAAGGG‑CAATCCCGCCAGTAGAGGCCAGGTTTGc                                                                                                                                                                                                            <  1:373563/144‑1 (MQ=255)                                                                                                                                              gcagcGGCGGGCGTAAGCGTAAGAACCACC‑AGGACTTGGGTCGAGGAACGACATT‑TTG‑GCGT‑ATACCATCGAGAGCCCCAGTGCACCGATAATCGTC‑ACCAACCATCCATAAATGG‑CAATCCCGCCAGTAGAGGCCAGGTTTGc                                                                                                                                                                                                            >  2:373563/1‑144 (MQ=255)                                                                                                                                                         gcgtagccgaaagaccccccgggACTGGGGTCGAGAAACGCAATT‑TGG‑GCTT‑ATCCAACCGAGACCCCCATGGCACCAATAATCGCC‑ACCACCCACCAATAAATGA‑CAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCg                                                                                       <  1:87277/230‑1 (MQ=255)                                                                                                                                                                  aagAACCACC‑AGGACTTGGGTCGAGGAACGACATT‑TTG‑GCGT‑ATACCATCGAGAGCCCCAGTGCACCGATAATCGTC‑ACCAACCATCCATAAATGG‑CAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTg                                                                              >  2:587004/1‑250 (MQ=255)                                                                                                                                                                   agAACCACC‑AGGACTTGGGTCGAGGAACGACATT‑TTG‑GCGT‑ATACCATCGAGAGCCCCAGTGCACCGATAATCGTC‑ACCAACCATCCATAAATGG‑CAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATATTTAAACTCCTGCGAAGGCGAGCTTCGTGAc                                                                            >  2:199002/1‑251 (MQ=255)                                                                                                                                                                   aaaaccccccgggcctggggTCGAGAAACGACATT‑TTG‑GCTT‑ATCCCACCGAGACCCCCAGGCCCCCGAAAACCGCC‑ACCAACCACCCAAAAATGC‑CAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGa                                                                             <  1:139216/234‑1 (MQ=255)                                                                                                                                                                              ggACTGGGGTCGAGGAACGAC‑TT‑TTGTGCGT‑ACCCAATCGAGAGCCCCAGTGCACCGATAATCGTC‑CCCAACCCTCCAAAAAAGA‑CAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTAc                                                                  <  1:210814/251‑1 (MQ=255)                                                                                                                                                                                      gTCGAGGAACGACATT‑TTG‑GCGT‑ATACCATCGAGAGCCCCAGTGCACCGATAATCGTC‑ACCAACCATCCATAAATGG‑CAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGtcatatc                                                          >  1:384176/1‑251 (MQ=255)                                                                                                                                                                                        ggAGGAACGACTTT‑TTG‑GCGT‑ATCCCATCGAGACCCCCAGTGCCCCAATAATCGCC‑CCAAACCATCAATAAATGG‑CAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGATCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGtcatatcat                                                        <  1:596227/250‑1 (MQ=255)                                                                                                                                                                                          aGGAACGACATT‑TTG‑GCGT‑ATACCATCGAACGCCCCGGTGCACCGATAATCGTC‑ACCACCCACCCATAAATGG‑CAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATAc                                                      <  1:311697/251‑1 (MQ=255)                                                                                                                                                                                           ggAACGACATT‑TTG‑GCGT‑ATACCATCGAGAGCCCCAGTGCACCGATAATCGTC‑ACCAACCATCCATAAATGG‑CAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACg                                                     >  1:157526/1‑251 (MQ=255)                                                                                                                                                                                            gAACGACATT‑TTG‑GCGT‑ATACCATCGAGAGCCCCAGTGCACCGATAATCGTC‑ACCAACCATCCATAAATGG‑CAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGc                                                    >  2:78969/1‑251 (MQ=255)                                                                                                                                                                                             ggcaacgatttttgggcttaccccATCGAGAGCCCCAGTGCACCGTAAATCGCC‑ACCACCCACCCAAAAAGGG‑CAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATAc                                                      <  1:175754/229‑1 (MQ=255)                                                                                                                                                                                                  cATT‑GTG‑GCGT‑AACCCATCGAGAGCCCCAGTGCCCCAAAAATCATC‑CCCAACCACCCAAAAATGC‑CAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCAt                                              <  1:448465/251‑1 (MQ=255)                                                                                                                                                                                                              t‑ataCCATCGAGAGCCCCAGTGCACCGATAATCGTC‑ACCAACCATCCATAAATGG‑CAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACCCCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATAAGACTGCATTCGATAGTTg                                    >  2:569816/1‑251 (MQ=255)                                                                                                                                                                                                               atacccATGGAGAGCCCCAGTCCCCCGATACGCGTC‑CCCACCCACCCAAAAATGG‑CACCCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTg                                    <  2:286761/247‑1 (MQ=255)                                                                                                                                                                                                                  cccacggagaGCCCCGGTGCCCCCATAATCCAC‑ACCAACCACCCAAAAAGGG‑CAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGaa                                  <  1:199002/244‑1 (MQ=255)                                                                                                                                                                                                                     aTCGAGAGCCCCAGTGCACCGATAATCGAC‑ACCAACCATCCATAAATGG‑CAATCCCGCCAGTAGAGGCCAGGTTTGCCGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGAGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAattat                              >  1:55616/1‑251 (MQ=255)                                                                                                                                                                                                                     aTCGAGACCCCCAGGGCACCGATAATCGTC‑ACCAACCACCCAAAAATGC‑CAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAattat                              <  1:40492/251‑1 (MQ=255)                                                                                                                                                                                                                              cccAGTGCACCGATAATCGTC‑ACCAACCATCCATAAATGG‑CAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCg                     >  1:26026/1‑251 (MQ=255)                                                                                                                                                                                                                                             aTCGTC‑ACCAACCATCCATAAATGG‑CAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGc      >  2:13269/1‑251 (MQ=255)                                                                                                                                                                                                                                                 cc‑accacCCACCCATAAAGGG‑CAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGtt  <  1:277091/250‑1 (MQ=255)                                                                                                                                                                                                                                                      caaccaTCCATAAATGG‑CAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAattat                              >  2:381235/1‑219 (MQ=255)                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                              TGGTTAATACCAATGGATCTTTCAGAATCGGGAAGAAGTAACTTAAATATCCTACGCCAATGACCACCATGGCGATATTGCCGATCC‑AGCAGGCCAGCCAGTAGAGGA‑CGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACC‑AGGACTTGGGTCGAGGAACGACATT‑TTG‑GCGT‑ATACCATCGAGAGCCCCAGTGCACCGATAATCGTC‑ACCAACCATCCATAAATGG‑CAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGTT  >  NC_000913/4336653‑4337137

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