Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,173,363	Δ2 bp	intergenic (‑490/+918)	gatD ← / ← gatB	galactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,173,363	0	C	.	89.4%	164.7 / 17.6	47	intergenic (‑490/+919)	gatD/gatB	galactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand):  ref base C (0/5);  new base . (29/13);  total (29/18)
Fisher's exact test for biased strand distribution p-value = 5.59e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.20e-01
*	NC_000913	2,173,364	0	C	.	89.4%	169.1 / 13.1	47	intergenic (‑491/+918)	gatD/gatB	galactitol‑1‑phosphate 5‑dehydrogenase/galactitol‑specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand):  ref base C (0/5);  new base . (29/13);  total (29/18)
Fisher's exact test for biased strand distribution p-value = 5.59e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCA‑TCGCCACGAAGAAGCCGA‑TGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCG  >  NC_000913/2173120‑2173613                                                                                                                                                                                                                                                       ||                                                                                                                                                                                                                                                          ccAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCA‑TCGCCACGAAGAAGCCGA‑TGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCCGATACCACC‑‑GCCGTATg                                                                                                                                                                                                                                                   >  2:36142/1‑251 (MQ=255) ccAGCATTAGCCGCCAGTGGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCA‑TCGCCACGAAGAAGCCGA‑TGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCCGATACCACC‑‑GCCGTATg                                                                                                                                                                                                                                                   <  1:188041/251‑1 (MQ=255)  cAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCA‑TCGCCACGAAGAAGCCGA‑TGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCCGATACCACC‑‑GCCGTATGt                                                                                                                                                                                                                                                  >  1:277117/1‑251 (MQ=255)    ccatgagcggccggttggtggggCAACCCAATAGTTTGCGCCGCGATCCACAGGGTGATCCTCATAATGACACCACCCGAGATTAGGGTGCGGAACAGATTCCCCCGATGCCCGGCGCCCGCCA‑TCGCCACGAAGAAGCCGA‑TGGTGGCAAGATCGCCAAACGGCAGCGCCTGATTACCCGGCACACAAACAGCAATTAAAATGGTAAGTGGGATAAAAATCAGG‑CTTCCGGATACCACC‑‑GCCGTATGTcc                                                                                                                                                                                                                                                <  2:456280/229‑1 (MQ=255)      aTTAGCCGCCAGTTGGGTGTGCAAACCAATATTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCA‑TCGCCACGAAGAAGCCGA‑TGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCCGATACCACC‑‑GCCGTATGTCCca                                                                                                                                                                                                                                              <  1:470288/251‑1 (MQ=255)       ttAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGCTGCTCATAATGATGACACCCGAGATTAAGGAGCGGAACAGATTTCCACGATGCACGGCGACCGCCA‑TCGCCACGAAGAAGCCGA‑TGGAGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAACCAGCAATTAAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCCGATACCACC‑‑GCCGTATGTccc                                                                                                                                                                                                                                               >  2:427656/1‑249 (MQ=255)                aGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCA‑TCGCCACGAAGAAGCCGA‑TGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGc                                                                                                                                                                                                                                    >  2:343886/1‑251 (MQ=255)                                cAATAGTTTGCGTCGCGCTCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCA‑TCGCCACGAAGAAGCCGA‑TGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATc                                                                                                                                                                                                                    >  2:202312/1‑251 (MQ=255)                                 aaTAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCA‑TCGCCACGAAGAAGCCGA‑TGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATc                                                                                                                                                                                                                    >  2:163998/1‑250 (MQ=255)                                  aTAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACCCCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCA‑TCGCCACGAAGAAGACGA‑TGGTGGCAAGATCGCCACACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATc                                                                                                                                                                                                                    >  1:460742/1‑249 (MQ=255)                                     gTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCA‑TCGCCACGAAGAAGCCGA‑TGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCCGATACCACC‑‑GCCGTATTTCCCAGCAGCCACGCCGGATCCCGCCCAATCAGGa                                                                                                                                                                                                                >  2:138423/1‑250 (MQ=255)                                          cGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCA‑TCGCCACGAAGAAGCCGA‑TGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTccc                                                                                                                                                                                                          >  2:222083/1‑250 (MQ=255)                                              gcgATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCA‑TCGCCACGAAGAAGCCGA‑TGGTGGCAAGATCGCCAAACGGCAGCCCCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAcc                                                                                                                                                                                                      <  1:343886/251‑1 (MQ=255)                                              gcgATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCA‑TCGCCACGAAGAAGCCGA‑TGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAcc                                                                                                                                                                                                      >  1:435569/1‑251 (MQ=255)                                              cggATCCACAGGGTGATCCTCATAATGAGGACCCCCGAGATTAAGTTGCGAAACAGATTTCCACGATGCACGGCGCCCGCCA‑TCGCCACGAAGAAGCCGA‑TGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGAACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAACCAGG‑CTTGCCTATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAcc                                                                                                                                                                                                      <  1:323451/249‑1 (MQ=255)                                                            gATGCTCATAATGATGACCCCCGAGATTAAGGGGCGGAACAGATTCCCACGATGCCCGGCGACCGCCA‑TCGCCACGAAGAAGCCGA‑TGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCt                                                                                                                                                                                        <  1:143086/251‑1 (MQ=255)                                                                              ccccGAGATTAAGGTGCGGAACAGATTCCCAGGATGCACGGCGACCGCCA‑TCGCCACGAAGAAGCCGA‑TGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGGGGGATAAAAATCAGG‑CTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCt                                                                                                                                                                      <  1:202312/250‑1 (MQ=255)                                                                                  ggaaTTAGGGTGCGGAAAAGATTCCCACGCTGCACGCCACCCCCCA‑TCGCCCCAAAGAAGCCGA‑GGGGGGAAAGATCGCCAAAGGGCAGCCCCTGATTCCCCGGCACAAAAACAGAAATAAAAATGGTGAGTGGGATAAAAA‑CAGGGCTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCtt                                                                                                                                                                  <  2:91567/248‑1 (MQ=255)                                                                                         aGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCA‑TCGCCACGAAGAAGCCGA‑TGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGAGAGTGGGATAAAAATCAGG‑CTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATg                                                                                                                                                           >  1:466449/1‑251 (MQ=255)                                                                                             gCGGAACAGATTTCCACGATGCACGGCGACCGCCA‑TCGCCACGAAGAAGACGA‑TGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCCGATACCACC‑‑GCCGTATGTCCCCGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCg                                                                                                                                                       >  1:331799/1‑251 (MQ=255)                                                                                                  acatatgccccgaagccACGGCGACC‑CCAGTCCCCACGAAAGAGCCGA‑TGGGGGCAAGACCCACAACCGGCACCACCTGATTCCCGGGCACACAAACAGAAAATAAAATGGTGAGGGGGATAAAAACCAGG‑CTTGCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCGGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTaaa                                                                                                                                                  <  2:264373/235‑1 (MQ=255)                                                                                                           cACGATGCACGGCGACAGCCA‑TCGCCACGAAGAAGCCGA‑TCGTGGCAAGATCGCAAAACGGCAGCACCTGATTACCCGGCACACAAACAGCACTTAAAATGGTGAGTGGGATCAAAATCAGG‑CGTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTCAAccatccatt                                                                                                                                         >  2:463737/1‑250 (MQ=255)                                                                                                                          tggggccgcagcggaaagaagCCGA‑TGGTGGCAAAATCCCCAAACGGCAGCACCTGATTCCCCGCCACAAAAACAGCAATTAAAATGGTGGGTGGGAAAAAAACCAGG‑CTTGCCGATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCa                                                                                                                            <  1:174640/235‑1 (MQ=255)                                                                                                                           cGCCA‑TCGCCAAGAAGAAGTCGA‑TGGTGGCAAGATCGCCCAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAACTGGTGAGTGGGATAAAAATCAGG‑CTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCAccc                                                                                                                         >  2:309350/1‑251 (MQ=255)                                                                                                                            gCCA‑TCGCCAAGAAGAAGCCGA‑TGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCg                                                                                                                        >  2:55035/1‑251 (MQ=255)                                                                                                                                cccGCCACGAAAAAGC‑GAGTGGGGGCAAGACCCCCAACCGGAAGCCCCGGATTCCCCGCCACAAAAACAGAAATAAAAAGGGTGAGGGGAAAAAAAACCAGG‑CT‑GGCGAAAAC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGACCAAGGCCAATCAGGAACTCCGGACCGCCTAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGAGGGGTTTAATCACCCGTgg                                                                                                                     <  2:213905/249‑1 (MQ=255)                                                                                                                                    ccACGAAGAAGCCGA‑TGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCCATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCAtt                                                                                                                 >  2:45501/1‑251 (MQ=255)                                                                                                                                    ccACGAAGAAGCCGA‑TGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTATTCACCCGTGGCAtt                                                                                                                 >  2:468920/1‑251 (MQ=255)                                                                                                                                        gaagaaGCCGA‑TGGTGGCAAGATCGCCAAACGGCAGCAACTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTACACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTagc                                                                                                              >  1:47328/1‑250 (MQ=255)                                                                                                                                           gaagCCGA‑TGGTGGCAAGATCGCAAAACGGCAGCACCTGATTCCCCGGCACACAAACAGAAATTAAAATGGTGAGTGGGAAAAAAATCAGG‑CTTGCCGATACCACC‑‑GCCTTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAgcagca                                                                                                          <  1:246621/251‑1 (MQ=255)                                                                                                                                                        ggcaagtcccccAAAGGGCAGCACCTGATTCCCCGGCACAAAAACAAAACTTAAATTGGGGAGGGGAATAAAAAGCAGG‑CTTGCCGAT‑CC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGccgccg                                                                                              <  1:268617/241‑1 (MQ=255)                                                                                                                                                                    aaaCGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGcc                                                                                   >  1:239899/1‑250 (MQ=255)                                                                                                                                                                    aaaCGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTAAGTGGGATAAAAATCAGG‑CTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGcc                                                                                   >  1:321565/1‑250 (MQ=255)                                                                                                                                                                              ccctaaTTACCCGGCACACAAAAAGAAATAAAAAGGGTGAGTGGGATAAAAATCGGG‑CTTGCCGA‑GCC‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATac                                                                        <  1:36142/246‑1 (MQ=255)                                                                                                                                                                                cTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGCGGGATAAAAATCAGG‑CTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGACCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGAAAACGACTACGAGCCTGCTTCGCGATCGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATacac                                                                      >  1:12066/1‑251 (MQ=255)                                                                                                                                                                                 tGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccg                                                                                                                                                                                                     <  1:357515/123‑1 (MQ=255)                                                                                                                                                                                 tGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccg                                                                                                                                                                                                     >  2:357515/1‑123 (MQ=255)                                                                                                                                                                                      aCCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCACTACACCTTTg                                                                 >  1:336355/1‑250 (MQ=255)                                                                                                                                                                                           gCACACAAACAGCAATAAAAATGGTGAGTGGGAAAAAAATCAGG‑CTTGCCGATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATTCACCTTTGACATCg                                                           <  1:99638/251‑1 (MQ=255)                                                                                                                                                                                               acaAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCCGCTACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAAc                                                       >  1:76546/1‑251 (MQ=255)                                                                                                                                                                                                                            aaaaaaTCAGG‑CTTGCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATc                          <  1:222083/250‑1 (MQ=255)                                                                                                                                                                                                                             aaaaaTCAGG‑CTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGCAGCCGAtccaccttcaatggg                         >  2:225757/1‑237 (MQ=255)                                                                                                                                                                                                                              aaaaTCAGG‑CTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTg                                                                              >  2:274163/1‑197 (MQ=255)                                                                                                                                                                                                                              aaaaCCAGG‑CTTGCCGATCCCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTg                                                                              <  1:274163/197‑1 (MQ=255)                                                                                                                                                                                                                                    aGG‑CTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGCTAATCAGCCCCATCACAAAAcc                  >  1:41859/1‑251 (MQ=255)                                                                                                                                                                                                                                        cTTGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCAc               >  2:82528/1‑251 (MQ=255)                                                                                                                                                                                                                                         ttGCCGATACCACC‑‑GCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACg              >  2:312509/1‑251 (MQ=255)                                                                                                                                                                                                                                            ggcgaccccccccGCCGTATGCCCCAGCAGCAACGCCGAATAAAGGCCAATAAGGACCTCCGGACCGCCGAACTTCGCCTGGAAACGACTTCGAGCCTGCTTGGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCGGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACgg             <  2:214745/242‑1 (MQ=255)                                                                                                                                                                                                                                                     caccgccgTATGTCCCAGGAGAAACGCGGGATCAAGGCCAATCAGGAACTCCGAACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCAGCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACGGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGAGAATCAGCCGCATCACAAAACCCGCGGTGACAGGCt    <  2:47328/249‑1 (MQ=255)                                                                                                                                                                                                                                                     caccgccgTATGTCCCAGCAGCCACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGCGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAACCCCACGGTGACAGGCt    >  2:216727/3‑251 (MQ=255)                                                                                                                                                                                                                                                       ccgccgAATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCg  >  1:141369/1‑251 (MQ=255)                                                                                                                                                                                                                                                       ||                                                                                                                                                                                                                                                          CCAGCATTAGCCGCCAGTTGGGTGTGCAAACCAATAGTTTGCGTCGCGATCCACAGGGTGATGCTCATAATGATGACACCCGAGATTAAGGTGCGGAACAGATTTCCACGATGCACGGCGACCGCCA‑TCGCCACGAAGAAGCCGA‑TGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGG‑CTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTAGCAGCATCACTGCCGCCGTTTTTACCGCCAGCTGCAATACACCTTTGACATCGTAACCCGCGAGGATGCCGATAATCAGCCCCATCACAAAACCCACGGTGACAGGCTCG  >  NC_000913/2173120‑2173613

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