Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,623,578	C→T	S180S (TCG→TCA)	lplA ←	lipoate‑‑protein ligase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,623,578	0	C	T	96.9%	197.7 / NA	66	S180S (TCG→TCA)	lplA	lipoate‑‑protein ligase A
Reads supporting (aligned to +/- strand):  ref base C (0/1);  major base T (40/22);  minor base G (0/2);  total (41/25)
Fisher's exact test for biased strand distribution p-value = 1.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAATGGGCGAAAAAGGCCTCGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAACAGCTCGGTGAGGTTGGTCACGCGGGAACGTACCGACGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTGAGGTCGGCATTGAGT‑AGCAAGGTGCCGTGGTGGAAGCCGCGATCTTTGGTTTCGCGATAGGCCGAGCCTGAGACTTTGCGGTCGCCTTCGACGGTTTTCACCACCAGATCGTTACGTCCGGACGCTTCGGCGCTGACGCCGAGCGCGTTCAGCGCATTGAGCACAATCGACGTGGAGA  >  NC_000913/4623338‑4623819                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                    aaaGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AAt                                                                                                                                                                                                                                            <  1:246815/249‑1 (MQ=255) aaaGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGc                                                                                                                                                                                                                                          <  1:15805/251‑1 (MQ=255)      gtacctCCAGCAGATGCTAGAAGGCCGGAGCCTGACCGAAGTTCCATTCTCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTc                                                                                                                                                                                                                                     <  2:565100/246‑1 (MQ=255)         cATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTcgcc                                                                                                                                                                                                                                  <  2:240777/251‑1 (MQ=255)            ccAGCAGATGCGAGAATGCCGGAGCCTGACCGCAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTATCCGGGGAGATGATTTCCGATTATACACGCTATACATAATGGGCGCAAAAGGCCTTGGTTATGGCCTCGAATACCAGCTCATGGGTGATCCCCGGCAGTAGCTCGGTGAGACTGGTAACGCGGGAACGCACTG‑CATGAATGCCTTTCgccgcc                                                                                                                                                                                                                               >  2:28774/1‑251 (MQ=255)              agcagATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTTGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCGCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAg                                                                                                                                                                                                                             <  2:371130/251‑1 (MQ=255)                    ggCGAGAATGCCGGGGCCTGACCGAAGTTCCTTTCCCAGCTTCTCTGACGGGCAAAGGTTTCGGCGAAGTTGGGCAAGGCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACCCGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGAGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCtt                                                                                                                                                                                                                       <  2:66009/250‑1 (MQ=255)                        agaATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGCTTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCtttttt                                                                                                                                                                                                                   >  2:486906/1‑250 (MQ=255)                             gCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCg                                                                                                                                                                                                               <  1:615894/250‑1 (MQ=255)                               cGGAGCCTGACCGAAGGTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGAt                                                                                                                                                                                                            <  1:87375/251‑1 (MQ=255)                                ggAGCCTGACCGACGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGACGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCACGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGA‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGAtt                                                                                                                                                                                                           >  1:622251/1‑251 (MQ=255)                                ggAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGAtt                                                                                                                                                                                                           >  1:625020/1‑251 (MQ=255)                                    ccTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGCTGATTTACGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGATTATGGCCTAGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTAGCCGCCAGTTTCTTTTTATcagggatgagg                                                                                                                                                                                                       >  2:640427/1‑241 (MQ=255)                                       gACCGAAGTACCATTCCCAGCTACTGTGAGGGGCAAAGGTTTCGGGTAAGTTGGGAAAGTCTGGCGTTTTGT‑CCGGGAGATGATTTCCGCTTCTACACTCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTAGCCGGATTGAGATAg                                                                                                                                                                                                    <  1:496789/250‑1 (MQ=255)                                            aaGTTCCATTCCCAGCTACTCTGGCGGGCAAAGGTTGCGGCGAAGTTGGGCAAGTCTGGCGTTTTGTCCGGGTAGATGATTTCCGCTGCGACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGGGAGGTTTGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTTAGATAGTTTGc                                                                                                                                                                                               <  2:622251/251‑1 (MQ=255)                                                                  gacgtgCAAAGGTTTCGGCGAAGTTTGGCAAGTTGGGTGTGTTGGCCGGGGAGATGATTGCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCCCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCTGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCAtt                                                                                                                                                                          <  2:625020/245‑1 (MQ=255)                                                                             gTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATACTGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAAGGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGAGAGTTTGCCAAGCGGCTAAGGGCGGCCTTGACC‑AGCAAggg                                                                                                                                                             >  2:322347/1‑250 (MQ=255)                                                                                  ggCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGC‑AGCAAGGTGCCgt                                                                                                                                                        >  1:584754/1‑251 (MQ=255)                                                                                  ggCGAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCCGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGC‑AGCAAGGTGCCgt                                                                                                                                                        >  1:83547/1‑251 (MQ=255)                                                                                  ggCGAAGTTTGGCAAGACTGGCGTTTTGTCCGGGGAGATGCTTTCAGCTTCTACACGCTCGCCATACTGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTAGGTGAGGTTGGTACCGCGGGAAAGCACTGAGGT‑AATGCCTTTCGCCGCCAGGTTCTTTTTATCCGGATTGAGATAGTTGGCCCGGCGGCTAAGGTCGGCATT‑TGTCAGCAAGGTGCCgg                                                                                                                                                        >  2:538409/1‑250 (MQ=255)                                                                                     gAAGTTTGGCAAGTCTGGCGTTTTGTCCGGGGAAATGATTTCCGCTTCTACACTCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCGGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGC‑AGCAAGGTGCCgtggt                                                                                                                                                     <  2:641912/251‑1 (MQ=255)                                                                                            ggCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTg                                                                                                                                                                         >  1:369649/1‑225 (MQ=255)                                                                                             gCAAGTCTGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGC‑AGCAAGGTGCCGTGGTGGAAGcca                                                                                                                                             <  1:619667/251‑2 (MQ=255)                                                                                              cAAGTCTGTGGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCCCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGC‑AGCAAGGTGCCGTGGTGGAAGccac                                                                                                                                            <  2:545451/251‑3 (MQ=255)                                                                                              cAAGTCTGGCGTTTTGCCCGGGGCGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCTGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGC‑AGCAAGGTGCCGTGGTGGAAGctcc                                                                                                                                            >  2:382994/1‑248 (MQ=255)                                                                                                    tGGCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGC‑AGCAAGGTGCCGTGGTGGAAGCCACGATCtt                                                                                                                                      >  1:363176/1‑251 (MQ=255)                                                                                                    gggCGTTTGGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTGGGTCACGCGGGAACGCACTGATGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGTCATTGAGC‑AGCAAGGTGCCGTGGTGGAAGCCACGATCtt                                                                                                                                      <  1:157567/250‑1 (MQ=255)                                                                                                      gCGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTTGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGC‑AGCAAGGTGCCGTGGTGGAAGCCACGATCTTTg                                                                                                                                    >  2:620229/1‑251 (MQ=255)                                                                                                       cGTTTTGTCCGGGGAGATGATTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGC‑AGCAAGGTGCCGTGGTGGAAGCCACGATCTTTgg                                                                                                                                   >  1:390881/1‑251 (MQ=255)                                                                                                            ggtcgggggggATGATTTCCGCTTCTACACTGTCGACATAATGGGCGAAAAAGGCCTTGGTTAGGGCCTCGCAATCCTGCTCATGGGTGATCCCCGGCAGCAGCGCGGTGAGGTTGGTCACGCGGGAACGCTCGGAGGT‑AATTCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTGGCCAGGCGGCTAAGGTCGGCATTGAGC‑AGCAAGGTGCCGTGGTGGAAGCCACGATCTTTGGTTTc                                                                                                                               <  2:584754/245‑1 (MQ=255)                                                                                                             gtccGGGGAGATGATTTCAGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGCTTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGCGC‑AGCAAGGTGCCGTGGTGGAAGCCACGATCTTTGGTTTcgc                                                                                                                             >  2:109943/4‑251 (MQ=255)                                                                                                                  gggggttgatTTCCGCGTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTGATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGC‑AGCAAGGTGCCGGGGTGGAAGCCACGATCTTTGGGTTCTCGATAg                                                                                                                        <  2:286843/245‑1 (MQ=255)                                                                                                                      gatgatTTCCGCTTCTACACGCTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGC‑AGCAAGGTGCCGTGGTGGAAGCCACGATCTTTGGTTTCGCGATAGGCCg                                                                                                                    >  1:666201/1‑251 (MQ=255)                                                                                                                          atTTCAGCTTCTACACGCGCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACAGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGGTTGAGATAGTTTGCCAGGCGGCTACTGTCGGCATTGAGC‑AGCAAGGTGCCGTGGTGGACGCCACGATCTTTGCGTTCGCGATAGGCCGaacc                                                                                                                >  1:73672/1‑248 (MQ=255)                                                                                                                                 cttctacacgcTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGC‑AGCAAGGTGACGAGGTGGAAGCCACGATTTTTGGTTTCGCGATAGGCCGAACCTGAGACt                                                                                                         >  2:117934/9‑251 (MQ=255)                                                                                                                                   tctacacgcTCGCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGC‑AGCAAGGGGCCGTGGTGGAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGGCttt                                                                                                       >  2:339185/7‑251 (MQ=255)                                                                                                                                              gCCATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AGTGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGC‑AGCAAGGTGCCGTGGTGGAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCtt                                                                                            >  1:295172/1‑251 (MQ=255)                                                                                                                                              gCCATAACGGGCGAAAAAGGCCTTGGTTATGGAATAGCAAACCTGCTCATGGGGGATCACCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGACCGCACTGAGGT‑AATGCCTTTCGCGGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGC‑AGCCAGGGGCCGAGGTGGATGCTACGTTCTTTGGCTTCGCGCTAGGCCGAGCCTGAGACTTTGCGGtggactt                                                                                            >  2:565540/1‑245 (MQ=255)                                                                                                                                               ccATAATGGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGC‑AGCAAGGTGCCGTGGTGGAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTc                                                                                           <  1:91461/251‑1 (MQ=255)                                                                                                                                               ccATAATAGGCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCCGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGC‑AGCACGGTGCCGTGGTGGAAGCCCCGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTTCGGTCGCCTTc                                                                                           >  2:71047/1‑251 (MQ=255)                                                                                                                                                      gggCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGC‑AGCAAGGTGCCGTGGTGGAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGtt                                                                                    >  1:135074/1‑251 (MQ=255)                                                                                                                                                      gggCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGC‑AGCAAGGTGCCGTGGTGGAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGtt                                                                                    >  1:403983/1‑251 (MQ=255)                                                                                                                                                       tgCGAAAAAGGCCTTGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGCTTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGC‑AGCACGGCGCCGTGGTGGACGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGAATTGAGGGTCGCCTTCggagggtt                                                                                   >  2:547248/2‑244 (MQ=255)                                                                                                                                                            aaaaGGCCTCGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAACAGCTCGGTGAGGTTGGTCACGCGGGAACGTACCGACGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTGAGGTCGGCATTGAGT‑AGCAAGGTGCCGTGGTGGAAGCCGCGATCTTTGGTTTCGCGATAGGCCGAGCCTGAGACTTTGCGGTCGCCTTCGACGGTTTTcacc                                                                              <  1:21174/251‑1 (MQ=255)                                                                                                                                                                  ccttggTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGGGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCCGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGC‑AGCAAGGTGCCGTGGTGGAAGCCCCGATCTTTGGTTTCGCGATAGGCCGCTCCTGAGACTTTGCGGTCGCCTTCGGCGGTTTTcaccccc                                                                           >  1:268086/5‑245 (MQ=255)                                                                                                                                                                     tggTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGC‑AGCAAGGTGCCGTGGTGGAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGTTTTCACCACCAGAt                                                                       <  1:224560/248‑1 (MQ=255)                                                                                                                                                                      ggTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGC‑AGCAAGGTGCCGTGGTGGAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGgtttacacaccaaaccct                                                                    >  2:661649/1‑237 (MQ=255)                                                                                                                                                                         tATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGTACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGC‑AGCAAGGTGCCGTGGTGGAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGTTTTCACCACCAGATCGTTa                                                                  <  1:356535/250‑1 (MQ=255)                                                                                                                                                                         tATGGCATCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGC‑AGCAAGGTGCCGTGGTGGAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGTTTTCACCAACAGATCGTTAc                                                                 >  1:113962/1‑251 (MQ=255)                                                                                                                                                                          aTGGCCTCGCAAACCTGCTCATGGGTGACCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGC‑AGCAAGGTGCCGTGGTGGAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGTTTTCACCACCAGATCGTTAc                                                                 <  1:491125/250‑1 (MQ=255)                                                                                                                                                                            ggCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGC‑AGCAAGGTGCCGTGGTGGAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGGTTTCACCACCAGATCGTTACgcc                                                              >  1:652859/1‑249 (MQ=255)                                                                                                                                                                              ccTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGC‑AGCAAGGTGCCGTGGTGGAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGTTTTCACCACCAGATCGTTaagcccg                                                            >  1:345520/1‑245 (MQ=255)                                                                                                                                                                               cTCGCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGCGGTTGGTCACGCGGGACCGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCCTTGAGC‑AGCAAGGTGCCGTGGGGGAAGCAACGATCTTTGGTTTCGCGCTCGGGCGCTCCTGAGCCTTTGCGGTCGCCTTCGGCGTGTTTCACCACCGGATCGtatccccccg                                                           >  2:657649/1‑242 (MQ=255)                                                                                                                                                                                  gCAAACCTGCTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGC‑AGCAAGGTGCCGTGGTGGAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGCCTTTGCGGTCGCCTTCGGCGGTTTTTACCCCCCGAGCGGTACGCCCGGACg                                                        >  2:170548/1‑251 (MQ=255)                                                                                                                                                                                          gcccatggggtaTCCCCGGCGTCAGCTCGGTGAGGTTGGTCACGCGGGAACGTACGGAGGT‑AATGCCTTTCGCCGCCAGTTTCGTTTTATCCCGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCAGTGAGC‑AGCAAGGGGCCGTTGTGGAAGCCACGATCTTTGGTTTCGCGATATGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGTTTTCACCACCAGATCGTTACGCCCGGACGCTTCgg                                                  <  2:577174/238‑1 (MQ=255)                                                                                                                                                                                           cTCATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGC‑AGCAAGGTGCCGTGGTGGAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGTTTTCACCAACAGATCGGTACGCCCGGACGCTTCGgcgc                                               >  1:506825/1‑251 (MQ=255)                                                                                                                                                                                             cATGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGC‑AGCAAGGTGCCGTGGTGGAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGTTTTCACCACCAGATCGTTACGCCCGGACGCTTCGGCGCTg                                             >  1:498080/1‑251 (MQ=255)                                                                                                                                                                                               tGGGTGATCCCCGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGC‑AGCAAGGTGCCGTGGTGGAAGCCACGATCTTTGGTTTCGCGATAGGCCGatcc                                                                                                                >  1:13621/1‑179 (MQ=255)                                                                                                                                                                                                   tGATCCCCGGCAGCAGCTCGGTGAGATTGGTCACGAGGGACCGCACTGAGGT‑AATGCCTTGCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGTCTAAGGTCGGCATTGAGC‑AGCAAGGTGCCGTGGTGGGAGCCACGATCTTTGGTTTCTCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGTTTTACCCCCCAACTCGTTCCGCCCGGACGCTTCTGCTCTCACGccc                                       >  2:336000/1‑250 (MQ=255)                                                                                                                                                                                                        cccGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGC‑AGCAAGGTGCCGTGGTGGAAGCCACGCTCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGTTTTCACCACCAGGTCGTTACGCCCCGGCGCTTCGgcggttgcgccaagccc                                  >  1:248254/1‑237 (MQ=255)                                                                                                                                                                                                        cccGGCAGCAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGC‑AGCAAGGTGCCGTGGTGGAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGTTTTCACCACCAGATCGTTACGCCCGGACGCTTCGGCGCTGACGCCAAgcgc                                  >  1:469651/1‑251 (MQ=255)                                                                                                                                                                                                            gcagcAGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGC‑AGCAAGGTGCCGTGGTGGAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGTTTTCACCACCAGATCGTTACGCCCGGACGCTTCGGCGCTGACGCCAAGCGCGTTc                              <  1:486906/247‑1 (MQ=255)                                                                                                                                                                                                                 aGCTCGGTGAGGTTGGTCACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGC‑AGCAAGGTGCCGTGGTGGAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGTTTTCACCCCCAGATCGTTACGCCCGGACGCTTCCGCGCCGACGCCAAGCGCGTTCAgccc                         >  1:644171/1‑249 (MQ=255)                                                                                                                                                                                                                                  cACGCGGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGC‑AGCAAGGTGCCGTGGTGGAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGTTTTCACCACCAGATCGTTACGCCCGGACGCTTCGGCGCTGACGCCAAGCGCGTTCAGCGCATTGAGCACAATCGAc         <  1:135981/250‑1 (MQ=255)                                                                                                                                                                                                                                    cgcgGGAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGC‑AGAAAGGTGCCGTGGTGGAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGATTTCACCACCAAATCGGTACGCACGGACGATTAGGCCCTGACGCCAAGCGCGCTACGCGCATTGAGCACAAtagccgtg      >  1:124203/1‑244 (MQ=255)                                                                                                                                                                                                                                        ggAACGCACTGAGGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTAAGGTCGGCATTGAGC‑AGCAAGGTGCCGTGGTGGAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGTTTTCACCACCAGGACGTTACACCCGGCCGCTTCGGCGCTGCCGCCACGCGCGTTCAGCGCagtgcgctcaacctacgttgagg  >  2:224981/1‑229 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                    AAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAATGGGCGAAAAAGGCCTCGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAACAGCTCGGTGAGGTTGGTCACGCGGGAACGTACCGACGT‑AATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTGAGGTCGGCATTGAGT‑AGCAAGGTGCCGTGGTGGAAGCCGCGATCTTTGGTTTCGCGATAGGCCGAGCCTGAGACTTTGCGGTCGCCTTCGACGGTTTTCACCACCAGATCGTTACGTCCGGACGCTTCGGCGCTGACGCCGAGCGCGTTCAGCGCATTGAGCACAATCGACGTGGAGA  >  NC_000913/4623338‑4623819

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: ‑