Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,120,989	G→A	S97S (AGC→AGT)	glcA ←	glycolate/lactate:H(+) symporter GlcA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,120,989	0	G	A	93.4%	171.0 / ‑6.3	65	S97S (AGC→AGT)	glcA	glycolate/lactate:H(+) symporter GlcA
Reads supporting (aligned to +/- strand):  ref base G (0/1);  major base A (32/25);  minor base C (1/3);  total (33/32)
Fisher's exact test for biased strand distribution p-value = 3.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCAGAATCGGCACGCCCAACGCACCAAACGCCACCGGCGCAGTATTGGCAATCAGACACAGCCCCGCCGCGTATAACGGTTTGAAGCCCAGGCCCACCAGCA‑GCGCACCGGTAATCGCCACCGGCGCACC‑AAAG‑CC‑AGCCGCTCCTTCCAGCAACGCACCAAAGGAGAAACCAAT‑CAGTAACACCTGCAAACGCTGATCGTCGGTGATGGAGATAACCGAGCTGCGGATAATATCGAACTGCCCGCTGGCAACGGTTAATTTATACAGGAACACCGCCGCGACAATAATCCACGCTATTGGCCATAATCCATAAATAAAGCCATAGCCCGCAGCAGCAAATGCCATATCAATCGGCATTTTAAAGGCGAATATTGCAATCAGGATAGATAATATAAGGGTTATTGCTCCAGCGACATGTCCTTTCAGACGTAATACCGCGAGTGCAACGAAGAAGAATATTATCGGGATCAGGGCGACCAGAGCGG  >  NC_000913/3120745‑3121230                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                  ccAGAATCGGCACGCCCAACGCACCAAACGCCACAGGCGCAGTATTGGCAATCAGACACAGCCCCGCCGCGTATAACGGTTTGAAGCCCAGGCCCACCAGCA‑GCGCACCGGTAATCGCCACCGGCGCACC‑AAAG‑CC‑AGCCGCTCCTTCCAGCAACGCACCAAAGGAGAAACCAAT‑CAGTAACACCTGCAAACGCTGATCGTCGGTGATGGAGATAACCGAGCTGCGGATAATATCGAACTGCCCACTGGc                                                                                                                                                                                                                                              >  1:575609/1‑250 (MQ=255)               cccAACGGTCCAAACGCGACCGGGGCAGTATTGCCGATCAGACCCAACCCCGCCGTGTATAACGGTTTGAAGCCCATGCCCACCAGCA‑GCGCACCGGTAATCGCCACCGGCTCACC‑AAAG‑CC‑AGCCGCTCCTTCCAGCAACGCACCAAAGGAGAAACCAAT‑CAGTAACGCCTGCAAACGCGGATCGTCGGTGATGGAGATAACCGAGCTGCGGATAATATCGAACTGCGCATTGGCAACGGTTAATTtata                                                                                                                                                                                                                               <  2:418791/251‑1 (MQ=255)               cccAACGCACCAAACGCCAACGGCGCAGTATTGGCAATCAGACACAGCCCCGCCGCGTATAACGGTTTGAAGCCCAGGCCCACCAGCA‑GCGCACCGGTAATCGCCACCGGCGCACC‑AAAG‑CC‑AGCCGCTCCTTCCAGCAACGCACCAAAGGAGAAACCAAT‑CAGTAACACCTGCAAACGCTGATCGTCGGTGATGGAGATAACCGAGCTGCGGATAATATCGAACTGCCCACTGGCAACGGTTAATTtat                                                                                                                                                                                                                                >  1:516708/1‑250 (MQ=255)                   aCGCACCAAACGCCACCGGCGCAGTATTGGCAATCAGACACAGCCCAGCCGCGTATAACGGTTTGAAGCCCAGGCCCACCAGCA‑GCGCACCGGTAATCGCCACCGGCGCACC‑AAAG‑CC‑AGCCGCTCCTTCCAGCAACGCACCAAAGGAGAAACCAAT‑CAGTAACACCTGCAAACGCTGATCGTCGGTGATGGAGATAACCGAGCTGCGGATAATATCGAACTGCCCACTGGCAACGGTTAATTTATACAgg                                                                                                                                                                                                                           <  1:318307/251‑1 (MQ=255)                     gCACCAAACGCCACCGGCGCAGTATTGGCAATCAGACACAGCCCCGCCGCGTATAACGGTTTGAAGCCCAGGCCCACCAGCA‑GCGCACCGGTAATCGCCACCGGCGCACC‑AAAG‑CC‑AGCCGCTCCTTCCAGCAACGCACCAAAGGAGAAACCAAT‑CAGTAACACCTGCAAACGCTGATCGTCGGTGATGGAGATAACCGAGCTGCGGATAATATCGAACTGCCCACTGGc                                                                                                                                                                                                                                              <  2:464287/230‑1 (MQ=255)                     gCACCAAACGCCACCGGCGCAGTATTGGCAATCAGACACAGCCCCGCCGCGTATAACGGTTTGAAGCCCAGGCCCACCAGCA‑GCGCACCGGTAATCGCCACCGGCGCACC‑AAAG‑CC‑AGCCGCTCCTTCCAGCAACGCACCAAAGGAGAAACCAAT‑CAGTAACACCTGCAAACGCTGATCGTCGGTGATGGAGATAACCGAGCTGCGGATAATATCGAACTGCCCACTGGc                                                                                                                                                                                                                                              >  1:464287/1‑230 (MQ=255)                     gCACCAAACGCCACCGGCGCAGTATTGGCAATCAGACACAGCACCGCAGCGTCTACCGGCTTGAAGCCCAGGCCAACCAGCA‑GCGCACCGGTAATCGCCACCGGCGCAAC‑AAAG‑CC‑AGCAGCTCCTTCCAGAAACGCACCAAAGGCGAAACCAAT‑CAGTAACACCTGCAAACGCTGACCGACGGTGAAGGAGATACCCGAGCTGCGGACAATATCGAACTGCCCCCCAGCAACGGTTAATTTATCCAGGa                                                                                                                                                                                                                          >  1:58422/1‑250 (MQ=255)                          aaaCGACACAGGCGAAGTATTGCGAATCAGACACAACCCCGCCGCGTATAACGGTTTGAAGCCCAGGCCCACCAGAA‑GCGCACCGGTAATCGCCACCCGCGCACC‑AACG‑AC‑AGCCGCTCCTTCCAGCAACGCACCAAAGGAGAAAAAAAT‑CAGTAACACCTGCAAACGCTGATCGTCACTGAGGGAGATAACCGAACCGCGGATAATCTCGAACTCCACACCGGCAACTGTTAATTAATACaaaaacacac                                                                                                                                                                                                                    >  2:98513/1‑242 (MQ=255)                             cccctacccgaagcccttccaatccAGACAAAGCCCCGCCGCGTATAAGGGTTTGAAGCCCAGGCCCCCCAGCG‑GGCCTCCTGTACTCGCTACCGGCGC‑CC‑AAAGGCC‑AGCCGCTTCTGCAAGCAGCGCCCTAAAGGAGAACCCAAT‑CAGTAACACCTGAAAACGTTGCTCGTCGGTGATTGAGAGAACCGAGCTGCGGATAATATCGAACTGCCCACTGCCAACGGTTAATTTATACAGGACCAccgccg                                                                                                                                                                                                                 <  2:162670/227‑1 (MQ=255)                              gCCACCGGCGCAGTATTGGCAATCAGACACAGCCCCGCCGCGTATAACGGTTTGAAGCCCAGGCCCACCAGCA‑GCGCACCGGTAATCGCCACCGGCGCACC‑AAAG‑CC‑AGCCGCTCCTTCCAGCAACGCACCAAAGGAGAAACCAAT‑CAGTAACACCTGCAAACGCTGATCGTCGGTGATGGAGATAACCGAGCTGCGGATAATATCGAACTGCCCACTGGCAACGGTTAATTTATACAGGAACACcgccgc                                                                                                                                                                                                                >  2:487574/1‑251 (MQ=255)                               ccACCGGCGCAGTATTGGCAATCAGACACAGCCCCGCCGCGTATAACGGTTTGAAGCCCAGGCCCACCAGCA‑GCGCACCGGTAATCGCCACAGGCGCACC‑AAAG‑CC‑AGCCGCTCCTTCCAGCAACGCACCAAAGGAGAAACCAAT‑CAGTAACACCTGCAAACGCTGATCGTCGGTGATGGAGATAACCGAGCTGCGGATAATATCGAACTGCCCACTGGCAACGGTTAATTTATACAGGAACACcgccgcc                                                                                                                                                                                                               >  2:218143/1‑250 (MQ=255)                                             ttGGCAATCAGACACAGCCCCGCCGCGTATAACGGTTTGAAGCCCAGGCCAACCAGCA‑GCGCACCGGTAATCGCCACCGGCGCACC‑AAAG‑CC‑AGCCGCTCCTTCCAGCAACGCACCAAAGGAGAAACCAAT‑CAGTAACACCTGCAAACGCTGATCGTCGGTGATGGAGATAACCGAGCTGCGGATAATATCGAACTGCCCACTGGCAACGGTTAATTTATACAGGAACACCGCCGCGACAATAATCCACGc                                                                                                                                                                                                 >  1:381729/1‑251 (MQ=255)                                             gtggtaacaacccaaaGCCCCACCTCGTATAACGGTTTGACGCCCAGGCCCACCACCA‑GCGCACCGGTAAGCGCCCCCGGCACACC‑ACAG‑CC‑CGCCGCTCCTTCCAGCAACTCACCAAAGGAGAAACCAAT‑CAGTACCACCTGCAAACGCTGCTCGTCGGTGATGTAGATAACCGAGCTGCGGATAATATCGACCTGCCCTCTGGCAACGGTTAATTTATACAGGAACACCGCCGCCACAATAATCCACg                                                                                                                                                                                                  <  1:305406/235‑1 (MQ=255)                                                acacccaaccccccccccGCCGCTTAAAACGGTTTGAAGCCCAGGCCCCCCAGAA‑GCGCCCCGGTAATCCCCACCGGCGCCCC‑AAAG‑CC‑AGCCGCTCCTTCCAGCAACGCACAAAAGGAAAAACCAAT‑CAGTAACCCCTGCAAACGCTGATCGTCGGTGATGAAGATAACCGAGCTGCGGATAATATCGAACTGCCCACGGGCAACGGTTAATTTATAAAGGAACACCGCCGCGACAATAATCCACGCTAt                                                                                                                                                                                              <  2:258554/237‑1 (MQ=255)                                                 cAATCAGACACAGCCCCGCCGCGTATAACGGTTTGAAGCCCAGGCCCACCAGCA‑GCGCACCGGTAATCGCCACCGGCGCACC‑AAAG‑CC‑AGCCGCTCCTTCCAGCAACGCACCAAAGGAGAAACCAAT‑CAGTAACACCTGCAAACGCTGATCGTCGGTGATGGAGATAACCGAGCTGCGGATAATATCGAACTGCCCACTGGCAACGGTTAATTTATACAGGAACACCGCCGCGACAATAATCCACGCTAtt                                                                                                                                                                                             >  2:191374/1‑251 (MQ=255)                                                  aaTCAGACACAGCCCCGCCGCGTATAACGGTTTGAAGCCCAGGCCCACCAGCA‑GCGCACCGGTAATCGCCACCGGCGCACC‑AAAG‑CC‑AGCCGCTCCTTCCAGCAACGCACCAAAGGAGAAACCAAT‑CAGTAACACCTGCAAACGCTGATCGTCGGTGATGGAGATAACCGAGCTGCGGATAATATCGAACTGCCCACTGGCAACGGTTAATTTATACAGGAACACCGCCGCGACAATAATCCACGCTATTg                                                                                                                                                                                            >  1:288099/1‑251 (MQ=255)                                                   cacAGACACACCCCCGCCGCGTATAACGGTTTGAAGCCCAGGCCCCCCAGCA‑GCCCACCGTTAATCGCCCCCGGCGCACC‑AAAG‑CC‑AGCCGCTCCTTCCAGCAACGCACCAAAGGAGAAACCAAT‑CAGTAACACCTGCAAACGCTGATCGTCGGGGATGGAGATAACCGAGCTGCGGATAATATCGAACTGCCCACTGGCAACGGTTAATTTAAACAGGAACACCGCCGCGACAATAATCCACGCTATTgg                                                                                                                                                                                           <  2:504574/249‑1 (MQ=255)                                                   aTCAGACACAGCCCCGCCGCGTATAACGGTTTGAAGCCCAGGCCCACCAGCA‑GCGCCCCGGTAATCGCCCCCGGCGCACC‑AAAG‑CC‑AGCCGCTCCTTCCAGCAACGCACCAAAGGAGAAACCAAT‑CAGTAACCCCTGCAAACGCTGATCGTCGGTGATGGAGATAACCGAGCTGCGGATAATATCGAACTGCCCACTGGCAACGGTTAATTTATACAGGAACACCGCCGCGACAATAATCCACGCTATTgg                                                                                                                                                                                           <  1:191374/251‑1 (MQ=255)                                                      aGACACAGCCCCGCCACGTATAACGGTTTGAAGCCCAGGCCCACCAGCA‑GCGCACCGGTAATCGCCAACGGCGCACC‑AAAG‑CC‑AGCCGCTCCTTCCAGCAACGCACCAAAGGAGAAACCAAT‑CAGTAACACCTGCAAACGCTGATCGTCGGGGATGGAGATAACCGAGCTGCGGATAATATCGAACTGCCCACTGGCAACGGCTAATTTATACAGGAACAccgccg                                                                                                                                                                                                                 >  1:162670/1‑226 (MQ=255)                                                       gaaacaGCCCGGCGGCGTATAACGGTTTGAAGCCCAGGCATACCAGCA‑GCCCACCGGTCACCGCCACCGGCGCACC‑AAAG‑CC‑AGCCGCTCCTTCCAGCAACGAAACACAGGAGAAACCAAT‑CAGTAACACCTGCAAACGCTGAACATCGGTGATGGAGATAACCGAGCTGCGGATAACATAGAACTGACCACTGGCAACGGTTTATCTCTACAGGAACACCGCCGCGACACTAATACACTCTATTGGNCAt                                                                                                                                                                                       >  2:605909/4‑251 (MQ=255)                                                         cccaGCCCCGCCGCGGATACCGGTTTGAAGCCAAGGCCCCCCA‑CACGCGCACCGGTAACCGCC‑CCAGCCCACCCAAAG‑CC‑ACCCGCTCCTTCCAGCAACCCACCAAAGGAGAAACCAGT‑AAGAAACACCTGAAAACGCTGATCGTCGGTGATGGAGATAACCGAGCTGCGGATAATATCGAACTGCCACCTGGCAACGTTTAATTTATACAGGACCCCCGCCGCGACAGTAATCCTCGCTATTGGCCATaa                                                                                                                                                                                     <  2:455293/249‑1 (MQ=255)                                                         cacaGCCCCGCCGCGTATAACGGTTTGAAGCCCAGGCCCACCAGCA‑GCGCACCGGTAATCGCCACCGGCGCACC‑AAAG‑CC‑AGCCGCTCCTTCCAGCAACGCACCAAAGGAGAAACCAAT‑CAGTAACACCTGCAAACGCTGATCGTCGGTGATGGAGATAACCGAGCTGCGGATAATATCGAACTGCCCACTGGCAACGGTTAATTTATACAGGAACACCGCCGCGACAATAATCCACGCTATTGGCCATa                                                                                                                                                                                      >  2:129672/1‑250 (MQ=255)                                                           caGCCCCCCCGCGTATAACGGTTTGAAGCCCAGGCCCCCCACCA‑GCGCACCGTTAACCCCCACCGGC‑CACCCAACC‑CC‑AGCCGCTCCTTCCACCAACGCACCAAAGTAGAAACCAAT‑CAGTACCTCCTGCAAACGCGGACCGTCGGTGATGGAGATAACCGAGCTTCGGATAATATCGACCTTCCCACTGCCAACGTTTAATTTATACAGGAACACCGCCGCGACAATAATCCACGCTATTGGCCATAATc                                                                                                                                                                                   <  2:205867/251‑1 (MQ=255)                                                                 cgccgcGTATAACGGTTTGAAGCCCAGGCCCAACAGCA‑GCGCACCGGAAATCGACACCGGCGCACA‑AAAG‑CC‑AGCCGCTCCCTCAAGCAACGCACCAAAGGAGAAACCAAT‑CAGTAACACCTGCAAACGCTGATCGACGGTGATGGAGATAACCGAGCTGCGGATAATATCGAACTGCCCACTGGCAACGGTTAATTTATACAGGCACAACGCCGCGACAATAATCCACGCTATTGGCCATAATCCagaa                                                                                                                                                                              >  2:326204/1‑247 (MQ=255)                                                                      cgTATAACGGTTTGAAGCCCCGCCCCACAATCA‑GCGCACCGGAAAGCGCCACGGGCGCACC‑AAGG‑CC‑AGCCGCTCCTTCCACCAACGCACCAAAGGAGAAACCAAA‑CAGTAACACCTGCAACCGCTGATCGTCGTTGATGGAGATAACCGAGCTGAGGATAATATCGACCTGCCCACTGGCAACGGTTAATTTATACAGGAACAACGCCGCGACAATAATCCACGCTATTGGCCATAATCCATAAATAAAg                                                                                                                                                                        <  2:53466/251‑1 (MQ=255)                                                                          taccggttggacgcccggtcccACCAGCA‑GCGCACGGGTAATCGCCCCCGCCGCACC‑AAAC‑CC‑AGCCGCTCCTGCCAGCAACGCACCAAAGGA‑AAACCAAAACAGTAACACCTGCAACCGCTGATCGTGTGTGATGGAGATATCCGAGCTGGGGATAATATCGAACTGCCCTCTGGCAACGGTTAATTTATACAGGAACACCGCCGCGACAATAATCCACGCTATTGGCCATAATCCATAAATAAAGCCa                                                                                                                                                                     <  2:17133/231‑1 (MQ=255)                                                                             cGGTTTGAAGCCCAGGCCCACCAGCA‑GCGCACCGGTAATCGCCACCGGCGCACC‑AAAG‑CC‑AGCCGCTCCTTCCAGCAACGCACCAAAGGAGAAACCAAT‑CAGTAACACCTGCAAACGCTGATCGTCGGTGATGGAGATAACCGAGCTGCGGATAATATCGAACTGCCCACTGGCAACGGTTAATTTATACAGGAACACCGCCGCGACAATAATCCACGCTATTGGCCATAATCCATAAATAAAGCCATAGc                                                                                                                                                                 >  2:575803/1‑251 (MQ=255)                                                                             cGGTTTGAAGCCCAGCCCCCCCAGCA‑GCGCACCGGTAATCGCCACCGGCCCACC‑AAAG‑CC‑AGCCGCTCCTTCCAGCAACGCACCAAAGGAGAACCCAAT‑CAGTAACCCCTGCAAACGCTGATCGTCGGTGATGGAGATAACCGAGCTGCGGATAATATCGAACTGCCCACTGGCAACGGTTAATTTATACAGGAACACCGCCGCGACAATAATCCACGCTATTGGCCATAATCCATAAATAAAGCCATAGc                                                                                                                                                                 <  2:448593/251‑1 (MQ=255)                                                                               gTTTGACGCCCAGGCCCACCATCC‑GCCCACCGGTAATCGCCACCGGCGCACC‑AACG‑CC‑AGCCGCTCCTTCCAGCAACGCACCAAAGGAGAAACCAAT‑CAGTAACACGTGCAAACGCTGATCGTCGGTGATGGAGATAACCGAGCTGCGGATAATATCGAACTGCCCACTGGCAACGGTTAATTTATACAGGAACACCGCCGCGACAATAATCCACGCTATTGGCCATAATCCATAAATAAAGCCATAAccc                                                                                                                                                               >  2:106722/1‑251 (MQ=255)                                                                                 ttGAAGCCCAGGCCCACCAGCA‑GCGCACCGGTAATCGCCCCCGGCGCACA‑AAAG‑CC‑AGCCGCTCCTTCCAGCAACGCACCAAAGGAGAAACCAAT‑CAGTAACACCTGCAAACGCTGATCGTCGGTGATGGAGATAACCGAGCTGCGGATAATATCGAACTGCCCACTGGCAACGGTTAATTTATACAGGAACACCGCCGCGACAATAATCCACGCTATTGGCCATAATCCATAAATAAAGCCATAGCCCg                                                                                                                                                              <  1:23099/250‑1 (MQ=255)                                                                                   gAAGCCCAGCCCCCCCAGCA‑GCCCCCCGGTAATCGCCACCGGC‑CCCC‑AAAG‑CCCAGCCGCCCTTTCCAGCAACGCACCAAAGGAGAAACCAAG‑CAGTAACCCCTGCAAACGCGGATCGCCGGTGATGGAGATAACCGAGGTGCGGATAATATCGAACTGCCCCCTGGCAACGGTTAATTTATACAGGAACACCGCCGCGACAATAATCCACGCTATTGGCCATAATCCATAAATAAAGCCATAGCCCgcag                                                                                                                                                           <  2:475791/251‑1 (MQ=255)                                                                                            gCCCACCAGCA‑GCGCACCGGTAATCGCCACCGGCGCACC‑AAAG‑CC‑AGCCGCTCCTTCCAGCAACGCACCAAAGGAGAAACCAAT‑CAGTAACACCTGCAAACGCTGATCGTCGGTGATGGAGATAACCGAGCTGCGGATAATATCGAACTGCCCACTGGCAACGGTTAATTTATACAGGAACACCGCCGCGACAATAATCCACGCTATTGGCCATAATCCATAAATAAAGCCATAGCCCGCAGCAGCAAATg                                                                                                                                                  >  2:219107/1‑251 (MQ=255)                                                                                                  cagca‑gCGCCCCGATAATCGCCACCGGCGCACC‑AAAG‑CC‑AGCCGCTCCTTCCAGCACCGCACCAAAGGAGAAACCAAT‑CAGTAACACCTGCAAACGCTGATCGTCGGTGATGGAGATAACGGAGCTGCGGATAATATCGAACTGCCCACTGGCAACGGTTAATTTATACAGGAACACCGCCGCGACAATAATCCACGCTATTGGCCATAATCCATAAATAAAGCCATAGCCCGCAGCAGCAAATGCCatat                                                                                                                                            <  2:490285/251‑1 (MQ=255)                                                                                                  cagca‑gCGCACCGGTAATCGCCACCGGCGCACC‑AAAG‑CC‑AGCCGCTCCTTCCAGCAACGCACCAAAGGAGAAACCAAT‑CAGTAACACCTGCAAACGCTGATCGTCGGTGATGGAGATAACCGAGCTGCGGATAATATCGAACTGCCCACTGGCAACGGTTAATTTATACAGGAACACCGCCGCGACAATAATCCACGCTATTGGCCATAATCCATAAATAAAGCCATAGCCCGCAGCAGCAAATGCCatat                                                                                                                                            >  2:482184/1‑251 (MQ=255)                                                                                                     ca‑gcGAACCGGTAATCGCCACCGGCGCACC‑AAAG‑CC‑AGCCGCTCCTTCCAGCAACGCACCAAAGGAGAAACCAAT‑CAGTAACACCTGCAAACGCTGATCGTCGGTGATGGAGATAACCGAGCTGCGGATAATATCGAACTGCCAACTGGCAACGGATAATTTATACAGGAACACCGCCGCGACAATAATCCACGCTATTAGCAATAATCCATAAATAAAGCCATAGCCCGCAGCAGCAAATGCCATATCa                                                                                                                                          >  2:257744/1‑250 (MQ=255)                                                                                                        cccccccGGTAAACGCCACCGGCGCACC‑AAAG‑CC‑AGCCGCTCCTTCCAGCAACGCACCAAAGGAGAAACCAAT‑CAGTAACACCTGCAAGCGCTGATCGTCGGTGATGGAGATAACCGTGCTGCGGATAATTTCGAACTGCCCACTGGCAACGGTTAATTTATACAGTAACACCGCCGCGACAATAATCCACGCTATTGGCCATAATCCATAAATAAAGCCATAGCCCGCAGCAGCAAATGCCATATCAATc                                                                                                                                       <  1:69244/246‑1 (MQ=255)                                                                                                          gcACCGGTAATCGCCCCCGGCGCACC‑AAAG‑CC‑AGCCGCTCCTTCCAGCAACGCACCAAAGGAGAAACCAAT‑CAGTAACACCTGCAAACGCTGATCGTCGGTGATGGAGATAACCGAGCTGCGGATAATATCGAACGGCCCACTGGCAACGGTTAATTTATACAGGAACACCGCCGCGACAATAATCCACGCTATTGGCCATAATCCATAAATAAAGCCATAGCCCGCGGCAGCAAATGCCATATCAATCgg                                                                                                                                     <  1:309555/251‑1 (MQ=255)                                                                                                               ggTAATCGCCACCGGCGCCCC‑AAAG‑CC‑AGCCGCTCCTTCAAGCAACGCACCAAAGGAGAAACCAAT‑CAGTACCACCTGCAAACGCTGATCGTCGGTGATGGAGATAACCGAGCTGCGGATAATATCGAACTGCCCACTGGCAACGGTTAATTTATACAGGAACACCGCCGCGACAATAATCCACGCTATTGGCCATAATCCATAAATAAAGCCATAGCCCGCAGCAGCAAATGCCATATCAATCGGCAttt                                                                                                                                <  1:373190/251‑1 (MQ=255)                                                                                                               ggTAATCGCCACCGGCGCACC‑AAAG‑CC‑AGCCGCTCCTTCCAGCAACGCACCAAAGGAGAAACCAAT‑CAGTAACACCTGCAAACGCTGATCGTCGGTGATGGAGATAACCGAGCTGCGGATAATATCGAACTGCCCACTGGCAACGGTTAATTTATACAGGAACACCGCCGCGACAATAATCCACGCTATTGGCCATAATCCATAAATAAAGCCATAGCCCGCAGCAGCAAATGCCATATCAATCGGCAttt                                                                                                                                <  1:487574/251‑1 (MQ=255)                                                                                                                 tAATCGCCACCGGCGCACC‑AAAG‑CC‑AGCCGCTCCTTCCAGCAACGCACCAAAGGAGAAACAAAT‑CAGTAACACCTGCAAACGCTGATCGTCGGTGATGGAAATAACCGAGCTGCGGATAATATCGAACTGCCCACTGACAACGGTTAATTTATACAGGAACACCGCCGCGACAATAATCCACGATATTGGCCATAATCCATAAATAAAGCCATATCCCGCAGCAGCAAATGCAATATCAATCGGCAtttt                                                                                                                               >  2:355017/1‑250 (MQ=255)                                                                                                                      cccACCGGCGCCCC‑AAAG‑CC‑AGCCGCCCCTGCCAGCACCGCACCAAAGGAGAAACCAAT‑CAGTAACACCTGCAAACGCTGATCGTCGGTGATGGATATAACCGAGCTGCGGATAATATCGAACCGCCCACTGGCAACGGTTAATTTAGACAGGAACACCGCCGCGACAATAATCCACGCTATTGGCCATAATCCATAAATAAAGCCATAGCCCGCAGCAGCAAATGCCATATCAATCGGCATTTTAAAGGc                                                                                                                         <  2:645320/250‑1 (MQ=255)                                                                                                                       ccACCGGCGCCCC‑AAAG‑CC‑AGCCGCTCCTTCCAGCAACTCACAAAAGGAGAAACCAAT‑CAGTACCCCCTTCAAACGCTGATCGTCGGTGAGGGAGAAAACCGATCTGCGGATAATAGCGAACTGCCCACTGGCAACGGTTAATTTATACAGGAACACCGCCGCGACAATAATCCACGCTATTGGCCATAATCCATAAATAAAGCCATAGCCCGCAGCAGCAAATGCCATATCAATCGGCATTTTAAAGGCg                                                                                                                        <  1:456030/251‑1 (MQ=255)                                                                                                                       ccACCGGCGCACC‑AAAG‑CC‑AGCCGCTCCTTCCAGCAACGCACCAAAGGAGAAACCAAT‑CAGTAACACCTGCAAACGCTGATCGTCGGTGATGGAGATAACCGAGCTGCGGATAATATCGAACTGCCCACTGGCAACGGTTAATTTATACAGGAACACCGCCGCGACAATAATCCACGCTATTGGCCATAATCCATAAATAAAGCCATAGCCCGCAGCAGCAAATGCCATATCAATCGGCATTTTAAAGGCg                                                                                                                        >  2:489397/1‑251 (MQ=255)                                                                                                                         aCCGGCGCACC‑AAAG‑CC‑AGCCGATCCTTCCAGCAACGCACCACAGGAGAAACCAAT‑CAGTAACACCTGCAAACGCTGATCGTCGGTGATGGAGATAACCGAGCTGCGGCTAATATCGAACTGCCCACTGGCAACGGTTAATTTATACAGGAACACCGCCGCGCCAATAATCCACGCTATTGGCCATAATCCATAAATAAAGCCATAGCCCGCAGCAGCAAATGCCATATCAATCGGCATTTTAAAGGCGa                                                                                                                       >  2:288923/1‑250 (MQ=255)                                                                                                                           cGGCGCACC‑AAAG‑CC‑AGCCGCTCCTTCCAGCAACGCACCAAAGGAGAAACCAAT‑CAGTAACACCTGCAAACGCTGATCGTCGGTGATGGAGATAACCGAGCTGCGGATAATATCGAACTGCCCACTGGCAACGGTTAATTTATACAGGAACACCGCCGCGACAATAATCCACGCTATTGGCCATAATCCATAAATAAAGCCATAGCCCGCAGCAGCAAATGCCATATCAATCGGCATTTTAAAGGCGAat                                                                                                                     >  2:409463/1‑250 (MQ=255)                                                                                                                                  cacaaaaccc‑agccGCTCCTTCCAGCAACGCACAAAAGGAGAAACCAAT‑CAGAAACCCCTGCAAACGCTGATCGTCGGTGATGGAGATAACCGAGCGGCGGATAATATCGAACTGCCCACTGGCAACGGTTAATTTATACAGGAACACCGCCGCGACAATAATCCACGCTATTGGCCATAATCCATAAATAAAGCCATAGCCCGCAGCAGCAAATGCCATATCAATCGGCATTTTAAAGGCGAATATTGCa                                                                                                               <  1:345776/243‑1 (MQ=255)                                                                                                                                      aaG‑CC‑AGCCGCTCCTTCCAGCAACGCACCAAAGGAGAAACCAAT‑CAGTAACACCTGCAAACGCTGATCGACGGTGATGGAGATAACCGAGCTGCGGATAATATAGAACTGCCCACTGGCAACGGTTAATTTATACAGGAACACCGCAGCGACAATAATCCACGCTATTGGCCATAATCCATAAATAAAGCCATAGCCCGCAGCAGCAAATGCCATATCAATCGGCATTTTAAAGGCGAATATTGCAANc                                                                                                            >  2:601597/1‑249 (MQ=255)                                                                                                                                           c‑agccGCTCCTCCCAGCAACGCCCCAAAGGAGAAACCAAT‑AAGTAACACCTGCAAACGCTGATCGTCGGTGATGGAGATAACCGAGCTGCGGATAATCTCGAACTGCCCTCTGGCAACGTTTAATTTATTCAGGAACACCGCCGCGACAATAATCCACGCTATTGGCCATAATCCATAAATAAAGCCATAGCCCGCAGCAGCAAAGGCGATAGCAATCGGCATTTTAAAGGCGAATATTGCAATCAGgata                                                                                                      <  2:88248/251‑1 (MQ=255)                                                                                                                                                cGCTCCTTCCAGCAACGCACCAAAGGAGAAACCAAT‑CAGTAACACCTGCAAACGCTGATCGTCGGTGATGGAGATAACCGAGCTGCGGATAATATCGAACTGCCCACTGGCAACGGTTAATTTATACAGGAACACCGCCGCGACAATAATCCACGCTATTGGCCATAATCCATAAATAAAGCCATAGCCCGCAGCAGCAAATGCCATATCAATCGGCATTTTAAAGGCGAATATTGCAATCAGgatagat                                                                                                   >  2:30892/1‑250 (MQ=255)                                                                                                                                                cGCTCCTTCCAGCAACGCACCAAAGGAGAAACCAAT‑CAGTAACACCTGCAAACGCTGATCGTCGGTGATGGAGATAACCGAGCTGCGGATAATATCGAACTGCCCACTGGCAACGGTTAATTTATACAGGAACACCGCCGCGACAATAATCCACGCTATTGGCCATAATCCATAAATAAAGCCATAGCCCGCAGCAGCAAATGCCATATCAATCGGCATTTTAAAGGCGAATATTGCAATCAGgatagat                                                                                                   >  2:314210/1‑250 (MQ=255)                                                                                                                                                   cccTTCCAGCAACGCACCAAAGGAAAAACCAAT‑CAGTAACACCTGCAAACGCTGATCGTCGGTGATGGAGATAACCGAGCTGCGGATAATATCGAACTGCCCACTGGCAACGGTTAATTTATACAGGAACACCGCCGCGACAATAATCCACGCTATTGGCCATAATCCATAAATAAAGCCATAGCCCGCAGCAGCAAATGCCATATCAATCGGCATTTTAAAGGCGAATATTGCAATCAGGATAGataata                                                                                               <  1:575803/250‑1 (MQ=255)                                                                                                                                                         caccAACGCACCAAAGGTGAAACCAAT‑CATTAACACCTGCAAACGCTGAGCATCGGTTATGGAGAAAACCGGGCGCCGGATAATATCGAACTGCCCCCTGGCAACTGTTAATTTAAACAGAAACACCGCCGCGACAATAATCCACGCTATTGGCCATAATCCATAAATAAAGCCATAGCCCGCAGCAGCAAATGCCAGATCAATCGGCATTTTAAAGGCGAATATTGCAATCATGATAGATAATAGAAggg                                                                                         <  2:274000/248‑1 (MQ=255)                                                                                                                                                         cAGCAACGCACCAAAGGAGAAACCAAT‑CAGTAACACCTGCAAACGCTGATCGTCGGTGATGGAGATAACCGAGCTGCGGATAATATCGAACTGCCCACTGGCAACGGTTAATTTATACAGGAACACCGCCGCGACAATAATCCACGCTATTGGCCATAATCCata                                                                                                                                                                               >  2:334267/1‑165 (MQ=255)                                                                                                                                                         cAGCAACGCACCAAAGGAGAAACCAAT‑CAGTAACACCTGCAAACGCTGATCGTCGGTGATGGAGATAACCGAGCTGCGGATAATATCGAACTGCCCACTGGCAACGGTTAATTTATACAGGAACACCGCCGCGACAATAATCCACGCTATTGGCCATAATCCata                                                                                                                                                                               <  1:334267/165‑1 (MQ=255)                                                                                                                                                             aacgcccaaaaGGAGAAACCAAT‑CAGTAACACCTGCAAACATTAATCGTCGGTGATGGAGTTAACCAAGCTGCGAATAATATCTAACTGCCCACTGCCAACGATTAATTTATACAGGAACACCGCCGCGACAATAGTCCACGCTATTGGCCATAATCCATAAATAAAGCCATAGCCCGCAGCAGCAAATGCCATATCAATCGGCATTTTAAAGGCGTATATTGCAATCAGGATAGATAATATAAGGGTTAt                                                                                     <  1:249229/243‑1 (MQ=255)                                                                                                                                                                      aaGGAGAAACCAAT‑CAGTAACACCTGCAAACGCTGATCGTCGGTGATGGAGATAACCGAGCTGCGGATAATATCGAACTGCCCACTGGCAACGGTTAATTTATACAGGAACACCGCCGCGACAATAATCCACGCTATTGGCCATAATCCATAAATAAAGCCATAGCCCGCAGCAGCAAATGCCATATCAATCGGCATTTTAAAGGCGATTATTGCAATCAGGATAGATAATATAAGGGTTATTGCTCCAGc                                                                            >  1:459826/1‑251 (MQ=255)                                                                                                                                                                                                 aaaacgcgaacaggcggggagggAGATAACCGAGGTGCGGATAATATCAACCTGCCAACTGGCAACGGTTAATTTATACAGGAACACCGCCGCGACAATAATCCACGCTATGGGCCATAATCCATAAATAAAGCCATAGCCCGCAGCAGCAAATGCCATATCAATCGGCATTTTAAAGGCGAATATTGCAATCAGGATAGATAATATAAGGGTTATTGCTCCAGCGACATGTCCTTTCAGACGTAATACcg                                                  <  2:141302/230‑1 (MQ=255)                                                                                                                                                                                                     cGCTGACCGACGGCGATGGAGATAAACGAGCTGCGGCTAATATCGAACTGACCACTGCCAACGGTTAATTTATACAGAAACACCGCCGAGACCCTAATCCACGCTATTGGCCATAATCCATAAATAAAGCAATGGCCGGCAGCATAAAATGCCATATCAATCGGCATTTTAAACGCGAATATTGCAATCAGAATAGATAATATAAGGGTTATTGCTCCAGCGACATGTCCTTTCAGCGGTAATACcgcg                                                >  2:627324/1‑249 (MQ=255)                                                                                                                                                                                                             ttcgGTAATGGAGATAACCGAGCTGCGGATAATATCAACCTGCCCACTGGCAACGGTTAATTTATACAGAAACACCGCCGCGACAATAATCCACGCTATTGGCCATAATCCATAAATAAAGCCATAGCCCGCAGCAGCAAATGCCATATCAATCGGCATTTTAAAGGCGAATATTGCAATCAGGATAGATAATATAAGGGTTATTGCTCCAGCGACATGTCCTTTCAGACGTAATACcgcg                                                <  2:615502/240‑1 (MQ=255)                                                                                                                                                                                                             gtcgGTGATGGAGATAACCGAGCTGCGGATAATATCGAACTGCCCACTGGCAACGGTTAATTTATACAGGAACACCGCCGCGACAATAATCCACGCTATTGGCCATAATCCATAAATAAAGCCATAGCCCGCAGCAGCAAATGCCATATCAATCGGCATTTTAAAGGCGAATATTGCAATCAGGATAGATAATATAAGGGTTATTGCTCCAGCGACATGTCCTTTCAGACGTAATACcgcg                                                >  1:615502/1‑241 (MQ=255)                                                                                                                                                                                                                            aaCCGAGCTGCGGATAATATCGAACTGCCCACTGGCAACGGTTAATTTATACAGGAACACCGCCGCGACAATAATCCACGCTATTGGCCATAATCCATAAATAAAGCCATAGCCCGCAGCAGCAAATGCCATATCAATCGGCATTTTaaa                                                                                                                            >  1:429977/1‑150 (MQ=255)                                                                                                                                                                                                                            aaCCGAGCTGCGGATAATACCGAACTGCCCACGGGCAACGGTTAATTTATACAGGAACACCGCCGCGACAATAATCCACGCTATTGGCCATAATCCATAAATAAAGCCATAGCCCGCAGCAGCAAATGCCATATCAATCGGCATTTTaaa                                                                                                                            <  2:429977/150‑1 (MQ=255)                                                                                                                                                                                                                              ccGAGCGGCGGATAATCTGGACCTGCGCGGGGGCAACGTTTAATTTCAACAAGAACTCCGCCGCGACAATAATCCACGCGATTGGCCATAAGCCATAAATAAAGCCATAGCTCGCAGCAGCAAATGCCATATCAATCGGCATTTTAAAGGCGAATATTGCAATCAGTAGAGATAATATAAGGGTTATTGCTCCAGGGACATATCCTTTCAGACTTAATACCGCGATTGCAACGAAGAAGTATATTATCggg                     <  2:640520/251‑1 (MQ=255)                                                                                                                                                                                                                                            atatCGAACTGCCCACTGGCAACGGTTAATTTATACAGGAACACCGCCGCGACAATAATCCACGCTATTGGCCATAATCCATAAATAAAGCCATAGCCCGCAGCAGCAAATGCCATATCAATCGGCATTTTAAAGGCGAATATTGCAATCAGGATAGATAATATAAGGGTTATTGCTCCAGCGACATGTCCTTTCAGACGTAATACCGCGAGTGCAACGAAGAAGAATATTATCGGGATCAGGGCGAcc         >  1:444396/1‑249 (MQ=255)                                                                                                                                                                                                                                                 gAACTGCCCACTGGCAACGGTTAATTTATACAGGAACACCGCCGCGACAATAATCCACGCTATTGGCCATAATCCATAAATAAAGCCATAGCCCGCAGCAGCAAATGCCATATCAATCGGCATTTTAAAGGCGAATATTGCAATCAGGATAGATAATATAAGGGTTATTGCTCCAGCGACATGTCCTTTCAGACGTAATACCGCGAGTGCAACGAAGAAGAATATTATCGGGATCAGGGCGACCAGAGCgg  >  2:109082/1‑251 (MQ=255)                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                  CCAGAATCGGCACGCCCAACGCACCAAACGCCACCGGCGCAGTATTGGCAATCAGACACAGCCCCGCCGCGTATAACGGTTTGAAGCCCAGGCCCACCAGCA‑GCGCACCGGTAATCGCCACCGGCGCACC‑AAAG‑CC‑AGCCGCTCCTTCCAGCAACGCACCAAAGGAGAAACCAAT‑CAGTAACACCTGCAAACGCTGATCGTCGGTGATGGAGATAACCGAGCTGCGGATAATATCGAACTGCCCGCTGGCAACGGTTAATTTATACAGGAACACCGCCGCGACAATAATCCACGCTATTGGCCATAATCCATAAATAAAGCCATAGCCCGCAGCAGCAAATGCCATATCAATCGGCATTTTAAAGGCGAATATTGCAATCAGGATAGATAATATAAGGGTTATTGCTCCAGCGACATGTCCTTTCAGACGTAATACCGCGAGTGCAACGAAGAAGAATATTATCGGGATCAGGGCGACCAGAGCGG  >  NC_000913/3120745‑3121230

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