Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,173,363	Δ2 bp	intergenic (‑490/+918)	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,363	0	C	.	90.0%	179.7 / 12.7	50	intergenic (‑490/+919)	gatD/gatB	galactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand):  ref base C (2/3);  new base . (34/11);  total (36/14)
Fisher's exact test for biased strand distribution p-value = 1.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*	NC_000913	2,173,364	0	C	.	92.0%	188.3 / 8.6	50	intergenic (‑491/+918)	gatD/gatB	galactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand):  ref base C (1/3);  new base . (35/11);  total (36/14)
Fisher's exact test for biased strand distribution p-value = 6.12e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

GCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTA‑CCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCG  >  NC_000913/2173123‑2173613                                                                                                                                                                                                                                                  ||                                                                                                                                                                                                                                                          gCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTA‑CCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATACAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTcc                                                                                                                                                                                                                                                >  2:439050/1‑251 (MQ=255)     gaGCCGCCAGTGGGGTGTGCAAACCAATAGTTTGCGTCGCGAGCCACAGGGTGATGCTCATAAGGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTA‑CCCGGCACACAAACAGCAATTAAAATGGTGAGTTGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCcagc                                                                                                                                                                                                                                            <  1:590454/250‑1 (MQ=255)      aGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTA‑CCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCcagc                                                                                                                                                                                                                                            >  1:408501/1‑250 (MQ=255)              gTTGGGTGTGCAAAACAAACGTTAGCGTCGAGATCAACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGCGCGGCACAGATTTCCACGATGCACGGCGACCGCCCTCGAAACGAAGAAGCCGATGGTGGCCAGATCGCCACACGGCAGCACCTGATCA‑CCCGGCACACAAGCAGCCATTACAATAGTGAGTGGGATCAAAATCAGGCTAGCCGGTACCACC‑‑GCCGTATGTCCCAGCAGCAACGcc                                                                                                                                                                                                                                   >  2:619273/1‑251 (MQ=255)                tGGGTTTGAAAACCAATAGGTTGGCTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGCGATTAAGGGGCGGAACAGTTTTCCACAAGGCCCGGCGCCCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTC‑CCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGAAAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCg                                                                                                                                                                                                                                  <  2:126621/250‑1 (MQ=255)                             cAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGCTGCGGAACAGATTTCCACCATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTA‑CCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATc                                                                                                                                                                                                                    >  1:113918/1‑251 (MQ=255)                                     tGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTA‑CCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTc                                                                                                                                                                                                            >  1:9881/1‑251 (MQ=255)                                            cgATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTA‑CCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccg                                                                                                                                                                                                     <  1:64363/251‑1 (MQ=255)                                             gATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTA‑CCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccgc                                                                                                                                                                                                    >  2:68089/1‑251 (MQ=255)                                                   caGGGTGATGCTCATAATGATGACACCCGCGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCACAAAACGGCAGCACCTGATTA‑CCCGGCACACAAACAGCAATTAAAATGGTGAGTGAGCTAAAAATAAGGCTTGCCAATA‑‑ACACCGCCGTATGTCCCAGAAGCAACGCAGGATCAAGGCCAATCAGGAACTACTGACCGCCGAACt                                                                                                                                                                                              >  1:60211/1‑251 (MQ=255)                                                    aGGGTGCTGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTA‑CCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCCACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACtt                                                                                                                                                                                             >  2:80348/1‑251 (MQ=255)                                                            gCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGCTGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTA‑CCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGt                                                                                                                                                                                      >  1:410855/1‑250 (MQ=255)                                                                aTAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTA‑CCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACg                                                                                                                                                                                 >  2:22687/1‑251 (MQ=255)                                                                              cGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTA‑CCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCa                                                                                                                                                                                                                             >  1:262580/1‑193 (MQ=255)                                                                                 ccgctctcccatcccaggATTCCCACGAGGCACGGCGACCGCAATCGCCACGAAGAAGCCGAGGGTGGCAAGACCGCAAAACGGCAGCACCTGATTA‑CCCGGCACAAAAACAGCAATTAAAAGGGTGAGTGGGAAAAAAACCAGGCTTGCCGATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAg                                                                                                                                                                <  2:552563/234‑1 (MQ=255)                                                                                 atttaagtggCGGAACAGATTTCCACGATCCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTA‑CCCGGCCCACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAg                                                                                                                                                                <  1:264635/242‑1 (MQ=255)                                                                                   ttAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCAACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTA‑CCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCg                                                                                                                                                              >  1:476422/1‑251 (MQ=255)                                                                                    tAAGGCGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTA‑CCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCg                                                                                                                                                              >  2:356682/1‑250 (MQ=255)                                                                                     aaGGTGAGGAACAGATTTCCACGAAGCCCGGCTACAGCCATCGCCACGCAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTA‑CCCGGCACACAAACAGCAATTAAAATGGTGCGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTCCGCCTGTACACGCCTACGAGCCTGCTTAGcaat                                                                                                                                                            >  2:69730/1‑248 (MQ=255)                                                                                       gggggggAAAAGATTTCCAGGAGGCAGGGCGACCGCCACCGCCACGAAGAAGCCGATGGTGGAAAGATCGCCAAACGGCAGCACCTGATTC‑CCCGGCACACAAACAGCAATTAAAAGGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATgg                                                                                                                                                          <  2:181776/246‑1 (MQ=255)                                                                                        gTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTA‑CCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATggg                                                                                                                                                         >  1:15316/1‑251 (MQ=255)                                                                                         tGCGGAACAGATTTCCACGATGCCCGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTA‑CCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGACATTCGCCAGTAAACGACAACGAGCCTGCTTAGCGAGGGGc                                                                                                                                                        >  2:382533/1‑251 (MQ=255)                                                                                                cAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTA‑CCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAg                                                                                                                                                                >  1:552563/1‑236 (MQ=255)                                                                                                   aTTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTA‑CCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGAAGGGCGTTAAAccat                                                                                                                                              >  2:186919/1‑251 (MQ=255)                                                                                                      gccACGATCCACGGCCACCCCCAGCACAACGAAGAAGCCGATGGTGGAAAAACCGCCAAACGGCAGCACCTGATTA‑CCCGGCACACAACCAGCAATTAAAATGGGGAGTGGGATAAAAATCAGGTGTGCCGATAC‑‑CTCCGCCGTATGCCCCAGCAGCAACGCCGGATCAAGGCCAAGCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccatcca                                                                                                                                           <  2:82005/250‑1 (MQ=255)                                                                                                       ccACGATGCCCGGCGCCCGCCACCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTA‑CCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccatcca                                                                                                                                           <  1:35027/250‑1 (MQ=255)                                                                                                       ccACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTA‑CCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTACGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccatcca                                                                                                                                           <  1:162523/250‑1 (MQ=255)                                                                                                        cACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGAACGCCAAACGGCAGCACCTGATTA‑CCCGGCACACCGAAAGCACTTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACACCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatg                                                                                                                                         >  1:680300/1‑251 (MQ=255)                                                                                                                 cGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTA‑CCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGCACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGttt                                                                                                                                 >  1:85141/1‑250 (MQ=255)                                                                                                                  ggCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTA‑CCCGGAACAAAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACAACC‑‑GAAGTATGTCCCCGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTa                                                                                                                                >  1:423848/1‑250 (MQ=255)                                                                                                                      aCCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGAACCTGATTA‑CCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTACACCATCCATGATGGGTTTAATCAc                                                                                                                           >  1:82902/1‑251 (MQ=255)                                                                                                                          ccATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTA‑CCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGt                                                                                                                       >  1:247514/1‑251 (MQ=255)                                                                                                                                ccACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTA‑CCCGGCACACAAACAGCAATTAAAAGGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCAtt                                                                                                                 <  1:68089/251‑1 (MQ=255)                                                                                                                                 cACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTA‑CCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCAtt                                                                                                                 >  2:619544/1‑250 (MQ=255)                                                                                                                                    gaagaaGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTA‑CCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTagc                                                                                                              >  1:413259/1‑250 (MQ=255)                                                                                                                                      agaagCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTA‑CCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGCTACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAAATTCGCCTGTAAACGACTACGCGCCTGCTTAGCGATGGGCGTTAAACCCTCCATGATGGGTTTACTCACCCGAGGCATTagcagc                                                                                                           >  2:610884/1‑251 (MQ=255)                                                                                                                                      aaacgccattggtgaCAAGATCACCAAACGGCAGCCCCT‑ATTAACCCGCCCCAAAAAAAGAAATAAAAAGGGGGAGGGGGATAAAAATCAGGCTTGCCGATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTagcagc                                                                                                           <  2:496683/242‑1 (MQ=255)                                                                                                                                            cGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTA‑CCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACt                                                                                                     >  2:520759/1‑251 (MQ=255)                                                                                                                                                ggtgcCAAGATCGCCAAACGCCACCACCTGATTA‑CCCGGCACACAACCAGAATTTAAAAGGTTGAGTGGGATAAAAATCAGTCTTGCCGAT‑CC‑CCCCGCCGTATGTCCCAGCATCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCCCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTgccg                                                                                                 <  1:404742/251‑1 (MQ=255)                                                                                                                                                  gggCAAGACCGCCAAACGGCAGCACCTGATTC‑CCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCAGCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTgccgcc                                                                                               <  1:80348/250‑1 (MQ=255)                                                                                                                                                     cAAGATCGCCAAACGGCAGCACCTGATTA‑CCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGtt                                                                                            >  1:171293/1‑251 (MQ=255)                                                                                                                                                            gCCAAACGGCAGCACCTGATTA‑CCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCg                                                                                     >  1:150597/1‑251 (MQ=255)                                                                                                                                                               aaaCGGCAGCACCTGATTA‑CCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGcc                                                                                   >  1:185805/1‑250 (MQ=255)                                                                                                                                                                      agcaCCTGATTA‑CCCGGCACACAACCAGCAATTAAACTGGTGAGTGGGATAAAAATCAGGCTTGACGATACCACCC‑‑CCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGACTGTAAACGCCTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTCACCGCCAGCTGCa                                                                            >  2:374605/1‑250 (MQ=255)                                                                                                                                                                                        cacacacaCAGCAATTCAAATGGCGAGTGGGCTAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGt                                                          >  2:80160/1‑251 (MQ=255)                                                                                                                                                                                                                  gTGGGAAAAAAATCAGGCGTGCGGAT‑CC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGc                                <  1:154637/251‑1 (MQ=255)                                                                                                                                                                                                                          aaaaaacagggctgcccatacccaccgccgTATGTCCCAGCAGCAACGCCGGAACAAGGCCAACCAGGAACTCCTGACCGCCGAACTGCGCCTGTAAACTACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACAGCGTAACCCGCGAGTATGCCGATAATCAGCCCCATc                          <  2:408501/227‑1 (MQ=255)                                                                                                                                                                                                                                aGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGACCTTCGCCTGTAAACGACTACGAGCCTGCTTAACGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCAAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCATCCCCATCACaaaccc                  >  2:624211/1‑248 (MQ=255)                                                                                                                                                                                                                                aGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCAAATCAGGAACTCCTGACCGCAGAACTTCGCCTGTAAAAGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAAcc                  >  1:663195/1‑251 (MQ=255)                                                                                                                                                                                                                                 ggCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCCATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAAccc                 >  2:428253/1‑251 (MQ=255)                                                                                                                                                                                                                                  aggttgccgtacacccccGCCGTATGTCCCAGCAGCACCGCCGGATCAAGCCCAACCAGGACCTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGTTGGGCGTGAAACCATCCATGATGGGTTGAATCCCCCGTGGCATTATCAGTATCACTGCCGCCGTTTTTACCGCCAGCGGCAAGACACCTTTTACATCGTGACCCGCGAGGATGCCGATAATCTTCCCCAGCACAAAAcc                  <  2:423848/237‑1 (MQ=255)                                                                                                                                                                                                                                     tGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATAACTGCCGCCGTTTTTACCGCTCGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACgg             >  1:432930/1‑251 (MQ=255)                                                                                                                                                                                                                                      tgccctataccaccgccgTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGATCTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCGTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGTATGCCGATAATCAGCCCCATCACAAAACCCAc               <  1:610884/238‑1 (MQ=255)                                                                                                                                                                                                                                       tgccgatcccaccgccgTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCAc               <  1:22687/238‑1 (MQ=255)                                                                                                                                                                                                                                             taccaccgccgTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGAc         >  2:189121/6‑249 (MQ=255)                                                                                                                                                                                                                                              accaccgccgTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGCACTCCTGACCGCCGAACTTCGCCTGACAACGACTACGAGCCTGCTTAGCGATGGCCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAgg      >  1:80272/5‑251 (MQ=255)                                                                                                                                                                                                                                              accaccgccgTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAg       <  1:439050/246‑1 (MQ=255)                                                                                                                                                                                                                                                 cccGCCGTATGGCCCAGCTGCAACGCCGGATCAGGGCCAATCAGGAATGCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAAGCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAAGACGCCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTc   <  1:374605/251‑1 (MQ=255)                                                                                                                                                                                                                                                  ccgccgTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCg  <  2:185805/251‑1 (MQ=255)                                                                                                                                                                                                                                                  ||                                                                                                                                                                                                                                                          GCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTA‑CCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCG  >  NC_000913/2173123‑2173613

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