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	.	87.3%	199.6 / 22.6	55	intergenic (‑490/+919)	gatD/gatB	galactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand):  ref base C (1/6);  new base . (30/18);  total (31/24)
Fisher's exact test for biased strand distribution p-value = 3.52e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*	NC_000913	2,173,364	0	C	.	87.3%	204.0 / 18.3	55	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/6);  new base . (30/18);  total (31/24)
Fisher's exact test for biased strand distribution p-value = 3.52e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAG‑ATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGC‑ACCTGA‑TTACCCGGCACACAAACA‑GCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCG  >  NC_000913/2173125‑2173613                                                                                                                                                                                                                                                    ||                                                                                                                                                                                                                                                          aTTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAG‑ATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGC‑ACCTGA‑TTACCCGGCACACAAACA‑GCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGAATGTccc                                                                                                                                                                                                                                               >  2:468953/1‑250 (MQ=255)        gccAGTTGGGTGTGCAACCCAATAGTTTGCGTCGCGATCCCCAGGGTGATGCTCATAATGATGACACCCGAG‑ATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGC‑ACCTGA‑TTACCCGGCACACAAACA‑GCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTTCCGATACCACC‑‑GCCGTATGTCCCAGCAGCaa                                                                                                                                                                                                                                       <  1:1493/251‑1 (MQ=255)          cAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAG‑ATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGC‑ACCTGA‑TTACCCGGCACACAAACA‑GCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACg                                                                                                                                                                                                                                     >  1:222471/1‑251 (MQ=255)           aGTTGGGTGTGCAAACCAATAGTTGGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAG‑ATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGC‑ACCTGA‑TTACCCGGCACACAAACA‑GCAATTAAAAGGG‑TGAGTGGGATAAAAATCAGGCTTGCCGAGACCACC‑‑GCCGTATGTCCCAGCAGCAACGc                                                                                                                                                                                                                                    <  1:183325/251‑1 (MQ=255)                ggTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAG‑ATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGC‑ACCTGA‑TTACCCGGCACACAAACA‑GCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGAt                                                                                                                                                                                                                               >  2:286624/1‑251 (MQ=255)                   gtgCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAG‑ATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGC‑ACCTGA‑TTACCCGGCACACAAACA‑GCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCa                                                                                                                                                                                                                             >  2:557440/1‑250 (MQ=255)                   gtgCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAG‑ATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGC‑ACCTGA‑TTACCCGGCACACAAACA‑GCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCa                                                                                                                                                                                                                             >  1:215365/1‑250 (MQ=255)                   gtgCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAG‑ATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGC‑ACCTGA‑TTACCCGGCACACAAACA‑GCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCa                                                                                                                                                                                                                             >  1:165444/1‑250 (MQ=255)                      cAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAG‑ATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGC‑ACCTGA‑TTACCCGGCACACAAACA‑GCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGc                                                                                                                                                                                                                         >  2:433132/1‑251 (MQ=255)                            aaTAGTTTGCGTGGCGATCCACAGGGTTAGGCCCATAATGATGAC‑CCCGAGGATTAAGGTGCGAAACAAATTCCCACGATGCCCGGCGACCGCCATCGCCCGGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGC‑ACCTGG‑TTACCCGGCACACAAACA‑GCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCa                                                                                                                                                                                                                   <  2:295467/251‑1 (MQ=255)                             aTAGTTTGCGTGGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAG‑ATAAAGGTGCGGAACAGATTTCCACGCTGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGC‑ACCTGA‑TTACCCGGCACACAAACA‑GCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAg                                                                                                                                                                                                                  <  2:350124/251‑1 (MQ=255)                             aTAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACACGAG‑ATTAAGGGGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGC‑ACCTGA‑TTACCCGGCACACAAACA‑GCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATc                                                                                                                                                                                                                    >  1:584040/1‑249 (MQ=255)                                      gTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAG‑ATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGC‑ACCTGA‑TTACCCGGCACACAAACA‑GCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTg                                                                                                                                                                                                         >  2:316612/1‑251 (MQ=255)                                            ccccgacgcgcttccaatctAATGATGACCCCCGAA‑ATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAAATCGCCAAACGGCAGC‑ACCTGA‑TTACCCGGCACACAAACA‑GCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccgcc                                                                                                                                                                                                   <  2:341513/232‑1 (MQ=255)                                             gccACCGGGTGATGCTCATAATGATGACCCCCGAG‑ATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCAGATGGTGGCAAGATCGCCAAACGGCAGC‑ACCTGA‑TTACCCGGCACACAAACA‑GCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccgccg                                                                                                                                                                                                  >  2:112150/2‑251 (MQ=255)                                                ccaGGGTGATGCCCAAAAGGATGACCCCCGAA‑ATTAAGGGGCGGAACAGATTCCCACAATGCACGGCGCCCGCCATCGCCACGAAGAAGCCAATGGTGGCAAGATCCCCAAACGGCAGC‑ACCTGA‑TTACCCGGCACACAAACA‑GCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTGGCCGATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAAc                                                                                                                                                                                               <  2:588096/250‑1 (MQ=255)                                                   gggTGATGCTCATAAGGATGACACCCGAG‑ATTAAGGGGCGGAACAGTGTTCCACGATGCACGGCACCCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAACCGGCAGC‑ACCTGA‑TTACCCGGCACACAAACA‑GCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTGCCGTTACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTc                                                                                                                                                                                            <  1:320998/251‑1 (MQ=255)                                                    ggggATGCTAATAAGGATGCCCCCCGAG‑ATTAGGGGGGGGAAAAGATTTCCCGGAGGCACGGCGACCGCCAGCGCCACGAAGAAGCCGAGGGTGGCAAGATCGCCAAACGGCAGC‑ACCTGA‑TTACCCGGCACACAAACA‑GCAATTAAAATGG‑TGAGTGGGATAAAAAGCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCg                                                                                                                                                                                           <  2:222235/248‑1 (MQ=255)                                                       gATGCTCATAATGATGACACCCGAG‑ATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGC‑ACCTGA‑TTACCCGGCACACAAACA‑GCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCt                                                                                                                                                                                        <  1:583912/251‑1 (MQ=255)                                                                aaTGATGACACCCGAG‑ATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGC‑ACCTGA‑TTACCCGGCACACAAACA‑GCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccgcc                                                                                                                                                                                                   >  1:341513/1‑231 (MQ=255)                                                                    atgaCACCCGAG‑ATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGC‑ACCTGC‑TTACCCGGCACACAAACA‑GCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACg                                                                                                                                                                           >  2:49203/1‑251 (MQ=255)                                                                          cccGAG‑ATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCACCAAACGGCAGC‑ACCTGA‑TTACCCGGCACACAAACA‑GCAATTAAAATGG‑GGAGTGGGATAAAAATCAGGCTTGCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGCGCCTg                                                                                                                                                                     >  1:235853/1‑251 (MQ=255)                                                                                   tAAGGTCCGGAACAGATTTCCACAATGCCCGGCGACCGCCACCGCCACGAAGAACCCGATGGTGGCAAGATCGCCAAACGGCAGC‑CCCTGA‑TTCCCCGGCACAAAAACA‑GCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCg                                                                                                                                                              <  1:572424/250‑1 (MQ=255)                                                                                    aaGGTGCGGACCAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGC‑ACCTGA‑TTACCCGGCACACAACCA‑GCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGa                                                                                                                                                             <  1:510053/250‑1 (MQ=255)                                                                                          cggacaaGATTCCCACGATGCACGGGGACCGCAATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAACCGGCAGC‑ACCTGA‑TTCCCCGGCACACAAACA‑GAAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTTTAAACGACTACGAGCCTGCTTAGCGATGGGCGt                                                                                                                                                      <  1:468953/245‑1 (MQ=255)                                                                                           ggAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGC‑ACCTGA‑TTACCCGGCACACAAACA‑GCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTGCTGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATGAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGAGGGGCGtt                                                                                                                                                     >  2:358309/1‑251 (MQ=255)                                                                                                          gATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGC‑ACCTGA‑TTACCCGGCACACAAACA‑GCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatgatg                                                                                                                                      >  1:179964/1‑251 (MQ=255)                                                                                                           aTGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGC‑ACCTGA‑TTACCCGGCACACAAACA‑GCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATgg                                                                                                                                     >  1:166715/1‑251 (MQ=255)                                                                                                               gcgccgaccccctccGCCACGAAGAAGCCGATGGTGTCAAGATCGCCAAGCGG‑AGCCACCTGT‑TTACCCGGCACACAAACA‑GCAATTAAAAGGG‑TGATTGGCATAAAAACCAGGTTGGCCGATACCACC‑‑GCCGGATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGttt                                                                                                                                 <  1:316612/237‑1 (MQ=255)                                                                                                                    gACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGC‑ACCTGA‑TTACCCGGCACACAAACA‑GCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATc                                                                                                                             >  1:250844/1‑250 (MQ=255)                                                                                                                          cATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGC‑ACCTGA‑TTACCCGGCACACAAACA‑GCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTg                                                                                                                      >  1:263013/1‑251 (MQ=255)                                                                                                                              gCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGC‑ACCTGA‑TTACCCGGCACACAAACA‑GCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCAt                                                                                                                  >  2:391940/1‑251 (MQ=255)                                                                                                                                       aagCGGATGGGGGAAAGATCGCCAACCGGCAGC‑ACCTGA‑TTCCCCGCCACACAAACA‑GCATTTAAAATGG‑GGAGGGGGATAAAAATCAGGCTTGCCGATGCCACC‑‑TCCGTATGTCCCAGCAGCAACGCCGTATCAAGGCCAATCAGGAACTCCTGACCGCCGATCTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGTGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCAt                                                                                                         <  1:36313/251‑1 (MQ=255)                                                                                                                                               ggtggCAAGATCGCAAAACGGCA‑C‑AGCGGAATTCCCCGCCACCCAACCA‑GCATTAAAAATGG‑TGAGGGGGATAAAAATCAGCCGTGCCTATAC‑‑CCCCGCCGTATGTCCCAACAGCAACGCCGGATCAAGGCCAATCAGGAACTCCGGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCGCCCGTGGCATTATCAGCAGCACTgcc                                                                                                  <  2:554368/250‑1 (MQ=255)                                                                                                                                                     aagcccGCCAAACGGCAGC‑CCTTGA‑TTCCCCGGCCAACAAACA‑G‑AATTTAAATGGGTGAGGGGGAAAAAACTCAGCCTTGCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGttt                                                                                           <  2:348222/246‑1 (MQ=255)                                                                                                                                                         ccGCCAACCGGCAGC‑CCCTGA‑TTACCCGGCACACAAACA‑GCAATAAAAAGGG‑TGAGTGGGAAAAAAAGCAGGTTTGCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTAc                                                                                       <  1:49203/250‑1 (MQ=255)                                                                                                                                                              aaaccggaggccccTGA‑T‑ACCCGGGCCACAAACAAGAAATAAAAATGG‑TGAGTGGGAAAAAAAGAAGGTTGGCCGAT‑CC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGCCCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGcc                                                                                   <  2:250844/239‑1 (MQ=255)                                                                                                                                                              aaaCGGCAGC‑CCCTGA‑TTACCCGGCACACAAACA‑GCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCa                                                                                  <  1:360745/251‑1 (MQ=255)                                                                                                                                                                      gc‑aCGTGA‑TTACCCGCCACACAAACA‑GCAATAAAAAGGG‑TGGGGGGGAAAAAAACCAGGCTTGCCGAT‑CC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGAACCTCCTGACCGCCGAACCTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAAt                                                                          <  2:222471/251‑1 (MQ=255)                                                                                                                                                                         aCCTGA‑TTACCCGGCACACAAACA‑GCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATac                                                                        >  1:482583/1‑251 (MQ=255)                                                                                                                                                                           cTGA‑TTACCCGGCACACAAACA‑GCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATacac                                                                      >  1:169822/1‑251 (MQ=255)                                                                                                                                                                                ttACCCGGCACACAAACA‑GCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCCGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCttt                                                                  >  1:493380/1‑251 (MQ=255)                                                                                                                                                                                ttACCCGGCACACAAACA‑GCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCCGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCttt                                                                  >  1:43666/1‑251 (MQ=255)                                                                                                                                                                                 accccGGCACAAAAAAA‑CAAATTAAAAGGG‑TGAGTGGAAAAAAAACCGGGCTGGCCGATACAACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTg                                                                 <  1:557440/249‑1 (MQ=255)                                                                                                                                                                                  aCCCGGCACACAAACA‑GCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTACACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTg                                                                 >  2:520947/1‑250 (MQ=255)                                                                                                                                                                                      ggaacacaACCA‑GAAATAAAAAGGG‑TGGGTGGGATAAAAACCAGCCTTGCCGATACCACC‑‑GCCGTATGCCCCAGCAGCAACGCCGGATCAAGGCCAATCGGGAACTCCTGCCCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATc                                                            <  2:241603/248‑1 (MQ=255)                                                                                                                                                                                          cacaAACA‑GCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTa                                                         >  1:44181/1‑250 (MQ=255)                                                                                                                                                                                               aCA‑GCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCcgc                                                   >  2:590332/1‑251 (MQ=255)                                                                                                                                                                                                          aaaTGG‑CGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGa                                              >  2:286319/1‑246 (MQ=255)                                                                                                                                                                                                           aaagggcgAGTGGGATAAAAATCAGGCTGGCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGa                                              <  1:286319/239‑1 (MQ=255)                                                                                                                                                                                                            aTGG‑TGAGTGGGGTAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGAt                                       >  1:416521/1‑251 (MQ=255)                                                                                                                                                                                                                  gattggaaTAAAAATCAGGCTCGCCGATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACGGCCGCCGTTTTTACCGCCAGCTGc                                                                             <  1:183581/201‑1 (MQ=255)                                                                                                                                                                                                                  gAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGc                                                                             >  2:183581/1‑208 (MQ=255)                                                                                                                                                                                                                           aaaaaTCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCTCCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATc                          >  1:189668/1‑250 (MQ=255)                                                                                                                                                                                                                                      ttGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACg              >  1:33381/1‑251 (MQ=255)                                                                                                                                                                                                                                              atcccaccgccgTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGc                                           <  1:1328/210‑1 (MQ=255)                                                                                                                                                                                                                                              ataccaccgccgTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGc                                           >  2:1328/7‑216 (MQ=255)                                                                                                                                                                                                                                                  ccccGCCGTATGTCCCAGAACCACCGCCGGACCAAGCCCAATCAGGAACCCCTGCCCGCCGAACTTCCCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATTATGGGTTTAATCCCCCGGGGCATTAGCAGCATCCCTGCCGCCGTTTTTACCGCCAGCTGCAATCCACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCt    <  2:169822/251‑1 (MQ=255)                                                                                                                                                                                                                                                   cccGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGACCTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTc   <  2:33381/251‑1 (MQ=255)                                                                                                                                                                                                                                                   accgccgTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGc                                >  1:439188/2‑222 (MQ=255)                                                                                                                                                                                                                                                    ccgcggTATGTCCCAGCAGCAACCCTGTATCACGCCCAATCAGAAACCCCCCCCCGCCGAACTTCGCCTGAAAACCACTACGCTCCTGCTTAGCGAGGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGGTTTTACCGCCAGCTGTAATACAACTCTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAGCCCACGGGGCCAGGCTCg  <  2:235853/251‑1 (MQ=255)                                                                                                                                                                                                                                                    ||                                                                                                                                                                                                                                                          ATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAG‑ATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGC‑ACCTGA‑TTACCCGGCACACAAACA‑GCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCG  >  NC_000913/2173125‑2173613

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