Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,102	T→C	G123G (GGT→GGC)	thrC →	threonine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,102	0	T	C	95.4%	173.1 / 0.1	67	G123G (GGT→GGC)	thrC	threonine synthase
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base C (27/35);  total (31/36)
Fisher's exact test for biased strand distribution p-value = 5.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.94e-01

TGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCACAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCT‑GACCCATATTGCGGGTGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGC‑AGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGG‑TTATCCTCTATCCACGAGG‑CAAAATCAGTCCACTGCAAGAAAAAC‑TGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCG‑ATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAACTGAAAGTGGC  >  NC_000913/3856‑4347                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                           tgaAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCT‑GACCCATATTGCGGGCGAt                                                                                                                                                                                                                                                         >  1:122206/1‑250 (MQ=255) tgaAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGAAGTGCGACGATCCTCTCGGCGTTTATTGGTGCTGAAATCCCGCAGGAAATCCTGGAAGCGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATCGTTCCACGGGCCAACGCTGGCATTTACAGATTTCGGCGGTCGCTTTTTTGCACAAATGCT‑GACACATCATGCGGGCGAt                                                                                                                                                                                                                                                         >  1:252535/1‑250 (MQ=255)  gaAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCT‑GACCCATATTGCGGGCGATa                                                                                                                                                                                                                                                        <  1:149845/250‑1 (MQ=255)   aaaTTGATGAGATGCGTAAGCTGGATGTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTGATTGGTGATGAAACCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCT‑GACCCATATTGCGGGCGATAAg                                                                                                                                                                                                                                                      <  2:300655/251‑1 (MQ=255)    aaTTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGACGCTTTATGGCACAAATGTT‑GTCCCATATTTCGGGCGAtt                                                                                                                                                                                                                                                        >  2:561051/1‑247 (MQ=255)            agaTGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGGTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATGTAAAGATTTCGGCGGTCGGTTTTTGGGCCAACAGCT‑GACACAAATTGTGGGCGATAAGCCAGGGAcc                                                                                                                                                                                                                                             >  2:313348/1‑251 (MQ=255)              aTGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCT‑GACCCATATTGCGGGCGATAAGCCAGTGACCAt                                                                                                                                                                                                                                           <  1:531250/251‑1 (MQ=255)                gCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCT‑GACCCATATTGCGGGCGATAAGCCAGTGACCAtt                                                                                                                                                                                                                                          <  1:431656/250‑1 (MQ=255)                 cTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGCATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTTTGGCCAAAATACT‑GACCCCTATTGCGGGGGATAAACCAGTGaacctatt                                                                                                                                                                                                                                        >  2:344184/1‑244 (MQ=255)                     aGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACCAATGCT‑GACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGAcc                                                                                                                                                                                                                                    >  1:456319/1‑251 (MQ=255)                      gCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCT‑GACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACcg                                                                                                                                                                                                                                   <  1:232735/251‑1 (MQ=255)                       cTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGACTTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCT‑GACCCCTATTGCGGGCGATAAGCCCGTGACCATTATGACcgc                                                                                                                                                                                                                                  >  1:78186/1‑251 (MQ=255)                           aTTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGACGCTTTATGGCACAAAAGCTTGACC‑ATATTGCGGGCGGTAAGcccgagcacattctgcacccgccc                                                                                                                                                                                                                              >  2:98832/1‑229 (MQ=255)                                  cACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCT‑GACCCATATTGCGGGCGATAAGCCAGTGACCATTCTTACCGCGCCCTCCGGTg                                                                                                                                                                                                                       >  1:588811/1‑251 (MQ=255)                                         aGTGCGAAGATCCTCTCGGCGTTTATTGGTGAAGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCT‑GACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCgg                                                                                                                                                                                                                <  2:343832/251‑1 (MQ=255)                                          gTGCGAAGATCCTTTCGGAGTTTATGGGTGATGAAACCCCGCAGGAAATCCTGGAAGGGCGCGTGCGCGCGGCGGTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGAGGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCT‑GACCCAGATTGCGGGCGATAAGCCAGTGTCCATTCTGACCGCGACCTCCGGTGATACCGGa                                                                                                                                                                                                               <  2:602442/251‑1 (MQ=255)                                               aaGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGCAAGCGATGTCGGATGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCT‑GACCCATATTGCGGGCGATAAGCCAGGGACCATTCTGACCGCCACCTCCGGCGATAACGGAGCGGc                                                                                                                                                                                                          >  1:342248/1‑251 (MQ=255)                                                     ctctCGGCGTTTATTGGTGATGAAATCCCGCAAGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTAAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTAGCTTTATGGCGCAAATCCT‑GACCCATATTGCGGGCGATAAACCAGTGGCCATTATGACCGCGAGCTCCGGTGaatatcgggagagcacccg                                                                                                                                                                                                    >  1:350777/1‑233 (MQ=255)                                                         cGGCGTTTATTTGTGATGAAATCCCGCAGGAAATCCTGGAAGGGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATTGCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCT‑GACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGC‑AGTGGCTCAt                                                                                                                                                                                               <  1:339337/251‑1 (MQ=255)                                                          ggCGTTTATTGGTGATGAAATCCCACAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGACCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCT‑GACCCATATTGCGGGTGATAAGCCAGTAACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGC‑AGTGGCTCATg                                                                                                                                                                                              <  1:144349/251‑1 (MQ=255)                                                           gCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCCATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCT‑GACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGa                                                                                                                                                                                                               >  1:437188/1‑234 (MQ=255)                                                           gCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGGGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTGCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCT‑GACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCTACCTCCGGTGATACCGGa                                                                                                                                                                                                               <  2:437188/234‑1 (MQ=255)                                                                 aTTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGGGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCT‑GACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGC‑AGTGGCTCATGCTTTCTa                                                                                                                                                                                       <  2:508450/251‑1 (MQ=255)                                                                   tGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCT‑GACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGC‑AGTGGCTCATGCTTTCTAc                                                                                                                                                                                      <  1:116731/250‑1 (MQ=255)                                                                     gTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCT‑GACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGC‑AGTGGCTCATGCTTTCTACGGt                                                                                                                                                                                   <  1:647825/251‑1 (MQ=255)                                                                        atgaAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCT‑GACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGC‑AGTGGCTCATGCTTTCTACGGTTTa                                                                                                                                                                                <  1:561051/251‑1 (MQ=255)                                                                         ggaAATCCCGAAGGAAACCCTGGAAGAGCGCGGGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATGGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCT‑GACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGC‑AGTGGCTCATGCTTTCTACGGTTTa                                                                                                                                                                                <  2:177311/249‑1 (MQ=255)                                                                          gaAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCT‑GACCCATATTGCGGGCGATAAGCCAGt                                                                                                                                                                                                                                                 <  1:155235/185‑1 (MQ=255)                                                                                  gcaGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCT‑GACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGC‑AGTGGCTCATGCTTTCTACGGTTTACCGAATGTGa                                                                                                                                                                      <  1:178545/250‑1 (MQ=255)                                                                                      gAAAGCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCT‑GACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGC‑AGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGt                                                                                                                                                                  <  1:452767/251‑1 (MQ=255)                                                                                        aaTCCTAGAAGGGCGCGAGCGCGCGGCGTTTGCATTCCCGGCTCCGGTCGCCAATGTTGCAAGCGATGTCTGTTGTCTGGATTTGTTACACGGGCCAACGCTGGCATTTAAAGATTTCGGAGGTCGCATTATGGACCAAATGAT‑AACCCATATTTCGGGAGAAAAACCAGTGAACATTCTGACCGGGACCTCAGGTGAGACCGGGGCGGC‑AGTTGCTCAGGCtttgtaaaggttaccccatgtgtaaga                                                                                                                                                                  >  2:610749/1‑225 (MQ=255)                                                                                                  gagCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCT‑GACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGC‑AGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGG‑TTATCCTCTa                                                                                                                                                     <  1:194429/251‑1 (MQ=255)                                                                                                   agCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCT‑GACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGC‑AGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGG‑TTATCCTCTAt                                                                                                                                                    <  1:488484/251‑1 (MQ=255)                                                                                                      gcgTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCT‑GACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGC‑AGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGG‑TTATCCTCTATcc                                                                                                                                                  >  1:384063/1‑250 (MQ=255)                                                                                                                   gttgtctttCCCGGCTCCGGTCGCCAAGGTGGAAAGCGATGTCGGTTGTCTGGAAGTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCT‑GACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGC‑AGTGGCTCATGCTTTCGCCGGTTTGCCGAATGTGAAAGTGG‑TTATCCTCTATCCACGAGG‑CAAAATCa                                                                                                                                   <  2:298858/244‑1 (MQ=255)                                                                                                                   gTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCT‑GACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGC‑AGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGG‑TTATCCTCTATCCACGAGG‑CAAAATCa                                                                                                                                   <  1:286167/251‑1 (MQ=255)                                                                                                                         cTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCT‑GACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGC‑AGTGGCTCATGCTTTCTACGGTTTAACGAATGTGAAAGTGG‑TTATCCTCTATCCACGAGG‑CAAAATCAGTCCAc                                                                                                                             >  1:314301/1‑251 (MQ=255)                                                                                                                                         gCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCT‑GACCCATATTGCGGGTGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGC‑AGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGG‑TTATCCTCTATCCACGAGG‑CAAAATCAGTCCACTGCAAGAAAAAC‑TGtt                                                                                                            >  1:341423/1‑251 (MQ=255)                                                                                                                                           cAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCAACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCT‑GACCCATATTGAGGGCGATAACCCAGTGACCATTCTGACCGAGCACTCCGGTGCTAACGGAGCGGC‑AGGGGCTCATGCTTTTTACGGTTTACCGCATGTGAAAGTGG‑GTATCCTCTATCCCCGAGT‑CAAAATCAGTCCACTGCACGGAAAAC‑TATTCt                                                                                                          >  1:562232/1‑251 (MQ=255)                                                                                                                                                  gaaatctaTGTCGGTTGTCTGGAATTGTTCCACGGGCCAACTCTGTCATTTAATGATTTCGGCGTTCGCTTTATGTCACAAATGCT‑GCCCCAGATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGCCCTCCGGTGATTCCGGAGCGTC‑AGTGTCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTTG‑TTATTCTCTATCCACGAGG‑CAAAATCAGTCCACTGCAAGAAAAAT‑TGGTCTGTACAt                                                                                                    <  2:223436/243‑1 (MQ=255)                                                                                                                                                  gaaaggtatgacGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCTGCGGTCGCTTTATGGCACAAATGCT‑GACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGTAGCGGC‑AGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGG‑TTATCCTCTATCCACGAGG‑CAAAATCAGTCCACTGCAAGAAAAAC‑TGTTCTGTACAtt                                                                                                   <  2:429231/240‑1 (MQ=255)                                                                                                                                                                  gTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCT‑GACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGC‑AGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGG‑TTATCCTCTATCCACGAGG‑CAAAATCAGTCCACTGCAAGAAAAAC‑TGTTCTGTACATTGGGCGGCAATATCGaa                                                                                   <  1:344184/251‑1 (MQ=255)                                                                                                                                                                       gAATTGTTCCACGAGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTAAGGCCCAAATACT‑GACCCATATTGCGGGCGATAAGCCAGTGAACATTCTGACACCGACCTCCGGAGAAACCGGAGCGGA‑CGTGGCTCATGCTTTCTTCGGTTTACCGAATGTGAAAGTGG‑CTATCCTCTATCCACGAGG‑CAAAATCAGTCCACTGCAAGAAAAAC‑TGTTCTGTAAATTGTGCGGCAATATCGACACTGt                                                                              >  1:202049/1‑251 (MQ=255)                                                                                                                                                                         aTTGTTCCAAGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGGTTTATGGCACAAATGAT‑GACCCATATTGAAGGCGATAAGCCAGTGACCATTCTGAGCGCGACATCCGGTGATACCGGAGCGGC‑AGGGGGTCATGGTTTTTTAGGTTTTACGAATGTGAAAGGGGGTTACCCT‑TAACCCCGAGGGCAAATACA‑TCCACAGCAAGAAAAAT‑AGTTCTGTACATTGGGCGGGAATATCGaaaatgatt                                                                            >  2:21390/1‑245 (MQ=255)                                                                                                                                                                            gTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCT‑GACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGC‑AGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGG‑TTATCCTCTATCCACGAGG‑CAAAATCAGTCCACTGCAAGAAAAAC‑TGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGcc                                                                          >  1:137228/1‑250 (MQ=255)                                                                                                                                                                             ttCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGGCGGTTTATGGAACACAAGCT‑GAACCCTATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATTCCGGAGCGGCCAG‑GGCACATGCCTTCTACGGGTTACCGAAGGGTCGAGTGG‑TTTTTCTTTATCCAAGTGG‑CAAATTCAGTCCAATGCAAGAAAAAC‑TTTTATGTACATTGGGCGGAAAAATCggacatggtgccat                                                                        >  2:637925/1‑238 (MQ=255)                                                                                                                                                                               ccccGGGCCAACGCTGGAATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCT‑GACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGC‑AGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGG‑TTATCCTCTATCCACGAGG‑CAAAATCAGTCCACTGCAAGAAAAAC‑TGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCg                                                                      <  2:36548/248‑1 (MQ=255)                                                                                                                                                                               ccACGGGCCAACGCTGGCATTTAAAGATTGCGGCGGTCGCGTGATGGCACAAATGCT‑GACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGC‑AGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGG‑TTATCCTCTATCCACGAGG‑CAAAATCAGTCCACTGCAAGAAAAAC‑TGTTCTGTACAGTGGGCGGCAATAGCGAAACTGTTGCCATCg                                                                      <  2:78186/251‑1 (MQ=255)                                                                                                                                                                                  cGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCT‑GACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGC‑AGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGG‑TTATCCTCTATCCACGAGG‑CAAAATCAGTCCACTGCAAGAAAAAC‑TGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACg                                                                   >  1:638610/1‑251 (MQ=255)                                                                                                                                                                                    ggCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCT‑GACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGC‑AGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGG‑TTATCCTCTATCCACGAGG‑CAAAATCAGTCCACTGCAAGAAAAAC‑TGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGc                                                                 <  1:21390/251‑1 (MQ=255)                                                                                                                                                                                     gCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCT‑GACCCATATTGCGGGCGATAAGCCAGAGACCATTCTGACCGCGACCTCCGGTGAAACCGGAGCGGC‑AGCGGCTCATGCTTTCTACGGCTTACCGAATGTGAAAGTGG‑TTATCCTATATCCACGAGG‑CAAAATAAGTCCACTGCAAGAAAAAC‑TGGTCTGTACATTGGGCGGCAAAAACGAAACTGTTGCCATCGACGGCg                                                                >  1:657483/1‑251 (MQ=255)                                                                                                                                                                                      ccAACGCTGGCATTTAAAGATTTCGGCGGGCGCTTTATGGCACAAATGCT‑GACCCATATTGCGGGCGATAAACCAGGGACCATTCTGACCGCGACCCCCGGTGATACCCGGGCGGC‑AGTGGCCCAAGCTTTCTACGGTTTACCGACAGTTAAAATGG‑TAATTCCCTATCCCCCAGG‑CAAAATCAATCCACTACAAGAAAAACATGT‑CTGTACATATGGCGGCCATAACGAAACTGTTGCCATCGACggga                                                                >  2:479856/1‑248 (MQ=255)                                                                                                                                                                                              ggCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCT‑GACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGC‑AGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGG‑TTATCCTCTATCCACGAGG‑CAAAATCAGTCCACTGCAAGAAAAAC‑TGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCG‑ATg                                                      <  2:638610/251‑1 (MQ=255)                                                                                                                                                                                                             cggcggTCGCTTTATGGCACAAATGCT‑GACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGCGCGGC‑AGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGG‑TTATCCTCTATCCACGAGG‑CAAAATCAGTCCACTGCAAGAAAAAC‑TGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCG‑ATGCCTGTCAGGCGCTgg                                       <  1:238623/251‑1 (MQ=255)                                                                                                                                                                                                               gcggTCGCTTTATGGCACAAATGCT‑GACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGC‑AGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGG‑TTATCCTCTATCCACGAGG‑CAAAATCAGTCCACTGCAAGAAAAAC‑TGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCG‑ATGCCTGTCAGGCGCTGGTg                                     <  1:114723/251‑1 (MQ=255)                                                                                                                                                                                                                cggTCGCTTTATGGCACAAATGCT‑GACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGC‑AGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGG‑TTATCCTCTATCCACGAGG‑CAAAATCAGTCCACTGCAAGAAAAAC‑TGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCG‑ATGCCTGTCAGGCGCTGGTGa                                    <  1:264132/251‑1 (MQ=255)                                                                                                                                                                                                                  gTCGCTTTATGGCACAAATGCT‑GACCAATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGC‑AGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGG‑TTATCCTCTATCCACGAGG‑CAAAATCAGTCCACTGCAAGAAAAAC‑TGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCG‑ATGCCTGTCAGGCGCTGGTGa                                    >  1:426292/1‑249 (MQ=255)                                                                                                                                                                                                                     gCTTTATGGCACAAATGCT‑GACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGC‑AGTGGCTAATGCTTTCTACGGTTTACCGAATGTGAAAGTGG‑TTATCCTCTATCCACGAGG‑CCAAATCAGTCCACTGCAAGAAAAAC‑TGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGACATCGACGGCGATTTCGTAT‑CCTGTCAGGCGCTGGTGAAGccg                               >  2:104744/1‑249 (MQ=255)                                                                                                                                                                                                                       tttATGGCACAAATGCT‑GACCCATATTGCGGGCGATAAGCAAGTGACCATTCTGACCGCGACCTCCGGAGATACCGGAGCGGC‑AGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGG‑TTATCCTCTATCCACGAGG‑CAAAATCAGTCCACTGCAAGAAAAAC‑TGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCG‑ATGCCTGTCCGGCGCTGGTGAAGCAGGc                             >  1:589504/1‑251 (MQ=255)                                                                                                                                                                                                                                caAATGCT‑GACCCATATTGCGGGTGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGC‑AGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGG‑TTATCCTCTATCCACGAGG‑CAAAATCAGTCCACTGCAAGAAAAAC‑TGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCG‑ATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTtgatg                     >  1:51029/1‑250 (MQ=255)                                                                                                                                                                                                                                  cgcttgtccgatctatTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCACCGGTGATACCGGAGCGGC‑AGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTTG‑TTATCCTCTATCCACGAGG‑CAAAATCAGTCCACTGCAAGAAAAAC‑TGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCG‑ATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTtgatgat                   <  2:452279/238‑1 (MQ=255)                                                                                                                                                                                                                                   aTGCT‑GACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGC‑AGTGGCTCATGCTTTCTACGGTTTACCGACTGTGAAAGTGG‑TTATCCTCTATCCACGAGG‑CAAAATCAGTCCACTGCAAGAAAAAC‑TGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGCTTTCG‑ATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTtgatgatg                  >  1:655131/1‑250 (MQ=255)                                                                                                                                                                                                                                         gACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGC‑AGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGG‑TTATCCTCTATCCACGAGG‑CAAAATCAGTCCACTGCAAGAAAAAC‑TGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCG‑ATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATgaagag             <  1:33812/250‑2 (MQ=255)                                                                                                                                                                                                                                          ccccATCTTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGC‑ATTGTCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGG‑TTATCCTCTATCCACGAGG‑CAAAATCAGTCCACTGCAAGAAAAAC‑TGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCG‑ATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATgaagagct           <  2:229175/250‑4 (MQ=255)                                                                                                                                                                                                                                                atTGCGGGCGATAAGCCAGTGACCATTCTGCCCGCGACCTCCGGTGATACCGGAGCGGC‑AGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGG‑TTATCCTCTATCCACGAGG‑CAAAATCAGTCCACTGCAAGAAAAAC‑TGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCG‑ATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAgg     >  1:368800/1‑250 (MQ=255)                                                                                                                                                                                                                                                atTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGC‑AGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGG‑TTATCCTCTATCCACGAGG‑CAAAATCAGTCCACTGCAAGAAAAAC‑TGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCG‑ATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTtgatgat                   >  1:452279/1‑237 (MQ=255)                                                                                                                                                                                                                                                   gCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGC‑AGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGG‑TTATCCTCTATCCACGAGG‑CAAAATCAGTCCACTGCAAGAAAAAC‑TGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCG‑ATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAGCTGAAAGTGgc  <  1:98832/251‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                           TGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCACAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCT‑GACCCATATTGCGGGTGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGC‑AGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGG‑TTATCCTCTATCCACGAGG‑CAAAATCAGTCCACTGCAAGAAAAAC‑TGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCG‑ATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGAACTGAAAGTGGC  >  NC_000913/3856‑4347

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