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%	145.5 / NA	48	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 (28/19);  minor base A (1/0);  total (29/19)
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.99e-01

AGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAG‑CAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCAC‑GGTGACAGGCTCGCCAAATGGACCAAATTTGCGCTGAATATCGTCGGCGCTAAATTTAATTCGGTTAACGCCTGGGATTTTCTCGATGATAGC  >  NC_000913/2173202‑2173693                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       agaTTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGCTACCACC‑‑GCCGTATGTCCCAG‑CAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAgcttgctt                                                                                                                                                                                                                                                   >  1:630346/1‑250 (MQ=255)    ttAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAG‑CAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCg                                                                                                                                                                                                                                               >  1:396484/1‑251 (MQ=255)         gggcgaaaCAGATTTCCACGAGGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAG‑CAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATggg                                                                                                                                                                                                                                          <  1:17920/245‑1 (MQ=255)                  aGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACAACC‑‑ACCGTATGTCCCAG‑CAGAAACGCCGGATCAAGGCCAATCAGGAACACCTGACCGCCGAACTTCGCCTGTAAACGCCTACTAACATGCTTAACGATGGGCGTTAAAcc                                                                                                                                                                                                                                 >  2:374757/1‑251 (MQ=255)                        ccACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAG‑CAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTACACGACTACGAGCTTGCTTAGCGATGGGCGTTAAAccatccat                                                                                                                                                                                                                           >  1:124139/1‑251 (MQ=255)                                 aCGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAG‑CAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGttt                                                                                                                                                                                                                  >  1:567982/1‑251 (MQ=255)                                   gccGACCGCCACCGCCACGAAGAAGCCGAGGGGGGCAAGATCGCCAAACGGTAGCACCTGATTACCCGGCACAAAAACAGAAATTAAAATGGTGAGTGGGATAAAAATCAGCCTTGCCGATCCCACC‑‑GCCGTATGTCCCAG‑CAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTaa                                                                                                                                                                                                                <  1:524881/249‑1 (MQ=255)                                        ccGCCACCGCCCCAAAGAAGCCGATGGTGGAAAGAACGCCAAACGGAGGCCCCTGATTCCCCGGCACACAAACAGCAATTAAAATGGGGAGTGGGATAAAAACCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAG‑CAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCAcc                                                                                                                                                                                                           <  1:320052/251‑1 (MQ=255)                                                      aagaagCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTGCCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAG‑CAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTagc                                                                                                                                                                                               >  1:374425/1‑249 (MQ=255)                                                       agaagCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGGGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAG‑CAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGACCTTCGCCTGTAAACGACTACGAGCTTCCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTagaagc                                                                                                                                                                                            >  2:567921/1‑251 (MQ=255)                                                            ccGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAG‑CAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCAc                                                                                                                                                                                       >  1:553996/1‑251 (MQ=255)                                                                tggtggCAAGATCGCCAAACGGCAGCCCCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGGTACCACC‑‑GCCGTATGTCCCAG‑CAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGACCTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTgcc                                                                                                                                                                                   >  1:536743/1‑251 (MQ=255)                                                                   gggaaaGATCGCCAAACGGCAGCCCCTGATTCCCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAG‑CAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTgccgcc                                                                                                                                                                                <  1:219415/247‑1 (MQ=255)                                                                      cAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAG‑CAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGtt                                                                                                                                                                             >  1:527392/1‑251 (MQ=255)                                                                                     gcacaccttgtttccccGGCACACAAACACAAATAAAAAGGGTGATGGGGAAAAAAATCAGGCTGTCGGATACCACC‑‑GCCGTATGTCCCAG‑CAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTGAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGc                                                                                                                                                              <  1:344005/237‑1 (MQ=255)                                                                                                     cGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAG‑CAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACAt                                                                                                                                              >  1:267125/1‑251 (MQ=255)                                                                                                     cGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTACCGATACCAAC‑‑GAAGTATGTACCAG‑CAGCAACGCCGGATCAAGGCCCATCAGGCACTCCTAACCGCCGAACTTCGCCTGTAAACGACTAAGCTCTTACTAAGCGATGGGCTTTACACCATCCATGATGGGTTTAATCAACCGTGTCATTAGCATCATCTCTGCCGCCGTTTATACCGCCAGCCGCAATACAACtttaaccc                                                                                                                                              >  2:213239/1‑246 (MQ=255)                                                                                                       gCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAG‑CAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCg                                                                                                                                            >  1:49253/1‑251 (MQ=255)                                                                                                            aaaccaccAATAAAAAGGGTGAGTGGAATAAAAACCAGGCTTGCCAATCCCACC‑‑GCCGTATGTCCCAG‑CAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTTTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAAcc                                                                                                                                       <  1:566180/244‑1 (MQ=255)                                                                                                             aaaCAGCAATTACAATGGTGAGAGGGATAAAAATCACGCTTGCCGCTACCACC‑‑GCCGTATGTCCCAG‑CAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTCCGAGCTTGCTCAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATAAGCAGCATCACTCCCGCCGTTTTTACAGCCAGCTGCCATACACCTTTGACATCGTAAccc                                                                                                                                      >  2:168179/1‑251 (MQ=255)                                                                                                                cAGCAATTACAATGGTGAGTGGGATCCAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAG‑CAGCAAAGCCGGATCAAGGACAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGCTGGTCGTTAAACCATCCATGACGGGTTTAATCACCCTGGGCATTAGCCGAATCACTGCCGCCGTATTTACCGCAAACTGCCATACACATTTGACATCGTAACAcgcg                                                                                                                                   >  2:656878/1‑251 (MQ=255)                                                                                                                   tacttAAAATGGCGAGTGGGATAAAAATCAGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAG‑CAGAAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTATACCATCCATGCTGGGTTTAATCACCCGCGGCATTAGCAGCATCCCTGCCGCCGTTTTTCCCGCCAGCTGCAATACACCTTTGACATCGTAACCcgcgcgg                                                                                                                                >  2:27663/4‑248 (MQ=255)                                                                                                                        aaaaGGGTGAGTGGAATAAAAACCAGGCTGGCCGATCCCACC‑‑GCCGTATGTCCCAG‑CAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGcc                                                                                                                           <  1:439507/251‑1 (MQ=255)                                                                                                                                gAGTGGTATAAAAATCAGGCTTGCCGATAC‑‑CCCCGCCGTATGTCCCAG‑CAGCAACGCCGGATCAAGGCCAATCAGGACCCCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCa                                                                                                                   <  1:12577/251‑1 (MQ=255)                                                                                                                                        aaaaaacaggggtgcgcacaaccaccCC‑CCGTTTGTCCCAGACA‑CACCGCCGGAGTAAGGCCAATAAGGAACCCCGGACCTCCGAACTTCGCCTTTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGTTGTGTTTAATCACCAGGTGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTACCCCGCGAGGATGCCGATAATCAGCCCCAt                                                                                                            <  2:195678/231‑1 (MQ=255)                                                                                                                                             aTCAGGCTTGCCGCTA‑‑ACACCGCCGAATGTCCCCA‑CAGCAACGCCGGATCAAGGCCCATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAAAGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCCTCCATGATGGGTTTAATCACCCGTGGCATCAGCAGCATCCCAGCCGCCGTTTTTACCGCCAGCTGCAATAACCCTTTGACATCGTAACCCGCGCAGATGCCGCTAATCAGCCCACTCACaa                                                                                                       >  2:310437/1‑250 (MQ=255)                                                                                                                                              tccggCTCGCCGATACCACC‑‑GCCGTATGTCCCAG‑CAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGCTGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCCGCTGCAATACACCTTTGACATCGTACCCCGCGAGGATGCCGATAATAAGCCCCCTCACaaa                                                                                                      >  2:100421/4‑250 (MQ=255)                                                                                                                                                aGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAG‑CAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAAcc                                                                                                   >  1:367696/1‑251 (MQ=255)                                                                                                                                                       gcccagtacagccgccgTATGCCCCAG‑AATCACCGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACCTCGCCTGTAAACGACTACGATCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGTTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTACCCCGCGAGGATGCCGATAATCAGCCCCATCTCAAAACCCAc                                                                                                <  1:374757/238‑1 (MQ=255)                                                                                                                                                            ataccaccgccgTATGTCCCAG‑CAGCAACGCCGGATCAAGGCCAATCAGGCACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAAccatccat                                                                                                                                                                                                                           >  2:271328/7‑121 (MQ=255)                                                                                                                                                            ataccaccgccgTATGTCCCAG‑CAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAAccatccat                                                                                                                                                                                                                           <  1:271328/115‑1 (MQ=255)                                                                                                                                                              accaccgccgTATGTCCCAG‑CAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCAC‑GGTGACAgg                                                                                      >  1:501981/5‑251 (MQ=255)                                                                                                                                                                        tATGTCCCAG‑CAGAAACGCCGGATCAAGGCCCATCAGGAACTCCTGAACGCCGAACTTCGCCTGTAAACGAATACGAGCTTGCTTAGCGATGGGCGTTAAACCAACCATGATGGGTTTACTCACACGTGGCATTAGCAGCCTCCCTGCCGCCGTTTTTAACGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGGGGATGCCGATAATCAGCCACAT‑AAAAAACCCACCGGAGACAGGCCAGCCaaa                                                                             >  2:663418/1‑250 (MQ=255)                                                                                                                                                                         aTGTCCCAG‑CAGCAACGCCGGATCAAGGCCAATCAGTAACTCCTGACCGCCGAACTTCGCATGTAAACGACTACGAGCTTGCTTAGCGATGAGAGTTAAACCATCCATGATGGGTTTAATCACCCGAGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGAAACCCGCGAGGATGCCGAAAATCAGCCCCATCACAAAACCCAC‑GGTGACAGGACCGCCAAATg                                                                           >  2:493645/1‑251 (MQ=255)                                                                                                                                                                             cccAG‑CACCAACGCCGAAACAGGACCAACCAGAAACCCCTGCCCGCCGACCTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAACCCATCCATGATGGGTGTAATCACCCGTGGCATTAGCAGCACCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCAC‑GGTGACAGGCTCGCCAAATGGAcc                                                                       <  2:659475/251‑1 (MQ=255)                                                                                                                                                                              ccAG‑CAGCAACGCCGTATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCGGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCCCCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCAC‑GGTGACAGGCTCGCCAAATGGACCa                                                                      <  1:656878/251‑1 (MQ=255)                                                                                                                                                                                         cGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCAC‑GGTGACAGGCTCGCCAAATGGACCAAATTTGCGCt                                                            >  1:37404/1‑251 (MQ=255)                                                                                                                                                                                                 aaaGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCAC‑GGTGACAGGCTCGCCAAATGGACCAAATTTGCGCTGAATAtc                                                     <  1:100421/249‑1 (MQ=255)                                                                                                                                                                                                  aaGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCAC‑GGTGACAGGCTCGCCAAATGGACCAAATTTGCGCTGAATAtcgt                                                   >  1:477765/1‑251 (MQ=255)                                                                                                                                                                                                    gcccattaagaaaccccTGACCGCCGAACCTCGCCTGTAAACGACTCCGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCCCAAAACCCAC‑GGTGACAGGCTCGCCAAATGGACCAAATTTGCGCTGAATAtcgtcg                                                 <  2:545653/236‑1 (MQ=255)                                                                                                                                                                                                      ccAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCAC‑GGTGACAGGCTCGCCAAATGGACCAAATTTGCGCTGAATATCGTCGgc                                               <  1:567921/251‑1 (MQ=255)                                                                                                                                                                                                      caaaccgggaccTCCGGACCTCCGAACTTCCCCTGTAACCGACTACAAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGAGGGGTTTAACCACCCGTGGCATTAGCAGCATCCCTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCAC‑GGTGACAGGCTCGCCAAATGGACCAAATTTGCGCTGAATATCGTCGgc                                               <  2:553996/240‑1 (MQ=255)                                                                                                                                                                                                                  tCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCCGCCCCATCACAAAACCCAC‑GGTGACAGGCTCGCCAAATGGACCAAATTTGCGCTGAATATCGGCGGCGCTAAATTTAAt                                   >  2:223399/1‑251 (MQ=255)                                                                                                                                                                                                                     tGACCGCCGAACTTCGCCTGTAAACGACTACGAGCTTGCTTAGCGATGGGCGTTCAACCATCCATGATGGGATTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACAGCCAGCTGCAATACAACTTTGAAATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCAC‑GGTGACAGGCTCGCCAAATGGACCAACTTTGCGCTGAATCTCGTCGGCGCTAAATTTAATTCg                                >  2:427942/1‑251 (MQ=255)                                                                                                                                                                                                                              caagtCGCCTGTAACCGACTGCGAGCTTGCCTAGCGCTGGGCCTTAAACCATCCCTGATGGGGTTAATCACCCGTGGCATTAGCCGCATCACTACCGCCGTTTTTACCGCCAGCTGCAATCAACCTTCGACAACGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACCAAACCCAC‑GGTGAAAAGCTCACCCACTGGACCAAATTTGCGCTGACTATAGTCGGCGCATACTTTAATTCgaataccgcc                       >  2:74333/5‑242 (MQ=255)                                                                                                                                                                                                                                  tCGCCTGTAAACGCCTACGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTCCCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCAC‑GGTGACAGGCTCGCCAAATGTACCAAATTTGCGCTGAATATCGTCGGCGCTAAATTTAATTCGGTTAACGCCTgg                    <  1:49128/250‑1 (MQ=255)                                                                                                                                                                                                                                         tAAACGACTACGAGATTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGTCCCATCACAAAACCCAC‑GGTGACAGGCTCGCCAAATGGACCAAATTTGCGCTGACTATCGTCGGCGCTAAATTTAATTAGATGAACGCCTGGGCTTTtct            >  2:261230/1‑251 (MQ=255)                                                                                                                                                                                                                                                  ccGAGCTTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGGGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGACAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCAC‑GGTGACAGGCTCGCCAAATGGACCAAATTTGCGCTGAATATCGTCGGCGCTAAATTTAATTCGGTTAACGCCTGGGATTTTCTCGATGATAg   <  1:105954/250‑1 (MQ=255)                                                                                                                                                                                                                                                   cGAGCTTGCTTAGCGATGGGAGTTAAACCATCCATGATGGGTTTAATCACACGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGAAATACACCTTTGACATCGTAACCCGCGCGGATGCCGATAATCAGCCCCATAACAAAACCCAC‑GGTGACAGGCCAGCCAAATGGACCAAATTTGCGCTGAATATCGTCGGCGCCAAATTTAATTCGGTTAACGCCTGGGATTTTATCGATGGTAgg  >  2:580273/1‑250 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       AGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAG‑CAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCAC‑GGTGACAGGCTCGCCAAATGGACCAAATTTGCGCTGAATATCGTCGGCGCTAAATTTAATTCGGTTAACGCCTGGGATTTTCTCGATGATAGC  >  NC_000913/2173202‑2173693

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