Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	26,648	C→T	H124H (CAC→CAT)	ispH →	4‑hydroxy‑3‑methylbut‑2‑enyl diphosphate reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	26,648	0	C	T	97.1%	203.3 / ‑3.8	69	H124H (CAC→CAT)	ispH	4‑hydroxy‑3‑methylbut‑2‑enyl diphosphate reductase
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (32/34);  total (33/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-01

CCGCTATGTGGTCGATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGTCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑TC‑GCCGTGGCGAAGAATCTATTCTCATCGGTCACGCCGGGCACCCGGAAGTGGAAGGGACAATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTCCCGAAAATTGTCGGTCCGCGCAAAGATGACATCTGCTACGCCACGACTAACCGTCAGGAAGCG  >  NC_000913/26402‑26894                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                       ccGCTATGAGGACGATAGCTCGCGTGAGCGTGGGGCTATCTTTCTTGAGCAGATTAGCGAAGTACAGGACGGCGCGATCCTGATTGTCTTCGCACACGGGGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGACGCGATTTGACGGAGGTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGGTGCGAAAAGTCGATTCTCATCGGCCACGCCGGGCATCCgg                                                                                                                                                                                                                                                    >  2:137683/1‑251 (MQ=255)      attgggTCGATAGCTTTCGTGAGCGTGGGGTTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTGCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGTCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑TC‑GCCGTGGCGAAGAATCTATTCTCATCGGTCACGCCGGGCCCCCGGAAg                                                                                                                                                                                                                                                 <  2:479571/245‑1 (MQ=255)        gtgGTCGAGATCTTGCTTTAGCGTGGGGCTATCTTTATTGAGCAGATTATCGAAGTTCCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTTCGTAACGAAGCAAAAAGCCTCGATTTGACGGTGTGTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAATTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGa                                                                                                                                                                                                                                             <  2:541023/251‑1 (MQ=255)               aTAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGCAGTACCGGCCGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGCGACCAAAGGGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGCGGTGGAAggtacg                                                                                                                                                                                                                                      >  1:651350/1‑247 (MQ=255)                 aGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAggtacgat                                                                                                                                                                                                                                    <  1:504640/251‑7 (MQ=255)                         ggAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCcagt                                                                                                                                                                                                                            <  1:290244/250‑1 (MQ=255)                           ggCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCcagtac                                                                                                                                                                                                                          <  1:1879/250‑1 (MQ=255)                                ggggCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCcagtacagta                                                                                                                                                                                                                      >  1:96062/1‑250 (MQ=255)                                 gggCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAAc                                                                                                                                                                                                                    >  1:58327/1‑251 (MQ=255)                                  ggCTATCTTTATTGGGCAGATTAGCGAAGTACCGGCCGGCGCGATCCTGATTTTCCCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGGCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAAcc                                                                                                                                                                                                                   <  2:43184/251‑1 (MQ=255)                                                 gCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACTGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGt                                                                                                                                                                                                    <  1:51101/251‑1 (MQ=255)                                                    gATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATc                                                                                                                                                                                                 <  1:233382/251‑1 (MQ=255)                                                     tttAGAGAAGAACCGGACGGCGCGATCCTGCTTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAAAGAAGCAAAAAGCCGCGATTTGACGGTGTTTGCTGCCACCTGTCCGCTGGTGACCAACGTGAATCTGGAAGTAGCACGCGCCAG‑CC‑GTCGTGGCGGAGAGTCGCTTCTCATAGGCCACGCCGGGCATCCGGAAGTGGAAGGTACAATGGGCCAGCACAGTAACCCGGCAggggggctgattct                                                                                                                                                                                                >  2:239633/2‑242 (MQ=255)                                                          cGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCg                                                                                                                                                                                           <  1:204726/251‑1 (MQ=255)                                                              gTACCGGACGGCGCGATCCTTATTTTCTCCGCACACGGTGTTTCGCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATGTGACGGTGTTTGATGCCCCCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATc                                                                                                                                                                                       <  2:194581/251‑1 (MQ=255)                                                               gaCCGGACGGCGAGATTCTGTTTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGTTGCGTTGTGGCGGTGTTTGATGCCACCGGTCCGCTGGGGTCCAAAGGGCATATGGAAGCCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCAGCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTCCAGTAACCCGGAAGGGGGAATGTATCTGGTAGaagcg                                                                                                                                                                                      <  2:376735/250‑4 (MQ=255)                                                                aCCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCg                                                                                                                                                                                      <  1:17025/250‑1 (MQ=255)                                                                  cGGACAGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGACGAATCGCCg                                                                                                                                                                                   >  1:489675/1‑251 (MQ=255)                                                                    gACGGCGCGATCCGTATTTCTTCCGCACCCGGTGTTTCGCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGAGTCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCTCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCCAGGCCGGGCATCCGGAAGTGGAAGGTACGATGTGCCAGTACAGTAGCCCGGAAGGGGGAATTTATCTGGTCGAATTGCCGGa                                                                                                                                                                                 <  2:640303/251‑1 (MQ=255)                                                                     aCGGCGCGATCCTGATTTTCTCCGCTCACGTTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGGTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGAc                                                                                                                                                                                <  1:524642/251‑1 (MQ=255)                                                                     aCGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGAc                                                                                                                                                                                >  1:511134/1‑251 (MQ=255)                                                                         cgcgATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGAtg                                                                                                                                                                            <  1:618366/251‑1 (MQ=255)                                                                          gcgATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGAtgt                                                                                                                                                                           >  1:194208/1‑251 (MQ=255)                                                                            gATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGCAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGGATGTATCTGGTCGAATCGCCGGACGAtgtgt                                                                                                                                                                         >  1:567684/1‑251 (MQ=255)                                                                                   atttctgccGCACACGGTGTTTCTCAGGCGGTCCGTAACGAAGCAAAAAGCCGAGATTTGACGGTGTTTGATGCCCCCTGTCCGCTGGTGACCAAAGTGAATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGGGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGGATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAAca                                                                                                                                                                  <  2:33752/244‑2 (MQ=255)                                                                                    atttatgccgactaggttgttCTTAGGCGGGACGTAAAGAAGCAAAAAGCAGCGATTTGACGGGGTTTTATGCAACCGGCCTGCTGGTGCCCAAAGTGCATATGGAAGGCGCCCGCGCCAGCTT‑GCAG‑GGCGAAGAGTCGATTCCCAACGGACACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCTGGAAGGGGGAATGTATCTGGTCGAATCGCCTGACGATGTGTGGAAACTg                                                                                                                                                                 <  2:126036/232‑1 (MQ=255)                                                                                     tttCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCTCCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTg                                                                                                                                                                 <  1:16671/250‑1 (MQ=255)                                                                                           cGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTc                                                                                                                                                           >  1:589189/1‑250 (MQ=255)                                                                                           cGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCa                                                                                                                                                          <  1:324194/251‑1 (MQ=255)                                                                                             cacaCGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGTAGTGGTAGGTACGATGGGCCAGTACAGTAACCCGGGAGGGGGAATGTATCTGGTCGAATCGCCGGACGAGTTGGGGGAACTAACGGTCaa                                                                                                                                                         >  2:253061/1‑250 (MQ=255)                                                                                                acGGTGGTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGAGTTGACGGTGTGTGAGGCCCCCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAAc                                                                                                                                                     <  2:58327/251‑1 (MQ=255)                                                                                                    tgtTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGa                                                                                                                                                   >  1:78523/1‑249 (MQ=255)                                                                                                            aGGCGGTACGTAACGAAGCAAAAAGACGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGACGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGACGAGAAGctc                                                                                                                                         >  1:483803/1‑251 (MQ=255)                                                                                                             ggCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAATTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGctct                                                                                                                                        <  1:382809/251‑1 (MQ=255)                                                                                                               gtgtatgTAACCAAGCAAAAAGCCGTGAGTTGAAGTTGTTTGAGGCACCCTGTCCGCTGGTGCCCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGTAGTTACGATGGGCCAGTACAGGAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCTGACGATGGGTGGAAACTGACGGTCAAAATCGAAGAGAAGctctac                                                                                                                                      <  2:479171/245‑3 (MQ=255)                                                                                                                   accgaaCGAAGCAAAAAGCCGTGATTTGACGGTTTTTGATGCCGCCTTGCCGCTGGTGTCCAAATTGCATTTGGAAGTCGCCCGCGCCAG‑CC‑GTCGGGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCttt                                                                                                                                   <  2:430311/246‑1 (MQ=255)                                                                                                                   aCGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTa                                                                                                                                  <  1:146492/251‑1 (MQ=255)                                                                                                                       aaCGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGAc                                                                                                                              >  1:397084/1‑251 (MQ=255)                                                                                                                                aaaaGCCGCGATTGGACTGTGTTTGTTGCTACCTGGCCGCTGGTGACCAAAGTGCTTATGGAAGTCGCCCGCGCCAG‑CC‑GGCGTGGCGAAGAGGCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTATCCCGGAAGGGGGAATGTATCTGGTAGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCAc                                                                                                                     <  2:356409/251‑1 (MQ=255)                                                                                                                                 aaaGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCAACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGAcaa                                                                                                                      >  1:389067/1‑247 (MQ=255)                                                                                                                                         gATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTc                                                                                                               <  1:103076/248‑1 (MQ=255)                                                                                                                                         gATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCCCTCCTTTCTGACCCAGACCAcactttc                                                                                                               >  2:103076/1‑242 (MQ=255)                                                                                                                                              gACGGTGCTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCGACGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGCAGGTACGATGGGCCAGTACCGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCAGGACGATGCGGGGAAACTGACGGTCAAAAACGAAGAAAAGCTCTCCTTTATGACCCAGACCACGCTTTCGGTGGATg                                                                                                       >  2:252639/1‑251 (MQ=255)                                                                                                                                                cGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGac                                                                                                     >  1:645707/1‑251 (MQ=255)                                                                                                                                                  gtgtTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGaa                                                                                                                                                                                                              <  2:158486/144‑1 (MQ=255)                                                                                                                                                  gtgtTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGaa                                                                                                                                                                                                              >  1:158486/1‑144 (MQ=255)                                                                                                                                                       ttaTGCCCCCTGTCCGGTGGTGACCAAATTGATTATGTAAGCCGCCCCCGCC‑G‑TCCGTCGTGGCGAAGAGCCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAATGTGCGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGTCGATGTGGGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACTCGTCt                                                                                               <  2:122761/248‑1 (MQ=255)                                                                                                                                                        gATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGa                                                                                             <  1:480145/251‑1 (MQ=255)                                                                                                                                                                      gCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAag                                                                                                                                                 >  1:126007/1‑185 (MQ=255)                                                                                                                                                                      gCTGGTGACCAAAGTGAATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGACACGCCGGGCATCCGGACGTGGAAGGAACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAAAGAAGAGAAGCTCTCCTTTATGAACCAGACCACGCTGTCGGTGGAAGACACGtaggttgtgaacgccgggc                                                                               >  1:669822/1‑233 (MQ=255)                                                                                                                                                                                 aaGTGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGCAGTGGAAGGTACGCTGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGGGTGGAAACTGACGGTCAAAAACGAAGAGACGCTCTCCTTTTTGACCCAGAACACGCTGTCGGTGGATGACACGTCTGATGTGATAGACGCGCTGCGGAAACGc                                                                    >  2:257378/1‑251 (MQ=255)                                                                                                                                                                                  aGTGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCt                                                                   <  1:38673/251‑1 (MQ=255)                                                                                                                                                                                   gTGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGGTGGGCCAGTACAGGAACCCGGAAGGGGGAATGGATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCATTTATGAACCAGACCACGCGATCGGAGGGTGACCCTTCTGGTGTGATCGAAGGGCGGCGTAGACGCtt                                                                  >  2:465199/1‑251 (MQ=255)                                                                                                                                                                                           ggAAGTCGGCCCCGCCAG‑CC‑GTCGGTGCGGAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGTTGGGGCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCTGACGATGTGTTGAAACTGACGGTCAAAAACGAAGAGAATCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTCCCGaaaa                                                          <  2:558736/251‑1 (MQ=255)                                                                                                                                                                                                            cagccgtcGTGGCGAAGCGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGAACGATGGGCCAGTACAGTCAACCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCAAGCAGTCGGTGGATGACACGTCTGATGTGATCGCAGCACTGCGTTAACGCTTCCCGAAACTTGTCGGACCgcgc                                            >  2:576564/8‑250 (MQ=255)                                                                                                                                                                                                            cagccgtcGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTCCCGAAAATTGTCGGTCCgcgc                                            >  1:546223/8‑250 (MQ=255)                                                                                                                                                                                                             ggccgtcggggCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGGGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTCCCGAAAATTGTCGGTCCGCGCaa                                          <  2:546223/242‑1 (MQ=255)                                                                                                                                                                                                                cgtcGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGCGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGAAACGTCTGATGGGATCGACGCGCTGCGTAAACGCTTCCCGCAAATTGTAGGTCCGCGCaa                                          >  2:337601/4‑248 (MQ=255)                                                                                                                                                                                                                   cGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGGCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCGTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTCCCGAAAATTGTCGGTCCGCGCAAAGATg                                     <  2:397084/250‑1 (MQ=255)                                                                                                                                                                                                                     gggCGAAGAGTCGATTCTCATCGGCCACGCCGGGCAGCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGTAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTCCCGAAAATTGTCGGTCCGCGCAAAGATGACa                                  <  2:496297/250‑1 (MQ=255)                                                                                                                                                                                                                         gaagagTCGATTCTCAGCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCGGATGTGATCGACGCGCTGCGTAAACGCTTGCCGAAAATTGTCGGTCCTCGCAATGATGGCATCTg                              <  2:287409/251‑1 (MQ=255)                                                                                                                                                                                                                           agagtcgaTTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGCTCGACGCGCTGCGTAAACGCTTCCCGAAAATTGTCGGTCCGCGCAAAGATGACATCTGCt                            >  1:100858/8‑251 (MQ=255)                                                                                                                                                                                                                                 gaTTCTCATCGGCCACGCCGGGAATCCGGAAGAGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTCCCGAAAATTGTCGGTACGCGCAAAGATGACATCTGCTACGccc                      >  2:455550/2‑250 (MQ=255)                                                                                                                                                                                                                                       cATCGGCCACGCCGGGCATCCGGAAGTGGCAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTCCCGAAAATTGTCGGTCCGCGCAAAGATGACATCTGCTACGCCACGACt                 >  1:529457/1‑250 (MQ=255)                                                                                                                                                                                                                                       cATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTCCCGAAAATTGTCGGTCCGCGCAAAGATGACATCTGCTACGCCACGACt                 <  1:337601/250‑1 (MQ=255)                                                                                                                                                                                                                                        aTCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTCCCGAAAATTGTCGGTCCGCGCAAAGATGACATCTGCTACGCCACGACTa                >  1:398895/1‑250 (MQ=255)                                                                                                                                                                                                                                            gTCACGCCGGGCACCCGGAAGTGGAAGGGACAATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTCCCGAAAATTGTCGGTCCGCGCAAAGATGACTTCTTCTCCGCCACTACTAAccat           >  2:215164/1‑249 (MQ=255)                                                                                                                                                                                                                                                 gCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTCCCGAAAATTGTCGGTCCGCGCAAAGATGACATCTGCTACGCCACGACTAACCGTCAgg       >  1:571500/1‑250 (MQ=255)                                                                                                                                                                                                                                                     ggCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTCCCGAAAATTGTCGGTCCGCGCAAAGATGACATCTGCTACGCCACGACTAACCGTCAGGAAGCg  >  1:89839/1‑251 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                       CCGCTATGTGGTCGATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGTCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑TC‑GCCGTGGCGAAGAATCTATTCTCATCGGTCACGCCGGGCACCCGGAAGTGGAAGGGACAATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTCCCGAAAATTGTCGGTCCGCGCAAAGATGACATCTGCTACGCCACGACTAACCGTCAGGAAGCG  >  NC_000913/26402‑26894

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