Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,336,901	0	T	A	28.9%	62.0 / 13.2	46	Y44F (TAT→TTT)	adiC	arginine:agmatine antiporter
Reads supporting (aligned to +/- strand):  ref base T (22/10);  new base A (0/13);  total (22/24)
Fisher's exact test for biased strand distribution p-value = 2.25e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ATGGTTAATACCAATGGATCTTTCAGAATCGGGAAGAAGTAACTTAAATATCCTACGCCAATGACCACCATGGCGATATTGCCGATCCAGCAGGCCAGCCAGTAGAGGACGTTGGTTTGATAACCGAGAAA‑CGGGCCAAA‑GCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGA‑CTTGGGTCGAGGAACGACA‑TTTTGGCGTATAC‑CATCGAGAGCCCCAGTG‑CACCG‑ATAATCGTCACCAACCA‑TCCATAAAT‑GGCAA‑TCCC‑GCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGA  >  NC_000913/4336652‑4337130                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                         aTGGTTAATAACAATGGATCTTTCAGCATCGGGACGAAGTAACTTACATCTCCTAAGCCAATAACCCCCATGGCGATATTGCCGATCCAGCAGGCCAGACAGTAGAGGACGTTGGTTTGATAACCGAGAAA‑CGGGCCAAA‑GCACCGGCGGGCGTAAGCGTAAGAACCACCCGGC‑CATCGGTCGAGGCACGACA‑TCTTGGCGTATAC‑CCTCGAGAGACCCAGTG‑GACCG‑ATAATCGACACCAAcca‑tccat                                                                                                                                                                                                                                          >  2:194552/1‑250 (MQ=255)             aaTGGATCTGTCAGAATCGGGAAGAAGTAACTTAAATTTCCTACGCCAATGACCACCATGGCGATATTGCCGATCCAGCAGGCCAGCCAGTAGAGGACGTTGGTTTGATAACCGAGAAA‑CGGGCCAAA‑GCAGCGGCGTGCGTAAGCGTAAGAACCACCAGGA‑CTTGGGTCGAGGAACGACA‑TTTTGGCGTATAC‑CATCGAGAGCCCCAGTG‑CACCG‑ATAATCGTCACCAACCA‑TCCATAAAT‑GGCAA‑Tccc                                                                                                                                                                                                                           <  2:156596/251‑1 (MQ=255)                      ttCAGACTCGGGAAGAAGCAACTTAAATATCCTACGCCAATGCCCACCCTGGCGATATTGCCGATCCAGCAGGCCAGCCAGTAGAGGACGTTGGTTTGATAAACGAGAAA‑CGGGCCAAA‑GCAGCGGCGGGCGTAAGCGTAAGAAACAACAGGC‑CTTGGGTCGGGGAACGACA‑TTTTGGCGTATAC‑CATCGAGAGCCCCAGGG‑CGCCG‑ACAATCGGCACCAACCA‑ACCATAAAT‑GGCAA‑TCCC‑GCCAGTagg                                                                                                                                                                                                                 >  2:184188/1‑250 (MQ=255)                         aaaaTCGGAAAGAAGTATCTTAAATATCCTACCCCAATGCCCACCGTGGCGATTTTGCCGATCCAGGAGGCCAGCCGGTAGAGGACGTTGGTTTAAAAACCGAGAAA‑CGGGCCAAA‑GCAGCGGCGGGGGTAAGCGTAAAACCCCCCGGGA‑CTTGGGTGGAGGAACGACA‑TTTTGGCGTATCC‑CATCGAGAGCCCCAGTG‑CACCT‑ATAATCGTCACCACCCA‑CCCATATTT‑GGCAA‑GCCC‑GCCATTAGAGGc                                                                                                                                                                                                              <  1:359107/249‑1 (MQ=255)                         aGAATCGGGAAGAAGTAACTTAAATATCCTACGCCAATGACCACCATGGCGATATTGCCGATCCAGCAGGCCAGCCAGTAGAGGACGTTGGTTTGATAACCGAGAAA‑CGGGCCAAA‑GCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGG‑CTTGGGTCGAGGAACGACA‑TTTTGGCGTATAC‑CATCGAGAGCCCCAGTG‑CACCG‑ATAATCGTCCCCAACCA‑TCCATAAAT‑GGCAA‑TCCC‑GCCCGTAGAGGc                                                                                                                                                                                                              >  2:145692/1‑251 (MQ=255)                                     aagTAACTTAAATATCCTACGCCAATGACCACCATGGCGATATTGCCGATCCAGCAGGCCAGCCAGTAGAGGACGTTGGTTTGATAACCGAGAAA‑CGGGCCAAA‑GCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGA‑CTTGGGTCGAGGAACGACA‑TTTTGGCGTATAC‑CATCGAGAGCCCCAGTG‑CACCG‑ATAATCGTCACCCACCA‑TCCATAAAT‑GGCAA‑TCCC‑GCCAGTAGAGGCCAGGTTTGCAgg                                                                                                                                                                                                  >  1:362650/1‑251 (MQ=255)                                       gTAACTTAAATATCCTACGCCAATGACCACCATGGCGATATTGCCGATCCAGCAGGCCAGCCAGTAGAGGACGTTGGTTTGATAACCGAGAAA‑CGGGCCAAA‑GCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGA‑CTTGGGTCGAGGAACGACA‑TTTTGGCGTATAC‑CATCGAGAGCCCCAGTG‑CACCG‑ATAATCGTCACCAACCA‑TCCATAAAT‑GGCAA‑TCCC‑GCCAGTAGAGGCCAGGTTTGCAGGt                                                                                                                                                                                                 >  1:16421/1‑250 (MQ=255)                                                  atCCTACGCCAATGACCACCATGGCGATATTGCCGATCCAGCAGGCCAGCCAGTAGAGGACGTTGGTTTGATAACCGCGAAA‑CGGGCCAAA‑GCAGCGGCGGGCGTAAGCGTAAAAACCACCAGGA‑CTTGGGTCGAGGAAAGACA‑ATTTGGCGTCTAC‑CATCGCGAGCCCCAGTG‑AACCG‑CTAATCGTCACCAACCA‑TCCCTaaa                                                                                                                                                                                                                                       >  2:155919/1‑204 (MQ=255)                                                  atCCTACGCCAATGACCACCATGGCGATATTGCCGATCCAGCAGGCCAGCCAGTAGAGGACGTTGGTTTGATAACCGATAAA‑CTGGCCAAA‑GCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGA‑CTTGGGTCGAGGAACGACA‑TTTTGGCGTATAC‑CATCGAGAGCCCCAGTG‑CACCG‑ATAATCGTCACCAACCA‑TCCATaaa                                                                                                                                                                                                                                       <  1:155919/204‑1 (MQ=255)                                                               gACCACCATGGCGCTATTGCCGATCCAGCAGGCCAGCCAGTAGAGGACGTTGGTTTGATAACCGAGAAA‑CGGGCCAAA‑GCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGA‑CTTGGGTCGAGGAACGACA‑TTTTGGCGTATAC‑CATCGAGAGCCCCCGTG‑CACCG‑ATAATCGTCACCAACCA‑TCCATAAAT‑GGCAA‑TCCC‑GCCAGTAGAAGACCGGTTTTCAGGTAACCGGAAAACACCTGACCCCata                                                                                                                                                                         >  2:473421/1‑250 (MQ=255)                                                               gACCACCATGGCGATATTGCCGATCCAGCAGGCCAGCCAGTAGAGGCCGTTGGTTTGATAACCGAGACA‑CGGGCCAAA‑GCCGCGGCGGGCGGAAGCGTAAGAACCACCCTGA‑CTTGGGTCCAGGAACGACA‑TTTTGGCGTATAC‑CATCGAGAGCCCCCGTG‑CACCG‑ATCATCGTTACCAACCA‑TCCATAAAT‑GGCAA‑TCCC‑GCCAGGAGTGGCCAGGATTGCAAGGAAAAGAAAAAAACCTGAacccctt                                                                                                                                                                         >  2:91259/1‑247 (MQ=255)                                                                    ccaTGGCGATATTGCCGATCCAGCAGGCCAGCCAGTAGAGGACGTTGGTTTGATAACCGAGAAA‑CGGGCCAAA‑GCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGA‑CTTGGGTCGAGGAACGAAA‑TTTTGGCGTATAC‑CATCGAGAGCCCCCGTG‑CACCG‑ATAATCGTCACCAACCA‑TCCCTAAAT‑GGCAA‑TCCC‑GCCAGTAGAGGCAAGGGTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTc                                                                                                                                                                   >  1:266972/1‑251 (MQ=255)                                                                              atgtgctacccagtgggCCAGCCAGAAGAGTACGTAGTTTGGAGAACC‑AGAAAAGGGGCCAAA‑GCAGCGGCGGGCTAAAGGTTACGAACCTCCAGGT‑CTGGGGTCGAGGAACGACA‑TTTTGGCTTATAC‑CATCGAGAGCCCCGGTG‑C‑CCGTATACTCTCCTCCACCCA‑‑CCAAAAAA‑TGTAA‑TCCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACAcc                                                                                                                                                          <  2:368531/235‑1 (MQ=255)                                                                                        cagcagGCCAGCCAGTAGAGGACGTTGGTTTGATAACCGAGAAA‑CGGGCCAAA‑GCAGCGGCGGGCGTAAACGTAAGAACCACCCGGA‑CTTGAGTCGAGGAACGACC‑TTTTAGCGTATAC‑CATCGCGAGCCCCAGTG‑CACCG‑CTAATCGTCACCCACCA‑TCCATAAAT‑GGCAA‑TCCC‑GACAGAAGAGGCCAGGTTTCCAGGGAACAGAAAAACAACTGCCCCCATACTCTTCCCCGACACCATCAGGGTGAc                                                                                                                                               >  1:189488/1‑251 (MQ=255)                                                                                               ccaccaagTAGAGGACGTTGGTTTGATACCCGAGAAC‑CGGCCAAAA‑GCGGGGGCGGGCGAAAGCGAAAGACCCACCAGGA‑CTGGGGTCAAGGAACGACA‑TTTTGGCGTATAC‑CATCGAGACCCCCAGTG‑CCCCG‑AAAAGCGTCACCAACCA‑CCCATAAAT‑GGCAA‑TCCC‑GCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGAtt                                                                                                                                        <  1:42006/245‑1 (MQ=255)                                                                                                       tAGAGGACGTTGGTTTGATAACCGAGAAA‑CGGGCCAAA‑GCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGA‑CTTGGGTCGAGGAACGACA‑TTTTGGCGTATAC‑CATCGAGAGCCCCAGTG‑CACCG‑ATAATTGTCACCAACCA‑TCCATAAAT‑GGCAA‑TCCC‑GCCAGTAGAGGCCCGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCCTCAGGGGGACGGGGAt                                                                                                                                         >  2:132651/1‑242 (MQ=255)                                                                                                       tAGAGGACGTTGGTTTGATAACCGAGAAA‑CGGGCCAAA‑GCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGA‑CTTGGGTCGAGGAACGACA‑TTTTGGCGTATAC‑CATCGAGAGCCCCAGTG‑CACCG‑ATAATCGTCACCAACCA‑TCCATAAAT‑GGCAA‑TCCC‑GCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCAc                                                                                                                                >  1:171746/1‑251 (MQ=255)                                                                                                       tAGAGGACGTTGGTTTGAAAACCAAGAAC‑CGGGCCAAA‑GCAGGGGGGGGCGTAAGCGTAAGACCCCCCAGGA‑CTTGGGTCGAGAAACGACA‑TTTTGGCGTATCC‑CATCGAGACCCCCAGTG‑CACCG‑ATAATCGTCCCCAACCA‑CCCAAAAAT‑GGCAA‑TCCC‑GCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGAt                                                                                                                                         <  1:132651/242‑1 (MQ=255)                                                                                                         gagGACGTTGGTTTGATAACCGAGAAA‑CGGGCCAAA‑GCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGA‑CTTGGGTCGAGGAACGACA‑TTTTGGCGTATAC‑CATCGAGAGCCCCAGTG‑CACCG‑ATAATCGTCACCAACCA‑TCCATAAAT‑GGCAA‑TCCC‑GCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTcccc                                                                                                                                                                >  2:78385/1‑217 (MQ=255)                                                                                                         gagGACGTTGGTTTGATAACCGAGAAA‑CGGGCCAAA‑GCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGA‑CTTGGGTCGAGGAACGACA‑TTTTGGCGTATAC‑CATCGAGAGCCCCAGTG‑CACCG‑ATAATCGTCACCAACCA‑TCCATAAAT‑GGCAA‑TCCC‑GCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGGAGGGGATTAAGCCCACtt                                                                                                                              >  2:194248/1‑251 (MQ=255)                                                                                                         gagGACGTTGGTTTGAAACCCGAGAAA‑CGGGCAAAA‑GCAGCGGCGGGCGTAAGCGTAAGACCCACCAGGA‑CTTGGGTCGAGGAACGACA‑TTTTGGCGTATCC‑CATCGAGAGCCCCAGTG‑CACCG‑ATAATCGTCCCCACCCA‑CCCATAAAT‑GGCAA‑TCCC‑GCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTcccc                                                                                                                                                                <  1:78385/217‑1 (MQ=255)                                                                                                           gggcGTTGGGTAGATACCCGAGAAA‑CGGGCAAAG‑GCAGGGGTGGGCGTACGCGAAAGAACCCCCAGGC‑CTGGGGTCGAGGAACGACA‑TGTGGGCGTATCC‑CATCGAGAGCCCCTGCG‑CACCG‑ATAATCTTCACCACCCA‑CCCAAAAAT‑GACAA‑TCCC‑GCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGTATTAAGCCCACTTtg                                                                                                                            <  2:446268/248‑1 (MQ=255)                                                                                                                ttggtttaacaccccagaAC‑CGGGCCACA‑GCAGCGGCGGTGGTATGCTAAAGACCCACAAGGT‑CTTGGGGCGAGGAACAAAA‑TTTTGGCTTAACC‑CAGCGAGCGCCCGAGTG‑CCCCA‑ATAACCGTCCCCACCCA‑TCCATAAAT‑GGCAA‑TCCC‑GCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACGCCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGag                                                                                                                        <  2:324437/235‑1 (MQ=255)                                                                                                                        aTAACCGAGAAA‑CGGGCCAAA‑GCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGA‑CTTGGGTCGAGGAACGACA‑TTTTGGCGTATAC‑CATCGAGAGCCCCAGTG‑CACCG‑ATAATCGTCACCAACCA‑TCCATAAAT‑GGCAA‑TCCC‑GCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGGGACGGGGATTAAGCCCACTTTGTGagcataagcct                                                                                                               >  2:302653/1‑249 (MQ=255)                                                                                                                          aaCCGAGAAA‑CGGGCCAAA‑GCAGCGGCGGGCGTAAGCGTAAGAACCACCAAGA‑CTTGGGTCGAGGAACGACA‑TTTTGGCGTATAC‑CATCGAGAGCCCCAGTG‑CACCG‑ATAATCGTCACCAACCA‑TCCATAAAT‑GGCAA‑TCCC‑GCCAGTAGAGGCCCGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATcc                                                                                                             >  1:298968/1‑251 (MQ=255)                                                                                                                           cccGAGAA‑‑CGGGCCAAAACCAGCGGGGGGCGTAGGCG‑AAAAACCCCCCGGAACTTGGGTCGAGAACCGCCT‑TTTTGGCGTATCC‑CACCGAGGGCCCCATGG‑CCCCA‑ATAACCGCCCCCACCCA‑‑CCATAAATTGGCAG‑CCCC‑GCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCg                                                                                                            <  1:184188/250‑1 (MQ=255)                                                                                                                                                cagcGGGGGGCGTAAGCGAAAAAACCACCAGGA‑TTTGGGTCGAGGAACGACT‑TTTGGGCGTATAC‑CACCGAGAGCCCCAGTG‑CACCG‑ATAATCGCCACCAACCA‑TCCAAAAAT‑GGCAA‑CCCC‑GCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAGCTCCt                                                                                         <  1:145692/251‑1 (MQ=255)                                                                                                                                                      ggggctaaagcttaagAACCCCAAGGA‑CTTGGGTCGAGGAACGACA‑TTTGGGCGTATCC‑CATCGAGAGCCCAAGGG‑CCCCG‑ATAATCTCCACCACCCA‑TCCATAAAT‑GGCAA‑TCCC‑GCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAg                                                                                   <  1:196901/244‑1 (MQ=255)                                                                                                                                                        ggcgTAAGCGTAAGAACCACCAGGA‑CTTGGGTCGAGGAACGACA‑TTTTGGCGTATAC‑CATCGAGAGCCCCAGTG‑CACCG‑ATAATCGTCACCAACCA‑TCCATAAAT‑GGCAA‑TCCC‑GCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAggg                                                                                 >  2:111989/1‑250 (MQ=255)                                                                                                                                                          ggtaagcgaaagACCCACCAGAA‑CTTGGGTCGAGGAACGACA‑TTTTGGCGTATCC‑CATCGAGACCCCCAGTG‑CACCG‑ATAATCGTCCCCACCCA‑CCCAAAAAG‑GGCAA‑TCCC‑GCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGa                                                                               <  1:302653/250‑1 (MQ=255)                                                                                                                                                              ggcgaaagAACCCCCAGGA‑TGGGGGTGGAGGAAC‑ACAATTTTGGCTTATACGCA‑CGAGAGCCCCGGTG‑CACCG‑ATACCCGCCACCACCCA‑TCAAAAAAT‑GGAAA‑TCCC‑GCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCtt                                                                           <  2:298968/246‑1 (MQ=255)                                                                                                                                                              agcgtaagAACCCCCAGGA‑CTTGGGTCGAGAAGCGATA‑TTTTGTCGTATCC‑CACCGTGCGCCCCAGTG‑CCCCT‑ATAATCGCCCCCAACCT‑TCAAAAAAG‑GGCAA‑TCCC‑GCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTTACCCCATAATATTCCCCGACACCATCAGGGTGACGGTGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCtt                                                                           <  1:306288/251‑1 (MQ=255)                                                                                                                                                                 gaaagACCCACCAGGA‑CTTGGGTCGAGGATCTATA‑TTTTGGTGTAGAC‑CATCGAGACCCCCGGTG‑CCCCG‑ATAATCGCCGCCACCCA‑CCCATAAAT‑GGCAA‑ACCC‑GCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGt                                                                        <  1:205238/249‑1 (MQ=255)                                                                                                                                                                    agAACCACCAGGA‑CTTGGGTCGAGGAACGACA‑TTTTAGCGTATAC‑CATCGGGAGCCCCAGTG‑CACCG‑ATAATCGTCACCAACCA‑TCCATAAAT‑GGCAA‑TCCC‑GCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAAATCCTGCGAAGGCGAGCTTCGTGAc                                                                     >  2:148955/1‑251 (MQ=255)                                                                                                                                                                                          gAGGAACGACA‑TTTTGGCGTATAC‑CATCGAGAGCCCCAGTG‑CACCG‑ATAATCGTCACCAACCA‑TCCATAAAT‑GGCAA‑TCCC‑GCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGtcatatcat                                                 >  1:435021/1‑250 (MQ=255)                                                                                                                                                                                            ggAACGACA‑TTGTGGCGTATAC‑CAACGAGTGCCC‑AGTGGCACCG‑ATAAGCGCCCCAAACCA‑TCCAAAAAT‑GCAAA‑CCCC‑GCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACg                                              <  2:10974/251‑1 (MQ=255)                                                                                                                                                                                                       tttGGCGTATAC‑CATCGAGAGCCCCAGTG‑CACCG‑ATAATCGTCACCAACCA‑TCCATAAAT‑GGCAA‑TCCC‑GCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCg                                    >  1:488727/1‑251 (MQ=255)                                                                                                                                                                                                                 aC‑CATCGAGAGCCCCAGTG‑CACCG‑ATAATCGTCACCAACCA‑TCAATAAAA‑GGCAA‑TCCC‑GCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAACACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGAAGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCATGCATTCGATAGTTGAat                          >  1:479825/1‑251 (MQ=255)                                                                                                                                                                                                                         agagCCCCAGGG‑CCCCG‑AAAAGCGCCACCAACCA‑CCAAAAAAG‑GGAAA‑TCCC‑GCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATAGCATACGCCTGCATTCGATAGTTGAATTATTATc                   <  2:171746/251‑1 (MQ=255)                                                                                                                                                                                                                          gagCCCCGGGG‑CACCT‑ATAATCGCCACCAACCC‑CCAAAAAAT‑GGCAA‑TCCC‑GCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCg                  <  1:221933/251‑1 (MQ=255)                                                                                                                                                                                                                          agacccccaggcccccG‑ATAATCATCACAAACCT‑TCCA‑AAAA‑GGCAACTCCC‑GCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTTACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATc                   <  2:16421/237‑1 (MQ=255)                                                                                                                                                                                                                              cccAGTG‑CACCG‑ATAATCGACACCAACAA‑TCCATAAAT‑GGCAA‑TCCC‑GCCAGTAGAGGCAAGGCTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTcccc                                                                                                                                                                >  2:422031/1‑105 (MQ=255)                                                                                                                                                                                                                                cccccccttcccgatcgcctccccaacca‑CCAATAAAG‑GGCAC‑CCCC‑GCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCg              <  2:120884/229‑1 (MQ=255)                                                                                                                                                                                                                                  gTG‑CCCCG‑ATAATCGTCCCCACCCA‑CCAAAAAAT‑GACAA‑TCCC‑GCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATGAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATGACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGtt          <  1:22262/251‑1 (MQ=255)                                                                                                                                                                                                                                   gtccACCG‑ATAATCGCCACCAACCA‑CCCAAAAAG‑GGCAA‑TCCC‑GCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGtt          <  1:184280/248‑1 (MQ=255)                                                                                                                                                                                                                                            aaaaTCGTCACCA‑CCACACCAAAAAT‑GGCAA‑CCCC‑GCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGa  <  1:473421/249‑1 (MQ=255)                                                                                                                                                                                                                                                 cGTCACCAACCA‑TCCATAAAT‑GGCAA‑TCCC‑GCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCg              >  1:120884/1‑234 (MQ=255)                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                         ATGGTTAATACCAATGGATCTTTCAGAATCGGGAAGAAGTAACTTAAATATCCTACGCCAATGACCACCATGGCGATATTGCCGATCCAGCAGGCCAGCCAGTAGAGGACGTTGGTTTGATAACCGAGAAA‑CGGGCCAAA‑GCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGA‑CTTGGGTCGAGGAACGACA‑TTTTGGCGTATAC‑CATCGAGAGCCCCAGTG‑CACCG‑ATAATCGTCACCAACCA‑TCCATAAAT‑GGCAA‑TCCC‑GCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGA  >  NC_000913/4336652‑4337130

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