Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,637,426	Δ1 bp	intergenic (‑16/‑216)	yhiN ← / → pitA	putative oxidoreductase YhiN/metal phosphate:H(+) symporter PitA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,637,426	0	G	.	100.0%	373.1 / NA	77	intergenic (‑16/‑216)	yhiN/pitA	putative oxidoreductase YhiN/metal phosphate:H(+) symporter PitA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (32/45);  total (32/45)

GAGTGCAGACTTACAAAAATGCGGATTCTGGCTCAGATAAGCGCCTGGTTCGACATAAAGGTTGGTAAAGTTGCAGCGCCCACCGCCAGACATAAGGATTTTGCGCCCTGGTTTTTTACCATTATCGATCAGCAGAACCCGGCGTCCTGCCTGACCTGCCAGCGCAGAACAGAACATACCCGCCGCACCAGCGCCTATAATAATGGCATCAAACCTTTCCACGTTGCGCTCCTCTTAGAAAAAACGGGCGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTACATTTGTTTGCTGGCCTGGAT  >  NC_000913/3637189‑3637668                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                   atgtgcacatcaaccacgacgactgtctggaaatcgccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTg                                                                                                                                                                                                                                      >  1:809135‑M2/223‑251 (MQ=255) atgtgcacatcaaccacgacgactgtctggaaatcgccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTg                                                                                                                                                                                                                                      >  1:217860‑M2/223‑251 (MQ=255)     gcacatcaaccacgacgactgtctggaaatcgccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggacgattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctacataccctatccgaGTTGCGCTCCTCTTA‑AAACAACGGACGTGAAtt                                                                                                                                                                                                                                  >  1:333082‑M2/219‑251 (MQ=255)       acatcaaccacgacgactgtctggaaatcgccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGt                                                                                                                                                                                                                                <  2:294419‑M2/35‑1 (MQ=255)         atcaaccacgacgactgtctggaaagcgccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtccttgggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTaa                                                                                                                                                                                                                              <  2:319411‑M2/37‑1 (MQ=255)          cacaccccgcagtctgtcttgaaattgccgtgttgaaaggtgacatgggtgacgtgaagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggtatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacttgattactctctaaataccctatccgaGTTGCGCTCCTTTTA‑AAAAAACGGTCGGGAATTGTaaa                                                                                                                                                                                                                             <  2:214807‑M2/38‑1 (MQ=255)                  gacgactgtctggaaatcgccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTc                                                                                                                                                                                                                     >  1:839515‑M2/206‑251 (MQ=255)                       ctgtctggaaatcgccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTgg                                                                                                                                                                                                                >  1:160981‑M2/201‑251 (MQ=255)                         gtctggaaatcgccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGACTCCTCAGTGGt                                                                                                                                                                                                               <  1:123540‑M2/52‑1 (MQ=255)                          tctggaaatcgccgtgttgaccggtgacatgggtgacgtgccgcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgcGTTGCGCTCCTCTTA‑AACAAACGGACGTGAATTGAAAAGATTCCTCAGTGGtcc                                                                                                                                                                                                             >  2:439679‑M2/198‑250 (MQ=255)                            tggaaatcgccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGc                                                                                                                                                                                                            >  1:562854‑M2/196‑250 (MQ=255)                               aaatcgccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCAcc                                                                                                                                                                                                         >  2:771741‑M2/193‑250 (MQ=255)                                aatcgccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAg                                                                                                                                                                                                       <  2:562854‑M2/60‑1 (MQ=255)                                 atcgccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGc                                                                                                                                                                                                      >  2:451158‑M2/191‑251 (MQ=255)                                 atcgccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGc                                                                                                                                                                                                      <  1:534938‑M2/61‑1 (MQ=255)                                   cgccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgccgttatcgcacagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAAGTGTAAAGATTCCTCAGTGGTCGCACAAGCAt                                                                                                                                                                                                    <  2:572253‑M2/63‑1 (MQ=255)                                       gtgttgaaaggtgacatgggtgtcgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCGCAGTGGTCGCACCAGCATCAat                                                                                                                                                                                                <  2:733987‑M2/67‑1 (MQ=255)                                          ttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATAtt                                                                                                                                                                                             >  1:537433‑M2/182‑251 (MQ=255)                                                ggtgacatgggtgacgtgcagcattttgccgaggacgttatcgcccagcgcggcgtgcggcacgggtatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctagccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTTAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTaaa                                                                                                                                                                                       <  2:689716‑M2/76‑1 (MQ=255)                                                        gggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTa                                                                                                                                                                               <  1:289085‑M2/84‑1 (MQ=255)                                                         ggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTAt                                                                                                                                                                              <  2:537433‑M2/85‑1 (MQ=255)                                                         ggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTAt                                                                                                                                                                              >  1:450519‑M2/167‑251 (MQ=255)                                                            gacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCGCTTA‑AAAAAACGGACGTTAATTGAAAAGACTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTg                                                                                                                                                                           <  2:189766‑M2/88‑1 (MQ=255)                                                               gtgcagcatgttgccgatgacgttatcgcccagcgaggcgtgcggaccgggcatttgcagtgcttgccgaaggaagattgagtcgattattgggtcctttggtgaactatcaacaaagggccactagcacacctgattcctctctaaatacccgatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAAGTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAGAAGGAAGTATTTGCCt                                                                                                                                                                        <  2:37646‑M2/91‑1 (MQ=255)                                                                 gcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGa                                                                                                                                                                      <  2:774213‑M2/93‑1 (MQ=255)                                                                   agcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAAt                                                                                                                                                                    <  1:889056‑M2/95‑1 (MQ=255)                                                                      attttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATtat                                                                                                                                                                 <  1:738884‑M2/98‑1 (MQ=255)                                                                        tttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATtata                                                                                                                                                                <  1:489705‑M2/99‑1 (MQ=255)                                                                         ttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAAGATTACTAAAAGGAAGTATTTGCCTGAATtatat                                                                                                                                                               <  2:37612‑M2/100‑1 (MQ=255)                                                                           gccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaatcccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAACGGAAGTATTTGCCTGAATtatata                                                                                                                                                              >  2:461935‑M2/149‑249 (MQ=255)                                                                                 gacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggcctctagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTTAATTATATAAGATAatt                                                                                                                                                      <  2:78438‑M2/109‑1 (MQ=255)                                                                                 gacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAatt                                                                                                                                                      <  1:677230‑M2/109‑1 (MQ=255)                                                                                 gacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAatt                                                                                                                                                      <  2:473204‑M2/109‑1 (MQ=255)                                                                                 gacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAatt                                                                                                                                                      >  1:78404‑M2/143‑251 (MQ=255)                                                                                 gacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgatgcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAatt                                                                                                                                                      <  2:887572‑M2/109‑1 (MQ=255)                                                                                   cgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCCATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAattat                                                                                                                                                    >  2:412194‑M2/141‑251 (MQ=255)                                                                                   cgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctagtattgggtcctttggtgaacgatcaacaaagggctactagcacacctgatttctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATGCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGgaagtatttgcctgaatgagagaagagaaggat                                                                                                                                                    <  1:276848‑M2/111‑17 (MQ=255)                                                                                       atcgcccagcgcggcgtgcggcacggtcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAATATTCCTCAGTTGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTAttttt                                                                                                                                                <  1:48489‑M2/115‑1 (MQ=255)                                                                                         cgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTg                                                                                                                                              <  1:615570‑M2/117‑1 (MQ=255)                                                                                         cgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCCCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTg                                                                                                                                              <  2:83824‑M2/117‑1 (MQ=255)                                                                                          gcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGa                                                                                                                                             <  2:171981‑M2/118‑1 (MQ=255)                                                                                                 gcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAAt                                                                                                                                      <  2:246511‑M2/125‑1 (MQ=255)                                                                                                 gcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAAGTATTTTTTGAGTGAAAt                                                                                                                                      <  1:199469‑M2/125‑1 (MQ=255)                                                                                                   ggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATcc                                                                                                                                    <  1:245218‑M2/127‑1 (MQ=255)                                                                                                    gcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCa                                                                                                                                   <  2:217860‑M2/128‑1 (MQ=255)                                                                                                      gtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATa                                                                                                                                 <  1:439679‑M2/130‑1 (MQ=255)                                                                                                        gcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaatgccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATAc                                                                                                                                <  1:221825‑M2/131‑1 (MQ=255)                                                                                                        gcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATAc                                                                                                                                >  2:653778‑M2/120‑250 (MQ=255)                                                                                                        gcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATAc                                                                                                                                >  2:800071‑M2/120‑250 (MQ=255)                                                                                                                   atttgcagtgcttgccgaaggaagcttgagtctattattgggtcctttgatgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑CAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATc                                                                                                                    >  1:300308‑M2/109‑251 (MQ=255)                                                                                                                   atttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATc                                                                                                                    <  1:436064‑M2/143‑1 (MQ=255)                                                                                                                   atttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATc                                                                                                                    <  2:145796‑M2/143‑1 (MQ=255)                                                                                                                                gccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCtat                                                                                                       >  1:722948‑M2/96‑251 (MQ=255)                                                                                                                                  cgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTatat                                                                                                     <  2:429950‑M2/158‑1 (MQ=255)                                                                                                                                   gaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATAtt                                                                                                    <  1:451158‑M2/159‑1 (MQ=255)                                                                                                                                       gaagattaagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATAAAAAAAAGTCTATATTTCAc                                                                                                >  2:443908‑M2/89‑251 (MQ=255)                                                                                                                                        aagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACt                                                                                               >  1:58368‑M2/88‑251 (MQ=255)                                                                                                                                                gtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCccgc                                                                                       >  1:566177‑M2/80‑251 (MQ=255)                                                                                                                                                gtctattattgggtcctttggtgaacgatcaacaaagggccactagcacaccggattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGAAAATCAAAAAAAGTCTATATTTCACTTTGCccgc                                                                                       <  1:213399‑M2/172‑1 (MQ=255)                                                                                                                                                           ggtccgttggtgaacgatcaacaaagggcccctagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGACTTGTAAAGATTCCTCAGTGGTCGCACCAGCATCACTATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGt                                                                            >  2:580249‑M2/69‑251 (MQ=255)                                                                                                                                                            gtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTc                                                                           <  1:321962‑M2/184‑1 (MQ=255)                                                                                                                                                                   ggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCATCTTA‑AAAAAACGGCCGTGAATTGTAAAGATTCATCAGTGGTCGCACCAGCATCAATATTACTACAAGGACGTATTTGCCTGAATTATATAAGCTAATTATATTTTGAGTTAAATCCATACAGGGGGCAAATAAAAAAAAGTCTATATTTCACTTTGCCCGCGACGCGAAAGTCACTGAt                                                                     >  2:148505‑M2/61‑250 (MQ=255)                                                                                                                                                                         cgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATgcgc                                                              <  2:419575‑M2/197‑1 (MQ=255)                                                                                                                                                                          gatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGcgcc                                                             >  1:55410‑M2/54‑251 (MQ=255)                                                                                                                                                                                caaagggccactagcacacctgagtcctctccaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGT‑AAGATTCCTCAGTGGTCGCACCAGCATCAATATTGCTAAAAGGAAGGATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTGGACTTTGCCCGCGCCGCGAAAGTCACTTATGTGGCGCCGCGTTCa                                                      <  2:720263‑M2/204‑1 (MQ=255)                                                                                                                                                                                 aaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTc                                                       >  1:804056‑M2/47‑250 (MQ=255)                                                                                                                                                                                        cactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAACACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTc                                               >  1:252280‑M2/40‑251 (MQ=255)                                                                                                                                                                                         actagcacacctgattcctctctaaataccctatacgaGTTGCGCTACTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAACAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTc                                               >  1:204625‑M2/39‑250 (MQ=255)                                                                                                                                                                                             gcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGACCTCAAAAt                                          <  1:580249‑M2/217‑1 (MQ=255)                                                                                                                                                                                              cacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATg                                         >  2:781997‑M2/34‑251 (MQ=255)                                                                                                                                                                                                    tgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAACACGGACGAGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTa                                    >  2:512412‑M2/28‑250 (MQ=255)                                                                                                                                                                                                     gattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAAc                                  >  2:882228‑M2/27‑251 (MQ=255)                                                                                                                                                                                                              ctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTAt                           <  1:559621‑M2/232‑1 (MQ=255)                                                                                                                                                                                                                aaataccctagccgaGGTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCATCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTa                       <  2:655810‑M2/236‑1 (MQ=255)                                                                                                                                                                                                                         atccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGATTTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTACAtttgttt              >  1:171918‑M2/7‑251 (MQ=255)                                                                                                                                                                                                                                     ctcctcTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTACATTTGTTTGCTGGCCTGGAt  >  1:631268/1‑251 (MQ=255)                                                                                                                                                                                                                                       cctcTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTACATTTGTTTGCTGGCCTGGAt  <  2:839515/249‑1 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                   GAGTGCAGACTTACAAAAATGCGGATTCTGGCTCAGATAAGCGCCTGGTTCGACATAAAGGTTGGTAAAGTTGCAGCGCCCACCGCCAGACATAAGGATTTTGCGCCCTGGTTTTTTACCATTATCGATCAGCAGAACCCGGCGTCCTGCCTGACCTGCCAGCGCAGAACAGAACATACCCGCCGCACCAGCGCCTATAATAATGGCATCAAACCTTTCCACGTTGCGCTCCTCTTAGAAAAAACGGGCGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTACATTTGTTTGCTGGCCTGGAT  >  NC_000913/3637189‑3637668

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