Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,455,940	0	A	C	31.5%	69.3 / 10.9	57	T6P (ACC→CCC)	paaD	phenylacetate degradation protein
Reads supporting (aligned to +/- strand):  ref base A (24/13);  new base C (0/17);  total (24/33)
Fisher's exact test for biased strand distribution p-value = 3.02e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GCACTGAGTGAAGAGGGTATTGCGGTTGATCCACGCACTTTACGCGCAGCGTG‑GGAAGCCGAAGTTTTTGCCGGGATCAACGAAGCCACATTGAACGTACCGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAA‑GGACTGCATACCGAACACC‑TTGGACCCATGCTGGCAGAAATGCAGTATCTCCAGCGTGT‑CTTGCCCGGTCAG‑CAATGGTAA‑CAGA‑GGAGATGGGTATGCAACGTCTGG‑CCACCATTGCA‑CC‑GCCACAGGTTCATGAGATAT‑GGGCGCTACTCAGCCAGATCCCGGACCCGGAGATCCCGGTGCTGACCATTACTGATTTAGGCATGGTGCGTAATGTGACACAGATGGGAGAAGGATGGGTGATCGGCTTTACGCCGACATATTCCGGTTGTCCGGCAACGGAACATTTGATTGGCGCGATACGTGAGGCAATGACAACCAACGGGTTT‑ACCCCCGTTCAGGTTGTGCTGCAACTCGACCC  >  NC_000913/1455695‑1456190                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                               gtaCTGAGTGAAGAGGGTATCGCGGTTGATCCACGCACTTTACGCGTCTC‑TGTGGAAGCCGAAGTTTTTGCCGGGATCAACGAAGCCACATTGAACGTACCGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAA‑GGACTGCATACCGAACACC‑TTGGACCCATGCTGGCAGAAATGCAGTATCTCCCGCGTGT‑CTTTTCCGGTCAG‑CAATGGTAA‑CAGA‑GGAGATGGGTATGCAACGTCTGG‑CCACCAtt                                                                                                                                                                                                                                                           <  1:583826/249‑1 (MQ=255)       aGTGAAGAGGGTATTGCGGTTGATCCACGCACTTTACGCGCAGCGTG‑GGAAGCCGAAGTTTTTGCCGGGATCAACGAAGCCACATTGAACGTACCGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAA‑GGACTGCATTCCGAACACC‑TTGGACCCATGCTGGCAGAAATGCAGTATCTCCAGCGTGT‑CTTGCCCGGTCAG‑CAATGGTAA‑CAGA‑GGAGATGGGTATGCAACGTCTGG‑CCACCATTGCA‑CC‑g                                                                                                                                                                                                                                                   <  2:554201/251‑1 (MQ=255)          gAAGAGGGTATTGCGGTTGATCCACGCACTTTACGCGCAGCGTG‑GGAAGCCGAAGTTTTTGCCGGGATCAACGAAGCCACATTGAACGTACCGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAA‑GGACTGCATACCGAACACC‑TTGGACCCATGCTGGCAGAAATGCAGTATCTCCAGCGTGT‑CTTGCCCGGTCAG‑CAATGGTAA‑CAGA‑GGAGATGGGTATGCAACGTCTGG‑CCACCATTGCA‑CC‑Gccc                                                                                                                                                                                                                                                >  1:244866/1‑250 (MQ=255)                     tGCGGTTGATCCACGCACTTTACGCGCAGCGTG‑GGAAGCCGAAGTTTTTGCCGGGATCAACGAAGCCACATTGAACGTACCGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAA‑GGACTGCATACCGAACACC‑TTGGACCCATGCTGGCAGAAATGCAGTATCTCCAGCGTGT‑CTTGCCCGGTCAG‑CAATGGGAA‑CAGA‑GGAGATGGGGATGCAACGTCTGG‑CCACCATTGCA‑CC‑GCCACAGGGTCATg                                                                                                                                                                                                                                      >  1:146225/1‑250 (MQ=255)                       gggTTGATTCACTCACTTTACGCGCAGCGTG‑GGAAGCCGAAGTTTTTGCCGGGATCAACGAAGCCACATTGAACGTACCGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAA‑GGACTGCATACCGAACACC‑TTGGACCCATGCTGGCAGAAATGCAGTATCTCCAGCGTGT‑CTTGCCCGGTCAG‑CAATGGTAA‑CAGG‑GGAGATGGGTATGAAACGTCGGG‑CCACCATTGCA‑CC‑GCCACAGGTTCATgaga                                                                                                                                                                                                                                   <  1:672350/250‑1 (MQ=255)                        ggTTGATCCACGCACTTTACGCGCAGCGTG‑GGAAGCCGAAGTTTTTGCCGGGATCAACGAAGCCACATTGAACGTACCGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAA‑GGACTGCATACCGAAAACC‑TTGGACCCATGCTTGCAGAAATGCAATATCTCCAGCGTGT‑CTTGCCCGGTCAG‑CAATGGTAA‑CAGC‑GGTGATGGGAATGAAACGTCTGGACCAC‑ATTGCA‑CC‑GCCACAAGGTCATgaagt                                                                                                                                                                                                                                  >  1:731325/1‑248 (MQ=255)                              tCCACGCACTTTACGCGCAGCGTG‑GGAAGCCGAAGTTTTTGCCGGGATCAACGAAGCCACATTGAACGTACCGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAA‑GGACTGCATACCGAACACC‑TTGGACCCATGCTGGCAGAAATGCAGTATCTCCAGCGTGT‑ATTTCCCGGTCCG‑CCATGGTAA‑CAGA‑GGGGATGGGTATGCAAAGTCTGG‑CCACCATTTCG‑CC‑GCCACAGGTTCCTGAGATGT‑gggg                                                                                                                                                                                                                           >  2:231952/1‑250 (MQ=255)                                                   gTG‑GGAAGCCGAAGTTTTTGCCGGGATCAACGAAGCCACATTGAACGTACCGCCAGAGCAGGCGTATCGCACTGGCGGTAAAAAA‑GAACTGCATAACGAACACC‑TTGGACCCATGCTGGCAGAAATGCAGGATCTCCAGCGGGT‑CTTGCCCTCTCCG‑CAATGGTAA‑CAGG‑GGAGATGGGTATGCAACGTCTGG‑CCACCCATGCA‑CC‑CCCCCAGGTTACTGAGATAT‑GGGCGGTAATCCGCCGGATcccggc                                                                                                                                                                                                      >  2:537259/1‑250 (MQ=255)                                                    gg‑ggAAGCCGAAGTTTTGGCCGGGAACACCGAAGCCACATTGACCGTACCGAAAGAGCAGGCGTACCGCACTGGCGGAAAAAAA‑GGACTGAATACCGAACACC‑GTGGACCCATGCTGGAAGAAATGCAGTATCTCCAGCGTGT‑CTTGCCCGGTCAG‑CAAGGGTAA‑CAGA‑GGAGAGGGGTATTCAACGTCTGG‑CCCCCATTGCA‑CC‑GCCACAGGTTCATGAGATAT‑GGGCGCTACTCAGCCAGATcccggac                                                                                                                                                                                                     <  2:106170/250‑1 (MQ=255)                                                           gcagaattgtgtgccggaaaaaaCGAAGCCACATGTAACGTACCGAAAGAGAAGGCGAACCGCCCTGGCGGTAAAAAA‑GTACTGCATCCCGAACCCC‑TGGGACCCATGTTGGAAGAAATCAAGTATCCCAAGCGTGT‑TTCGCCCGGTCAG‑CAATGGTAA‑CAGA‑GAAGATGGGTATGCAACGTCTGG‑CCCCAATTGCA‑CC‑GCCACAGGTTCATGAGATAT‑GGGCGCTACTCAGCCAGATCCCGGACCCGgag                                                                                                                                                                                               <  2:577811/230‑1 (MQ=255)                                                              gaattttttGCCGGGATCAACGAAGCCACATTGAACGTACCGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAA‑GGACGGCATACCGAACACC‑TTGGACCCATGCTGGCAGAATTGCAGTATCTCCAGCGTGT‑CTTGCCCGGTCAG‑CAATGGTAA‑CAGA‑GGAGATGGGTATGCAATTTCTGG‑CCACCATTGCA‑CC‑GCCACAGGTTCATGAGATAT‑GGGCGCTACTCAGCCAGATCCCGGACCCGGAGATc                                                                                                                                                                                            <  1:286417/247‑1 (MQ=255)                                                                       gCCGGGATCAACGAAGCCACATTGAACGTACCGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAA‑GGACTGCATACCGAACACC‑TTGGACCCATGCTGGCAGAAATGCAGTATCTCCAGCGTGT‑CTTGCCCGGTCAG‑CAATGGTAA‑CAGA‑GGAGATGGGTATGCAACGTCTGG‑CCACCATTGCA‑CC‑GCCACAGGTTCATGAGATAT‑GGGCGCTACTCAGCCAGATCCCGGACCCGGAGATCCCGGTGCTg                                                                                                                                                                                   >  1:638336/1‑251 (MQ=255)                                                                         ggggatcacaaaaGCCACATTGAACGTACCGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAA‑GGACTGCATACCGAACCCC‑TTGGACCCATGCTGGAAGAAATGCAGTATCTCCGGCGTGT‑CTTGCCCGGTCAG‑CAATGTTAA‑CAGA‑GGAGATGGGTATGCAACGTCTGG‑CCACCATTGCA‑CC‑GCCACAGGTTCATGAGATAT‑GGGCGCTACTCAGCCAGATCCCGGACCCGGAGATCCCGGTGCTGAc                                                                                                                                                                                 <  1:537259/240‑1 (MQ=255)                                                                         cGGGATCAACGAAGCCACCTTGAACGTACCGCAAGAGCAGGCGTATCGCACTGGCGGAAAAAAA‑AGACAGCCTACCGAACACC‑TTGGACCCTTGATGGCAGAAATGCAGTATTTCCAGCGTGT‑CTTGCCCGGCCCG‑CCATGGTAA‑AAGA‑GGAGATGGGAATGCAACGACTGG‑CCACCATTGCA‑CC‑GCCACAGGTTCATGGGAACT‑GGGGCCGACTCAGCCAGATCCCGGAACCGAGGATCCCGGTGCTGAc                                                                                                                                                                                 >  2:294891/1‑251 (MQ=255)                                                                              ccaaggAAGCCCCATTAAACGTACCCCAAGAGCAGGCTTACCGCACGGGCGGTAAAAAA‑GGACTGCATACCGACCACC‑TTGGACCCAGGCTGGCAAAAATGCAGTATCTCCGGCGTGT‑CTTGCCCGGTCAG‑CAAGGGTAA‑CAGG‑GGAGAGGGGATTGCAACGTCTGG‑CCTCCATTGCA‑CC‑GCCACAGGTTCATGAGATAT‑GGGCGCTACTCAGCCAGATCCCGGACCCGGAGATCCCGGTGCTGACCATTa                                                                                                                                                                            <  1:49518/246‑1 (MQ=255)                                                                                accgaaccaacatggaacgtccccaaaGAGCGGGCGTATCGCACGGGGGGAAAAAAA‑GGACTGCATACCGAACCCC‑TGGGCCCCCGGCGGGAAGAAAGGCATTATCTCCAGTGTGC‑CTTGCC‑GGGCAGGAAAGGGAAA‑AAGA‑GGAGAGGGGTATGCAACTTCTGG‑CCCCCATTGCA‑CC‑GCCACAGGTTCATGAGATAT‑GGGCGCTACTCAGCCAGATCCCGGACCCGGAGATCCCGGTGCTGACCATTACt                                                                                                                                                                          <  2:219997/226‑1 (MQ=255)                                                                                  cgacgCCACATTGAACGTACCGCAAGAGCAGGCTTATCGCCCTGGCGGTAAAAAA‑GGACTGCATCCCGAACCCC‑TTGGCCCCATGCTGGCAGAAATGCAGTATCTCCAGCGTGT‑CTTGCCCGGTCAG‑CATTGTTAA‑CAGA‑GGAGATGGGTATGCACCTTCTGG‑CCCCCATTGCA‑CC‑GCCACAGGTTCATGAGATAT‑GGGCGCTACTCAGCCAGATCCCGGACCCGGAGATCCCGGTGCTGACCATTACTGa                                                                                                                                                                        <  1:596243/247‑1 (MQ=255)                                                                                   gAAGCCACATTGAACGTACCGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAA‑GGACTGCATACCGAACACC‑TTGGACCCATGCTGGCAGAAATGCAGTATCTCCAGCGTGT‑CTTGCCCGGTCAG‑CAATGGTAA‑CAGA‑GGAGATGGGTATGCAACGTCTGG‑CCACCATTGCA‑CC‑GCCACAGGTTCATGAGATAT‑GGGCGCTACTCAGCCAGATCCCGGACCCGGAGATCCCGGTGCTGACCATTACTGAt                                                                                                                                                                       >  1:233315/1‑251 (MQ=255)                                                                                       ccACATTGAACGTACCGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAA‑GGACTGCATACCGAAAACC‑TTGGACCCATGCTGGCCGAAATCCAGTATCTCCAGCGTGT‑CTTTCCCGGTCAG‑CAATGGTAA‑CAGA‑GGAGATGGGTATGCAACGTCTGG‑CCACCATTGCA‑CC‑GCCACAGGTTCATGAGATAT‑GGGCGCTACTCCGCCAGATCCCCGACCCGGAGATCCCGGTGCTGACCATTACTGAttttg                                                                                                                                                                   >  1:743375/1‑249 (MQ=255)                                                                                          caTTGAACGTCCCGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAA‑AGACAGCATACCGAACACC‑TTTGACCCCTGCTGGCAGAAATGCAGTATCTCCACCGTGT‑CTTGCCCGGACCG‑CCATGGTAA‑CAAA‑GGAAATGGGAATGCCAAGTCTGG‑CCACCAATGCA‑CC‑GCCACCGGTGAATTAG‑TATTGGGAGCGAATCACCCCAATCCCAGACCCGGAGATCCCGGtccttacccattctggtttaagc                                                                                                                                                                 >  2:569430/1‑228 (MQ=255)                                                                                           aTTGAACGTACCGCAAGAGCAGGCGTATCGCACTGGGGGTAAAAAA‑GGACTGCATACCGAACACC‑TTGGACCCATGCTGGCAGAAATGCAGTATCTCCAGGGTGT‑CTTGCCCGGTCAG‑CAATGGTAA‑CAGA‑GGAGATGGGTATGCAACGTCTGG‑CCACCATTGCA‑CC‑GCCACAGGTTCATGAGATAT‑GGGCGCTACTCAGCCAGATCCCGGACCCGGAGATCCCGGTGCTGACCATTACTGATTTAGGCAt                                                                                                                                                               <  1:110988/251‑1 (MQ=255)                                                                                              gaaggTACCGCAAGAGCAGGCGTACCGCACTGGCGGTAAAAAA‑GGACTGCATACCGAACACC‑TTGGACCCATGCTGGCAGAAATGCAGTATCTCCAGCGTGT‑CTTGCCCGGTCAG‑CAATGGTAA‑CAGA‑GGAGATGGGTAGGCAGCGTCTGG‑CCACCATTGCA‑CC‑GCCACAGGTTCATGAGATAT‑GGGCGCTACTCAGCCAGATCCCGGACCCGGAGATCCCGGTGCTGACCATTACTGATTTAGGCATGGt                                                                                                                                                            <  1:525305/247‑1 (MQ=255)                                                                                              aaaCGTACCGCAAGAGCAGGCTGATCGCACTGGGGGTAAAAAA‑GGCCTGCAGCCCGAACCCC‑TTGGCCCCATGCGGGCAGAAAGGAAGAACCTCCAGCGTGC‑CTTGCCCGGTAAG‑AAAGGGAAA‑CAGA‑GGAGATGGGTATGAAACGTCTGG‑CCCCCATTGCA‑CC‑GCCACAGGTTCATGAGATAT‑GGGCGCTACTCAGCCAGATCCCGGACCCGGAGATCCCGGTGCTGACCATTACTGATTTAGGCATGGt                                                                                                                                                            <  2:83367/250‑1 (MQ=255)                                                                                                     ccGCAAGGGCAGGTGTATCGCACTGGCGGTAAAAAA‑GGACTGCATAC‑GAACCCCCTTGACCCCATGCGGCCAGAAATGCGGTATCCCCAGCGTGC‑TTTGCCCGGTCAG‑CAATGGTAA‑AGGG‑GGAGATGGGTATGCAACGTTTGG‑CCACCATTGCA‑CC‑GCCACAGGTTCATGAGATAT‑GGGCGCTACTCAGCCAGATCCCGGACCCGGAGATCCCGGTGCTGACCATTACTGATTTAGGCATGGTGCGTAAt                                                                                                                                                     <  1:292496/251‑1 (MQ=255)                                                                                                      cGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAA‑GGACTGCATACCGAACACC‑TTGGACCCATGCTGGCAGAAATGCAGTATCTCCAGCGTGT‑CTTGCCCGGTCAG‑CAATGGTAA‑CAGA‑GGAGATGGGTATGCAACGTCTGG‑CCACCATTGCA‑CC‑GCCACAGGTTCATGAGATAT‑GGGCGCTACTCAGCCAGATCCCGGACCCGGAGATCCCGGTGCTGACCATTACTGATTTAGGCATGGTGCGTAAtg                                                                                                                                                    >  1:622220/1‑251 (MQ=255)                                                                                                       gCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAA‑GGACGGCATACCGAACCCC‑TTGGACCCATGCTGGCAGAAATGCAGTATCCCCAGCGTGT‑CTTGCCCGGTCAG‑AAATGGTAA‑CAGA‑GGAGATGGTTATGCAACGTCTGG‑CCCCCATTGCA‑CC‑GCCACAGGTTCATGAGATAT‑GGGCGCTACTCAGCCAGATCCCGGACCCGGAGATCCCGGTGCTGACCATTACTGATTTAGGCATGGTGCGTAAtgt                                                                                                                                                   <  1:403724/251‑1 (MQ=255)                                                                                                       gCAAGAGCAGGCGTACCGCACTGGCGGTAAAAAA‑GGACTGCATACCGAACCCC‑TTGGACCCATGCTGGCAAAAATGAAGTATCTCCAGCGTGT‑CTTGCCCGGTCAG‑CAATGGTAA‑CAGG‑GGAGAGGGGTATGCAACGTCTGG‑CCCCCATTCCA‑CC‑GCCACAGGTTCATGAGATAT‑GGGCGCTACTCAGCCAGATCCCGGACCCGGAGATCCCGGTGCTGACCATTACTGATTTAGGCATGGTGCGTAAtgt                                                                                                                                                   <  1:475866/251‑1 (MQ=255)                                                                                                          aggccAGGTGTAGCGCACTGG‑GGTAAAAAAATGTCTGCATACCAAACCCC‑GTGGCCCCATGCGGGAAGAAATGCAGTATCTCCAGAGGTT‑CTGTCCCGGTAAG‑CAAGGGTAA‑AAGA‑GGAGATGGGTATGCAACGTCGGG‑CCACCATGGCA‑CC‑GCCACAGGTTCATGAGATAT‑GGGCGCTACTCAGCCAGATCCCGGACCCGGAGATCCCGGTGCTGACCATTACTGATTTAGGCATGGTGCGTAATGTGac                                                                                                                                                <  2:680840/247‑1 (MQ=255)                                                                                                           gagCAGGCGTATCGCACTGGCGGTAAAAAA‑GGACTGCATACCGAACACC‑TTGGGCCCATGCTGGCAGAAAAGCAGTATCTCCAGCGTGT‑CTTGCCCGGTCAG‑CAATGGTAA‑CAGA‑GGAGATGGGTATGCAACGTCTGG‑CCACCATTGCA‑CC‑GCCACAGGTTCATGAGATAT‑GGGCGCTACTCAGCCAGATCCCGGACCCGGAGATCCCGGTGCTGAACATTACTGATTTAGGCATGGGGCGTAATGTGac                                                                                                                                                >  1:455232/1‑250 (MQ=255)                                                                                                                    tATCGCACTGGCGGTAAAAAA‑GGACTGCATACCGAACACC‑TTGGACCCATGCTGGCAGAAATGCAGTATCTCCAGCGTGT‑CTTACCCGGTCAG‑CAATGGTAA‑CAGA‑GGAGATGGGTATGCAACGTCTGG‑CCACCATTGCA‑CC‑GCCACAGGGTCATGAGATAT‑GGGCGCTACTCAGCCAGATCCCGGACCCGGAGATCCCGGTGCTGACCATTACTGATTTAGGCATGGTGCGGAATATGACACAGAttggg                                                                                                                                      >  1:418691/1‑247 (MQ=255)                                                                                                                                   aaaaaa‑GGACTGCATACCGAACCCC‑TTGGCCCCAGGCTGGCAGAAATGCAGTATCCCCCGCGTGT‑CTTGCCCGGTCAG‑CAATGGTAA‑CAGA‑GGAGATGGGTAGGCAACGTCTGG‑CCCCCATTGCA‑CC‑GCCACAGGTTCATGAGATAT‑GGGCGCTACTCAGCCAGATCCCGGACCCGGAGATCCCGGTGCTGACCATTACTGATTTAGGCATGGTGCGTAATGTGACACAGATGGGAGAAGGATGGGTGATc                                                                                                                       <  1:231952/251‑1 (MQ=255)                                                                                                                                    aaaaa‑GGACTGCATACCGAACACC‑TTGGACCCATGCTGGCAGAAATGCAGTATCTCCAGCGTGT‑CTTGCCCGGTCAG‑CAATGGTAA‑CAGA‑GGAGATGGGTATGCAACGTCTGG‑CCACCATTGCA‑CC‑GCCACAGGTTCATGAGATAT‑GGGCGCTACTCAGCCAGATCCCGGACCCGGAGATCCCGGTGCTGACCATTACTGATTTAGGCATGGTGCGTAATGTGACACAGATGGGAGAAGGATGGGTGATCg                                                                                                                      >  1:664249/1‑251 (MQ=255)                                                                                                                                                        acacC‑TTGGACCCATGCTGGCAGAAATGCAGTATCTCCAGCGTGT‑CTTGCCCGGTCAG‑CAATGGTAA‑CAGA‑GGAGATGGGTATGCAACGTCTGG‑CCACCATTGCA‑CC‑GCCACAGGTTCATGAGATAT‑GGGCGCTACTCAGCCAGATCCCGGACCCGGAGATCCCGGTGCTGACCATTACTGATTTAGGCATGGTGCGTAATGTGACACAGATGGGAGAAGGATGGGTGATCGGCTTTACGCCGACATATTc                                                                                                   >  1:413661/1‑251 (MQ=255)                                                                                                                                                          ccc‑TTGGCCCCAGGCTGGAAGAAATGCTGTATTCCCAGTGTGT‑TTTGCCCGGTCAG‑CAATGGTAA‑CAGA‑GGAGATGGGTTTGCAACGTCTGG‑CCCCCATTGCA‑CC‑GCCACAGGTTCATGAGATAT‑GGGCGCTACTCAGCCAGATCCCGGACCCGGAGATCCCGGTGCTGACCATTACTGATTTAGGCATGGTGCGTAATGTGACACAGATGGGAGAAGGATGGGTGATCGGCTTTACGCCGACATATTCCg                                                                                                 <  1:73510/250‑1 (MQ=255)                                                                                                                                                                ggACCCATGCTGGCAGAAATGCAGTATCTCCAGCGTGT‑CTTGCCCGGTCAG‑CAATGGTAA‑CAGA‑GGAGATGGGTATGCAACGTCTGG‑CCACCATTGCA‑CC‑GCCACAGGTTCATGAGATAT‑GGGCGCTACTCAGCCAGATCCCGGACCCGGAGATCCCGGTGCTGACCATTACTGATTTAGGCATGGTGCGTAATGTGACACAGATGGGAGAAGGATGGGTGATCGGCTTTACGCCGACATATTCCGGTTGt                                                                                            >  1:760080/1‑251 (MQ=255)                                                                                                                                                                  aCCCATGTTGGCAAAAATGCAGTATCTCCAGC‑TGTTCTTGCCCGGTCAG‑AACTGGTAA‑CAGG‑GGATATGGGTATGCAACTTTTGG‑CC‑CCCTTGCC‑CCCGCCACAGGTTCATGAGATAT‑GGGCGCTACTCAGCCAGATCCCGGACCCGGAGATCCCGGTGCTGACCATTACTGATTTAGGCATGGTGCGTAATGTGACACAGATGGGAGAAGGATGGTTGAACTGCTTTACGCCGACATATTCCGGTTGTcc                                                                                          <  1:664266/251‑1 (MQ=255)                                                                                                                                                                    ccaggcgggCAGATAGGCAGTATCTCCAGCGGTT‑CTGGCCCTGCCAG‑CAA‑GGTAAACGGA‑GGTTATGTGTATCCACCGTCGGC‑CCCCCATTGCG‑CC‑TCCACAGGTTCATGAGATAT‑GGGCGCTACTCAGCCAGATCCCGGACCCGGAGATCCCGGTGCTGACCATTACTGAGTTAGGCATGGTGCGTAATGTGACACAGATGGGAGAAGGATGGGTGATCGGCTTTACGCCGACATATTCCGGTTGTCCgg                                                                                        <  2:797420/244‑1 (MQ=255)                                                                                                                                                                        gCTGGCAGAAATGCAGTATCTCCGGCGTGT‑CTGGCCCGGTCAG‑CAATGGTAA‑CAGG‑GGAGAGGGGTATGCACCGTCTGG‑CCACCATTGCA‑CC‑GCCACAGGTTCATGAGATAT‑GGGCGCTACTCAGCCAGATCCCGGACCCGGAGATCCCGGTGCTGACCATTACTGATTTAGGCATGGTGCGTAATGTGACACAGATGGGAGAAGGATGGGTGATCGGCTTTACGCCGACATATTCCGGTTGTCCGGCAAc                                                                                    <  1:690168/251‑1 (MQ=255)                                                                                                                                                                        gCTGGCAAAAATCCCGTATCTCCCGCTTGT‑CTTGCCCGTTCAG‑CAATGCTAA‑CAGA‑GGATATGGGTATGCAACTTCTGC‑CCCCCTTTCCA‑CC‑GCCACAGGTTCATGAGATAT‑GGGCGCTACTCAGCCAGATCCCGGACCCGGAGATCCCGGTGCTGACCATTACTGATTTAGGCATGGTGCGTAATGTTACACAGATGGGAGAAGGATGGGTGATCGGCTTTACGCCGACATATTCCGGTTGTCCGGCAAc                                                                                    <  1:343103/251‑1 (MQ=255)                                                                                                                                                                            gCAGAAATGCAGTATCTCCAGCGTGT‑CTTGCCCGGTCAG‑CAATGGTAA‑CAGA‑GGAGATGGGTATGCAACGTCTGG‑CCACCATTGCA‑CC‑GCCACAGGTTCATGAGATAT‑GGGCGCTACTCAGCCAGATCCCGGACCCGGAGATCCCGGTGCTGACCATTACTGATTTAGGCATGGTGCGTAATGTGACACAGATGGGAGAAGGATGGGTGATCGGCTTTACGCCGACATATTCCGGTTGTCCGGCAACGGa                                                                                 >  1:431241/1‑250 (MQ=255)                                                                                                                                                                                aaaTGCAGTATCTCCAGCGTGT‑CTTGACCGGTCAG‑CAATGGTAA‑CAGA‑GGAGATGGGTCTGCCAAGTCTGG‑CCACCATTGCA‑CC‑GCCACAGGTTCATAAGATAT‑GGGCGCTACTCAGCCAGATCCCGGACCCGGAGATCCCGGTGCTGAACATTACTGATTTGGGCATGGTGCGTAATGTGAAACAGATGGGCGAAGGCTGGGTGCACGGATTTACGCCGCCATATTACGGTTGTCCGGCACCGGCACAtt                                                                            >  2:782708/1‑251 (MQ=255)                                                                                                                                                                                 aaggCAGTATCGCCAGTGTGT‑TTGGCCCGGTCAG‑CAATGTAAA‑CAGG‑GGAGATGGGTATGCAACGTCTGG‑CCCCCATGGCA‑CC‑GCCACAGGTTCATGAGATAT‑GGGCGCTACTCAGCCAGATCCCGGACCCGGAGATCCCGGTGCTGACCATTACTGATTTAGGCATGGTGCGTAATGTGACACAGATGGGAGAAGGATGGGTGATCGGCTTTACGCCGACATATTCCGGTTGTCCGGCAACGGAACAttt                                                                           <  1:442324/248‑1 (MQ=255)                                                                                                                                                                                     gtgTATCCCCAGCGGTT‑CTTGCCTGGTCAG‑AAATGTTAA‑C‑GATGGAGATGGTTAGTCACCGGCTGG‑CCCCCATTGCA‑CC‑GCCACAGGTTCATGAGAGAT‑GGGCGCTACTCAGCCAGATCCCGGACCCGGAGATCCCGGTGCTGACCATTACTGATTTAGGCATGGTGCGTAATGTGACACAGATGGGAGAAGGATGGGTGATCGGCTTTACGCCGACATATTCCGGTTGTCCGGCAACGGAACATTTGAtt                                                                       <  2:492017/249‑1 (MQ=255)                                                                                                                                                                                          tctcCAGCGTGT‑CTAGCCCGGGCAG‑CAATGGTAA‑CAGA‑GGCGATGGGTATGCAACGTCTGG‑CCACCATTGCA‑CC‑GCCACAGGTTCATGAGATAT‑GGGCGCTCCTCAGCCAGATCCCGGCCCCGGAGATCCCGGTGCTGACCATTACTGATTTAGGCATGGTGCGTAATGTGACACAGATGGGCGAAGGATGGGGGATCGGCTTTACGCCGACATATTCCGGTTGTCGGGCAACGGAACATTTGATTGgcgg                                                                  >  2:780203/1‑250 (MQ=255)                                                                                                                                                                                              ccatcgtgcctgtcgaggTCAG‑AAATGGAAA‑CAGA‑GGAGTGGGGGAGAAAACGTCGGG‑CCACCATTGCA‑CC‑GCCACAGGTTCATGAGATAT‑GGGCGCTACTCAGCCAGATCCCGGACCCGGAGATCCCGGTGCTGACCATTACTGATTTAGGCATGGTGCGTAATGTGACACAGATGGGAGAAGGATGGGTTATCGGCTTTACGCCGACATATTCCGGTTGTCCGGCAACGGAACATTTGATTGgcg                                                                   <  1:298954/231‑1 (MQ=255)                                                                                                                                                                                              cAGCGTGT‑CTTGCCCGGTCAG‑CAATGGTAA‑CAGC‑GGAGATGGGTATGCAACGTCTGG‑CCACCATTGCA‑CC‑GCCACAGGTTCATGAGATAT‑GGGCGCTACTCAGCCAGATCCCGGACCCGGAGATCCCGGTGCTGACCATTACTGATTTAGGCATGGTGCGTAATGTGAAACAGATGGGAGAAGGATGGGTGATCGGCTTTACGCCGACATATTCCGGTTGTCCGGCAACGGaaggtcggaggagcgtcgtg                                                              >  2:762319/1‑233 (MQ=255)                                                                                                                                                                                               aGCGTGT‑CTTGCCCGGTCAG‑CAATGGTAA‑CAGA‑GGAGATGGGTATGCAACGTCTGG‑CCACCATTGCA‑CC‑GCCACAGGTTCATGAGATAT‑GGGCGCTACTCAGCCAGATCCCGGACCCGGAGATCCCGGTGCTGACCATTACTGATTTAGGCATGGTGCGTAATGTGACACAGATGGGAGAAGGATGGGTGATCGGCTTTACGCCGACATATTCCGGTTGTCCGGCAACGGAACATTTGATTGGCGCGATAc                                                             >  2:220933/1‑251 (MQ=255)                                                                                                                                                                                                    tt‑CTTGCCCGGTAAG‑CAAGGGTAA‑CAGA‑GGAGATGGGTATGCACCGTCTGT‑CCCCCATTGCA‑CC‑GCCACAGGTTCATGAGATAT‑GGGCGCTACTCAGCCAGATCCCGGACCCGGAGATCCCGGTGCTGACCATTACTGATTTAGGCATGGTGCGTAATGTGACACAGATGGGAGAAGGATGGGTGATCGGCTTTACGCCGACATATTCCGGTTGTCCGGCAACGGAACATTTGATTGGCGCGATACGTGa                                                         <  1:753839/249‑1 (MQ=255)                                                                                                                                                                                                               ccggcgaacaagggcaaagggggggaTGGGTAGGCAACGGCGGC‑CCCCCATGGCA‑CC‑GCCACAGGTTCATGAGATAT‑GGGCGCTACTCAGCCAGATCCCGGACCCGGAGATCCCGGTGCTGACCATTACTGATTTAGGCATGGTGCGTAATGTGACACAGATGGGAGAAGGATGGGTGATCGGCTTTACGCCGACATATTCCGGTTGTCCGGCAACGGAACATTTGATTGGCGCGATACGTGAGGCAATGa                                                 <  2:49282/227‑1 (MQ=255)                                                                                                                                                                                                                      aagggaaa‑CAGA‑GGAGATGGGTAGTCAAGGTCGGG‑CCGCCAT‑GCAACC‑GCCACAGGTTCATGAGATAT‑GGGCGCTACTCAGCCAGATCCCGGACCCGGAGATCCCGGTGCTGACCATTACTGATTTAGGCATGGTGCGTAATGTGACACAGATGGGAGAAGGATGGGTGATCGGCTTTACGCCGACATAGTCCGGTTGTCCGGCAACGGAACATTTGATTGGCGCGATACGTGAGGCAATGACAACCAACg                                        <  1:351820/245‑1 (MQ=255)                                                                                                                                                                                                                                    ggagtgggttatccactggcTGG‑CCCCCTTTGCA‑CC‑GCCACAGGTTCATGAGATAT‑GGGCGCTACTCAGCCAGATCCCGGACCCGGAGATCCCGGTGCTGACCATTACTGATTTAGGCATGGTGCGTAATGTGACACAGATGGGAGAAGGATGGGTGATCGGCTTTACGCCGACATATTCCGGTTGTCCGGCAACGGAACATTTGATTGGCGCGATACGTGAGGCAATGACAACCAACGGGTTT‑ACCCCCg                           <  2:244866/232‑1 (MQ=255)                                                                                                                                                                                                                                     gagaTGGGTATGCAACGTCTGG‑CCACCATTGCA‑CC‑GCCACAGGTTCATGAGATAT‑GGGCGCTACTCAGCCAGATCCCGGACCCGGAGATCCCGGTGCTGACCATTACTGATTTAGGCATGGTGCGTAATGTGACACAGATGGGAGAAGGATGGGTGATCGGCTTTACGCCGACATATTCCGGTTGTCCGGCAACGGAACATTTGATTGGCGCGATACGTGAGGCAATGACAACCAACGGGTTT‑ACCCCCGt                          >  2:165764/1‑251 (MQ=255)                                                                                                                                                                                                                                             gaggcatcGTCTGG‑CCACCATTGCA‑CC‑GCCACAGGTTCATGAGATAT‑GGGCGCTACTCAGCCAGATCCCGGACCCGGAGATCCCGGTGCTGACCATTACTGATTTAGGCATGGTGCGTAATGTGACACAGATGGGAGAAGGATGGGTGATCGGCTTTACGCCGACATATTCCGGTTGTCCGGCAACGGAACATTTGATTGGCGCGATACGTGAGGCAATGACAACCAACGGGTTT‑ACCCCCGTTCAGGTtg                  <  1:717652/244‑1 (MQ=255)                                                                                                                                                                                                                                              agacaactccggg‑CCCCCATGGCC‑CC‑GCCACAGGTTCATGAGATAT‑GGGCGCTACTCAGCCAGATCCCGGACCCGGAGATCCCGGTGCTGACCATTACTGATTTAGGCATGGTGCGTAATGTGACACAGATGGGAGAAGGATGGGTGATCGGCTTTACGCCGACATATTCCGGTTGTCCGGCAACGGAACATTTGATTGGCGCGATACGTGAGGCAATGACAACCAACGGGTTT‑ACCCCCGTTCAGGTtg                  <  2:622220/239‑1 (MQ=255)                                                                                                                                                                                                                                                gCAACGTCTGG‑CCACCATTGCA‑CC‑GCCACAGGTTCATGAGATAT‑GGGCGCTACTCAGCCAGATCCCGGACCCGGAGATCCCGGCGCTGACCCTTACTGCTTTAGGCATGGTGCGTAATGTGACACAGATGGGAGAAGGATGGGTGATCGGCTTTACGCCGACATATTCCGGTTGTCCGGCAACGGAACATTTGATTGGCGCGATACGTGAGGCAATGACAACCAACGGGTTT‑ACCCCCGTTCAGGTTGtgc               >  1:750764/1‑251 (MQ=255)                                                                                                                                                                                                                                                     gTCTGG‑CCACCATTGCA‑CC‑GCCACAGGTTCATGAGATAT‑AGGCGCTACTCAGCCAGCTCCCGGACCCGGGGCTCCCGGTGCTGACCATTACTGATTTAGGCATGGTGCGTAATGTACCACAGCTGGGAGAAGGATGCGTGATCGGCTTTACGCCGACATATTCCGGTTGTCCTGCATCGTAACATTTAATTGGCGCTATCCGTGAGGCAGTGACACCCAACGGGTTTCACACCC‑TTCAGGTTGTGCTGCac          >  2:788717/1‑250 (MQ=255)                                                                                                                                                                                                                                                         gg‑CCACCATTGCA‑CC‑GCCACAGGTTCATGAGATAT‑GGGCGCTACTCAGCCAGATCCCGGACCCGGAGATCCCGGTGCTGACCATTACTGATTTAGGCATGGTGCGTAATGTGACACAGATGGGAGAAGGATGGGTGATCGGCTTTACGCCGACATATTCCGGTTGTCCGGCAACGGAACATTTGATTGGCGCGATACGTGAGGCAATGACAACCAACGGGTTT‑ACCCCCGTTCAGGTTGTGCTGCAACTCg      <  1:632716/251‑1 (MQ=255)                                                                                                                                                                                                                                                              accaTTGCA‑CC‑GCCACAGGTTCATGAGATAT‑GGGCGCTACTCAGCCAGATCCCGGACCCGGAGATCCCGGTGCTGACCATTACTGATTTAGGCATGGTGCGTAATGTGACACAGATGGGAGAAGGATGGGTGATCGGCTTTACGCCGACATATTCCGGTTGTCCGGCAACGGAACATTTGATTGGCGCGATACGTGAGGCAATGACAACCAACGGGTTT‑ACCCCCGTTCAGGTTGTGCTGCAACTCGAccc  >  1:222127/1‑251 (MQ=255)                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                               GCACTGAGTGAAGAGGGTATTGCGGTTGATCCACGCACTTTACGCGCAGCGTG‑GGAAGCCGAAGTTTTTGCCGGGATCAACGAAGCCACATTGAACGTACCGCAAGAGCAGGCGTATCGCACTGGCGGTAAAAAA‑GGACTGCATACCGAACACC‑TTGGACCCATGCTGGCAGAAATGCAGTATCTCCAGCGTGT‑CTTGCCCGGTCAG‑CAATGGTAA‑CAGA‑GGAGATGGGTATGCAACGTCTGG‑CCACCATTGCA‑CC‑GCCACAGGTTCATGAGATAT‑GGGCGCTACTCAGCCAGATCCCGGACCCGGAGATCCCGGTGCTGACCATTACTGATTTAGGCATGGTGCGTAATGTGACACAGATGGGAGAAGGATGGGTGATCGGCTTTACGCCGACATATTCCGGTTGTCCGGCAACGGAACATTTGATTGGCGCGATACGTGAGGCAATGACAACCAACGGGTTT‑ACCCCCGTTCAGGTTGTGCTGCAACTCGACCC  >  NC_000913/1455695‑1456190

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