Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,364,402	G→A	17.0%	R533K (AGA→AAA)	yhcD →	putative outer membrane fimbrial subunit usher protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,364,402	0	G	A	17.0%	84.7 / 3.1	45	R533K (AGA→AAA)	yhcD	putative outer membrane fimbrial subunit usher protein
Reads supporting (aligned to +/- strand):  ref base G (32/5);  new base A (3/5);  total (35/10)
Fisher's exact test for biased strand distribution p-value = 8.25e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.59e-01

GTCTAACTCCATTCTGGGTGGTTTCAACATCAGTGGCTATAAAACATATAGTTACGACGGTGAC‑AATGACTCAAGCAACCTTAATATTAACTGG‑AACAAAGCGTTCAAA‑CACGCCACGGTTTCCGTCAACTGGCAGCATCAAC‑TTAGTG‑CTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTC‑AACAT‑CAGTATTCCATTTGGC‑AGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTAC  >  NC_000913/3364190‑3364622                                                                                                                                                                                                                             |                                                                                                                                                                                                                             gTCTAACTCCATTCTGGGTGGTTTCAACATCAGTGGCTATAAAACATATAGTTACGACGGTGAC‑AATGACTCAAGCAACCTTAATATTAACTGG‑AACAAAGCGTTCAAA‑AACGCCACGGTTTCCGTCAACTGGCAGCATCAAC‑TTAGTG‑CTTCAGAAAATAATGACGACGATGGTGAGCTGTTCTACCAC‑AAACT‑CAGTATTCCATTTGGC‑AAAACAAAAACCGCcac                                                                                                                                                                                                               >  1:328480/1‑228 (MQ=255)  tCTAACTCCATTCTGGGTGGTTTCAACATCAGTGGCTATAAATCATATAGTTACGACGGTGAC‑AATGACTCAAGCAACCTTAATATTAACTGG‑AACAAAGCGTTCAAA‑CACGCCAAGGTTTCCGTTAACTTGCAACATCAAC‑ATCGTG‑CTTCAGAACAAAATTAATAAGATGGACAACTGTTTTACGTC‑AAATT‑CAGCATTCATTTTGGC‑AAATAAAACAAGGCcaa                                                                                                                                                                                                               >  1:192720/1‑226 (MQ=255)  tCTAACTCCATTCTGGGTGGTTTCAACATCAGTGGCTATAAAACATATAGTTACGACGGTGAC‑AATGACTCAAGCAACCTTAATATTAACTGG‑AACAAAGCGTTCAAA‑CACGCCACGGTTTCCGTCAACTGGCAGCATCAAC‑TTAGTG‑CTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTC‑AACAT‑CAGTATTCCATTTGGC‑AGATCAAACACCGCcac                                                                                                                                                                                                               >  2:556666/1‑227 (MQ=255)  tCTAACTCCATTCTGGGTGGTTTCAACATCAGTGGCTATAAAACATATAGTTACGACGGTGAC‑AATGACTCAAGCAACCTTAAAATTAACTGG‑AACAAAGCGTTCAAA‑CACGCCACGGTTTCCGTCAACTGGCAGCATCAAC‑TTAGTG‑CTTCAGCAAATAATGAAGACGAGGGTGATCTGTTCTACGTC‑AACAT‑CAGAATTCCATTTGGC‑AGATCAAAAACCGCcac                                                                                                                                                                                                               >  2:643129/1‑227 (MQ=255)      aCTCCATTCTGGGTGGTTTCAACATCAGTGGCTATAAAACATATAGTTACGACGGTGAC‑AATGACTCAAGCAACCTTAATATTAACTGG‑AACAAAGCGTTCAAA‑CACGCCACGGTTTCCGTCAACTGGCAGCATCAAC‑TTAGTG‑CTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTC‑AACAT‑CAGAATTCCATTTGGC‑AGATCAAACACCGCCACACtca                                                                                                                                                                                                          >  1:198617/1‑226 (MQ=255)           aTTCTGGGTGGTTTCAACATCAGTGGCTATAAAACATATAGTTACGACGGTGAC‑AATGACTCAAGCAACCTTAATATTAACTGG‑AACAAAGCGTTCAAA‑CACGCCACGGTTTCCGTCAACTGGCAGCATCAAC‑TTATAG‑CTTCAGAAAATAATTAAGACGATGCTGATCGGTTCTACGTC‑AACAT‑CAGAATTCCATTTGTC‑AGATCAAAAACCCACACCCTgaaaaa                                                                                                                                                                                                      >  1:161122/1‑222 (MQ=255)           aTTCTGGGTGGTTTCAACATCAGTGGCTATAAAACATATAGTTACGACGGTGAC‑AATGACTCAAGCAACCTTAATATTAACTGG‑AACAAAGCGTTCAAA‑CACGCCACGGTTTCCGTCAACTGGCAGCATCAAC‑TTAGTG‑CTTCAGAAAATAACGAAGACGATGGGGATCTGTTCTACGTC‑AACAT‑AAGTATTCAATTTGGC‑AGCACAAAAACCGCCACACCGTAAACt                                                                                                                                                                                                     >  1:37218/1‑228 (MQ=255)                 ggtggtTTCAACATCAGTGGCTATAAAACATATAGTTACGACGGTGAC‑AATGACTCAAGCAACCTTAATATTAACTGG‑AACAAAGCGTTCAAA‑CACGCCCCGGTTTCCGTCAACTGGCAGCATCAAC‑TTAGTG‑CTTCAGAAAATAATGAAGACGATGGTGATATGTTCTACGTC‑AACAT‑CAATATTCCATTTGGC‑AGATAAAACACCGCCTCACGGaataacagccca                                                                                                                                                                                               >  1:587026/1‑220 (MQ=255)                 ggtggtTTCAACATCAGTGGCTATAAAACATATAGTTACGACGGTGAC‑AATGACTCAAGCAACCTTAATATTAACTGG‑AACAAAGCGTTCAAA‑CACGCCACGGTTTCCGTCAACTGGCAGCATCAAC‑TTAGTG‑CTTCAGAAAATAATGAAGACGATGGGTATCTGTTCTACGTC‑AACAT‑CAGTATTCCATTTGGC‑AGATCAAAAACCCCCACACTGTAAAcccacaaa                                                                                                                                                                                               >  1:681770/1‑221 (MQ=255)                  gtggtTTCAACATCAGGGGCTATAAAACATATAGTTACGACGGTGAC‑AATGACTCCAGCAACCTTAATATTAACTGG‑AACAAAGCGTTCAAA‑CACGCCACGGTTTCCGTCAACTGGCAGCATCAAC‑TTAGTG‑CTTCAGAACATAATGCAGACGATGGTTATCTGTTCTACGTC‑AAACT‑CAGGATTCCATTTTGC‑AGATAAGACACGGCCACACTGAATACCCGCCaca                                                                                                                                                                                              >  1:88281/1‑226 (MQ=255)                      tttCAACATCAGTGGCTATAAAACATATAGTTACGACGGTGAC‑AATGACTCAAGCAACCTTAATATTAACTGG‑AACAAAGCGTTCAAA‑CACGCCACGGTTTCCGTCAACTGGCAGCATCAAC‑TTATTG‑CTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTC‑AAAAT‑CAGCATTCAATTTGGC‑AGATCAAACACCGCCAACCTGAATACTCGCCATGACGa                                                                                                                                                                                          >  2:312931/1‑228 (MQ=255)                      tttCAACATCAGTGGCTATAAAACATATAGTTACGACGGGGAC‑AATGACTCAAGCAACCTTAATATTAACTGG‑AACAAAGCGTTCAAA‑CACGCCACGGTTTCCGTCAACTGGCAGCATCAAC‑TTAGTG‑CTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTC‑AACAT‑CAGTATTCCATTTGGC‑AGATCAAACACCGCCACACTGTATACTCGCCATGACg                                                                                                                                                                                           >  2:487397/1‑227 (MQ=255)                       ttCAACATCAGTGGCTATAAAACATATAGTTACGACGGTGAC‑AATGACTCAAGCAACCTTAATATTAACTGG‑AACAAAGCGTTCAAA‑CACGCCACGGTTTCCGTCAACTGGCAGCATCAAC‑TTAGTA‑CTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTC‑AACAT‑CAGTATTCCATTTGGC‑AGATCAAACACCGCCACACTGTATACTCGCCaatacaat                                                                                                                                                                                         >  1:637123/1‑221 (MQ=255)                           aCATCAGTGGCTATAAAACATATAGTTACGACGGTGAC‑AATGACTCAAGCAACCTTAATATTAACTGG‑AACAAAGCGTTCAAA‑CACGCCACGGTTTCCGTCAACTGGCAGCATCAAC‑TTAGTG‑CTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTC‑AACAT‑CAGTATTCCATTTGGC‑AGATCAAACACCGCCACACTGTATCCTCGCCATGAAGATCATa                                                                                                                                                                                     >  2:410940/1‑228 (MQ=255)                                     taaaaaaacattaatttacaggggTAAC‑AATGACCCAACCAACTTTAATATTAACTGA‑AAAAAAGCGTTCAAA‑CACCCCAGGTTTTCCGTAACCTGGAACAACCACT‑TTATTT‑CTTAAAAAAATAAGAAAGACGATGGTATTCTTTTCTACGTC‑AACAT‑CAGTTTTCCTTTTGGA‑AAATAAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTa                                                                                                                                                                           <  1:556666/206‑1 (MQ=255)                                          aaaCATAGAGTTCCAACGGTGAC‑AGTAGATCAAGCAACCATATTATT‑ACTGA‑GAAAAAGTCTTAAAACCACGACATGGTTTCGGCAATCTGCTAGCTACACT‑TTAGTG‑ATCAAAAACAAAAAGAAGCCGATATTGCTCTTTTCTACGCA‑AACAC‑CAGTATACCTTTTGGA‑AGAGCAACCACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTa                                                                                                                                                                      <  2:192720/228‑1 (MQ=255)                                           aaccataaattaacgaggTGACAAATGACGCAAGCACCTTTAAT‑TTAACTGGAAACAAAGCGTTAAAC‑CACGCCAGGGTTTCCGTCAATGGGCAG‑ATCCACTTTAGTG‑CTTCAGAAAATAAGGAAGAGTATGGTAATCTGTTTTACGTC‑AACAT‑CATTATTCTATTTGGC‑AGATCAACCACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTa                                                                                                                                                                      <  2:82173/212‑1 (MQ=255)                                                agagtcacggcgtgaaC‑AAGGACCCAAGAAACTTTAATATAACAGGG‑AAAAAGGGGTGCAAA‑CCCGGCATGGTTCCCGTAAATTGGCAGCATCAAT‑TTAGGG‑CTTCAAAAAATAAGAAGGACGAGGAGGATCGGTTCTACGCA‑AACAT‑AAGTATTCCTTTTGGC‑AAATCAAAAACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGgt                                                                                                                                                                 <  2:161122/212‑1 (MQ=255)                                                       cgacgGTGAC‑AATGACTCAAGCAACCTTAATATTAACTGG‑AACAAAGCGTTCAAA‑CACGCCACGGTTTCCGTCAACTGGCAGCATCAAC‑TTAGTG‑CTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTC‑AACAT‑CAGTATTCCATTTGGC‑AGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGAACTGGGGTCATggg                                                                                                                                                         >  2:581331/1‑228 (MQ=255)                                                       cgacgGTGAC‑AATGACTCAAGCAACCTTAATATTAACTGG‑AACAAAGCGTTCAAA‑CACGCCACGGTTTCCGTCAACTGGCAGCATCAAC‑TTAGTG‑CTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTC‑AACAT‑CAGAATTCACTTTGGC‑AGATCAAACCCCGCCCCCCTGTATACTCGCCATACCGATCATAAAACCCCCCAAGGACCTGGTGTTATgga                                                                                                                                                         >  1:611443/1‑227 (MQ=255)                                                              taaaaataccaccagacacattAATAGTAACAGG‑AAAAAGACGTCAAAA‑CCAGCCACGGTTAGCGTAAACTGGCAGCATCAAC‑TTGGTG‑CTTAAAAAAATAATGAAGACGATGGTGCACTGTTCCACTCC‑AACAT‑CAGTATCCCATTTGGA‑AAAACAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAgtcgt                                                                                                                                                   <  1:568684/208‑1 (MQ=255)                                                              gAC‑AATGACTCAAGCAACCTTAATATTAACTGG‑AACAAAGCGTTCAAA‑CACGCCACGGTTTCCGTCAACTGGCAGCATCAAC‑TTAGTG‑CTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTC‑AACAT‑CAGTATTCCATTTGGC‑AGATCAAACACCGcc                                                                                                                                                                                                                 >  2:158255/1‑165 (MQ=255)                                                                         aaaGCAATCTAAAATTTACCCGG‑AACAAAGCGTAAAAA‑CCCGCCCTGTTTTCCGACAACTGCGAACACCAAC‑‑TAGCGCCTAAAGAAAATAATGAAGACGATGGTGATCTTTCCTACGTC‑‑ACATACAGTTTTCCTTTTGGC‑AAATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCagatgag                                                                                                                                         <  2:62083/226‑1 (MQ=255)                                                                           aGCAACCTTAATATTAACTGG‑AACAAAGCGTTCAAA‑CACGCCACGGTTTCCGTCAACTGGCAGCATCAAC‑TTAGTG‑CTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTC‑AACAT‑CAGTATTCCATTTGGC‑AGATCAAACACCGCCACACTCTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGTTGTCATGGGAGTCGTCTCagataaaagg                                                                                                                                      >  1:31256/1‑222 (MQ=255)                                                                              aaCTTTAATATAAACGGA‑AAAAAACCTTTCAAA‑CACGCCAGTTTCTCAGTAAACTGGCAGCATCAGC‑TTAGGT‑CTTAAGAAAAATATGAAAGCGTTGGTACTGTTTTCTACGCA‑AACAT‑CAGTATTCCATTTGGA‑AAATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTcc                                                                                                                                   <  2:328480/228‑1 (MQ=255)                                                                                          aaCTGG‑AACAAAGCGTTCAAA‑CACGCCACGGTTTCCGTCAACTGGCAGCATCAAC‑TTAGTG‑CTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTC‑AACAT‑CAGTATTCCATTTGGC‑AGATCAAACACCGCCACACTGTATACACGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGGTGCGTTTCATATGAGATGTCCTACTATGGGaag                                                                                                                       >  1:498274/1‑227 (MQ=255)                                                                                           aCTGG‑AACAAAGCGTTCAAA‑CACGCCACGGTTTCCGTCAACTGGCAGCATCAAC‑TTAGTG‑CTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTC‑AACAT‑CAGTATTCCATTTGGC‑AGATc                                                                                                                                                                                                                           >  1:72793/1‑127 (MQ=255)                                                                                           aCTGG‑AACAAAGCGTTCAAA‑CACGCCACGGTTTCCGTCAACTGGCAGCATCAAC‑TTAGTG‑CTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTC‑AACAT‑CAGTATTCCATTTGGC‑AGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAAt                                                                                                                       >  2:497882/1‑227 (MQ=255)                                                                                           aCTGG‑AACAAAGCGTTCAAA‑CACGCCACGGTTTCCGTCAACTGGCAGCATCAAC‑TTAGTG‑CTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTC‑AACAT‑CAGTATTCCATTTGGC‑AGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGGGTCATGGGAg                                                                                                                                                       >  2:208249/1‑195 (MQ=255)                                                                                           aCTGG‑AACAAAGCGTTCAAA‑CACGCCACGGTTTCCGTCAACTGGCAGCATCAAC‑TTAGTG‑CTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTC‑AACAT‑CAGTATTCCATTTGGC‑AGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACAATGCACCTGGTGCCATGGGAGTCGCCCCACATGAAAAGTCCAACAATTTAAAta                                                                                                                      >  1:669769/1‑227 (MQ=255)                                                                                           aCTGG‑AACAAAGCGTTCAAA‑CACGCCACGGGTTCCGTCAACTGGCAGCATCAAC‑TTAGTG‑CTTCAGAAAATAATGAAGACGATGAAGATCTGTTCTACGTC‑AACAT‑CAGTATTCCATTTGGC‑AAATCAAACACCGGCACCCTGTATACTCGTCATGACGATCATAAAACCCACTATGGAACTGGTGACATTGGAGTCCTCTAagatgcatagtcctcataagtaaaag                                                                                                                      >  1:54846/1‑207 (MQ=255)                                                                                                 aaCAAAGCGTTCAAA‑CACGCCACGGTTTCCGTCAACTGGCAGCATCAAC‑TTAGTG‑CTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTC‑AACAT‑CAGTATTCCATTTGGC‑AGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCagatgagacggtctcttatacccatctccg                                                                                                                  >  2:259054/1‑205 (MQ=255)                                                                                                 aaCAAAGCGTTCAAA‑CACGCCACGGTTTCCGTCAACTGGCAGCATCAAC‑TTAGTG‑CTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTC‑AACAT‑CAGAATTCCATTTGGC‑AGATCAAACACCGCCAAACTGTATACTCGCCATGACGATCATAAAACCCACTATGGCACTGGTGTCAGGGGCGTCGTCTAAGAGGAAATACCCTCCTATGTGACTGCTGa                                                                                                                  >  1:131697/1‑227 (MQ=255)                                                                                                          ttCAAA‑CACGCCACGGTTTCCGTCAACTGGCAGCATGAAA‑TTAGTG‑CTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTC‑AACAT‑CAGTATTACATTTGGC‑AGATCAAACACCGCCACACTGTATACACGTCATGACGACCATAAAACCCACAATGGTACTGGAGTACTGGGGGACGTGTCAGATGAGATGGCCTACTAtggaaaagcggtttcttacaca                                                                                                        >  1:468369/1‑208 (MQ=255)                                                                                                          ttCAAA‑CACGCCACGGTTTCCGTCAACTGGCAGCATCAAC‑TTAGTG‑CTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTC‑AACAT‑CAGTATTCCATTTGGC‑AGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCagatga                                                                                                                                          >  2:447534/1‑194 (MQ=255)                                                                                                          ttCAAA‑CACGCCACGGTTTCCGTCAACTGGCAGCATCAAC‑TTAGTG‑CTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTC‑AACAT‑CAGTATTCCATTTGGC‑AGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACCGGTGTCATGGGAGTCGTCTCAGATGAGCTGTACTACTATCTGAATGCTCACCGCGATCat                                                                                                        >  2:314579/1‑227 (MQ=255)                                                                                                           tCAAA‑CACGCCACGGTTTCCGTCAACTGGCAGCATCAAC‑TTAGTG‑CTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTC‑AACAT‑CAGTATTCCATTTGGC‑AGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTTTCAGATGAGAGGTCCTACTATGTGAATGCTGAACGAGATCacg                                                                                                       >  2:205266/1‑228 (MQ=255)                                                                                                               aaccgcccAGGGTTTCCGGCACCGGGAAGAATACAG‑TTATGG‑CGTAAGAAAATAAAGAAGACGAGGGGGACCTG‑TCTTAGTCAAAAAC‑CAGTTTGCCTTTTGCC‑AGATAAACCTCCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCacgacg                                                                                                    <  2:37218/222‑1 (MQ=255)                                                                                                                                           ccaacAAT‑TTAGTC‑CTTAAAAAAATAATAAAGCGGTTGGTGACCTTTTCACCGAC‑AAGAC‑CA‑TATTCCTTTTTGCAAGACAAAACCCCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTAt                                                                        <  2:418302/224‑1 (MQ=255)                                                                                                                                                        tG‑CTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTC‑AACAT‑CAGTATTCCATTTGGC‑AGTTCAAAAACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACAATGGACCTGGTGTCCCTGGCGTCGACTCACAGAAGATGTCCTACTCTGTGAATGCTGACGCGGAGCACAACGAACGAGAAACAACCTTTAAGGCCAGCAGACGTTCCaaccc                                                            >  1:331955/1‑225 (MQ=255)                                                                                                                                                              cAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTC‑AACAT‑CAGTATTCCATTTGGC‑AGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTgg                                                                                                                                                                  >  2:91720/1‑121 (MQ=255)                                                                                                                                                              cAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTC‑AACAT‑CAGTATTCCATTTGGC‑AGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGGGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGGGAAACGAGCTTGACAGGCAGTATCAGTTCCacttccca                                                         >  1:518373/1‑219 (MQ=255)                                                                                                                                                               aGAAAAAAAGAAAGACGAGGGTGATCTGTTCTACGAC‑AACAT‑CAGTTTGCCATTTGGC‑AGACAAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGGGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAAt                                                              <  2:518373/220‑1 (MQ=255)                                                                                                                                                                             tcGATGGGGATCTGCTCTACGCC‑AACAT‑CAGCATTCCATTTGGC‑AGATCTAACACCGCCACACTGTATACTCGCCATGACGATAATAAAACCCACTATGGTACTGGTGTCATCGGACTCGTCCCCGATTATATTCCCTACTCTGTGAATGCTGTACCAGTTCACGACGAGCGCGCAACGAGCTTGACCGGCAGTATCCGTTCCAATACACATTTCACCCaggcaaaaa                                        >  1:331772/2‑220 (MQ=255)                                                                                                                                                                                                                     tttGGC‑AGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGGGAAACGAGCTTGAACGGACGTATCAGTTCCAATCTCCATTACACCCAAGCCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAAccaaccttac  >  1:553019/1‑220 (MQ=255)                                                                                                                                                                                                                             |                                                                                                                                                                                                                             GTCTAACTCCATTCTGGGTGGTTTCAACATCAGTGGCTATAAAACATATAGTTACGACGGTGAC‑AATGACTCAAGCAACCTTAATATTAACTGG‑AACAAAGCGTTCAAA‑CACGCCACGGTTTCCGTCAACTGGCAGCATCAAC‑TTAGTG‑CTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTC‑AACAT‑CAGTATTCCATTTGGC‑AGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTAC  >  NC_000913/3364190‑3364622

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