Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	12,484	C→A	19.0%	Q108K (CAG→AAG)	dnaK →	chaperone Hsp70, with co‑chaperone DnaJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	12,484	0	C	A	19.0%	105.6 / 2.8	68	Q108K (CAG→AAG)	dnaK	chaperone Hsp70, with co‑chaperone DnaJ
Reads supporting (aligned to +/- strand):  ref base C (44/10);  new base A (5/9);  total (49/19)
Fisher's exact test for biased strand distribution p-value = 1.60e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTAGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTT‑T‑CC‑A‑TCA‑TGCCGTTCAAAATTATTGCTGCTGATAACGGCG‑AC‑GCAT‑GGGTCGAAGTTAAAG‑GCCAGAAAATGG‑CA‑CCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGAC  >  NC_000913/12259‑12708                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                   gATCGCACCACGCCTTCTATCATTGCCTATGCCCAGGATGGTGAAACTCTAGTTGGTCAGCCGGCTAATCGTCAGGCAGTGACGAACCCGCAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTT‑T‑CC‑A‑TCA‑TGCCGTTCAAAATTATTGCTGCTGATAACGGCG‑AC‑GCAT‑GGGTCGAAGTTAAAG‑GCCa                                                                                                                                                                                                                                   <  2:627922/227‑1 (MQ=255)    cGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAACCTCTAGTGGGTCAGCCGGCTAAACGTCAGTCAGTGACGAACCCACAAAACACTCACTTTGCGATTAAACGCCTGGTTGGTAGACGCTTCCAGGCCGAAGCAGCACATCGTAATCTT‑T‑CC‑C‑TCT‑TGCCGTTCCAGAGTTTTGCTGCTGAAAACAGCG‑AC‑CCAT‑TGGACCACGTGAAAG‑GCCAGaaa                                                                                                                                                                                                                               >  1:190188/1‑228 (MQ=255)      caccacGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTAGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTT‑T‑CC‑A‑TCA‑TGCCGTTCAAAATTATTGCTGCTGATAACGGCG‑AC‑GCAT‑GGGTCGAAGTTAAAG‑GCCAGAAAAt                                                                                                                                                                                                                             >  2:457403/1‑228 (MQ=255)         cacGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTAGTTGGTCAGCCGGCTAAACGTCAGCCAGTGACGAACCCGCAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGTAGAAGTACAGCGTGAGGTT‑T‑CC‑A‑TCA‑TGCCGTTCCAAATTATTGCTGCTGATAACGGCG‑AC‑CCCT‑GGGGCGAAGTTAAAG‑GCCGAAAAATGG‑c                                                                                                                                                                                                                         >  1:251376/1‑228 (MQ=255)          acGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTAGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTCCAGCGTGATGTT‑G‑CC‑C‑TCA‑TGCCGTTCAAACTTATTGCTGCTGATAACGGCG‑AC‑CCAT‑GGGGCGAAGTTAAAG‑GCAaaaaaacagga                                                                                                                                                                                                                         >  1:693806/1‑223 (MQ=255)                     catcgccgtaaCCCAGGTTGGTGAATCTATAGTTGGCGAGCCGCATAAACGCCAGGCAGTGACGCCCCCGAAAACGCCGCTGTTGGCGATAAAACGCCTGTTGGGTGGCCTCTTCAAGGACGAAGAGGACCGGCGTG‑TGTT‑T‑TC‑A‑AGACGGCGGTAAAAAATTATTGCGGCTAATACCGGCG‑AC‑GTTT‑GGGTCAAATTTAAGG‑GCCAGAAAATGG‑CA‑CCGCCGCAGAt                                                                                                                                                                                                            <  1:444985/218‑1 (MQ=255)                       ttGCCTATACCCAGGATGGTGAAACTCTAGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTT‑T‑CC‑A‑TCA‑TGCCGTTCAAAATTATTGCTGCTGATAACGGCG‑AC‑GCAT‑GGGTCGAAGTTAAAG‑GCCGGAAAATGG‑CA‑CCGCCGCAGAttt                                                                                                                                                                                                          >  2:560096/1‑228 (MQ=255)                        tGCCTATACCCAGGATGGTGAAACTCTAGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTT‑T‑CC‑A‑TCA‑TGCCGTTCAAAATTATTGCTGCTGATAACGGCG‑AC‑GCAT‑GGGTCGAAGTTAAAG‑GCCAGAAAATGG‑CA‑CCGCCGCAGAtt                                                                                                                                                                                                           >  2:467535/1‑226 (MQ=255)                          ccTATACCCAGGATGGTGAAACTCTAGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACAATCGGTTTGCGATTAAACGCCTGATTGGTCGCAGCTTCCAGGACGAAGAAGTACAGCGTGATGTT‑T‑CC‑A‑TCA‑TTCCGTTCAAAATTATTGGTGCTGATAACGCCG‑AC‑GTCT‑GGGTAGACGTGAAAG‑GCCAGGAAATGG‑CA‑CCCCCGCAGATTTttg                                                                                                                                                                                                       >  2:329439/1‑228 (MQ=255)                                    ggATGGTGAAACTCTAGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTT‑T‑CC‑T‑TCA‑TGCCGTTCAAAATTATTGATGCTGATAACGGCG‑AA‑GCAA‑GGGCCGAAGTAAAAG‑GCAGGAAAATGG‑AA‑CAGCCGCAGGTTTCTGCTGAAGTGCt                                                                                                                                                                                             >  2:678853/1‑228 (MQ=255)                                    ggATGGTGAAACTCTAGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTT‑T‑CC‑A‑TCA‑TGCCGTTCAAAATTATTGCTGCTGATAACGGCG‑AC‑GCAT‑GGGTCGAAGTTAAAG‑GCCAGAAAATGG‑CA‑CCGCCGCAGATTTCTGCTGAAGCGCt                                                                                                                                                                                             >  2:464580/1‑228 (MQ=255)                                    ggATGGTGAAACTCTAGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACTCTGTTTGCGATTAAACGCATGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTT‑T‑CC‑A‑TCA‑TGACGTTCAAAATTATTGCAGCTGATAACGGCA‑CC‑GCCT‑GGGTCGAAGTTTAAA‑GCCAGAAAATGG‑CA‑CCGCCGCAGATTTCTGGAGAAGTGcg                                                                                                                                                                                             >  1:330720/1‑227 (MQ=255)                                     gATGGTGAAACTCTAGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTT‑T‑CC‑A‑TCA‑TGCCGTTCAAAATTATTGCTGCTGATAACGGCG‑AC‑GCAT‑GGGTCGAAGTTAAAG‑GCCTGTAAATGG‑CA‑CCGCCGCAGTTTTCTGCTGACgtggtg                                                                                                                                                                                            >  2:645859/1‑228 (MQ=255)                                        ggTGAAACTCTAGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTT‑T‑CC‑A‑TCA‑TGCCGTTCAAAATTATTGCTGCTGATAACGGCG‑AC‑GCAT‑GGGTCGAAGTTAAAG‑GCCAGAAAATGG‑CA‑CCGCCGCAGATTTCTGCTGAAGTGCTGaaa                                                                                                                                                                                         >  2:307973/1‑228 (MQ=255)                                        ggTGAAACTCTAGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACACTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTT‑T‑CC‑A‑TCA‑TGCCGTTCAAAATTATTGCTGCTGATAACGGCG‑AC‑GCAT‑GGGTCGCAGTTAAAG‑TCCGGAAAATGG‑CA‑CCGCCGCAGATTTCTGCTGAAGTGCTGaaa                                                                                                                                                                                         >  2:259710/1‑228 (MQ=255)                                         gTGAAACTCTAGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTT‑T‑CC‑A‑TCA‑TGCCGTTCAAAATTATTGCTGCTGATAACGGCG‑AC‑GCAT‑GGGTCGAAGTTAAAG‑GCAAGAAAATGG‑CA‑CCGCCGCAGATTTCTGCTGAAGGGCGGagaa                                                                                                                                                                                        >  2:51295/1‑228 (MQ=255)                                           gAAACTCTAGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTT‑T‑CC‑A‑TCA‑TGCCGTTCAAAATTATTGCTGCTGATAACGGCG‑AC‑GCAT‑GGGTCGAAGTTAAAG‑GCCAGAAAATGG‑CA‑CCGCCGCAGATTTCTGCTGAAGTGCTGaaaaa                                                                                                                                                                                       >  2:351287/1‑227 (MQ=255)                                               ctctAGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGACGAAGTACAGCGTGATGTT‑T‑CC‑A‑TCA‑TTCCGTTCCAAATTATTGCTGCTGATAACGGCG‑AC‑GCCT‑GGGTCGTAGTTAAAA‑GCCAGAAAATTG‑CA‑ACGCCGCAGACTTCTGCCGAAGGGCTGaaaaaaaaaa                                                                                                                                                                                  >  1:236397/1‑225 (MQ=255)                                               ctctAGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTT‑T‑CC‑A‑TCA‑TGCCGTTCAAAATTATTGCTGCTGATAACGGCG‑AC‑GCAT‑GGGTCGAAGTTAAAG‑GCCAGAAAATGG‑CA‑CCGCCGCAGATTTCTGCTGAAGGGCTGAAAAAAATga                                                                                                                                                                                  >  2:289296/1‑228 (MQ=255)                                                   aGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACTCTGTTTGCGATTAACCGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTT‑T‑CC‑A‑TCA‑TGCCGTTCAAAATTATTGCTGCTGAAAACGGCG‑AC‑GCAT‑GGGACGAAGTTAAAG‑GCCAGAAAATGG‑CA‑CCGCCGCAGATTTTTGCTGGAGGGCGGAAAAAAAAgaagaa                                                                                                                                                                              >  2:328794/1‑228 (MQ=255)                                                   aGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGGTGTT‑C‑CC‑A‑TCA‑TGCCGTTCAAAATTATTGCTGCTGATAACGGCG‑AC‑GCAT‑GGGTCGACGTTAAAG‑GCCAGACAATGG‑CA‑CCCCCGCAGTTTTCTGCTGAAGTGCTGAACAAACTcaagaa                                                                                                                                                                              >  2:166320/1‑228 (MQ=255)                                                          cAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTT‑T‑CC‑A‑TCA‑TGCCGTTCAAAATTATTGCTGCTGATAACGGCG‑AC‑GCAT‑GGGTCGAAGTTAAAG‑GCCAGAAAATGG‑CA‑CCGCCGCAGATTTCTGCTGAAGTGCTGGAAAAAATGAAGAAAAcc                                                                                                                                                                          >  2:114351/1‑225 (MQ=255)                                                            gCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTT‑T‑CC‑A‑TCA‑TGCCGTTCAAAATTATTGCTGCTGATAACGGCG‑AC‑GCAT‑GGGTCGAAGTTAAAG‑GCCAGAAAATGG‑CA‑CCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGACGAAAACCGCTGa                                                                                                                                                                     >  2:238690/1‑228 (MQ=255)                                                            gCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTT‑T‑CC‑A‑TCA‑TGCCGTTCAAAATTATTGCTGCTGATAACGGCG‑AC‑GCAT‑GGGTCGAAGTTAAAG‑GCCAGAAAATGG‑CA‑CCGCCCCAGATTTCTGCTGAAGTGCTGAAAAAAATAAAGAAAACCACTGa                                                                                                                                                                     >  1:346924/1‑228 (MQ=255)                                                             ccGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTT‑T‑CC‑A‑TCA‑TGCCGTTCAAAATTATTGCTGCTGATAACGGCG‑AC‑GCAT‑GGGTCGAAGTTAAAG‑GCCGGAAAATGG‑CA‑CCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGaa                                                                                                                                                                    >  2:619972/1‑228 (MQ=255)                                                                         ggggCTGTGCCGACCCCAAAAATCACGTTTTTTGCGTTTAAACGTTTGTGTGTCCGCCGCTCCCAGGCAATAGATGTACAGTG‑GAAGTTGT‑CC‑A‑TCA‑TGCCTTTAAAACTTGTGCCTACTAATAACGCCG‑CC‑TCTG‑GGGTAGAA‑TTAAAGAGACAAAAAAGTG‑CA‑CCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGt                                                                                                                                                        <  1:560096/226‑1 (MQ=255)                                                                                         cacaaaccgctgtttttGCTATTAACGGCCTGTTGGGGCTCCTTTCCAAGTACGAAGAAGCACAGCGGGAGTTT‑C‑CC‑A‑TAT‑GGGCGTTAAAAATAATTGCTCTTGA‑AACCGCG‑GCAGAAG‑GGGGCGAAGTAAAGG‑AACAAAAAATGG‑CA‑CCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAg                                                                                                                                        <  1:418333/214‑1 (MQ=255)                                                                                          gCAAAATACCCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGCACAGAGTACTATT‑T‑CC‑A‑TTA‑TGCCGTTCAAAATTATTGCCGCTGTTAACGGCG‑AC‑GCAT‑GGGTCGACGTTAAAG‑GCCGAACATTGG‑CC‑CCCCCGCAGGTTCCTGCCGACGTCCTGAAGACAATGTACAAAACCCCTGCCGCTTTCCAGGGGGAACCGGAAAATAAAgg                                                                                                                                       >  1:124879/1‑227 (MQ=255)                                                                                          gCAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTT‑T‑CC‑A‑TCA‑TGCCGTTCAAAATTATTGCTGCTGATAACGGCG‑AC‑CCAT‑GGGTCGAAGTTAAAG‑TCCAGAAAATGG‑CA‑CCGCCCCAGATTTCTGCTGAGTTGCTGAAAAAAAGGACAAAAAACGCTGAAGATACCCTGGGGGAACGGGaaaaaaaaga                                                                                                                                       >  1:602336/1‑221 (MQ=255)                                                                                                cattctTTTTGCGATTAAAGGCTTGTTGGGCCCCCGTTTCAAGGACGAAGAAGTAAGGCTTGATGTT‑T‑CA‑A‑TAA‑GGCGTTTAAAAATAATTGTTGCTGATAACGGCG‑AC‑GAAG‑GGGTCAAAGTAAAAG‑GCAAGAAAATGC‑CA‑CCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTAt                                                                                                                                 <  2:342195/225‑1 (MQ=255)                                                                                                     tGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTT‑T‑CC‑A‑TCA‑TGCCGTTCAAAATTATTGCTGCTGATAACGGCG‑AC‑GCAT‑GGGTCGCAGTTAAAG‑GCCAGAACATGG‑CA‑CCGCCGCAGATTTCTGCTGAAGTGCTGGAGAAAATAAAGAAAACCGCTGAAGATTACCATGGTGAACCGGAAACAGaaacagtaaaccc                                                                                                                              >  1:333146/1‑215 (MQ=255)                                                                                                     tGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTT‑T‑CC‑A‑TCA‑TGCCGTTCAAAATTATTGCTGCTGATAACGGCG‑AC‑GCAT‑GGGTCGAAGTTAAAA‑GCCATGAAAAGG‑AA‑CCGCCGCAGATTTCTGCTGAAGGGCTGGAAAAAATGAAGAAAACCGCGGACGATTTCCTGGGTGACCCGGAAACAGAAGCGGTAATccc                                                                                                                              >  1:493885/1‑224 (MQ=255)                                                                                                       ttgccGATTAACCGCTGGATTGGCCGCCGTTTCCAGGACAAAGAATTACAGC‑TGGTTTT‑T‑CCAA‑TAA‑GGCCGTCAAAAATAATTGCTGCTAATAACGGCG‑AC‑GAAT‑GGGTCAAATTAAAAG‑GCAAAAAAAGGG‑CA‑CCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTa                                                                                                                          <  2:346924/224‑1 (MQ=255)                                                                                                             tttaaaccgttgtttGGTCGCCCTTACCAGGACGAAGAACGACAGTGGGATTTT‑T‑CC‑T‑GCA‑TGCGATTAAAAATGATTGCTGCTGATAACGGGG‑AC‑GCAT‑GGGTGAAAGTTAAAA‑ACAAAAAAAAGA‑AA‑CCGCCGAAGATTTCTGCTGAAGTGCCGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCa                                                                                                                    <  1:51295/215‑1 (MQ=255)                                                                                                             ttaaaactggtgttgtgccgccgtTCCCAGGAAGAAAAAGAAAAGCG‑GAGATT‑T‑CC‑ACGCA‑TGCATTAAAAAATTATTGTGGCTGACAACGGCG‑AG‑GTAG‑GGGGCAAAGTTAAAG‑GCAAAAAAGTGG‑CA‑CCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCa                                                                                                                    <  1:307973/210‑1 (MQ=255)                                                                                                             aTTAGACGCATGATTGGTGGCCGTTCTCAGGGAGAGGAAGTACCGCGTTATTTT‑T‑CC‑A‑GCA‑TGCCGTTCAAAATTATTGCTGCTGATAACGGCG‑AC‑GAAT‑GGTGCAAATTAAAAG‑GCAAGAAAATGG‑CA‑CCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCa                                                                                                                    <  2:134929/228‑1 (MQ=255)                                                                                                                 aaCGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTT‑T‑CC‑A‑TCA‑TGCCGTTCAAAATTATTGCTGCTGATAACGGCG‑AC‑GCAT‑GGGTCGAAGTTAAAG‑GCCAGAAAATGG‑CA‑CCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTAACGGCATACt                                                                                                                >  2:330445/1‑228 (MQ=255)                                                                                                                        gATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTT‑T‑CC‑A‑TCA‑TGCCGTTCAAAATTATTGCTGCTGATAACGGCG‑AC‑GCAT‑GGGTCGAAGTTAAAG‑GCAAGAAAATGG‑CA‑CCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGGAAAAGAAGCGGATATCACCGAACCGGACAACTTTAAAGa                                                                                                         >  1:287868/1‑228 (MQ=255)                                                                                                                            ggTCGCCGCTGCCAGGAGAAAGAAGTAAGGCGTG‑TGTT‑TGCC‑C‑ATA‑GGCCGTAAAAAATTGTTGCTGCTGAAACCGGCG‑AC‑GATT‑GGGTCGAAGTAAAAG‑GCAAAAAAAGGG‑CA‑CCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCt                                                                                                     <  1:467535/228‑1 (MQ=255)                                                                                                                              tCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTT‑T‑CC‑A‑TCA‑TGCCGTTCAAAATTATTGCTGCTGATAACGGCG‑AC‑GCAT‑GGGTCGAAGTTAAAG‑GCCAGAAAATGG‑CA‑CCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGATGTTATAACCGTACCGGCCTACTTTAACGAtcccc                                                                                                    >  2:690237/1‑223 (MQ=255)                                                                                                                                gccgcTTCCAGGACGAAGAAGTACAGCGTGATGTT‑T‑CC‑A‑TCA‑TGCCGTTCAAAATTATTGCTGCTGATAACGGCG‑AC‑GCAT‑GGGTCGAAGTTAAAG‑GCCAGAAAATGG‑CA‑CCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGc                                                                                                 >  2:47424/1‑228 (MQ=255)                                                                                                                                      tCCAGGACGAAAAAGTACAGCTGGATTTT‑T‑CA‑A‑ACA‑GGCCTGTAAAAATTATTGCTGCGGATAACGGCG‑AC‑GCAT‑GGGTCAAAGTTAAAG‑GCAAAAAAATGG‑CA‑CCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAgg                                                                                           <  2:389705/228‑1 (MQ=255)                                                                                                                                                     tACAGCGTGATGTT‑T‑CC‑A‑TCA‑TGCCGTTCAAAATTATTGCTGCTGATAACGGCG‑AC‑GCAT‑GGGTCGAAGTTAAAG‑GCCAGAAAATGG‑CA‑CCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATCACCTGGGTGAACGGGTACCTGAAGCTGTTATAACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGAAACCAAAgcag                                                                               >  2:459437/1‑222 (MQ=255)                                                                                                                                                         gCGTGATGTT‑T‑CC‑A‑TCA‑TGCCGTTCAAAATTATTGCTGCTGATAACGGCG‑AC‑GCAT‑GGGTCGAAGTTAAAG‑GCCAGAAAATGG‑CA‑CCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGc                                                                          >  2:567809/1‑226 (MQ=255)                                                                                                                                                          cGTGATGTT‑T‑CC‑A‑TCA‑TGCCGTTCAAAATTATTGCTGCTGATAACGGCG‑AC‑GCAT‑GGGTCGAAGTTAAAG‑GCCAGAAAATGG‑CA‑CCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACGTAAGCGGTTATAACCGTACCGGCACACTCTAACGATGCTCAGCGTCAGGCACCCAAAAACCCCGGc                                                                          >  2:319162/1‑225 (MQ=255)                                                                                                                                                            tGATGTT‑T‑CC‑A‑TCA‑TGCCGTTCAAAATTATTGCTGCTGATAACGGCG‑AC‑GCAT‑GGGTCGAAGTTAAAG‑GCCAGAAAATGG‑CA‑CCGCCGCAGATTTCTGCTGAAGTGCTGACAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTAt                                                                     >  2:510116/1‑228 (MQ=255)                                                                                                                                                            tGATGTT‑T‑CC‑A‑TCA‑TGCCGTTCAAAATTATTGCTGCTGATAACGGCG‑AC‑GCAT‑GGGTCGAAGTTAAAG‑GCCAGAAAATGG‑CA‑CCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGGAt                                                                     >  2:519902/1‑228 (MQ=255)                                                                                                                                                            cagcgatatgtgtcatgtgtGCCTTTAAAAATTATTGCAGAGGAGAA‑GGCGTAC‑GCAT‑GGTATAAAGTAAAGC‑GCAAAAAAAGGG‑CA‑CCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCa                                                                             <  2:330720/206‑1 (MQ=255)                                                                                                                                                               tGTT‑T‑CC‑A‑TCA‑TGCCGTTCAAAATTATTGCTGCTGATAACGGCG‑AC‑GCAT‑GGGTCGAAGTTAAAG‑GCCAGAAAATGG‑CA‑CCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGAATCGc                                                                  >  2:115912/1‑228 (MQ=255)                                                                                                                                                                    t‑CC‑A‑TCA‑TGCCGTTCAAAATTATTGCTGCTGATAACGGCG‑AC‑GCAT‑GGGTCGAAGTTAAAG‑GCCAGAAAATGG‑CA‑CCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGAACCGGCATAATTTAACGATGCTCAGCGTAAGGCAACCAAAGAGCCAGGCCGTATCCCTggc                                                              >  1:459763/1‑227 (MQ=255)                                                                                                                                                                      cc‑A‑TCA‑TGCCGTTCAAAATTATTGCTGCTGATAACGGCG‑AC‑GCAT‑GGGTCGAAGTTAAAG‑GCCAGAAAATGG‑CA‑CCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCCACCAAAGAAGCAGGCCGGATCGCTGGTc                                                             >  2:603817/1‑228 (MQ=255)                                                                                                                                                                                gcagTTCAAAATTATTGCTGCTGATAACAGCG‑AC‑GCAT‑GGGTCGAAGTTAAAG‑GCAAGAAAATGG‑CA‑CCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGGGAACGGGTAACTGAACCGGTTATAACCGTACCGCCATACTTTAACGATGCGCGCAGCCAGGCCACCAAAACCACAGGCCGTATCGCTGATCTGaaaaa                                                      >  1:44115/4‑226 (MQ=255)                                                                                                                                                                                gCCGTTCAAAATTATTGCTGCTGATAACGGCG‑AC‑GCAT‑GGGTCGAAGTTAAAG‑GCCAGAAAATGG‑CA‑CCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGAATACTTTAACGTAGCTACGCGTCAGGCAACAAAAGACGCAGGCCGTCTTGTTGGTCGGGaaag                                                      >  1:587861/1‑226 (MQ=255)                                                                                                                                                                                 ccGTTCAAAATTATTGCTGCTGATAACGGCG‑AC‑GCAT‑GGGTCGAAGTTAAAG‑GCCAGAAAATGG‑CA‑CCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGGACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAAACGCAGGCCACATCGCgggcttgacagta                                                     >  1:66104/1‑218 (MQ=255)                                                                                                                                                                                  tgtaaaaaaTTATTGGTGGTGAAAGCGGCG‑GC‑GCGG‑GGGTAGAAGTAAAAG‑GCCA‑AAAATGGACA‑CCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTaa                                                    <  2:2107/223‑1 (MQ=255)                                                                                                                                                                                  ctataaaaaTTATTCGTG‑TGAAAACGCGG‑GC‑GCATGGGGCCGAATTAAAAG‑ACAAAAAAAGGG‑CA‑CCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAAGTAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTaa                                                    <  1:510116/223‑1 (MQ=255)                                                                                                                                                                                    ttCAAAATTATTGCTGCTGATAACGGCG‑AC‑GCAT‑GGGTCGAAGTTAAAG‑GCCAGAAAATGG‑CA‑CCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCCCAGCGTCAGGCAACCAAAACCCCAGGCCGTACGGCGGGTCTGGAAGAaaaa                                                  >  1:515693/1‑228 (MQ=255)                                                                                                                                                                                    ttCAAAATTATTGCTGCTGATAACGGCG‑AC‑GCAT‑GGGTCGAAGTTAAAG‑GCCAGAAAATGG‑CA‑CCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACAGCTGAAGATTACCTGGGAGAACCGGTAACTGAAGCTGTTATCACCGTACCGCCATACTTGAACGATGGCCATCGTCAGCCTACCAAAGACGCAGGCCGCAACGCTGGTCTGGAGGTaaa                                                   >  2:494945/1‑227 (MQ=255)                                                                                                                                                                                         aattatggatggagaaaggggCG‑AA‑GCAT‑GGGGCGAAATTAAAG‑GACAAAAAGTGG‑AA‑CCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTat                                             <  2:251376/209‑1 (MQ=255)                                                                                                                                                                                                   gctgATAACGGCG‑AC‑GCAT‑GGGTCGAAGTTAAAG‑GCCAGAAAATGG‑CA‑CCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAACTTCAGGCAACCAAAGACGCAGGCCGTATCGATGGTCTGGACGTAAAACGAAACATCAACGaa                                   >  2:704930/1‑228 (MQ=255)                                                                                                                                                                                                           cGGCG‑AC‑GTAG‑GGGTCAAAGATAAGG‑GCCAGAAATGAG‑CC‑CCGGCGAAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGc                           <  2:287868/228‑1 (MQ=255)                                                                                                                                                                                                                 aC‑GCAT‑GGGTCGAAGTTAAAG‑GCCAGAAAATGG‑CA‑CCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTCTCACTGGTCCGGAAGTAAAACGTATCATCAACGAACCAACCGccactg                      >  1:143243/1‑223 (MQ=255)                                                                                                                                                                                                                         gggTCGAAGTTAAAG‑GCCAGAAAATGG‑CA‑CCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAACCGCAGGCCGTATCGCTGGTCTGGAAGTAAAAGGAAGACCCAACGAACAGACCGCaagcccggtgg                >  1:459202/1‑218 (MQ=255)                                                                                                                                                                                                                             gaaattaaaaG‑GACAGAAAGTGG‑CA‑CCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTCCCTGGGTGAACCGGTAACTGAAGCTGTTATCGCCGGACCGGCATACTTTAACGATGCTCAGCGGCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCtt             <  2:66104/220‑1 (MQ=255)                                                                                                                                                                                                                                aGTTAAAG‑GCCAGAAAATGG‑CA‑CCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACgg         >  2:688748/1‑228 (MQ=255)                                                                                                                                                                                                                                  ttAAAG‑GCCAAAAAAGGG‑CA‑CCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGt        <  2:459763/227‑1 (MQ=255)                                                                                                                                                                                                                                     aaG‑GCCAGAAAATGG‑CA‑CCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTa                                              >  2:705027/1‑186 (MQ=255)                                                                                                                                                                                                                                         gCCAAAAAATGC‑CACCCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTCCCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGAc  <  1:457403/228‑1 (MQ=255)                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                   GATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTAGTTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTT‑T‑CC‑A‑TCA‑TGCCGTTCAAAATTATTGCTGCTGATAACGGCG‑AC‑GCAT‑GGGTCGAAGTTAAAG‑GCCAGAAAATGG‑CA‑CCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGAC  >  NC_000913/12259‑12708

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