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

GTGATGGTTAATACCAATGGATCTTTCAGAATCGGGAAGAAGTAACTTAAATATCCTACGCCAATGACCACCATGGCGATATTGCCGATCC‑AGCAGGCCAGCCAGTAGAGGA‑CGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACC‑AGGACTTGGGTCGAGGAACGACATT‑TTG‑GCGT‑ATACCATCGAGAGCCCCAGTGCACCGATAATCGTC‑ACCAACCATCCATAAATGG‑CAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGTT  >  NC_000913/4336649‑4337137                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                  gTGATGGTTAATACCAATGGATCTTTCAGAATCGGGAAGAAGTAACTTAAATATCCTACGCCAATGACCACCATGGCGATATTGCCGATCC‑AGCAGGCCAGCCAGTAGAGGA‑CGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACC‑AGGACTTGGGTCGAGGAACGACATT‑TTG‑GCGT‑ATACCATCGAGAGCCCCAGTGCACCGATAATCGTC‑ACCAAccatcc                                                                                                                                                                                                                                                 <  2:311842/251‑1 (MQ=255)     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)                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                  GTGATGGTTAATACCAATGGATCTTTCAGAATCGGGAAGAAGTAACTTAAATATCCTACGCCAATGACCACCATGGCGATATTGCCGATCC‑AGCAGGCCAGCCAGTAGAGGA‑CGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACC‑AGGACTTGGGTCGAGGAACGACATT‑TTG‑GCGT‑ATACCATCGAGAGCCCCAGTGCACCGATAATCGTC‑ACCAACCATCCATAAATGG‑CAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGTT  >  NC_000913/4336649‑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: ‑