Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,804	T→G	A357A (GCT→GCG)	thrC →	threonine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,804	0	T	G	92.8%	188.9 / 3.6	69	A357A (GCT→GCG)	thrC	threonine synthase
Reads supporting (aligned to +/- strand):  ref base T (3/2);  new base G (24/40);  total (27/42)
Fisher's exact test for biased strand distribution p-value = 3.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.90e-01

AACGATACCGTGCCACGTTTCCTGCACGACGGTCAGTGGTCACCCAAAGCGACTCAGGCGACGTTATCCAACGCGATGGACGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACTTCGGAGCCGCACGCTGCCGTAGCTTATCGTGCGCTGCGTGATCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGCAGAACGTGCTG‑ATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTT  >  NC_000913/4559‑5052                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                          aaCGATACCGTGCCACTTTTCCTGCACGACGGTCAGTGGTCACCCAAAGCGACTCAGGCGACTTTATCCAACGCGATGGACGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGTAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCg                                                                                                                                                                                                                                                      <  1:110312/251‑1 (MQ=255) aaCGATACCGTGCCACGTTTCCTGCACGACGGTCAGTGGTCACCCAAAGCGACTCAGGCGACGTTATCCAACGCGATGGACGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCg                                                                                                                                                                                                                                                      <  1:424058/251‑1 (MQ=255)  aCGATACCGTGCCACGTTTCCTGCACGACGGTCAGTGGTCACCCAAAGCGACTCAGGCGACGTTATCCAACGCGATGGACGTGATTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGt                                                                                                                                                                                                                                                     <  1:314330/251‑1 (MQ=255)       aCCGTGCCACGTTTCCTGCACGACGGTCAGTGGTCACCCAAAGCGACTCAGGCGACGTTATCCAACGCGATGGACGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAACCGCACGCTGCCGTAGCGTATCGTGCGCt                                                                                                                                                                                                                                                >  2:256979/1‑251 (MQ=255)           tGCCACGTTTCCTGCACGACGGTCAGTGGTCACCCAAAGCGACTCAGGCGACGTTATCCAACGCGATGGACGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACTTCGGAGCCGCACGCTGCCGTAGCTTATCGTGCGCTGCGt                                                                                                                                                                                                                                            >  2:73750/1‑251 (MQ=255)              cACGTTTCCTGCACGACGGTCAGTGGTCACCCAAAGCGACTCAGGCGACGTTATCCAACGCGATGGACGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCAGCCGTAGCGTATCGTGCGCTGCGTgcc                                                                                                                                                                                                                                         >  2:155834/1‑249 (MQ=255)                     ccTGCACGACGGTCAGTGGTCACCCAAAGCGACTCAGGCGACGTTATCCAACGCGATGGACGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTg                                                                                                                                                                                                                                   >  1:273863/1‑250 (MQ=255)                                gTCAGTGGTCACCCAAAGCGACTCAGGCGACGTTCTCCAACGCGATGGACGGGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGaa                                                                                                                                                                                                                       <  1:398352/251‑1 (MQ=255)                                gTCAGTGGTCACCCAAAGCGACTCAGGCGACGTTATCCAACGCGATGGACGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGTTGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGa                                                                                                                                                                                                                        >  2:116093/1‑250 (MQ=255)                                     tGGTCACCCAAAGCGACTCAGGCGACGTTATCCAACGCGATGGACGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGCAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATgg                                                                                                                                                                                                                  >  2:418795/1‑251 (MQ=255)                                         cACCCAAAGCGACTCAGGCGACGTTATCCAACGCGATGGACGTGAGTCAGCCGAACACCTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCtt                                                                                                                                                                                                               <  1:34561/250‑1 (MQ=255)                                          ccccAAAGCGACTCAGGCGTCGTTATCCAACGCGATGGCCGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCCCGCTGCCGTAGCGTATCGTGCGTTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGt                                                                                                                                                                                                             <  1:297593/250‑1 (MQ=255)                                            ccAAAGCGACTCAGGCGACGTTATCCACCGCGATGGACGTGAGCCGGCCGAACAACTGTCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTc                                                                                                                                                                                                           <  1:414253/251‑1 (MQ=255)                                               aaGCGACTCAGGCGACGTTATCCAACGCGATGGACGTGAGTCAGCCGAACAACTGGCCGGGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTc                                                                                                                                                                                                        <  1:463741/251‑1 (MQ=255)                                                ggtgtctagggCGACGTTACACAACGTGAGGGACGGAAGTCAGCCGAACAACTGGCCGCGGGGGGAAGAGGTGTACCGCCGCAAAAACTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTTCACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGAGAATATGGCTTGTTCCTCg                                                                                                                                                                                                       <  2:284611/242‑1 (MQ=255)                                                  cGACTCAGGCGACGTTATCCAACGCGATGGACGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGGGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGc                                                                                                                                                                                                     >  2:368653/1‑251 (MQ=255)                                                    ccTCAGGCGACGTTATCCAACGCGATGGACGTGAGTCAGCCGAACAACTGGCCGCGGGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCAc                                                                                                                                                                                                   <  1:116093/250‑1 (MQ=255)                                                      tCAGGCGACGTTATCCAACGCGATGGACGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACcg                                                                                                                                                                                                 >  1:272649/1‑251 (MQ=255)                                                           cGACGTTATCCAACGCGATGGACGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCAt                                                                                                                                                                                            >  1:293130/1‑251 (MQ=255)                                                            gACGTTACCCAAGGCGTGGGAGGTGAGTCAGCCGAACAACTGGCCGGGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATc                                                                                                                                                                                           <  2:272649/251‑1 (MQ=255)                                                                   tCCAACGCGATGGACGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGa                                                                                                                                                                                     >  1:312740/1‑250 (MQ=255)                                                                          cgATGGACGTGCGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTTTTCCCCGGCCCCGCCCATCCGGCGAAATTTaa                                                                                                                                                                              >  2:57865/1‑250 (MQ=255)                                                                           gATGGACGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGt                                                                                                                                                                                                             >  1:300640/1‑218 (MQ=255)                                                                              ggACGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCAGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGACAttttagg                                                                                                                                                                            >  1:479725/1‑244 (MQ=255)                                                                               gACGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAagagag                                                                                                                                                                        <  1:478783/251‑1 (MQ=255)                                                                                 cGTGAGTCAGCCGAAAAGCTGGCCGCGTGGGGAAGAGTGATGCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCg                                                                                                                                                                      <  2:127693/251‑1 (MQ=255)                                                                                        cAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAg                                                                                                                                                                <  1:172918/250‑1 (MQ=255)                                                                                             gAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATtc                                                                                                                                                          <  1:394235/251‑1 (MQ=255)                                                                                               acaacTGGCCGCGTGTGGAAAAGTGGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATtctc                                                                                                                                                        <  2:293130/251‑1 (MQ=255)                                                                                                 agcgggCCGGGTGTGGAATGGTTGTTCCGCCGCAAAAGCTGGCAACTGAAAGGGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCGCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCgg                                                                                                                                                      <  2:241495/247‑1 (MQ=255)                                                                                                 agagggtggggtgtGGAAGAGTTGTCCGGCCGCAAAATCGGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCgg                                                                                                                                                      <  2:229784/241‑1 (MQ=255)                                                                                                 aacTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCgg                                                                                                                                                      <  1:416985/251‑1 (MQ=255)                                                                                                    gggCCGCTTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGa                                                                                                                                                   <  1:234080/250‑1 (MQ=255)                                                                                                     ggCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGa                                                                                                                                                   <  1:80541/250‑1 (MQ=255)                                                                                                          ggtgtGGACGAGGTGTCCCGCCGCAAAATCTGGCAACTGAAAGAGCTTGGTTACGCAGCCGTGGATGATGACACCACGCACCAGACAATGCGTGAGTTAACAGCCCTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTCAACCCAGGCGAATATGGCTTGTTCCTCGGCACCGCTCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGtt                                                                                                                                             <  2:107271/250‑1 (MQ=255)                                                                                                             gtgGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTgg                                                                                                                                           <  2:216008/250‑1 (MQ=255)                                                                                                              gggaggagTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGAt                                                                                                                                         <  2:72559/246‑1 (MQ=255)                                                                                                               ggAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACTTCGGAGCCGCACGCTGCCGTAGCTTATCGTGCGCTGCGTGATCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGAt                                                                                                                                         <  1:73750/250‑1 (MQ=255)                                                                                                                gAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTTCTCGGCACCGCGCCTCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTggtact                                                                                                                                       >  2:95219/1‑247 (MQ=255)                                                                                                                 aaGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCCTCCGGCGAAATTTAAAGAAAGCGTGGAAGCGCTTCTCGGTGAAACGTTGGATCTg                                                                                                                                      >  1:156937/1‑251 (MQ=255)                                                                                                                 aaGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTg                                                                                                                                      <  1:451975/251‑1 (MQ=255)                                                                                                                                   aaaTCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGGGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGCAGaa                                                                                                                    <  2:273863/251‑1 (MQ=255)                                                                                                                                         ggCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGCAGAACGTGc                                                                                                               <  1:474321/250‑1 (MQ=255)                                                                                                                                            aaCTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGCAGAACGTGCTG‑At                                                                                                          >  1:58289/1‑251 (MQ=255)                                                                                                                                                   agagCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGCAGAACGTGCTG‑ATTTACCCt                                                                                                   <  2:312740/251‑1 (MQ=255)                                                                                                                                                    gagCTGGGTTAGGCAAGCGTGGATGATGAAACCACGCAAAAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGCAGAACGTGCTG‑ATTTACCCtt                                                                                                  <  2:478511/251‑1 (MQ=255)                                                                                                                                                           gctATGCAGCCGGGGATGATGAAACCACGCAACAGACAATGCGTGCGTTAAGAGTACGGGGCTACACCTCGGGCCCGCACGCTGCCGTAGCGTATCGGGGGCTGCGTGACCAATTGACTCCAGCCGAGTATAGCTTTTTTCAGGGCACCGCCCAAAGCGCGACATTTTAAGAGAGCGTGGAGCCCGTTCTCGGGGCAAAGTTTGCACGTCCAAAAGAGATGGCAGAACGAGTCG‑ATTTACCCTTGTTTTc                                                                                            >  2:35770/3‑250 (MQ=255)                                                                                                                                                           gTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGCAGAACGTGCTG‑ATTTACCCTTGCTTTca                                                                                           <  2:442370/251‑1 (MQ=255)                                                                                                                                                             tATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGCAGAACGTGCTG‑ATTTACCCTTGCTTTcac                                                                                          <  1:95219/250‑1 (MQ=255)                                                                                                                                                                    ccGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGCAGAACGTGCTG‑ATTTACCCTTGCTTTCACATAATCTg                                                                                  <  1:256979/251‑1 (MQ=255)                                                                                                                                                                      gTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACTTCGGAGCCGCACGCTGCCGTAGCTTATCGTGCCCTGCGTGATCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAACTTGAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGCAGAACGTGCTG‑ATTTACCCTTGCTTTCACATAATCTGcc                                                                                >  1:156666/1‑251 (MQ=255)                                                                                                                                                                        ggATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGTCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGCAGAACGTGCTG‑ATTTACCCTTGCTTTCACATAATCTGCccg                                                                              <  2:25999/251‑1 (MQ=255)                                                                                                                                                                             atgaAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGCAGAACGTGCTG‑ATTTACCCTTGCTTTCACATAATCTGCCCGCCGAt                                                                         >  1:362427/1‑251 (MQ=255)                                                                                                                                                                                 aaCCACGCAACAGACAATGCGTGAGTTCAAAGAACTGGGCTACACCTCGGAACCGCACGCTGCCGTAGCGTATCGAGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCCCGGCACCGCGCCTCCGGCGAAATTTAAAGAGGGCGTGGAAGCGATTCTCGGAGAAACGGTGGATATGCCAAAAGAGGTGGCAGAACGAGATG‑ATTTTCCCTTGTTTTCACCTATACTGCCCCCCGAttttt                                                                     >  2:294296/1‑250 (MQ=255)                                                                                                                                                                                      cgcaacagacaagggggtatttAAAAGAACTGGTCTACCCCCCGGAGCCGCACGCGGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATAGGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACTTTGGATCTGCAAAAAGAGCTGGCAGATCGTGCTG‑ATTTCCCCTTGCTTTAACATAATCTGCCCGCCGATTTTGCtgcg                                                                <  2:362427/231‑1 (MQ=255)                                                                                                                                                                                        cAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATAGGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGCAGAACGGGCTG‑ATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCtgcgtt                                                              <  1:418795/251‑1 (MQ=255)                                                                                                                                                                                          aCAGACAATGAGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGAGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTTAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGCAGAACGTGCTG‑ATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCtgcgttgc                                                            >  1:83719/1‑251 (MQ=255)                                                                                                                                                                                                 aggCGTGAGTTAAAAGAACGGGGCTACACCCCTGAGCCGCACGCGCCCGGAGCGTTTCGTGCGCTGCGTGACCATTTGAATCCGGGCGAATATGGTTTGTTCCCCGCCACCGCGCATCCTGCGAAATTTAAAGAGAGCGGGGAAGCAATTCTCGGTGAAACGTTTGATCTGCCAAAAGAGCTGGCAGAACGTG‑TGTATTTCCCCTTGCTTTCACATATCCTGCCCGCCGATTTTGCTGCGTTGCGTAAAtt                                                     <  1:60698/249‑1 (MQ=255)                                                                                                                                                                                                 aTGCGTGAGTTACAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCg                                                                                                                                                                                                       >  2:478645/1‑106 (MQ=255)                                                                                                                                                                                                 aTGCGTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGGAGCGGATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCg                                                                                                                                                                                                       <  1:478645/106‑1 (MQ=255)                                                                                                                                                                                                     gTGAGTTAAAAGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGCAGAACGTGCTG‑ATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATtgatg                                                 >  1:320584/1‑251 (MQ=255)                                                                                                                                                                                                             aaaGAACTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGGGTATCGTGCTCTGCGTGACCAGTTGAATCCAGGAGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCCCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGCAGAACGTGCTG‑ATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAatca                                         <  1:23503/251‑1 (MQ=255)                                                                                                                                                                                                                   cTGGGCTACACCTCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTACTCGGCACCGCGCATCCGGCGAAAATTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAAAGTTGGATCTGCCAAAAAAGCTGGCAGAACGGGCTG‑ATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATAATCAGt                                    >  2:49592/1‑250 (MQ=255)                                                                                                                                                                                                                           cacctCGGAGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGGTCTGCCAAAAGAGCTGGCAGAACGGGCAG‑ATTTACCCTTGCTCTCACATAATCTGCCCGCCGATTTTGTTGCGGTGCGTAAATTTATAGTGAATAATCAGTAAAATCTAt                           >  2:407594/5‑251 (MQ=255)                                                                                                                                                                                                                                 gcaGCCGCACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCCGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCCCCGCGCATCCGGCGAAATTTAACGAGACCGTGGAAGCGATTCTCGGCGCAACGTTGGATCTGCCAAAAGAGCTGGCAGAACCCGCTG‑CTTTACCCCTGCCTTCCCATCATCCGCCCGCCCCCTTTGCTGCGTTGCCTACATTGCTGCTGaaccacccaccaaacctactccctc                     >  2:308322/3‑228 (MQ=255)                                                                                                                                                                                                                                  gAGCCGCACGCTGCCGTAGCGTAGCGTGCGCTGCGTGACCAGTTGAACCCAGGCGAATATGGCTTATTCCTCGGCCCCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGCAGAACGTGCTG‑ATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTAt                    <  1:244796/251‑1 (MQ=255)                                                                                                                                                                                                                                        cACGCTGCCGTAGCGTATCGTGCGCTGCGTGACCAGTTGAATCCAGGGGAATATGGCTTGTTCCTCGGCCCCGCGCAGCCGGCGAAATTTAAAGAGAGCGGGGAAGCGTTTCTCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGCAGAACGTGCTG‑ATTGACCCTTGCTTGCACAAAATCGGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAAt              <  1:57865/251‑1 (MQ=255)                                                                                                                                                                                                                                                cGTAGCTTATCGTGCGCTGCGTGATCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGACACGTTGGATCTTCCAAAAGAGCTGGCAGAACGTGCTG‑ATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCTTCAGTAAAATCTATTCATTATCTCAATCAGGCCgg      >  1:245834/1‑251 (MQ=255)                                                                                                                                                                                                                                                    gCTTATCGTGCGCTGCGTGACCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAACCGTTGGATCTGCCAAAAGAGCTGGCAGAACGTGCTG‑ATTTCCCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGttt  <  1:147435/251‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                          AACGATACCGTGCCACGTTTCCTGCACGACGGTCAGTGGTCACCCAAAGCGACTCAGGCGACGTTATCCAACGCGATGGACGTGAGTCAGCCGAACAACTGGCCGCGTGTGGAAGAGTTGTTCCGCCGCAAAATCTGGCAACTGAAAGAGCTGGGTTATGCAGCCGTGGATGATGAAACCACGCAACAGACAATGCGTGAGTTAAAAGAACTGGGCTACACTTCGGAGCCGCACGCTGCCGTAGCTTATCGTGCGCTGCGTGATCAGTTGAATCCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGCAGAACGTGCTG‑ATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTT  >  NC_000913/4559‑5052

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