Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,298,199	C→T	A390T (GCA→ACA)	fdhF ←	formate dehydrogenase H

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,298,199	0	C	T	95.2%	176.5 / NA	64	A390T (GCA→ACA)	fdhF	formate dehydrogenase H
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base T (37/22);  minor base G (0/3);  total (39/25)
Fisher's exact test for biased strand distribution p-value = 6.08e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTCAACCGCCTTGAAGAAACGCTGGAAGCCACGGTCAGCCGCAGTAAACACGCCTTCATGCTCGCCCCACGACGTTGACGGTAAAATAACATCCGCCGCCGACGCGGTTTTGGTCATAAAGATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGCGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCGAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGT  >  NC_000913/4297956‑4298445                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                       gTTCAACCGCCTTGAAGAAACGCTGGAAGCCACGGTCAGCCGCAGTAAACACGCCTTCATGCTCGCCCCACGACGTTGACGGTAAAATAACATCCGCCGCCGACGCGGTTTTGGTCATAAAGATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACACGGCACGCACTTCGCCATGCGCTGTGCGGTGc                                                                                                                                                                                                                                                 >  2:259160/1‑251 (MQ=255)     aaCCGCCTTGAAGAAACGCTGGAAGCCACGGTCAGCCGCAGTAACCACGCCTTCATGCTCGCCCCACGCCGTTGACGGTAAAATAACATCCGCCGCCGACGCGGTTTTGGTCATAAAGATGTCCTGAACGATAACCAGTTCCCGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGc                                                                                                                                                                                                                                              >  1:384174/1‑250 (MQ=255)          ccTTGAAGAAACGCTGGAAGCCACGGTCAGCCGCAGTAAACACGCCTTCATGCTCGCCCCACGACGTTGACGGTAAAATAACATCCGCCGCCGACGCGGTTTTGGTCATAAAGATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCg                                                                                                                                                                                                                                        <  2:380561/251‑1 (MQ=255)             tGAAGAAACGCTTGAAGCCACGGTCATCCGCAGTAAACACGCCTTCATGCTCGCCCCACGACGTTGACGGTAAAATAACATCCGCCGCCGACGCGGTTGTGGTCATAAAGATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTCCTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCt                                                                                                                                                                                                                                      <  2:306426/250‑1 (MQ=255)               aagaaACGCTGGAAGCCACGGTCAGCCGCAGTAAACACGCCTTCATGCTCGCCCCACGACGTTGACGGTAAAATAACATCCGCCGCCGACGCGGTTTTGGTCATAAAGATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGAt                                                                                                                                                                                                                                   >  1:87643/1‑251 (MQ=255)                      gCTGGAAGCCACGGTCAGCCGCAGTAAACACGCCTTCATGCTCGCCCCACGACGTTGACGGTAAAATAACATCCGCCGCCGACGCGGTTTTGGTCATAAAGATGTCCTGAACGATAACCAGTTCCAGAACTTCAAAGGCTTTACGTACGGCCGCCAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCACATAATGTACGCGGCACGCACTTCGCAATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGAt                                                                                                                                                                                                                              >  2:297726/1‑249 (MQ=255)                          gAAGCCACGGTCAGCCGCAGTAAACACGCCTTCATGCTCGCCCCACGACGTTGACGGTAAAATAACATCCGCCGCCGACGCGGTTTTGGTCATAAAGATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGGACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCTCTGATGCGATAGCCgg                                                                                                                                                                                                                        <  2:33105/251‑1 (MQ=255)                            aGCCACGGTCAGCCGCAGTAAACACGCCTTCATGCTCGCCCCACGACGTTGACGGTAAAATAACATCCGCCGCCGACGCGGTTTTGGTCATAAAGATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTa                                                                                                                                                                                                                      <  2:315966/251‑1 (MQ=255)                                aCGGTCAGCCGCAGTAAACACGCCTTCATGCTCGCCCCACGACGTTGACGGTAAAATAACATCCGCCGCCGACGCGGTTTTGGTCATAAAGATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATgcg                                                                                                                                                                                                                  >  1:208208/1‑251 (MQ=255)                                  ggTCAGCCGCAGTAAACACGCCTTCATGCTCGCCCCACGACGTTGACGGTAAAATAACATCCGCCGCCGACGCGGTTTTGGTCATAAAGATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCt                                                                                                                                                                                                                >  1:176102/1‑251 (MQ=255)                                         cGCAGTAAACACGCCTTCATGCTCGCCCCACGACGTTGACGGTAAAATAACATCCGCCGCCGACGCGGTTTTGGTCATAAAGATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGc                                                                                                                                                                                                         >  1:149335/1‑251 (MQ=255)                                         cGCAGTAAACACGCCTTCATGCTCGCCCCACGACGTTGACGGTAAAATAACATCCGCCGCCGACGCGGTTTTGGTCATAAAGATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGATCGCTGATGCGATAGCCGGTATGCGCTTGCAGGc                                                                                                                                                                                                         >  1:202833/1‑251 (MQ=255)                                             gTAAACAAGCGTTCGTGATCGCCCCACGCCGGTGACGGTAAAATAACATCCGCCGCCGACGCGATCTTGGTAATAAAGATATCCTGACAAATAACCACTACCAGCTCTTAGACGGCTTTACGTGCTGCAAGCAGATCCGAGTCAGCTTGTAGCGGATCTTAGCCCATAATGTAAGCGGCACGCACTTCGCCATGTGCTGTGCGGTCCGGCAGTTCGCTGCTTCGATAGCCGGTATCCCCTTGCGGGCTTTc                                                                                                                                                                                                     >  2:207878/1‑251 (MQ=255)                                                                   cccACGACGTTGACGGTAAAATAACATCCGCCGCCGACGCGGTTTTGGTCATAAAGATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACACCCCAGGCTTTGGCGAAt                                                                                                                                                                               >  2:262621/1‑251 (MQ=255)                                                                             tGACGGTAAAATAACATCCGCCGCCGACGCGGTTTTGGTCATAAAGATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGt                                                                                                                                                                     >  1:138655/1‑251 (MQ=255)                                                                                 ggaaaaaTAACCGCCGTCGCCGACGCGGTTTTTGTCAAAAAGATGTCCTTAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTGCTGCCGACAGCTCCGCGTCAGTTTGTGGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCTATGCGCTGTGCGGTGCGGCAGCTCGCTTATGCGATAGCCGGTATGCGCTGGCAGTCTTTCCACGCCTCAGGCTTTGGCGAATTTCTCGCGGTTAGc                                                                                                                                                                 <  2:310492/248‑1 (MQ=255)                                                                                   tAAAATAACATCCGCCGCCGACGCGGTTTTGGTCATAAAGATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGCTAGCCGGTATGCGCTGGCAGGCTTTCCCCGCCCCAGGCGTTGGCGAATTTCTCGCGGTTAGCCg                                                                                                                                                               >  2:104431/1‑251 (MQ=255)                                                                                    aaaaTAACATCCGCCGCCGACGCGGTTTTGGTCATAAAGATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGGCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCg                                                                                                                                                               <  1:204388/250‑1 (MQ=255)                                                                                        tAACATCCGCCGCAGACGCGGTTTTGGTCATAAAGATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCt                                                                                                                                                          >  1:308516/1‑251 (MQ=255)                                                                                         aaCATCCGCCGCCGACGCGGTTTTGGTCATAAAGATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGAGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCtt                                                                                                                                                         >  1:338236/1‑251 (MQ=255)                                                                                            aTCCGCCGCCGACGCGGTTTTGGTCATAAAGATGTCCTGAACGATAACCAGTTCCAGATCTTCAACGGCTTTCCGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCCCTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTc                                                                                                                                                       >  1:196841/1‑250 (MQ=255)                                                                                                     cgACGCGGTTTTGGTCATAAAGATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGt                                                                                                                                             >  1:60624/1‑251 (MQ=255)                                                                                                     cgACGCGGTTTTGGTCATAAAGATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTCTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGCTCTTTCACGTACTGGt                                                                                                                                             >  1:41355/1‑251 (MQ=255)                                                                                                           ggTTTTGGTCATAAAGATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCAGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCg                                                                                                                                       >  1:134436/1‑251 (MQ=255)                                                                                                             ttttGGTCATAAAGATGTCCTGAACGATAACCAGTTCCAGATCTTCACAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCgg                                                                                                                                      >  1:392939/1‑250 (MQ=255)                                                                                                                    cATAAAGATGTCCTGAACGATAACCAGTGACAGCTCTTCAAAGCCTATACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTCCGCGGCACGCACTTCGCCATGCGCTGCGCGGTGCGGCAGCTCGCTGATGCGCTCGCCGGCATGCGCTGGCAGGCTTTCCACGCCCCAGACTTTAGCGAATTTCTCGCGGTGAGCCGGCTCTTGCGCGACCTGGGCTCCCGGATACGTAt                                                                                                                              >  2:61655/1‑251 (MQ=255)                                                                                                                    aaTAAAGATGTCTTGAAAGATAACCAGGTACAGAGTGTCAAAGGCTGTACGTAATGCCGACAGCTCCGCGTCAGGTTGGAGCGGATCTTCGCCAAGAATGTACGCGGCACGCAAGTCGCCATGCGCTGTGAGGTGCGGCAGCGCGCTGATGCGATAGCCGGTATGCGCTGGAAGGCTTTCCACGCCCCAGGCTTTGGCGAAGTGCGCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTAt                                                                                                                              <  2:93109/250‑1 (MQ=255)                                                                                                                     atacaGATGTCCGTAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTTCGTCCGACAGCTCAGCGTCAGGTTGTGGCGGATCTTCGCCCATAATGTACGCGTCACGCACGTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACGGGTATCCCGGATACGTATc                                                                                                                             <  2:276270/247‑1 (MQ=255)                                                                                                                     aTAAAGATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTTCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATc                                                                                                                             <  1:345104/251‑1 (MQ=255)                                                                                                                     aTAAAGATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATc                                                                                                                             >  1:129504/1‑251 (MQ=255)                                                                                                                                 tGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGCATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGcc                                                                                                                 >  1:170962/1‑251 (MQ=255)                                                                                                                                  gAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATcc                                                                                                                            <  1:229866/239‑1 (MQ=255)                                                                                                                                  gAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATcc                                                                                                                            >  2:229866/1‑239 (MQ=255)                                                                                                                                    aCGATACCCAGTTCCAGCTCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTCCGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCat                                                                                                              <  2:308516/251‑1 (MQ=255)                                                                                                                                     cGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCat                                                                                                              >  1:131777/1‑250 (MQ=255)                                                                                                                                           ccAGTTCCAGATCGTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTGCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAg                                                                                                       <  1:262621/251‑1 (MQ=255)                                                                                                                                            cAGTTCCAGATCTTCAAAGGCTTTACGTAATGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAgg                                                                                                      >  1:293116/1‑251 (MQ=255)                                                                                                                                                  cAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCt                                                                                                <  1:340535/251‑1 (MQ=255)                                                                                                                                                      tCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATa                                                                                                                                                                                                                             <  1:361279/122‑1 (MQ=255)                                                                                                                                                      tCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATa                                                                                                                                                                                                                             >  2:361279/1‑122 (MQ=255)                                                                                                                                                          caatggCTTTCCGTACCGCCGCCAGCTCCGCGTCTGGTTGTAGCGGTCCCTCGCCCATAATTGACGCGGCACGCACTTCGTCATGCGCTGGGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCTAAGTTCTCGCGGTTAGCCGTATCTTTCACGTACTGGTATCCCGGATACGTAGCCGGCGGCGCGCCCATATCTCAGGCACCCGTAACgttg                                                                                        <  2:48497/247‑1 (MQ=255)                                                                                                                                                            aaGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACgttgtt                                                                                      >  1:313096/1‑251 (MQ=255)                                                                                                                                                                     cGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACg                                                                             >  1:312912/1‑251 (MQ=255)                                                                                                                                                                      gTACTGCCGACAGCTCCGCGGCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGACGAATTTCTCGCGGTTAGCCGGATCTTTCACGTAATGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCCCCCTGAACGTTGTTCTGACCACg                                                                             >  1:82206/1‑250 (MQ=255)                                                                                                                                                                               aCAGCTACGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGtt                                                                    >  1:125488/1‑250 (MQ=255)                                                                                                                                                                                     ccGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACAcca                                                             <  2:129504/251‑2 (MQ=255)                                                                                                                                                                                        cgTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGc                                                           >  1:253203/1‑250 (MQ=255)                                                                                                                                                                                           cAGTTTGGAGCGGATCTTCGCCCATAATGTACGCGGCACGCCCTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGGATGCGCTGGCTGGCTGTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTAGCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTAACCACGAACCGGGTTAACACCCGCATg                                                        <  2:205007/250‑1 (MQ=255)                                                                                                                                                                                                     cGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCgc                                             >  2:129991/1‑250 (MQ=255)                                                                                                                                                                                                        aTCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTCTGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCCTATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCGAGGt                                          >  1:187158/1‑251 (MQ=255)                                                                                                                                                                                                                 cATAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCGAGGTTACCGGTCa                                 <  2:313096/251‑1 (MQ=255)                                                                                                                                                                                                                  aTAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCGAGGTTACCGGTc                                  >  1:246754/1‑249 (MQ=255)                                                                                                                                                                                                                   tAATGTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGTTACGTATCCGGCAGCGCGCCCATATCGCAGCCCCCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATCCGCCTTACCGAGGTTACCGGTCAGc                               >  2:73085/1‑251 (MQ=255)                                                                                                                                                                                                                      ggTACGCGGCACGCACTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCGAGGTTACCGGTCAGCATc                            <  2:100452/250‑1 (MQ=255)                                                                                                                                                                                                                                 cgcccTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCGAGGTTACCGGTCAGCATCGCGAGGCTGGt                 <  2:250957/247‑1 (MQ=255)                                                                                                                                                                                                                                  gcacTTCGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCGAGGTTACCGGTCAGCATCGCGAGGCTGGTc                >  1:38985/1‑251 (MQ=255)                                                                                                                                                                                                                                    ccTTCGCCATGCGCTGGGCGGTGCGGCAGCTCGCTGATGCGTTAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCGAGGTTACCGGTCAGCATCGCGAGGCTGGGCag              <  2:297047/250‑1 (MQ=255)                                                                                                                                                                                                                                     cgtCGCCCTGCGCTGTGCGGTGCGTCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCCGGCTTTCCACGCCCCGGGCTTTGGCGAATTTCTCGCGGTTAGCCGGTTCTTTCACGTACTGGTATCCCGGATCCGTATCCGGCAGCGCGCTCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCTCATGCGGCTTACCGAGTTTTCCGGTCAGCAGCGCGAGGCTGTTCaga             <  2:296854/249‑1 (MQ=255)                                                                                                                                                                                                                                       tCGCCATGCGCTGGGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCGAGGTTACCGGTCAGCATCGCGAGGCTGGTCagag            <  1:207878/250‑1 (MQ=255)                                                                                                                                                                                                                                       tCGCCATGCGCGGTGCGGTGTGGCAGCTCGCTGATGCGATAGCCGGTAGTCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACAGAGGTTACCGGTCAGCATCGCGAGGCTGGTCag              <  2:279915/248‑1 (MQ=255)                                                                                                                                                                                                                                        cGCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCACGCGGCTTACCGAGGTTACCGGTCAGCATCGCGGGGCTGGTCagcga           >  2:170729/1‑247 (MQ=255)                                                                                                                                                                                                                                        cGCCATGCGCTGTGCGGGGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCGTACCGAGGTTACCGGTCAGCATCTCGAGTCTGGTCAgaga           <  2:394684/250‑1 (MQ=255)                                                                                                                                                                                                                                         gCCATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCGAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAAc         >  1:366245/1‑251 (MQ=255)                                                                                                                                                                                                                                          ccATGCGCTGTGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCGAGGTTACCGGTCAGCATCGCGAGGCTGGGCAGAGAACg        >  1:339636/1‑251 (MQ=255)                                                                                                                                                                                                                                               cgctgtgcGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCGAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACgg   >  1:45892/1‑251 (MQ=255)                                                                                                                                                                                                                                               cgctgtgcGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCGAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACgg   <  1:104431/251‑1 (MQ=255)                                                                                                                                                                                                                                               cgctgtgcGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCACCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCGAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACgg   <  2:9172/251‑1 (MQ=255)                                                                                                                                                                                                                                                gctgtgcGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCGAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGt  <  1:13725/251‑1 (MQ=255)                                                                                                                                                                                                                                                gctgtgcGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCGAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGt  >  1:122469/1‑251 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                       GTTCAACCGCCTTGAAGAAACGCTGGAAGCCACGGTCAGCCGCAGTAAACACGCCTTCATGCTCGCCCCACGACGTTGACGGTAAAATAACATCCGCCGCCGACGCGGTTTTGGTCATAAAGATGTCCTGAACGATAACCAGTTCCAGATCTTCAAAGGCTTTACGTACTGCCGACAGCTCCGCGTCAGTTTGTAGCGGATCTTCGCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGCGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCGAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGT  >  NC_000913/4297956‑4298445

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