Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,256,301	G→A	S117S (TCC→TCT)	cyuA ←	putative L‑cysteine desulfidase CyuA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,256,301	0	G	A	100.0%	241.4 / NA	73	S117S (TCC→TCT)	cyuA	putative L‑cysteine desulfidase CyuA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (29/44);  total (29/44)

CGCAAACGGGACTTCATTGACGAACTTCAGGATCT‑CAGCCAGCGTCGTTCTGGAA‑AGCACCGTCAGCGGAGACTC‑TTGCTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATG‑TTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCC‑GCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGC  >  NC_000913/3256055‑3256534                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                           cgcAGACGGGACTTCATTGACGAACTTCAGGATCT‑CAGCCAGCGTCGTTCTGGAA‑AGCACCGTCAGCGGAGACTC‑TTGCTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGGCTCGATATTCACAATG‑TTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTGTGGCTCGTGAGAAGAGGATTTCATCGCAAGGGTCCTGGATCTTAACAgaga                                                                                                                                                                                                                                        <  1:356361/251‑1 (MQ=255)       cGGGTCTTCATTGACGAACGTCAGGA‑CTGCAGCCAGCGTCGTTCTGGAA‑AGCACCGTCTGCGGAGACTC‑TGGCGCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATG‑TTGGTATTCCCGCCGACGATGGTGACACACGCCCACTTCTCTCCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCTCAAGGTTCCTGGATTTTAACAGAGACTTTcc                                                                                                                                                                                                                                  <  1:475549/251‑1 (MQ=255)       cGGGAATTCATTGACGAACTTAAGGTTAT‑CAGCCAGCGTCGTTCTGGAA‑AGCACCGTCAGCGGAGACTC‑TTGCTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACCATG‑TTGGTATGCCCGCCGACGCTGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGCGAAGCGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTcc                                                                                                                                                                                                                                  >  2:313868/1‑251 (MQ=255)              tCATTGACGAACTTCAGGATCT‑CAGCCAGCGTCGTTCTGGAA‑AGCACCGTCAGCGGAGACTC‑TTGCTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATG‑TTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCc‑gccgc                                                                                                                                                                                                                          >  2:8539/1‑251 (MQ=255)                 ttGACGAACTTCAGGATCT‑CAGCCAGCGTCGTTCTGGAA‑AGCACCGTCAGCGGAGACTC‑TTGCTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATG‑TTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCC‑gccgcc                                                                                                                                                                                                                         >  1:56577/1‑249 (MQ=255)                   gacaaaCTTCAGGATCT‑CAGCCAGCTTCGTTC‑GGAACACCACCTGCAGCGGAGCCTC‑TTGCCCGCCCCCTGCCACCCACGCCTGCTGGGTAAACACCACACCATCGGGCGTCCCGATATGCCCAATG‑TTGGTATCCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCACCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCC‑GCCGCcagc                                                                                                                                                                                                                      <  2:474568/246‑1 (MQ=255)                   gACGAACTGCTGGATCT‑CAGCCCCCGGCGTTCTGAAA‑AGCACCGCCAGCGGATACTC‑TCGCCCCCCCTCTGCCACCCACGCCTGCTGGTTAAACCCCACACCATCGGGCGTCTCTCTATGCACAATG‑TTGGTATGCCCGCCGACGATGGTGACACACCCCCACTTCTCCCCGTTCCAGACTTTGGCGCGTGAGAAGAGAATTTCATCGCAAGGTTCCGGGATCTTAACAGAGACTTTCCCC‑GCCGCcagca                                                                                                                                                                                                                     <  1:482108/251‑1 (MQ=255)                    aCGAACTTCAGGATCT‑CAGCCAGCGTCGTTCTGGAA‑AGCACCGTCAGCGGAGACTC‑TTGCTCGCCCTCTGCCACACACGCCTGCTGGGTAGACACCACACCATCGTGCGTCTCGATATGCACAATG‑TTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTTCTGGATCTTAACAGAGACTTTCCCC‑GCCGCcagc                                                                                                                                                                                                                      >  1:153704/1‑249 (MQ=255)                           tCAGGATCT‑CAGCCAGCGTCGTTCGGGAA‑AGCACCGTCAGCGGAG‑CTCCTTGCTCGCCCTCTGCCACACACGCCGGCTGGGTAACCACCACACCATCGTGCGTCTCGATATGCACAATG‑TGGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCC‑GCCGCCAGCAGTGCttt                                                                                                                                                                                                              <  2:517236/250‑1 (MQ=255)                             aGGATCT‑CAGCCAGCGTCGTTCTGGAA‑AGCACCGTCAGCGGAGACTC‑TTGCTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCCCCCCACCGTGCGTCTCGATATGCACAATG‑TTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTGAACAGAGACTTTCCCC‑GCCGCCAGCAGTGCTTTGGc                                                                                                                                                                                                           <  2:395717/251‑1 (MQ=255)                              ggtcct‑cAGCCAGCGTAGTTCTGGAA‑AGCACCGACACCGGAGACTC‑TTGCCCGCCCTCCGCCACACTCGCCTGCTGGTTAAACCCCACACCCTCGTGCGGCTCGATACGCACAATG‑TGGGTAGGCCCGCCGACGAGGGTGACACACGCCCACGTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTGCATCGCAAGGTTCCTGGAGCTTAACAGAGACTTTCCCC‑GCCGCCAGCAGTGCTTTGGCa                                                                                                                                                                                                          <  2:294615/247‑1 (MQ=255)                                     cagccagcGTCGTTCTGGAA‑AGCACCGCCAGCGGAGACTC‑TTGCTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATG‑TTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTGGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCGGGATCTTAACAGAGACTTTCCCC‑GCCGCCAGCAGTGCTTTGGCATc                                                                                                                                                                                                        <  1:383543/247‑1 (MQ=255)                                             gtcgtTCTGGAA‑AGCACCGTCAGCGGAGACTC‑TTGCTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATG‑TTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCC‑GCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTg                                                                                                                                                                                            >  1:13919/1‑251 (MQ=255)                                             ggcgtTCGGGAA‑AGCACCGCCAGCGGAGACTC‑TTGCTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCCCCCCATCGTGCGTCTCGATATGCACAATG‑TTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCC‑GCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTg                                                                                                                                                                                            <  2:167068/249‑1 (MQ=255)                                               cgtTCTGGAA‑AGCACCGTCAGCGGAGACTC‑GTGCTCGCCCTCAGCCACACACGCCTGCTGGGTAAACACCACACCATCGGGCGTCTCGATATGCACAAGG‑TTGGTAGGCCCGCCGACGATGGTGACACACGCCCGCTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCC‑GCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTgc                                                                                                                                                                                           <  1:181255/250‑1 (MQ=255)                                                    gggAA‑AGCACCGTAAGCGGAGCCAC‑TTGCTCGCCCTCTGCCACCCACGCCTGCTGGGTAAACACCACTCCACCGTGCGGCTCGATATGCACAATG‑TTGGTATGCCCGCCGACGATGGTGACACAGGCCCACTTCTCCCCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTGCCTGGATCTTAACAGAGGCTTCCCCC‑GCCGCCAGCAGTGCTTTGGCAGCGGCAATTGCCGGCGCTGtt                                                                                                                                                                                     <  2:471168/250‑1 (MQ=255)                                                            cACCGTCAGCGGAGACTC‑TTGCTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATG‑TTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCC‑GCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCtt                                                                                                                                                                              <  1:396257/251‑1 (MQ=255)                                                            cACCGTCAGCGGAGACTC‑TTGCTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATG‑TTGGTAGGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCC‑GCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCtt                                                                                                                                                                              <  1:474691/251‑1 (MQ=255)                                                               ccgccgccggagctctgtGCTCCCCCCCTGCCACCCCCGCCGGCTGGGAAAACCCCACCCCATCGTGCGTCTGGATTGGCA‑AATGGTGGGTATGCCCGCCGACGAGGGTGACACCCGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGTCTTTCCCC‑GCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTc                                                                                                                                                                            <  2:22466/234‑1 (MQ=255)                                                                 ccAGCGGGGACTC‑TTCCTCCCCCTCTGCCCCACACGCCTGCGGGGTAAACCCCCCCCCATCGGGCTCCTCGTAATGCACAATG‑TGGGTATGCCCGCCGACGATGGTGACACACGCCCCCTTCCCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCGGGATCTTAACAGAGACTTTCCCC‑GCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGc                                                                                                                                                                         <  2:207594/250‑1 (MQ=255)                                                                   aGCGGAGACTC‑TTGCTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACCCCATCGGGCGTCTCGATATGCACAATG‑TTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCC‑GCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCAc                                                                                                                                                                       <  1:139631/251‑1 (MQ=255)                                                                     cGGAGACTC‑TTGCTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACCCCATCGTGCGTCTCGATATGCCCAATG‑TTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCC‑GCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACtt                                                                                                                                                                     >  1:366328/1‑251 (MQ=255)                                                                      ggAGACTC‑TTGCTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATG‑TTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCC‑GCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTc                                                                                                                                                                    <  1:265299/251‑1 (MQ=255)                                                                              cgtGCTCGCCCTCTCCCACACCCGCCTGCTGGTTAACCCCCACCCCATCGTGCGTCCAGATATGCACAATG‑GTGTTATGCCCGCCGACGATTGTGACCCCCTCCCACTTCTCACCGTTCCAGCCTTTGGCGCGTGAGAAGAGGATTTCATCGCAATGTTCCTGGATCTTAACAGAGACTTGCCCC‑GCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCTCTGTTGCGTCTTTCAGCACTTCCAGCCCg                                                                                                                                                             <  2:355853/249‑1 (MQ=255)                                                                                ggCTCCCCCTCTGCCACCCACGCCTGCGGGGTAAACCCCACCCCATCGTGCGTCTCGATAGGCACAATG‑GTGGTAGGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCTGACTTTGGCGCGTGAGAAGAGGATTTCACCGCAAGGTTCCTGGATCTAAACAGAGATTTCCCCC‑GCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCTCTGTTGCGTCTGTCAGCACTTCCAGCCCgg                                                                                                                                                            <  2:199958/249‑1 (MQ=255)                                                                                 ccccccccccggCCACACACGCTGGCTGGGAAAACCCCCCCCCATCGGGCGTCTGGATATGCACAATG‑TTGGTATGCCCGCCGACGATGGTGACACCCGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCC‑GCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGc                                                                                                                                                           <  2:54615/239‑1 (MQ=255)                                                                                   ccGCCCTCGACCACGCACGCCTACTGGGTAAGCAACACCCCATCGTTCGTCTCGTTATGCACAATG‑TTGGTATGCCCGCCGACGATGGTAACGCACGCCCACTACTCACCGTACCAGACTTTGGCACGAAGGAAGAGAATTTCAACGCAAGGTTCCAGGATCTTAACAGAGTCTTTCCCC‑GCCGCCAGCAGAGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCGCTTCCAGCCCTGCGtt                                                                                                                                                        <  1:54012/250‑1 (MQ=255)                                                                                         cctGCCACCCACGCCTGCTGGGTAACCACCACACCATCGGGCGTCTCGATATGCACAATG‑TTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCC‑GCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCAtt                                                                                                                                                  <  1:193713/250‑1 (MQ=255)                                                                                          ctGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATG‑TTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCC‑GCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCAttt                                                                                                                                                 >  2:187587/1‑251 (MQ=255)                                                                                             ccACACACGCCGTCTGGATAAACTCCCCACCATCGTGCGTCTCGATTTACACAATT‑TGGGTTGGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCTTGAGAAGAGGATTTCATCGCAAGTTTCCGGGATCTTAACATAGCCTTTCCCC‑GCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACGTCCAGCCCGGCGTTGGCATTTCCa                                                                                                                                              <  2:391987/251‑1 (MQ=255)                                                                                                    cccctgtggggTAAACACCCCACCATCGTGCGTCTCGATATGCACAATG‑TGGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCGGGATCTTAACAGAGACTTTCCCC‑GCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCGCCTAACg                                                                                                                                       <  2:390199/243‑1 (MQ=255)                                                                                                       ctgctgGGTAAACACCCCACCATCGTGCGTCTCGATATGCACAATG‑TTGGTCTGCCCGCCGACGATGGTGACCCACGCCAACTTCTCACCGTTCCAGACTTTGGAGCGTGAGAAGAGGATTTCCTCGCAAGGCTCATGGATCTTACCAGAGCCTTTCCCC‑GCCGCCAGCAGTGCTTTGGAATCGGCAATTGCCTGCGCTGTTGCGTCTTTCCGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGccc                                                                                                                                    >  1:545582/1‑251 (MQ=255)                                                                                                         gctgGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATG‑TTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCC‑GCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGcccc                                                                                                                                   <  1:446293/250‑1 (MQ=255)                                                                                                            gggTAAACACCACCCCATCGTGCGTCTCGATAGGCACAATG‑TTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCC‑GCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCa                                                                                                                                 <  2:173747/249‑1 (MQ=255)                                                                                                                aaaCACGACACCATCGGGCGACTCACTATGCACAATG‑TTGGTTTGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGCGACGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCC‑GCCGCCAGCAGTGCTTTCGCATCGTCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCCTTTCCACCTAACGCCCCCAGcgccg                                                                                                                           >  2:262452/1‑251 (MQ=255)                                                                                                                aaaCACCACACCATCGTGCGTCTCGATATGCACAATG‑TTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCCCCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCC‑GCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGcgcc                                                                                                                            <  1:261173/250‑1 (MQ=255)                                                                                                                 aaCACCACACCCTCGTGCGTCTCGATATGCACAATG‑TTGGTATGCCCGCCGACGAAGGTGACACACGCCCACTTCTCACCGTTCCAGACATTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCC‑GCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTGGGCATTTCCACCTAACGCCCCCAGcgccgc                                                                                                                          >  2:545251/1‑251 (MQ=255)                                                                                                                            aTCGTGCGTCTCGATATGCCCAATG‑TTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCC‑GCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGc                                                                                                                >  2:51583/1‑250 (MQ=255)                                                                                                                             tttccgatcccGATATGCCCAATG‑TTGGTATGCCCGCCGAAGATGTTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTCCCCC‑GCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCa                                                                                                               <  2:451949/240‑1 (MQ=255)                                                                                                                                       cGATATGCACAATG‑TTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCC‑GCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGcgccgccg                                                                                                                        >  1:160222/1‑231 (MQ=255)                                                                                                                                       cGATATGCACAATG‑TTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCC‑GCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCa                                                                                                               >  1:451949/1‑240 (MQ=255)                                                                                                                                       cGATATGCACAATG‑TTGGTAGGCCCGCCGACGATGGTGACCCACGCCCACTTCTCACCGTCCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGCCTTTCCCC‑GCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGcgccgccg                                                                                                                        <  2:160222/231‑1 (MQ=255)                                                                                                                                          aatGCACAAGG‑TTGGTATGCCCGCCGACGATGGGGACACACGCCCACTTCTCACCGTTCCAGACTGTGGCGCGGGAGAAGAGGAGTTCATCGCAAGGTTCCGGGATCTTAACAGAGACTTTCCCC‑GCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTccc                                                                                                 <  1:106545/250‑1 (MQ=255)                                                                                                                                            tGCACAATG‑TTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCC‑GCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGt                                                                                               >  1:236187/1‑251 (MQ=255)                                                                                                                                                    gggtGGTATGCCCGCTGATGATGGGGACCCCCGCCCACTTCTCCCCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCACCGCAAGGTTCCTGGATCTTAACAGAGCCTTCCCCC‑GCCGCCAGCAGTGCGTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCCCCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCAGTCCCGTGCCggg                                                                                         <  2:415406/247‑1 (MQ=255)                                                                                                                                                      ttGGTATGCCCGCCGACGATGGTGACACACGCCCCCTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCC‑GCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGacc                                                                                      >  2:341686/1‑249 (MQ=255)                                                                                                                                                      ttGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCCTCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCC‑GCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCCCTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGacc                                                                                      >  2:438130/1‑249 (MQ=255)                                                                                                                                                                cgccgACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCC‑GCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACgacc                                                                            >  2:59382/1‑248 (MQ=255)                                                                                                                                                                       gaTGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGCGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCC‑GCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCCTTCCCGTGCCGGGCACGGTGACGCCCAGACCGt                                                                     >  1:540331/1‑251 (MQ=255)                                                                                                                                                                       gaTGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCC‑GCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGt                                                                     >  1:254822/1‑251 (MQ=255)                                                                                                                                                                          ggggacccacccccccTTCCCACCTTTCCAGATTTTGGCGCGTGAGAAGAGGATTCCACCGCAAGGTTCCTGGATCTTAACAGAGTCTTCCCCC‑GCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCCCCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGttct                                                                  <  2:5444/236‑1 (MQ=255)                                                                                                                                                                          ggTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCC‑GCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGttct                                                                  >  1:91087/1‑251 (MQ=255)                                                                                                                                                                            tGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCACGGTTCCTGGATCTTAACAGAGACTTTCCCC‑GCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCACTCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGa                                                                         >  2:545699/1‑242 (MQ=255)                                                                                                                                                                            tGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCC‑GCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGa                                                                         <  1:545699/242‑1 (MQ=255)                                                                                                                                                                            tGACACACGCCAACTTCTCACAGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCC‑GCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCGTTTCCACCTATCGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGttcttc                                                                >  1:435319/1‑251 (MQ=255)                                                                                                                                                                             gACACACGCCCCCTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATGTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCC‑GCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTtca                                                               <  1:478412/251‑1 (MQ=255)                                                                                                                                                                                    gCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCC‑GCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGAtt                                                        <  1:8539/251‑1 (MQ=255)                                                                                                                                                                                      ccACTTCTCACCGGTCCAGACTTTGGCTCGTGAGAAGAGGATTTCAGCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCC‑GCCGCCAGCAGTGCTTTTGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCCCCTAACGCCCCCAGCTCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTg                                                      <  1:436261/251‑1 (MQ=255)                                                                                                                                                                                            ctcACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCC‑GCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGaa                                                 >  2:157848/1‑250 (MQ=255)                                                                                                                                                                                                     ccAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCC‑GCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCtt                                       >  1:27543/1‑251 (MQ=255)                                                                                                                                                                                                      cAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCC‑GCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGAAGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTc                                      >  1:128529/1‑251 (MQ=255)                                                                                                                                                                                                            ttGGCGGGTGAGAAGAGGTTTTCATCGCAAGGTTTCTGGATCTTAACAGAG‑CTTTCCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCTCCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTaca                                  <  2:153704/249‑1 (MQ=255)                                                                                                                                                                                                            tgggCGCGTGAGAAAAGGATGTCATCGCAAGGGTCCTGGATCTAAACAGAGACTTTCCCC‑GCCGCCAGCAGTGCTTGGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGGTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGGGCCGGGAACGGTGAGGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTaca                                  <  1:217073/247‑1 (MQ=255)                                                                                                                                                                                                              ggctggtgggaagaGGAGTTCATCGCAGGGTTCCGGAATCTTAACAGAGACTTCCCCC‑GCCGCCAGCAGTGCGTTGGCACCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGTCATGTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCAGTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGTCGAAACCCAGGCTTCTACACGtt                              <  1:313868/242‑1 (MQ=255)                                                                                                                                                                                                                   gTGAGAAGAGGTTTTCATCGCAAGGTTCCTGGATTTTACAAGAGACTTTCCCC‑GCCGCCAGCAGTTCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCTAAACCCAGGCTTCTACACGTTCAAcc                         <  2:250690/251‑1 (MQ=255)                                                                                                                                                                                                                          gagGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCC‑GCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAgc                                                                                                                               >  2:460473/1‑142 (MQ=255)                                                                                                                                                                                                                          gagGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCC‑GCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAgc                                                                                                                               <  1:460473/142‑1 (MQ=255)                                                                                                                                                                                                                                ttCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCC‑GCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAgc                                                                                                                               >  2:152321/1‑136 (MQ=255)                                                                                                                                                                                                                                ttCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCC‑GCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAgc                                                                                                                               <  1:152321/136‑1 (MQ=255)                                                                                                                                                                                                                                     cGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCC‑GCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGAACGTTCTTCATCAGCTTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGcc       >  1:94213/1‑251 (MQ=255)                                                                                                                                                                                                                                      caaaGGTTCCTGAATCTAACCAGAGACTTTCCCC‑GCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCg      <  2:75608/249‑1 (MQ=255)                                                                                                                                                                                                                                          ggTTCCGGGATCTAAACAGAGATTTGCCCC‑GCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCGTTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAAGCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATGTGGCGAAACCCAGGCTTCTTCACGTTCAACCGGACCGTCCATTGCTGCGgcagc  <  2:526662/251‑1 (MQ=255)                                                                                                                                                                                                                                                      ttAACAGAGACTTTCCCC‑GCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCa                                            >  1:418232/1‑197 (MQ=255)                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                           CGCAAACGGGACTTCATTGACGAACTTCAGGATCT‑CAGCCAGCGTCGTTCTGGAA‑AGCACCGTCAGCGGAGACTC‑TTGCTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATG‑TTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCC‑GCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGC  >  NC_000913/3256055‑3256534

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