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	97.9%	164.7 / NA	47	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);  major base T (25/21);  minor base A (1/0);  total (26/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

AAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCGCCAAATGGACCAAATTTGCGCTGAATATCGTCGGCGCTAAATTTAATTCGGTTAACGCCTGGGATTTTCTCGATG  >  NC_000913/2173207‑2173688                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                 aaGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGAt                                                                                                                                                                                                                                       >  2:215834/1‑251 (MQ=255)  aGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATg                                                                                                                                                                                                                                      <  1:30904/251‑1 (MQ=255)     ggCGGAACAGATTTCCCCGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGc                                                                                                                                                                                                                                   <  1:400794/250‑1 (MQ=255)           aCAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTaaa                                                                                                                                                                                                                             <  1:215834/251‑1 (MQ=255)                  tCCACGAGGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGG‑TGAGTGGGATAAAAACCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAAccatcc                                                                                                                                                                                                                       <  1:327878/250‑1 (MQ=255)                   ccACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCACACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGATTGCTTAGCGATGAGCGTTAAAccatccat                                                                                                                                                                                                                     >  2:252758/1‑251 (MQ=255)                      cGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCatgat                                                                                                                                                                                                                  >  1:266318/1‑251 (MQ=255)                                        aTCGCAACGAAGACGCCGATGACGGCACGGTAGCCAAACGGCAGCACCTGACTAACCGGCACACACACAGCAATTCAAATGC‑TGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACACCGGATAAAGGACAATCAGGAACTTCTGCCAGCCGCACTTCGCCTGTAACAGACTACGAGCTTTCTTAGCGCGGGTCGTTAAACCATTCATGATGGGTTTAATCACCCTTgg                                                                                                                                                                                                >  2:426292/1‑251 (MQ=255)                                            ccACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCAtt                                                                                                                                                                                            >  1:373966/1‑251 (MQ=255)                                                 aagaagCCGAGGGGGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAAAAGAAATTAAAAGGG‑TGAGTGGGATAAAAACCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTagcag                                                                                                                                                                                       <  2:263951/251‑1 (MQ=255)                                                   gaagCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAgcagca                                                                                                                                                                                     <  1:351805/251‑1 (MQ=255)                                                         gagggtgcaaaGATCGCCAAACGGCAGCACCGGATTACCCGGCACCCAAACAGCAATTAAAAGGG‑TGATTGGGATAAAACTCGGTCTTGCCGATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCTAATCAGGAACTCCTGCCCGGCGAACGTCTCCTGTAAACGTCTAGGAGTTTGCTTAGCTATGTGCTTTAAGCGTTCCTGGATGGGTTTAATCCCCCGTGGCATTAGCATCATCACTg                                                                                                                                                                               <  1:39930/242‑1 (MQ=255)                                                               ggaaaGATCCCAAAAGGGCAGCACCAGATTACCCGGCACACAAACAG‑AAATAAAATGGGTGAGTGGGAAAAAAACCAGGCTGGCCGAT‑CC‑CGCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGccgccg                                                                                                                                                                         <  2:184061/248‑1 (MQ=255)                                                                 aaaGATCGCCAACCGGCAGCCCCTGATTCCCCGGCACAAAAACAGAAATTAAAATGG‑TGAGTGGGAAAAAAATCAGGCTGCCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGtt                                                                                                                                                                       <  1:355877/250‑1 (MQ=255)                                                                                  agcaCCTGATTACCCGGCACACAAACAGCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAGCTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCa                                                                                                                                                       >  2:441778/1‑250 (MQ=255)                                                                                   gcaCCTGATTCCCCGCCACACAAACAGCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTGCCGATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAAt                                                                                                                                                     <  1:236795/251‑1 (MQ=255)                                                                                   gcaCCTGATTACCCGGCACACAAACAGCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAAt                                                                                                                                                     >  1:325364/1‑251 (MQ=255)                                                                                        tGATTACCCGGCACACAAACAGCAATTAAAATGG‑TGCGTGGGCTAACAATCAGGCTTGCCGATACAACC‑‑GCCGTATGTCCCAACAACAACGCCGGATCAAGACCAAACCGGAACTCCTGACCGCCGAACTTCGACTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATACATGATGGGTTTAATCACCCGTGGCATTAGCAGCAACACTGCCGCCGTTTTTACCGCCAGCTGCAATACAcc                                                                                                                                                >  1:440733/1‑251 (MQ=255)                                                                                                 gggacccgaatcgCAATCAAAATGG‑CGAGCGGGATAACAATCAGGCTTGCCAATACCACC‑‑GCCGTATGTCCCAGCCGCAACGCCGGATCCAGGCCAATCAGGAACTCCTGACCGCCGACCTTCGCCTGTAACCGACTAAGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGCTTAATCACCCCTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACACTTGACATc                                                                                                                                       >  2:67827/13‑251 (MQ=255)                                                                                                  gCACACAAACAGCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCg                                                                                                                                      >  1:115684/1‑251 (MQ=255)                                                                                                                 ttAAAATGG‑TGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATg                                                                                                                       >  2:140338/1‑251 (MQ=255)                                                                                                                   aaaaatgggtgttggaaaaaaaaTCAGGCTTGCC‑ACACC‑CCCTGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTAGCCGCCGAACTTCGCCTGTAAACGACTCCTAGCTTGCTTAGCGATGGGCTTTAAACCTTCCATGATGGGTTTAATCACCCGTGGCATTTGCAGCATCACTGCCGCCGTTTTTACCGCCATCTGCTATACACCTTTTATATCGTAACCCGCGAGGATg                                                                                                                       <  2:226752/232‑1 (MQ=255)                                                                                                                   aaaaagggtgggggggaaaaaaacgaGGCTTGCGGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGAGGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATg                                                                                                                       <  2:120285/225‑1 (MQ=255)                                                                                                                    aaaTGG‑TGAGTGGGATAAAAATCAGGCTTGCCGCTCCCAAC‑‑GCCGACTGTCACAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATAGTAACCCGCGAGGCTGCCg                                                                                                                    >  1:421167/1‑251 (MQ=255)                                                                                                                    aaaTGG‑TGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCg                                                                                                                    >  1:82858/1‑251 (MQ=255)                                                                                                                     aaTGG‑TGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCg                                                                                                                    >  1:235071/1‑250 (MQ=255)                                                                                                                      aTGG‑TGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGAt                                                                                                                  >  1:2727/1‑251 (MQ=255)                                                                                                                      aTGG‑TGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGAt                                                                                                                  >  1:255997/1‑251 (MQ=255)                                                                                                                       tGG‑TGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGAt                                                                                                                  >  1:412358/1‑250 (MQ=255)                                                                                                                            gAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAAc                                                                                                                                                                                                                            <  1:287389/140‑1 (MQ=255)                                                                                                                            gAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAAc                                                                                                                                                                                                                            >  2:287389/1‑140 (MQ=255)                                                                                                                                       aaaTCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCACGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAgcttgctt                                                                                                                                                                                                                                             >  2:188420/1‑112 (MQ=255)                                                                                                                                       aaaTCAGGCTTGCCGAGAC‑‑CGCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAgcttgctt                                                                                                                                                                                                                                             <  1:188420/112‑1 (MQ=255)                                                                                                                                                             accgccgactatCCCAGCACAAACGCAGGATCAAGGACAAACAGGAACTCCTGACCGCAGAACTTCGCCTTTAAACGACTACGAGCTTGCTTAGCGATAGACGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCATCTGCAATACACCTTTGACATCGTCACCCGCGAGGCTGCCAATAATCAGCCCCATCACACAACCCACGGTGACAGGCTc                                                                              >  2:285717/12‑251 (MQ=255)                                                                                                                                                              ccgccgTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCg                                                                             >  1:187671/1‑251 (MQ=255)                                                                                                                                                                            agcagcAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCGCCAAATGGACCaa                                                                >  2:339438/1‑250 (MQ=255)                                                                                                                                                                                           acAAGCCCAATCAGGAACTCCTGCCCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGTTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAGCCCACGGTGACAGGCTCGCCAAATGGACCAAATTTGCGCTGAATAtc                                                <  2:266318/250‑1 (MQ=255)                                                                                                                                                                                            cAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCGCCAAATGGACCAAATTTGCGCTGAATAtc                                                <  1:204134/250‑1 (MQ=255)                                                                                                                                                                                                   aaTCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCGCCAAATGGACCAAATTTGCGCTGAATATCGTCGgcgc                                        >  1:335058/1‑251 (MQ=255)                                                                                                                                                                                                       aggagccccttaCCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGTGTTTAATCCCCCGTGGCATTAGCAGCATCACTTCCGCCGTTTTTACCTCCAGCTGCAATACTCCTTTGACATCGTAACCCGCGAGGATGCCGATAATCTGCCCCATCACAAAACCCACGGTGTCAGGCTCGCCAAATGGACCAAATTTGCGCTGAATATCGTCGGCGCTaaa                                    <  1:434163/240‑1 (MQ=255)                                                                                                                                                                                                           ccTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCGCCAAATGGACCAAATTTGCGCTGAATATCGTCGGCGCTAAATTTa                                <  2:412358/250‑1 (MQ=255)                                                                                                                                                                                                           aCTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGAATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCGCCAAATGGACCAAATTTGCGCTGAATATCGTCGGCGCTACAttt                                 >  2:385908/1‑250 (MQ=255)                                                                                                                                                                                                             tCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCATGCTTAGCGCTGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACAGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCGCCAAATGGACCAAATTTGCGCTGAATATCGTCGGCGATAAATTTAAt                              >  1:282366/1‑251 (MQ=255)                                                                                                                                                                                                                 gACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCGCCAAATGGACCAAATTTGCGCTGAATATCGTCGGCGCTAAATTTAATTCgg                          >  1:148802/1‑251 (MQ=255)                                                                                                                                                                                                                   cggccgAACTGCGCCTGTAACCGACTACGAGCTTGCTTAGGGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCTCTGTGACAGGCTCGCCAAATGGACCAAATTTGCGCTGAATATCGTCGGCTCTAAATTTAATTCGGtt                        <  2:270233/249‑1 (MQ=255)                                                                                                                                                                                                                      ccgAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCGCCAAATGGACCAAATTTGCGCTGAATATCGTCGGCGCTAAATTTAATTCGGTTAAc                     >  1:161031/1‑251 (MQ=255)                                                                                                                                                                                                                              cGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCGCCAAATGGACCAAATTTGCGCTGAATATCGTCGGCGCTAAATTTAATTCGGTTAACGCCTggg              <  2:82858/250‑1 (MQ=255)                                                                                                                                                                                                                              cGCCTGGAACCGACTACGACCTTGCTTAGCGTGGGGCGTTAAACCAGCCATGATGGGTTTAATCACCCGTGGCATTATCAGCAGCACTGTCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGTTGACAGGCTCGCCAAATTGACCAAATTTGCGCTGAATATCGTCGGCGCTAAATTTAATTCGGTTAACGCCTGGGa             <  2:124154/251‑1 (MQ=255)                                                                                                                                                                                                                                         gACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTGTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCGCCAAATGGACCAAATTTGCGCTGAATATCGGCGGCGCTAAATTTAATTCGGTTAACGCCTGGGATTTTCTCgatg  <  1:232521/251‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                 AAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGG‑TGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCGCCAAATGGACCAAATTTGCGCTGAATATCGTCGGCGCTAAATTTAATTCGGTTAACGCCTGGGATTTTCTCGATG  >  NC_000913/2173207‑2173688

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