Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,637,436	G→A	intergenic (‑26/‑206)	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,436	0	G	A	100.0%	130.9 / NA	41	intergenic (‑26/‑206)	yhiN/pitA	putative oxidoreductase YhiN/metal phosphate:H(+) symporter PitA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (13/28);  total (13/28)

GAGTGCAGACTTACAAAAATGCGGATTCTGGCTCAGATAAGCGCCTGGTTCGACATAAAGGTTGGTAAAGTTGCAGCGCCCACCGCCAGACATAAGGATTTTGCGCCCTGGTTTTTTACCATTATCGATCAGCAGAACCCGGCGTCCTGCCTGACCTGCCAGCGCAGAACAGAACATACCCGCCGCACCAGCGCCTATAATAATGGCATCAAACCTTTCCACGTTGCGCTCCTCTTAGAAAAAACGGGCGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGC  >  NC_000913/3637189‑3637681                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      atgtgcacatcaaccacgacgactgtctggaaatcgccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaagaccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTg                                                                                                                                                                                                                                                   <  2:374224‑M2/29‑1 (MQ=255)   gtgcccatcaaccacgactactgtcttgaaatcgccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcgcacctgattcctctctaaatgccctatccgaGTTGTGCTCCTCTTA‑AAAAAACGGACGTGaa                                                                                                                                                                                                                                                 <  1:180796‑M2/31‑1 (MQ=255)     gcacatcaaccacgacgactgtctggaaatcgccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAAtt                                                                                                                                                                                                                                               <  1:179274‑M2/33‑1 (MQ=255)            aaccacgacgactgtctggaaatcgccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTaa                                                                                                                                                                                                                                           >  1:481621‑M2/212‑248 (MQ=255)                acgacgacggtctggaaatcgccgtgttgaaaggtgacatgggggacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaagaccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTcc                                                                                                                                                                                                                                    <  2:484429‑M2/44‑1 (MQ=255)                      actgtctggaaatcgccgtgttgaaaggtgacatgggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTg                                                                                                                                                                                                                              <  1:99864‑M2/50‑1 (MQ=255)                         gtcgtgaaatcgccgtgtggaaaggggacatgggtgacgggcagcattttgacgatgacgttatcggccagcgctgcgtgaggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacaccatattcctctctaaaggccctatccgaGTTGCGTTCCGCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCGTCAGTGGTc                                                                                                                                                                                                                           <  2:142904‑M2/53‑1 (MQ=255)                             ggaaatcgccgtgttgaaaggtgacatgggggacgtgcagcattctgccgattacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCAc                                                                                                                                                                                                                       <  2:430529‑M2/57‑1 (MQ=255)                                                ggtgacatgggttacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcactggtatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTaaa                                                                                                                                                                                                    <  1:339757‑M2/76‑1 (MQ=255)                                                         ggtgacgtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTAt                                                                                                                                                                                           <  1:231108‑M2/85‑1 (MQ=255)                                                               gtgcagcattttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCt                                                                                                                                                                                     <  1:261802‑M2/91‑1 (MQ=255)                                                                       ttttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCt                                                                                                                                                                                     >  1:133532‑M2/153‑243 (MQ=255)                                                                       ttttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCt                                                                                                                                                                                     <  2:133532‑M2/91‑1 (MQ=255)                                                                        tttgccgatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATtata                                                                                                                                                                             <  2:133598‑M2/99‑1 (MQ=255)                                                                           gccgatgacgttatcgcccatcgcggcgtgcggcacgggtatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATaa                                                                                                                                                                          <  1:303151‑M2/102‑1 (MQ=255)                                                                              gatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATa                                                                                                                                                                      <  1:22138‑M2/106‑1 (MQ=255)                                                                              gatgacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGAATTTGCCTGAATTATATAAGAt                                                                                                                                                                       >  2:463749‑M2/146‑250 (MQ=255)                                                                                 gacgttatcgcccagcgcggcgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAatt                                                                                                                                                                   >  1:450818‑M2/143‑251 (MQ=255)                                                                                         cgcctagcgcggcgtgcggcacgggcatttgtagtgcttgccgaagtaagattgagtctaatattgggtccgttggtgaacgatcaacaaagggctactagcactcctgattcctctctaaataccctgtccgaGTTTCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTTAGGGGTCGCACCAGCATCTATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATtatttttgg                                                                                                                                                           <  2:472355‑M2/117‑3 (MQ=255)                                                                                              agcgcggtgtgcggcacgggcatttgcagtgcttgccgaaggaagattgagtctattattgggtcctgtggtgaacgatcaacaaagggccactagcactcctgattcctctctaaataccctatccggGTTGCGCTCCTGTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAAGATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGa                                                                                                                                                      <  1:313224‑M2/122‑1 (MQ=255)                                                                                                                           ggcttgccgaaggaagattgagtctattattgggtcctttggtgaacgagcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCGCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATGACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATCCAGGGGGCAAATCAAAAAAAg                                                                                                                         <  2:467600‑M2/151‑1 (MQ=255)                                                                                                                               tgccgaaggaagattgagtatattattgggtcatttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTAAGTGGGCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATAAAAAAAAGTCt                                                                                                                      >  1:354381‑M2/97‑250 (MQ=255)                                                                                                                               tgccgaaagaagattgcgtctattattgggtcctttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaattcccgatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTc                                                                                                                       <  2:331073‑M2/153‑1 (MQ=255)                                                                                                                                 ccgaaggaagattgagtctattattgggtcctgtggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATTAATATTACTAAAAGGAAGTATTTGTCTGAATTATAGAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCtata                                                                                                                   <  1:217715‑M2/157‑1 (MQ=255)                                                                                                                                    aaggaagattgagtctattattgggtcctttggtgaacgatcaacaaagggcaactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGTGGGCAAATCAAAAAAAGTCTATAttt                                                                                                                >  1:332603‑M2/92‑251 (MQ=255)                                                                                                                                                                 ttggtgaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTg                                                                                    >  1:341539‑M2/63‑250 (MQ=255)                                                                                                                                                                    gtgaacgatcatcaacgggccactagcacacctgattcctctctaaataccctatccgaGGTGCGCTCATCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCACAGTGGTCGCCCCAGCATCAATCTTCCTAAAAGGAAGTATTTGCCTGAATTATATCAGATAATTATTTTTTGAAAGAAATCCATACGCGGGGCAAATCAAAAAAAGTCTATATTTNACTTTGCCCGCGAcgcgaaagtaccggata                                                                                 >  2:227984‑M2/60‑242 (MQ=255)                                                                                                                                                                      aaacgatcaacaaagggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTACGTGGTCGCACCAGCACCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAAACCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCACGCGCCGAGAAAGTCACTGATAATg                                                                              >  2:116633‑M2/58‑251 (MQ=255)                                                                                                                                                                             caacaaatggccactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCGGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGcgccgc                                                                        <  1:319424‑M2/200‑1 (MQ=255)                                                                                                                                                                                        cactagcacacctgattcctctctaaataccctatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCACACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTc                                                            <  1:486096‑M2/212‑1 (MQ=255)                                                                                                                                                                                             gcacgcctgattcctctctaaagtccctatccgaGGTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAAGTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGTGGGCAAATCAAAAAAAGTCTATATTTCACTTGGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAAt                                                       <  1:65393‑M2/217‑1 (MQ=255)                                                                                                                                                                                                                 aataccctatccgaGTTGCGCTACTCTTA‑AAAAAACGGACGTGAATTGAAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTCGAGTGAAATCCATAAAGGGGGCAAATCAAAAAAAGACTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGCTAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATgctcc                                   >  2:132535‑M2/15‑249 (MQ=255)                                                                                                                                                                                                                     ccctatccgaGGTGAGCTCCGCTTA‑AAAAAACGGACGTGAATTGTAAAGATTACTCAGTGGTCGCACCAGCATAAATATTACTAAAATGAAGTATTTGCCTGAAATATATAAGATAATTATTTTTTTAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATAATTCACTTTGCCCGCGCCGCGAAAGTCACTTATAATGCTCCGCGTTCATGTCCACAAAATGTCGTGACGTCCTATTCTCCAttt                               <  2:153204‑M2/241‑1 (MQ=255)                                                                                                                                                                                                                     cccgatccgaGTTGCGCTCCTCTTA‑AAAAAACGTAAGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTACTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCTCTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGTCGTAACGTCCTATGCTACAttt                               <  2:263712‑M2/241‑1 (MQ=255)                                                                                                                                                                                                                        tatccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTACAtttgtt                            >  1:435233‑M2/8‑251 (MQ=255)                                                                                                                                                                                                                         atccgaGTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTACAtttgttt                           <  1:207373‑M2/245‑1 (MQ=255)                                                                                                                                                                                                                               gTTGCGCTCCTCTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTACATTTGTTTGCTGGc                     >  1:62849/1‑251 (MQ=255)                                                                                                                                                                                                                                  tgcgctcacattaaaaaaaCTGACGTGACTTGTAAAGATTACTCAGTGGGCTCAACAGCATCAATATTACTAAAAGGAAGTAGTTGCCTGAATTATAAAAGATAATGATTTTTTGATTGAAATCCATACATGGGGCAAATCTAAAAAATGCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATTATTCTCCGCGTTCATGTCCGCATAATGGCGTATAGTCCTATGCTACATATGTTAGCTGGcc                    <  2:22834/237‑1 (MQ=255)                                                                                                                                                                                                                                       cctcTTA‑AAAAAACGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTACATTTGTTTGCTGGCCTGGAttt             <  1:200023/251‑1 (MQ=255)                                                                                                                                                                                                                                             taaaaaaaCGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTACATTTGTTTGCTGGCCTGGATTTGCAt         <  1:249943/248‑1 (MQ=255)                                                                                                                                                                                                                                                   aaCGGACGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTAACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGc  >  2:317032/1‑251 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      GAGTGCAGACTTACAAAAATGCGGATTCTGGCTCAGATAAGCGCCTGGTTCGACATAAAGGTTGGTAAAGTTGCAGCGCCCACCGCCAGACATAAGGATTTTGCGCCCTGGTTTTTTACCATTATCGATCAGCAGAACCCGGCGTCCTGCCTGACCTGCCAGCGCAGAACAGAACATACCCGCCGCACCAGCGCCTATAATAATGGCATCAAACCTTTCCACGTTGCGCTCCTCTTAGAAAAAACGGGCGTGAATTGTAAAGATTCCTCAGTGGTCGCACCAGCATCAATATTACTAAAAGGAAGTATTTGCCTGAATTATATAAGATAATTATTTTTTGAGTGAAATCCATACAGGGGGCAAATCAAAAAAAGTCTATATTTCACTTTGCCCGCGCCGCGAAAGTCACTGATAATGCGCCGCGTTCATGTCCTCAAAATGGCGTAACGTCCTATGCTACATTTGTTTGCTGGCCTGGATTTGCATACCGGGC  >  NC_000913/3637189‑3637681

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