Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,064,900	G→C	12.4%	S33S (TCG→TCC)	scpC →	propionyl‑CoA:succinate CoA transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,064,900	0	G	C	12.4%	94.8 / 2.3	50	S33S (TCG→TCC)	scpC	propionyl‑CoA:succinate CoA transferase
Reads supporting (aligned to +/- strand):  ref base G (27/14);  new base C (3/4);  total (32/18)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.62e-01

TATTCAGGGGATGCGCCGCGCGGTGTATGACAGCGAAGATTACCAGGAAGGGATGAACGCTTTCCTCGAAAAACGTAAACCTAATTTCGTT‑GGTCATTAATCCCTGCGAACGAAGGAGTAAAAATGGAAACTCAGTGGACAAGGATGACCGCCAATGAAGCGGCAGAAATTATCCAGCATAACGACATGGTGGCATTTAGCGGCTTTACCCCGGCGGGTTCGCCGAAAGCCCTACCCACCGCGATTGCCCGCAGAGCTAACGAACAGCATGAGGCCAAAAAGCCGTATCAAATTCGCCTTCTGACGGGTGCGTCAATCAGCGCCGCCGCTGACGATGTACTTTCTGACGCCGATGCTGTTTCCTGGCGTGCGCC‑ATATCAAACATCGTCCGGTTTACGTAAAAAGATCAATCAGGGCGCGGTGAGTTTCGTTGACCTGCATTTGAG  >  NC_000913/3064679‑3065124                                                                                                                                                                                                                               |                                                                                                                                                                                                                                  tATTCAGGGGATGCGCCGCGCGGCGTATGACAGCGAAGATTACCAGGAAGGGATGAACGCTTTCCTCGAAAAACGTAAACCTAATTTCGTT‑GGTCATTAATCCCTGCGAACGAAGGAGTAAAAATGGAAACTCAGTGGACAAGGATGACCGCCAATGAAGCGGCAGAAATTATCCAGCATAACGACATGGTGGCATTTAGCGGCTTTACCCCGGCGGGGTCGCCGaa                                                                                                                                                                                                                              >  1:2826/1‑227 (MQ=255)   gcaaGGAGAGGCGCCGCGCGGTGTAAGACAGCAAAAATTACGAGGAAGGGATGAACGTTTCCCTCGAAAAACGAAAACCTAATTGCGTT‑GGTCATTAATCCCGGCGAACGAAGGAGTAAAAATGAAAACTCAGTGAAAAAGGATGACCGCCAATGAAGCGGCAGAAATTATCCAGCATAACGACATGGTGGCATTTAGCGGCTTTACCCCGGCGGGTTCGCCGAAAGc                                                                                                                                                                                                                           <  1:173537/225‑1 (MQ=255)           aTGCGCCGCGCGGTGTATGACAGCGAAGATTACCAGGAAGGGATGAACGCTTTCCTCGAAAAACGTAAACCTAATTTCGTT‑GGTCATTAATCCCTGCGAACGAAGGAGTAAAAATGGAAACTCAGTGGACAAGGATGACCGCCAATGAAGCGGCAGAAATTATCCAGCATAACGACATGGTGGCATTTAGCGGCTTTACCCCGGCGGGTTCGCCCAAAGCCTTAccc                                                                                                                                                                                                                    >  1:611983/1‑227 (MQ=255)           aTGCGCCGCGCGGTGTATGACAGCGAAGATTACCAGGAAGGGATGAACGCTTTCCTCGAAAAACGTAAACCTAATTTCGTT‑GGTCATTAATCCCTGCGAACGAAGAAGTAAAAATGGAAACTCAGTGGACAAGGATGACCGCCAATGAAGCGGCAGAAATTATCCAGCATAACGACATGGAGGCATTTAGCGGCTTTACCCCGGCGGGGTCGCCGAAACCCCTAccc                                                                                                                                                                                                                    >  1:418051/1‑227 (MQ=255)            tGCGCCGCGCGGTGTATGACAGCGAAGATTACCAGGAAGGGATGAACGCTTTCCTCGAAAAACGTAAACCTAATTTCGTT‑GGTCATTAATCCCTGCGAACGAAGGAGTAAAAATGGAAACTCAGTGGACAAGGATGACCGCCAATGAAGCGGCAGCAATTATCCACCATAACGACATGGTGGCATTTAGCGGCTTTCCCCCGGCGGGTTACCCGAAACCCCTACCCAc                                                                                                                                                                                                                  >  2:552774/1‑228 (MQ=255)            tGCGCCGCGCGGTGGATGACAGCGAAGCTTACCAGGAAGGGATGAACGCTTTCCTCGAAAAACGTAAACCTAATTTCGTT‑GGACATTAATCCCTGCGAACGAAGGAGTCAAAATGGTAACTCAGTGGACAAGGATGACCGCCAATGAAGCGGCAGAAATTATCCAGCACAACGACATGGTGGCATTTAGCGGCTTTACCCCGTCAGGTTCTACCAAACCCCTCCCCa                                                                                                                                                                                                                   >  1:645252/1‑227 (MQ=255)                                  cGAAGATTACCAGGAAGGGATGAACGCTTTCCTCGAAAAACGTAAACCTAATTTCGTT‑GGTCATTAATCCCTGCGACCGAAGGAGTAAAAATGGAAACTCAGTGGACAAGGATGACCGCCAATGAAGCGGCAGAAATTATCCAGCATAACGACATGGTGGCATTTAGCGGCTTTACCCCGGCGGGGTGGCCGAAAGACCCACCCCACCCCATGTCCCCCAGAGCCAAc                                                                                                                                                                                            >  2:482675/1‑228 (MQ=255)                                  cGAAGATTACCAGGAAGGGATGAACGCTTTCCTCGAAAAACGTAAACCTAATTTCGTT‑GGTCATTAATCCCTGCGAACGAAGGCGTAAAAATGGAAACTCAGTGGACAAGGATGACCGCCAATGAAGCGGCAGAAATTATCCAGCATAACGACATGGTGGAATTTAGCGGCTTTACCCCGGCGGGTTCGCCGAAAGCCCTACCAACCGCAATTAACCGCAGGGCTAAc                                                                                                                                                                                            >  1:438734/1‑228 (MQ=255)                                  cGAAGATTACCAGGAAGGGATGAACGCTTTCCTCGAAAAACGTAAACCTAATTTCGTT‑GGTCATTAATCCCTGCGAACGAAGGAGTAAAAATGGAAACTCAGTGGACAAGGATGACCGCCAATGAAGCGGCAGAAATTATCCAGCATAACGACATGGTGGCATTTAGCGGCTTTACCCCGGCGGGGTCGCCGCAAGCCCTACCCACCGCGATTGCCCGCAGGGCTaag                                                                                                                                                                                            >  1:350706/1‑227 (MQ=255)                                  cGAAGATTACCAGGAAGGGATGAACGCTTTCCTCGAAAAACGTAAACCTAATTTCGTT‑GGTCATTAATCCCTGCGAACGAAGGAGTAAAAATGGAAACTCAGTGGACAAGGATGACCGCCAATGAAGCGGCAGAAATTATCCAGCATAACGACATGGTGGCATTTAGCGGCTTTACCACGGCGGGTTCTCCGAAAACCCTACCCACCGCGATTGCCCGCAAGGCTAAc                                                                                                                                                                                            >  1:336144/1‑228 (MQ=255)                                  cGAAGATTACCAGGAAGGGATGAACGCTTTCCTCGAAAAACGTAAACCTAATTTCGTT‑GGTCATTAATCCCTGCGAACGAAGGAGTAAAAATGGAAACTCAGTGGACAAGGATGACCGCCAATGAAGCGGCAGAAATTATCCAGCATAACGACATGGTGGAATTTAGCGGCTTTACCCCGGCGGGTTTGCCGAAAGCCCCACCCCCCGCGACTGCCCGCAGAGCTaa                                                                                                                                                                                             >  1:406500/1‑227 (MQ=255)                                  cGAAGATTACCAGGAAGGGATGAACGCTTTCCTCGAAAAACGTAAACCTAATTTCGTT‑GGTCATTAATCCCTGCGAACGAAGGAGTAAAAATGGAAACTCAGTGGACAAGGATGACCGCCAATGAAGCGGCAGAAATTATCAAGCATAACGACATGGTGGCATTTAGCGGCTTTACCCCGGCGGGTTCGCCGAAAGCCCTAACCACCCACATTGCCCGccagactaac                                                                                                                                                                                            >  2:181601/1‑219 (MQ=255)                                  cGAAGATTACCAGGAAGGGATGAACGCTTTCCTCGAAAAACGTAAACCTAATTTCGTT‑GGTCATTAATCCCTGCGAACGAAGGAGTAAAAATGGAAACTCAGTGGACAAGGATGACCGCCAATGAAGCGGCAGAAATTATACAGCATAACGACATGGGGGCATTTAGCGGCTTTACCCCAGCGGGTTCCCCGAAAGCCCTCCCCACCGCGATTGCCCGccacgcaaac                                                                                                                                                                                            >  1:321690/1‑219 (MQ=255)                                    aaGATTACCAGGAAGGGATGAACGCTTTCCTCGAAAAACGTAAACCTAATTTCGTT‑GGTCATTAATCCCTGCGAACGAAGGAGTAAAAATGGAAACTCAGTGGACAAGGATGACCGCCAATGAAGCGGCAGAAATTATCCAGCATAACGACATGGTGGCATTTAGCGGCTTTACCCCGGCGGGTTCGCCGAAAGCCCTACCCACCGCGATTGCCCGCAGAGCTAAccg                                                                                                                                                                                          >  1:226181/1‑226 (MQ=255)                                    aaGATTACCAGGAAGGGATGAACGCTTTCCTCGAAAAACGTAAACCTAATTTCGTT‑GGTCATTAATCCCTGCGAACGAAGGAGTAAAAATGGAAACTCAGTGGACAAGGATGACCGCCAATGAAGCGGCAGAAATTATCCAGCATAACGACATGGTGGCATTTAGCGGCTTTACCCCGGCGGGTTAGCCGCAAGCCCTACCCACCGCGGTTGCCCGCAGAGCTCACg                                                                                                                                                                                           >  1:574500/1‑227 (MQ=255)                                                 aGGGATGAACGCTTTCCTCGAAAAACGTAAACCTAATTTCGTT‑GGTCATTAATCCCTGCGAACGAAGGAGTAAAAATGGAAACTCAGTGGACAAGGATGACCGCCAATGAAGCGGCAGAAATTATCCAGCATAACGACATGGTGGCATTTAGCGGCTTTACCCCGGCGGGTTCGCCGAAAGCCCTACCCACCGCGCATGCCCGCCGAGCTAACGAACAGCATGAGGc                                                                                                                                                                              >  1:455797/1‑227 (MQ=255)                                                  ggatgtaaccttttCATCAAAAAACG‑AAACCATATTTCGTTGGGTCAGTACACCCTGAAAACAAGGAAGTAAAAATGGAACCACAGGGAAAAAGGAGGACCGCCAATGAAGCGGCAGAAATTATCCAGCATAACGACATGGTGGCATTTAGCGGCTTTACCCCGGCGGGTTCGCCGAAAGCCCTACCCACCGCGATTGCCCGCAGAGCTAACGAACAGCATGAGGCCa                                                                                                                                                                            <  2:604777/217‑1 (MQ=255)                                                                                  tAATTTCGTT‑GGTCATTAATCCCTGCGAACGAAGGAGTAAAAATGGAAACTCAGTGGACAAGGATGACCGCCAATGAAGCGGCAGAAATTATCCAGCATAACGACATGGTGGCATTTAGCGGCTTTACCCCGGCGGGTTCGCCGAAAGCCCTACCCACCGCGATTGCCCGCAGAGCTAACGAACAGCATGAGGCCAAAAAGCCGTATCAAATTCGCCTTCTGACggg                                                                                                                                             >  1:571549/1‑227 (MQ=255)                                                                                          tt‑GGTCATTAATCCCTGCGAACGAAGGAGTAAAAATGGAAACTCAGTGGACAAGGATGACCGCCAATGAAGCGGCAGAAATTATCCAGCATAACGACATGGTGGCATTTAGCGGCTTTACCCCGGCGGGTTCGCCGAAAGCCCTACCCACCGCGATTGCCCGCAGAGCTAACGAACAGCATGAGGCCAAAAAGCCGTATCAAATTCGCCTTCTGAAGGGTGCggaaa                                                                                                                                     >  1:410479/1‑223 (MQ=255)                                                                                           t‑GGTCATTAATCCCTGCGAACGAAGGAGTAAAAATGGAAACTCAGTGGACAAGGATGACCGCCAATGAAGCGGCAGAAATTATCCAGCATAACGACATGGTGGCATTTAGCGGCTTTACCCCGTCGGGTTCGCCGAAAGCCCTACCCACCGCGATTTCACGCAGGGCTCACGCACAGCATGAGGCCAAAAAGCCGTATCAAATTGCCCTTCTGACGGGTGAGTCAATc                                                                                                                                   >  1:169643/1‑228 (MQ=255)                                                                                           t‑GGTCATTAATCCCTGCGAACGAAGGAGTAAAAATGGAAACTCAGTGGACAAGGATGACCGCCAATGAAGCGGCAGAAATTATCCAGCATAACGACATGGTGGCATTTAGCGGCTTTACCCCGGCGGGTTCGCCGAAAGCCCTACCCACCGCGATTGCCCGCAGAGCTAACGAACAGCATGAGGCCAAAAAGCCGTATCAAATTCGCCTTCTGCCGGGGGCGTCAATc                                                                                                                                   >  1:371059/1‑228 (MQ=255)                                                                                                        ccTGCGAACGAAGGAGTAAAAATGGAAACTCAGTGGACAAGGATGACCGCCAATGAAGCGGCAGAAATTATCCAGCATAACGACATGGTGGCATTTAGCGGCTTTACCCCGGCGGGTTCGCCGAAAGCCCTACCCACCGCGATTTCCCGCAGAGCTAACGAACAGCATGAGGCCAAAAAGCCGTATCAAATTCGCCTTCTTACGGGTGCGTCAATCAGCCCCGCCGCt                                                                                                                       >  1:499532/1‑228 (MQ=255)                                                                                                           gCGAACGAAGGAGTAAAAATGGAAACTCAGTGGACAAGGATGACCGCCAATGAAGCGGCAGAAATTATCCAGCATAACGACATGGTGGCATTTAGCGGCTTTACCCCGGCGGGTTCGCCGAAAGCCCTACCCACCGCGATTGCCCGCAGAGCTAACGAACAGCATGAGGCCAAAAAGCCGGATCAAATTCGCCTTCTGACGGGTGCGCCACTCAGCTCCGCCCCTg                                                                                                                      >  1:327335/1‑226 (MQ=255)                                                                                                                 aaagaaaaaaaaaGGGAAACCCGGGGGACAGGGAGTGCCGCCAATGAAGCGGCAGAAATTATCCAGCATAACGACATGGTGGCATTTAGCGGCTTTACCCCGGCGGGTTCGCCGAAAGCCCTACCCACCGCGATTGCCCGCAGAGCTAACGAACAGCATGAGGCCAAAAAGCCGTATCAAATTCGCCTTCTGACGGGTGCGTCAATCAGCGCCGCCGCTGACGATGTa                                                                                                              <  1:482675/220‑1 (MQ=255)                                                                                                                                                    cccccccAATGAGGCGGCAAAAATTATCCAGCATACCGACATGGTGTCATTTAGCGGCTTTACCCCGGGGGGTTCGCCGAAAGCCCTACCCACCGCGATTGCCCGCAGAGCTAACGAACAGCATGAGGCCAAAAAGCCGTATCAAATTCGCCTTCTGACGGGTGCGTCAATCAGCGCCGCCGCTGACGATGTACTTTCTGACGCCGATGCTGTTTCCTGGCGTGCGcc                                                                           <  2:571549/223‑1 (MQ=255)                                                                                                                                                    cacccgaaaTGAAGGGGAAAAAATTTCCCATCATAGCGACATGGGGGCATTTGGCGGCCTGACCCCGGCGGTTTACCAGAAAGCCCCACCCCCCGCAATTGCCCCCCGAGTTAACGAACAGCATGGGGCCAAAACGCCGTACCAAATTCGCCTTCTGACGGGTGCGTCAATCAGCGCCGCCGCTGACGATGTACTTTCTGACGCCGATGCTGTTTCCTGGCGTGCGcc                                                                           <  2:227137/221‑1 (MQ=255)                                                                                                                                                     cccctcAATGAGGGGGGAGAAATGTCCCCGCATAAAGACAGGGTGCCATTTAGCGGCTTTCCCCCGGCGGGTTCGCCGAAAGCCCTACCCACCGGGATTGCCCGCAGAGCTAACGAACAGCATGAGGCCAAAAAGCCGTATCAAATTCGCCTTCTGACGGGTGCGTCAATCAGCGCCGCCGCTGACGATGTACTTTCTGACGCCGATGCTGTTTCCTGGCGTGCGcc                                                                           <  2:321690/222‑1 (MQ=255)                                                                                                                                                     cccGCCAATGAACGGGCAGAAATTCTCCGGCAGAACCACATGGTGGCCTTTGGCGGTTTTACCCCGGGCGGTTCGCCGAAAGCCCCACCCACCGCGATTGCCCGCAGCGCTAACGAACAGCATGAGGCCAAAAAGCCGTATCAAATTCGCCTTCTGACGGGTGCGTCAATCAGCGCCGCCGCTGACGATGTACTTTCTGACGCCGATGCTGTTTCCTGGCGTGCGcc                                                                           <  2:350706/226‑1 (MQ=255)                                                                                                                                                     cccGCCAAAGACGGGGCAAAAATGACCCAGCAGAACGAAAGGGGGGAATGTAGCGGCTTTACCCAGGCGGGTGCGCAAAAAGCCCTCCCCACCGCGATTGCCCGCAGAGCTAACGAACAGCATGAGGCCAAAAGGCCGTATCAAATTCGCCTTCTGACGGGTGCGTCAATCAGCTCCGCCGCTGACGATGTACTTTCTGACGCCGATGCTGTTTCCTGGCGTGCGcc                                                                           <  2:523668/226‑1 (MQ=255)                                                                                                                                                            aTGAAGCGGCAGAAATTATCCAGCATAACGACATGGTGGCATTTAGCGGCTTTACCCCGGCGGGGTCGCCGAAAGCCCTACCCCCCGCGATTGCCCGCAGAGCTAACGAACAGCATGAGGCCCAACAGCCGTATCCAATTCGCCTTCTCACGGGAGCGGCAATCAGCGCCGCCGCTGCCGCTGGTCTTTCTATCCCCGCGGCTGTTTCCCGGCGTGCGCC‑CAATCaaa                                                                  >  2:279588/1‑228 (MQ=255)                                                                                                                                                              taaggtggaaaaaTTAACCATGATAAGGGTATGGTGGCATTTAGCGGCTTTACCCCGGCGGGTTCGCCGAAAGCCCTACCCCCCGCGATTGCCCGCAGAGCTAACGAACAGCATGAGGCCAAAAAGCCGTATCAAATTCGCCTTCTGACGGGTGCGTCAATCAGCGCCGCCGCTGACGATGTACTTTCTGACGCCGATGCTGTTTCCTGGCGTGCGCC‑ATACCAAAc                                                                 <  2:226181/217‑1 (MQ=255)                                                                                                                                                                    gCAGAAATTATCCAGCATAACGACATGGTGGCATTTAGCGGCTTTACCCCGGCGGGTTCGCCGAAAGCCCTACCCACCGCGATTGCCCGCAGAGCTAACGAACAGCATGAGGCCAAAAAGCCGTATCAAATTCGCCTTCTGACGGGTGCGTCAATCAGCGCCGCCGCTGACGATGTACTTTCTGACGCCGAGGCTGTTTCCTGGCGTGCGCC‑ATATCAAACATCGTcc                                                          >  1:441637/1‑228 (MQ=255)                                                                                                                                                                         aaTTATCCCGCCTAACGTATTGGTGGCTTTTAGAGCTTTTCCCCCGGAGGGTTCCCCGAAACCCCTACCCACCGCGATTGCCCGCAGATCTAACGAACAGCATGAGGCCAAAAAGCCGTATCAAATTCGCCTTCTGACGGGTGCGTCAATCAGCGCCGCCGCTGACGATGTACTTTCTGACGCCGATGCTGTTTCCTGGCGTGCGCC‑ATATCAAACATCGTCCGGttt                                                     <  2:264494/228‑1 (MQ=255)                                                                                                                                                                                caagaaaaCGAGAGGGTGGCTTTAAGGGGCTTTCCCACGGGGGATGTCCCAAAAGCCCTCCCCACGGGGATTGCCCGCAGGGCTAACGAACAGCATGAGGCCAAAAAGCCTTATCAAATGCGCCTTCTGACGGGTGCGTCAATAAGCGCCGCCGCTGACGATGTACTTTCTGACGCCGATGCTGTTTCCTGGCGTGCGCC‑ATATCAAACATCGTCCGGTTTACGTaa                                               <  2:657681/221‑1 (MQ=255)                                                                                                                                                                                  gCATAACGACATGGTGGCATTTAGCGGCTTTACCCCGGCGGGTTCGCCGAAAGCCCTACCCACCGCGATTGCCCGCAGAGCTAACGAACAGCATGAGGCCAAAAAGCCGTATCAAATTCGCCTTCTGACGGGTGCGTCAATCAGAGCCGCCGCTGACGATGTACTTTCTGACGCCGATGCTGTTTCCTGGCGTGCGCC‑CTATCAAACATCGTCCGGTTTACGTaaaaa                                            >  1:490868/1‑228 (MQ=255)                                                                                                                                                                                   cATAACGACATGGTGGCATTTAGCGGCTTTACCCCGGCGGGTTCGCCGAAAGCCCTACCCACCGCGATTTCCCGCCGAGCTAACGCACAGCATGAGGCCAAAAAGCCGTATCAAATTCGCCTTCTGACGGGTGCGTCAATCCGCGCCCCCGCTGACGATGTACTTTCTGACCCCGATGCTGTTTCCTGGCGTGCGCC‑ATACCCAACATCGCCCGGTTTCCCTaaaaaa                                           >  2:539676/1‑227 (MQ=255)                                                                                                                                                                                   cATAACGACATGGTGGCATTTAGCGGCTTTACCCCGGCGGGTTCGCCGAAAGCCCTACCCACCGCGATTGCCCGCAGAGCTAACGAACAGCATGAGGCCAAAAAGCCGTATCAAATTCGCCTTCTGACGGGGGCGTCAATCAGCGCCGCCGCTGACCATGGACCTTCTTACGCCGGAGCCGTTTCCCGGCGGGCGCCAATTTCAAAAAAC‑TCCGGGTTAAGaaaaaaa                                           >  2:304522/1‑227 (MQ=255)                                                                                                                                                                                   cATAACGACATGGTGGCATTTAGCGGCTTTACCCCGGCGGGTTCCCCGAAAGCCCTACCCACCGCGATTGCCCGCAGAGCTAACGAACAGCATGAGGCCAAAAAGCCGTATCAAATTCGCCTTCTGACGGGTGCGTCAATCAGCGCCGCCGCTGACGATGTACTTTCTGACGCCGAt                                                                                               >  1:181090/1‑177 (MQ=255)                                                                                                                                                                                    aaaaCGACATGGTGGCTTTTAGCGGCTTTCCCCGGGCGGGTTCGCCGAAAGCCCTACCCACCGGGATTGCCCGAAGAGCTAACGAACAGCATGAGGCCAAAAAGCCGTATCAAATTCGCCTTCTGACGGGTGCGTCAATCAGCGCCGCCGCTGACGATGTACTTTCTGACGCCGATGCTGTTTCCTGGCGTGCGCC‑ATATCAAACATCGTCCGGTTTACGTaaaa                                             <  1:497935/223‑1 (MQ=255)                                                                                                                                                                                    aTAACGACATGGTGGCATTTAGCGGCTTTACCCCGGCGGGTTCCCCGAAAGCCCTACCCACCGCGATTGCCCGCAGAGCTAACGAACAGCATGAGGCCAAAAAGCCGTATCAAATTCGCCTTCTGACGGGTGCGTCAATCAGCGCCGCCGCTGACGATGTACTTTCTGACGCCGATg                                                                                              <  2:181090/177‑1 (MQ=255)                                                                                                                                                                                             tggtggCATTTAGCGGCTTTACCCCGGCGGGTTCGCCGAAAGCCCTACCCACCGCGATTGCCCGCAGAGCTAACGAACAGCATGAGGCCAAAAAGCCGTATCCAATTCGCCTTCTGACGGGTGCGTCAATCAGCGCCGCCGCTGACGATGTACTTTCTGACGCCGATGCTGTTTCCTGGCGGGCGCC‑ATATCCAACATCGTACGGGTTACGGAAAAAGAAAAATCAgg                                 >  2:175667/1‑228 (MQ=255)                                                                                                                                                                                             tggtggCATTTAGCGGCTTTACCCCGGCGGGTTCGCCGAAAGCCCTACCCACCGCGATTGCCCGCAGAGCTAACGAACAGCATGAGGCCAAAAAGCCGTATCAAATTCGc                                                                                                                                                        >  1:334713/1‑110 (MQ=255)                                                                                                                                                                                             tggtggCATTTAGCGGCTTTACCCCGGCGGGTTCGCCGAAAGCCCTACCCACCGCGATTGCCCGCAGAGCTAACGAACAGCATGAGGCCAAAAAGCCGTATCAAATTCGCCTTCTGACGGGTGCGTCAATCAGCGCCGCCGCTGACGATGTACTTTCTGACGCCGATGCTGTTTCCTGGCGTGCGCC‑ATATCAAACATCGTCCGGTTTACGTAAAAAAATCAATCAgg                                 >  2:390265/1‑228 (MQ=255)                                                                                                                                                                                              ggtggCATTTAGCGGCTTTACCCCGGCGGGTTCGCCGAAAGCCCTACCCACCGCGATTGCCCGCAGAGCTAACGAACAGCATGAGGCCAAAAAGCCGTATCAAATTCGcc                                                                                                                                                       <  2:334713/110‑1 (MQ=255)                                                                                                                                                                                               gtggCATTTAGGGGCTTTCCCCCGGGGGGTTAGCCGAAAGCCCTACCCACCGCGATTGCCCGCGGGGCTAACGAACAGCAGGAGGCAAAAAAGCCGTATCAAATTCGCCTTCTGACGGGTGCGTCAATCAGCGCCGCCGCTGACGATGTACTTTCTGACGCCGATGCTGTTTCCTGGCGTGCGCC‑ATATCAAACATCGTCCGGTTTACGTAAAAAGATCAATCAGGgc                               <  2:495463/228‑1 (MQ=255)                                                                                                                                                                                               gtggCATTTAGCGGCTTTACCCCGGCGGGTTCGCCGAAAGCCCTACCCACCGCGATTGCCCGCAGAGCTAACGAACAGCATGAGGCCAAAAAGCCGTATCAAATTCGCCTTCTGACGGGTGCGTCAATCAGCGCCGCCGCTGACGATGTACTTTCTGACGCCGATGCTGTTTCCTGGCGTGCGCC‑ATATCAAACATCGTCCGGTTTACGTAAAAAGATCAATCAggg                                <  1:390265/227‑1 (MQ=255)                                                                                                                                                                                                    ctttGGCGGCTTTCCCCCGGCGGGTTCGCAGAAAGCCCTCACCACCCGGATTGCCCGAAGGGCTAACGAACAGCAGGAGGCCAAAAAGCCGTATCAAATTCGCCTTCTGACGGGTGCGTCAATCAGCGCCGCCGCTGACGATGTACTTTCTGACGCCGATGCTGTTTCCTGGCGTGCGCC‑ATATCAAACATCGTCCGGTTTACGTAAAAAGATCAATCAGGGCGCGGt                          <  1:246507/227‑1 (MQ=255)                                                                                                                                                                                                          cGGCTTTACCCCGGCGGGTTCGCCGAAAGCCCTACCCACCGCGATTGCCCGCAGAGCTAACGAACAGCATGAGGCCAAAAAGCCGTATCAAATTCGCCTTCTGACGGGTGCGTCAATCAGCGCCGCCGCTGACGATGTACTTTCTGACGCCGATGCTGTTTCCTGGCGGGCGCC‑ATCTCAAACATCGTCCGGTTTACGTAAAAAGATCAATCAGGGCGCGGGGAGttt                    >  1:174983/1‑228 (MQ=255)                                                                                                                                                                                                                 cccccGGGGGGTTCGCCGAAAGCCCCACCCCCCGCGTTTGCCCGCAGGGCTAACGAACAGCATGAGGCCAAAACGCCCTATCAAATTCGCCTTCTGACGGGTGCGTCAATCAGCGCCGCCGCTGACGATGTACTTTCTGACGCCGATGCTGTTTCCTGGCGTGCGCC‑ATATCAAACATCGTCCGGTTTACGTAAAAAGATCAATCAGGGCGCGGTGAGTTTCGTTGAc             <  2:174983/227‑1 (MQ=255)                                                                                                                                                                                                                            ttCGCCGAAAGCCCTACCCACCGCGATTGCCCGCAGAGCTAACGAACAGCCTGAGGCCCAAACGCCGTATCCAATTCGCCTTCTGACGGGTGCGGCAATCAGCGCCGCCGCTTACGATGTACTTTCTGTCGCCGATTCCGTTTCCTGTGGGGCGCC‑ATATCAAACATCGTCCGGCTTTAGTAACAAGCGCAATCCGGGAGGGGCGAGATTCGTGGACCCCCATTTGAg  >  2:102691/1‑228 (MQ=255)                                                                                                                                                                                                                               |                                                                                                                                                                                                                                  TATTCAGGGGATGCGCCGCGCGGTGTATGACAGCGAAGATTACCAGGAAGGGATGAACGCTTTCCTCGAAAAACGTAAACCTAATTTCGTT‑GGTCATTAATCCCTGCGAACGAAGGAGTAAAAATGGAAACTCAGTGGACAAGGATGACCGCCAATGAAGCGGCAGAAATTATCCAGCATAACGACATGGTGGCATTTAGCGGCTTTACCCCGGCGGGTTCGCCGAAAGCCCTACCCACCGCGATTGCCCGCAGAGCTAACGAACAGCATGAGGCCAAAAAGCCGTATCAAATTCGCCTTCTGACGGGTGCGTCAATCAGCGCCGCCGCTGACGATGTACTTTCTGACGCCGATGCTGTTTCCTGGCGTGCGCC‑ATATCAAACATCGTCCGGTTTACGTAAAAAGATCAATCAGGGCGCGGTGAGTTTCGTTGACCTGCATTTGAG  >  NC_000913/3064679‑3065124

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