Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,173,448	C→T	intergenic (‑575/+834)	gatD ← / ← gatB	galactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,173,448	0	C	T	100.0%	138.1 / NA	41	intergenic (‑575/+834)	gatD/gatB	galactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (19/22);  total (19/22)

GAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGG‑CAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGT‑AAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCGCCAAATGGACCAAATTTGCGCTGAATATCGTCGGCGCTAAATTTAATTCGGTTAACGCCTGGGATTTTCTCGATGATA  >  NC_000913/2173201‑2173691                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                    gtgctctcccaatctaaCAGATTCCCACGATGCACGGCGCCCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGG‑CAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAACCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGT‑AAACGACTACGAgcttgctt                                                                                                                                                                                                                                                <  1:366234/236‑1 (MQ=255) gagaTTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGG‑CAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGT‑AAACGACTACGAgcttgctt                                                                                                                                                                                                                                                >  2:114906/1‑251 (MQ=255)   gaTTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGG‑CAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGT‑AAACGACTACGAgcttgctt                                                                                                                                                                                                                                                >  1:228165/1‑249 (MQ=255)            gcgaaaCAGATTCCCCCGATGCACGGCGCCCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGG‑CAGCACCTGATTCCCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGT‑AAACGACTACGAGCTTGCTTAGCGATGGGCg                                                                                                                                                                                                                                     <  1:11451/247‑1 (MQ=255)              ggAACAGTTTCCCACGAGGCACGGCGACCGCCACCGCCACGAAGAAGCCGATGGGGGCAAGATCGC‑AAACGGGCAGCCCCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGGGGGATAAAAACCAGGCTGGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGTCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGT‑AAACGACTACGAGCTTGCTTAGCGATGGGCGtt                                                                                                                                                                                                                                   <  2:350357/251‑1 (MQ=255)                aaCAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAACCGG‑CAGCACCTGATTACCCGGCACACACACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGT‑AAACGACTACGAgcttgctt                                                                                                                                                                                                                                                >  2:366234/1‑236 (MQ=255)                         ccACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGG‑CAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGT‑GAACGACTACGAGCTTGCTTAGCGATGGGCGTTCAAccatccat                                                                                                                                                                                                                        >  2:266553/1‑251 (MQ=255)                                            ccATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGG‑CAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGT‑AAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGt                                                                                                                                                                                                     <  1:91780/251‑1 (MQ=255)                                             caccGCCCCGAAGAGGCCGATGGGGGCAAGATCGCAAAAGGC‑CAGCCCCTGATTACCCGGCACACAAACAGAAATAAAAAGGGGGAGTGGGATAAAAATCAGGTTTGCCAATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGT‑AAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTg                                                                                                                                                                                                    <  2:334811/248‑1 (MQ=255)                                              aTCGCCACGCAGAAGCCGATGGTGGCAAGATCGCCAAACGG‑CAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGCGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGT‑AAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTgg                                                                                                                                                                                                   >  2:240965/1‑251 (MQ=255)                                                                          gATCGCCAAACGG‑CAGAACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGAGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGT‑CAACGACTACGAGCTTGCTTAGCGATGGTCGTTCAACCATCAATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGtttt                                                                                                                                                                        >  2:215698/1‑250 (MQ=255)                                                                            tCGCCAAACGG‑CAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGT‑AAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTAc                                                                                                                                                                     >  2:195057/1‑251 (MQ=255)                                                                               ccAAACGG‑CAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGT‑AAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGc                                                                                                                                                                  >  1:326752/1‑251 (MQ=255)                                                                                     gg‑CAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGT‑AAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTg                                                                                                                                                            >  1:89753/1‑251 (MQ=255)                                                                                     gg‑CAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGT‑AAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTg                                                                                                                                                            <  1:240965/251‑1 (MQ=255)                                                                                              cTGATTACCCGGCACACAAACAGCAATTCAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCACTCAGGAACTCCTGACCGCCGAACTTCGCCTGT‑AAACGCCTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATacac                                                                                                                                                    >  2:45203/1‑251 (MQ=255)                                                                                                              caAACAGCAATTAAAATGGTGAGTGGTATAAAAATCAGTCTTGCCTTTCCCACC‑‑GCCGTTTGTCCCAGCAGCAACGCCGGATCAATGCCAATTAGGAACTCCTGACCGCCGAACTTCGCCTGT‑AAACGACTACGAGCTTGCGTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCTTCACTGCCGCCGTTGTTACCGCCAGCTGCAATACACCTTTGACATCGTAAcc                                                                                                                                    <  1:364133/251‑1 (MQ=255)                                                                                                                aaCAGCAATTAAACTGGCGAGAGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGAATGTCACAGCAGCCACGACGGATCAAGGCCAACCAGGAACTCCAGACCGCCGAACTACGCCTGA‑AAACGACTACGAACTTGCTTAGCAATGGGCATTAAACCATCCATGATGAGTTTAATCCCCCGTTACATTAGCAGCATCACTGCCGCCTTTTTTACCGCCAGCTGCAATACACCTTTGACATCGAAAcacc                                                                                                                                  >  2:414545/1‑248 (MQ=255)                                                                                                                aaCAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGT‑AAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGttttt                                                                                                                                                                       <  1:37087/214‑1 (MQ=255)                                                                                                                aaCAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGT‑AAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGttttt                                                                                                                                                                       >  2:37087/1‑214 (MQ=255)                                                                                                                    gCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGT‑AAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCcgcg                                                                                                                                >  1:245830/1‑249 (MQ=255)                                                                                                                                       ggATAAAAACCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGAAACTCCTGACCGCCGAACTTCGCCTGT‑AAACGACTACGAGCTTGCTTAGCGATGGGCTTTAAACCAGCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACTCCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGcccc                                                                                                           <  1:190738/251‑1 (MQ=255)                                                                                                                                            aaaaTCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGT‑AAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATcac                                                                                                      >  1:78912/1‑251 (MQ=255)                                                                                                                                             aaaccAGGCTGGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGT‑AAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATTACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATcaca                                                                                                     <  1:114906/247‑1 (MQ=255)                                                                                                                                                        gCCGATACCACC‑‑GACGTCTGTACCAGCAGCACCGCCGGATCCAGGCCAATCAGGAACTCATGACCGCCGAACTTCGCCTGT‑AAACGCCTACGAGCTTGCTTCGCGATGGGCGTTAAACCATACATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACggg                                                                                          >  2:392419/1‑250 (MQ=255)                                                                                                                                                         ggcccaacccgccgccgTAGGTCCAAGCAGCAACGCCGGAACAAGGCCAAACAGGACCTCCTGACCGCCGAACTTCGCCTGT‑AAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGTTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCAc                                                                                             <  2:84828/238‑1 (MQ=255)                                                                                                                                                            gatccaccCC‑CCTTATCTCCCAGCAGCAACGCCGGTTCAAGGCCAATCATGACCACCTGACCGCCGAACTTCGCCTGT‑AAACGACTACGAGCTTGCTTAGCGATGGGCGTTAACCCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTg                                                                                         <  2:83565/246‑1 (MQ=255)                                                                                                                                                                        cgTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGT‑AAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCAcc                                                                                                                                                                                                        <  2:182701/127‑1 (MQ=255)                                                                                                                                                                        cgTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGT‑AAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCAcc                                                                                                                                                                                                        >  1:182701/1‑127 (MQ=255)                                                                                                                                                                             gTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGT‑AAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCGCCAAATgg                                                                        >  1:260020/1‑250 (MQ=255)                                                                                                                                                                             gTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGT‑AAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCGCCAAATgg                                                                        <  1:27402/250‑1 (MQ=255)                                                                                                                                                                               cccAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGT‑AAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATc                                                                                                        <  1:149079/216‑1 (MQ=255)                                                                                                                                                                               cccAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGT‑AAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATc                                                                                                        >  2:149079/1‑216 (MQ=255)                                                                                                                                                                                  agcagcACCGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCAGAACTTCGCCTGT‑CAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCGCCAAATGGAcaaa                                                                   >  2:421014/1‑247 (MQ=255)                                                                                                                                                                                     agcaACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGT‑AAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCGCCAAATGGACCAAAttt                                                               >  1:289481/1‑251 (MQ=255)                                                                                                                                                                                               gATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGT‑AAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCGCCAAATGGACCAAATTTGCGCTGAat                                                      >  2:144035/1‑250 (MQ=255)                                                                                                                                                                                                  aaaGGCCAATCATGAACTCCTGACCGCCGATCTTCGAC‑GTGAAACTACTACGGGCTGGTTTAGCGATGGGCGTTAAACCACCCAAGAGGGGTTTAAGCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCCTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCGCGCCAAATTGACCAAATTTGCGCTGAATAtcg                                                  <  2:162636/250‑1 (MQ=255)                                                                                                                                                                                                           ccAGGAACTCCTGACCGCCGAACTTCGCCTGT‑AAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCGCCAAATGGACCAAATGTGCGCTGAATATCGTCGGCGCTa                                         <  2:245830/250‑1 (MQ=255)                                                                                                                                                                                                              ggAACTCCTGACCGCCGAACTTCGCCGTT‑AAACGTCTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCGCCAAATGGACCAAATTTGCGCTGAATATCGTCGGCGCTAAAt                                      <  2:326752/251‑1 (MQ=255)                                                                                                                                                                                                                   cccTGACCGCCGAACTTCGCCTGT‑AAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCGTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCGCCAAATGGACCAAATTTGCGCTGAATATCGTCGGCGCTAAATTTAAt                                 <  2:401090/250‑1 (MQ=255)                                                                                                                                                                                                                                    cGCCTGT‑AACCGACTACGGGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCGCCAAATGGACCAAATTTGCGCTGAATATCGTCGGCGCTAAATTTAATTCGGTTAACGCCTGGGa                <  2:228165/251‑1 (MQ=255)                                                                                                                                                                                                                                       cTGT‑AAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCGCCAAATGGACCAAATTTGCGCTGAATATCGTCGGCGCTAAATTTAATTCGGTTAACGCCt                    <  1:161973/244‑1 (MQ=255)                                                                                                                                                                                                                                       cTGT‑AAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCGCCAAATGGACCAAATTTGCGCTGAATATCGTCGGCGCTAAATTTAATTCGGTTAACGCCt                    >  2:161973/1‑244 (MQ=255)                                                                                                                                                                                                                                         gT‑AAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCGCCAAATGGACCAAATTTGCGCTGAATATCGTCGGCGCTAAATTTAATTCGGTTAACGCCTGGGATTTtc           <  1:349861/251‑1 (MQ=255)                                                                                                                                                                                                                                                    aCGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCGCCAAATGGACCAAATTTGCGCTGAATATCGTCGGCGCTAAATTTAATTCGGTTAACGCCTGGGATTTTCTCGATGATa  <  2:161286/250‑1 (MQ=255)                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                    GAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGG‑CAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGT‑AAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCGCCAAATGGACCAAATTTGCGCTGAATATCGTCGGCGCTAAATTTAATTCGGTTAACGCCTGGGATTTTCTCGATGATA  >  NC_000913/2173201‑2173691

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