Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	157,585	C→T	D50N (GAT→AAT)	folK ←	2‑amino‑4‑hydroxy‑6‑ hydroxymethyldihydropteridine diphosphokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	157,585	0	C	T	94.9%	172.3 / ‑3.0	59	D50N (GAT→AAT)	folK	2‑amino‑4‑hydroxy‑6‑ hydroxymethyldihydropteridine diphosphokinase
Reads supporting (aligned to +/- strand):  ref base C (2/1);  new base T (33/23);  total (35/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

CAAACAGCGGCCACAGCATAAATCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATCTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGG‑GCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACGCCACTGTCATGCGGTACCCTCACGACGTGGTGCGCGTT‑TGCGTGGACGACGAGTACGACGACGCGGTGACGGTTCTTCATCCAGCTCGTTGAGCATCCCTTTTTGGTCT  >  NC_000913/157342‑157831                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                         cAAACAGCGGCCACAGCATAAATCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGGGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGc                                                                                                                                                                                                                                                   >  1:317081/1‑251 (MQ=255)    aCAGCGGCCACAGCATAAATCCACGATTCTTCATATCGTAGTGCGGCACGTTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTGGAGATCCAGCGAGAGTGGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGCGTAGCTCTTCAGGTGCAAGATAGGATTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTTTTCTTGCGGaccc                                                                                                                                                                                                                                                >  2:249987/1‑251 (MQ=255)     cAGCGGCCACAGCATAAATCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGCGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGcccc                                                                                                                                                                                                                                                >  1:215739/1‑250 (MQ=255)          gCCACAGCATAAATCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGt                                                                                                                                                                                                                                          <  1:448516/251‑1 (MQ=255)                        ccACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAAATCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGAACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGAGCGGTaaa                                                                                                                                                                                                                             >  1:334030/1‑250 (MQ=255)                                                           gCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGAAGCGACCTTGCTGCCATTCACTACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAAAGAAGAAACGGTAAGAATGTGGCTTTCAGGGttat                                                                                                                                                                                         >  2:361507/1‑251 (MQ=255)                                                            cGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATc                                                                                                                                                                                        >  1:321722/1‑251 (MQ=255)                                                             gTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGTTCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGCGATATCg                                                                                                                                                                                       <  2:30245/251‑1 (MQ=255)                                                                cAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCt                                                                                                                                                                                    >  1:185777/1‑251 (MQ=255)                                                                 aGTATTTATCACTTCATTCCCAAACAGCATGATGTCGAGATCCAGCGTGAGAGGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGAGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTCAGTAATCGGGTTGATTTTGAGGCCCCCGCGGTGGGGTGCGGCACAACGAAGAAACGGGAAGAAGGTGGCTTTCTGGGATATCgact                                                                                                                                                                                    >  2:525888/1‑247 (MQ=255)                                                                  gTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGATCTTCAGGTGCAAGAGAGGTTTCCGGCGCCACGGCTGCGTTTAAGTAATCGGGTAGATTTTGCGGCCCCAGCGGAGGGGTGCGGTAAAACGAAGAACCGGAAAGAATGGGGCTTTCAGGGATATCGCCtt                                                                                                                                                                                   >  2:509599/1‑249 (MQ=255)                                                                  gTATTTATCACTTCATTAACAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATCTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAGGAATGCGGCTTTCAGGGGTGTCGCCTa                                                                                                                                                                                   >  2:95139/1‑250 (MQ=255)                                                                    attaatcacctaaTTACCAAACAGCATGATGTCGAGAACCAGCGGGCATGGTCCCCAGCGTTCAGGTGTGCGGACGCGACCTTGCTGCAATTCAATACTCTGTGTGTGATTGAGTAGCTCTTCAGGGGCAAGAGAGGGTGCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTTGCTTTCAGGGATAGCGCCTAATg                                                                                                                                                                                <  2:22365/239‑1 (MQ=255)                                                                     tttATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGc                                                                                                                                                                               >  1:293323/1‑251 (MQ=255)                                                                            cTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAgg                                                                                                                                                                        >  1:356566/1‑251 (MQ=255)                                                                              tCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCC‑AGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGG‑gc                                                                                                                                                                     >  1:430787/1‑250 (MQ=255)                                                                                aTTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCa                                                                                                                                                                                                                                               <  1:184789/176‑1 (MQ=255)                                                                                          aGCATGATGTCGATATCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGG‑GCAGCATTGACCTg                                                                                                                                                         >  1:9249/1‑251 (MQ=255)                                                                                            cATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGG‑GCAGCATTGACCTGCt                                                                                                                                                       >  1:495499/1‑251 (MQ=255)                                                                                            cATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGAGTGATTGAGTAGCTCTGCAGGTGCAAGAGAGGTTTACAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGAGTTTGCAGCCACAGCGGTGGGGTGCGGTAAATCGAAGAACCGGTAAGAATGGGGCTTTCGGGGATAACGCCAACTGCTTT‑AGGAGCAGCATTgcactggc                                                                                                                                                       >  2:163484/1‑244 (MQ=255)                                                                                                 tgTCGAGATCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGG‑GCAGCATTTAACTGCTCCAGc                                                                                                                                                  >  2:256282/1‑251 (MQ=255)                                                                                                     gagaTCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGG‑GCAGCATTGACCTGCTCCAGCGgag                                                                                                                                              <  2:408523/251‑1 (MQ=255)                                                                                                          ccAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGG‑GCAGCATTGACCTGCTCCAGCGGAGAGGcc                                                                                                                                         <  1:531033/251‑1 (MQ=255)                                                                                                                       tCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGCGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGTGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTACGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGG‑GCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCt                                                                                                                             >  2:457683/1‑250 (MQ=255)                                                                                                                                ttCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTGGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTc                                                                                                                                                                                                 <  1:372749/174‑1 (MQ=255)                                                                                                                                  cAGCGTTGCGGACGCTACCGTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCGCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCGGCGTTTAAGTAATCTGGTGGAGTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATAGCGCCTAATGCTTACATG‑GCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATAc                                                                                                                 <  2:215479/251‑1 (MQ=255)                                                                                                                                          cGGACGCGACCTTGCTGCAATTCACTACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGAGCGGTAAAACGACGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCATG‑GCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACGCCACTGt                                                                                                         >  1:88007/1‑251 (MQ=255)                                                                                                                                           ggACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTACGTAATCGGGGTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAAAGCTTTCAGG‑GCAGCATTGACCTGCTCCAGCGGAGAGGCCAGTTTGCTGCCTATGGCCCTAAACGCCCCTGTc                                                                                                        >  2:116818/1‑251 (MQ=255)                                                                                                                                            gACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGG‑GCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACGCCACTGTc                                                                                                        >  2:406274/1‑250 (MQ=255)                                                                                                                                                cgACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGG‑GCAGCATTGACCTGCTCCAGCGGCGAGGCCAGATTGCTGCCTATGGCAATATACGCCACTGTCATGCg                                                                                                   >  2:158023/1‑251 (MQ=255)                                                                                                                                                   ccTTGCTGCAATTCAATACTCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGG‑GCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACGCCACTGTCATGCGGTa                                                                                                <  2:321722/251‑1 (MQ=255)                                                                                                                                                      tgctgcAATTCAATCCGCGGTGTGTGATTGAGTAGCTTTTCAGTTGTAAGAGAGGTTTCCAGCGCCACGGTTTCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCCGCGTTGTGGTGTGGTAAAACGAAGAAACGGTAAGAATGTGGCGTGCAGGGATACCGCCTAATGCTTTCAGG‑GCAGCATTGACCTGCTCCAGCGTAGAGTCCAGATTTCTTCCTATGGCAATATACGCCACGGTCATGCGGTAcc                                                                                              <  2:522147/250‑1 (MQ=255)                                                                                                                                                      tgctgcAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTAAGG‑GCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATAAGCCACTGTcctgaggttccc                                                                                             >  2:455875/1‑240 (MQ=255)                                                                                                                                                      tgctgcAATTCAATACGCTGTGTGTGATTGAGTAGCTCCTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATCTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGG‑GCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACGCCACTGTCATGCGGTAcc                                                                                              <  1:173744/250‑1 (MQ=255)                                                                                                                                                           caagtCAATACGCTGTGTGTGATTGATTAGCTCGTCAGGTGCAAGAGAGGTTCCCAGCGCCACGGCTGCTTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCGTTCAGG‑GCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACGCCACTGTCATGCGGTACCCTctcg                                                                                        <  2:381700/247‑4 (MQ=255)                                                                                                                                                            aaTTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGG‑GCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACGCCACTGTCATGCGGTACCCTCacga                                                                                       <  1:406274/251‑1 (MQ=255)                                                                                                                                                                            tgtgATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAAtgt                                                                                                                                                                                                        <  1:19830/123‑1 (MQ=255)                                                                                                                                                                            tgtgATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAAtgt                                                                                                                                                                                                        >  2:19830/1‑123 (MQ=255)                                                                                                                                                                            tgtgATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGG‑GCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCaa                                                                                                                      >  1:404424/1‑204 (MQ=255)                                                                                                                                                                            tgggATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGG‑GCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCaa                                                                                                                      <  2:404424/201‑1 (MQ=255)                                                                                                                                                                                aTTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGG‑GCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACGCCACTGTCATGCGGTACCCTCACGACGTGGTGCGCGTT‑TGCGTgg                                                                  >  1:159996/1‑251 (MQ=255)                                                                                                                                                                                       aGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGG‑GCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAAtata                                                                                                                  >  2:278295/1‑197 (MQ=255)                                                                                                                                                                                       aGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGG‑GCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAAtata                                                                                                                  <  1:278295/197‑1 (MQ=255)                                                                                                                                                                                          tctTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGG‑gcag                                                                                                                                                                   >  1:63775/1‑145 (MQ=255)                                                                                                                                                                                          tctTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGG‑gcag                                                                                                                                                                   <  2:63775/145‑1 (MQ=255)                                                                                                                                                                                                ggTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGG‑gcag                                                                                                                                                                   <  2:239953/139‑1 (MQ=255)                                                                                                                                                                                                ggTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGG‑gcag                                                                                                                                                                   >  1:239953/1‑139 (MQ=255)                                                                                                                                                                                                      agagagGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGG‑GCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACGCCACTGTCATGCGGTACCCTCACGACGTGGTGCGCGTT‑TGCGTGGACGACGAGTACGACGACGCGGt                                            <  1:260643/251‑1 (MQ=255)                                                                                                                                                                                                      agagagGTTTCCAGCGCCACGGCTGCGGTTAAGTAAGCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAAAGCTTTCAGG‑GCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGTCTATGGCAATATACGCCACTGTCATGCGGTACCCTCACGACGTGGTGCGCGTT‑TGCGTGGACGACGACTACGACGACGCGGt                                            <  2:127033/251‑1 (MQ=255)                                                                                                                                                                                                           ggTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGG‑GCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACGCCACTGTCATGCGGTACCCTCACGACGTGGTGCGCGTT‑TGCGTGGACGACGAGTACGACGACGCGGTGACgg                                       <  2:116259/251‑1 (MQ=255)                                                                                                                                                                                                             tttCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGG‑GCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACGCCACTGTCATGCGGTACCCTCACGACGTGGTGCGCGTT‑TGCGTGGACGACGAGTACGACGACGCGGTGACGGtt                                     >  1:351042/1‑251 (MQ=255)                                                                                                                                                                                                                cctgcgcCTCGGCTGCTTTTAAGTAATCGGGTTGATTTTTCTGCCGCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTGCAGGGATACCGCCTACTGTTTTCAGG‑GCAGCATTGACCCGCTCCAGCTGAGAGGCCAGATTGCTGCCTATGGCAATATACGCCACTGTCATGCGGGACCCTCCCGACGTGGTTCGCGTT‑TGCGTGGACGACGAGTACGACGACGCGGTGACgtgtct                                   <  2:336501/247‑6 (MQ=255)                                                                                                                                                                                                                     gccaagccTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAGGAAACGGTAAGAATGGGGCTTTCAGGGATATCGCATAATGCTTTCAGG‑GCAGCATTGAACTGCTCCAGCGGAGAGGCAAGATTGCTGCCTCTGGCCATATACTCCCCTGTCCTGCGGAACCCTCACGACGGGGTGCGCGGT‑TGAGTAGACGACGAGTACGACGACGCGGTGGCGGTTCTTCATcc                             >  2:332644/8‑251 (MQ=255)                                                                                                                                                                                                                      ccACGGCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGG‑GCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACGCCACTGTCATGCGGTACCCTCACGACGTGGTGCGCGTT‑TGCGTGGGCGACGAGTACGACGACGCGGTGACGGTTCTTCATcc                             >  1:67495/1‑250 (MQ=255)                                                                                                                                                                                                                      caaCGGCTGCGTTTCAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAGACGATAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGG‑GCAGCATTGACCTGCTCCAGCGGAGAGGCGAGATTGCTGCCTATGGCAGTATACGCCACTGTCATGGGGTGCCCTCACGCCGTGGTGCGAGTG‑TGCGTGGACGACGCGGAGGACGACGCGGTGCAGGGTCTTCATccc                            >  2:235771/3‑250 (MQ=255)                                                                                                                                                                                                                       cACGGCTGCGTTTAAGTAATCGGGTTGAGTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGG‑GCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACGCCACTGTCATGCGGTACCCTCACGACGTGGTGCGCGTT‑TGCGTGGACGACGAGTACGACGACGCGGTGACGGTTCTTCATCCAg                           <  1:361507/251‑1 (MQ=255)                                                                                                                                                                                                                          ggCTGCGTTTAAGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGG‑GCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACGCCACTGTCATGCGGTACCCTCACGACGTGGTGCGCGTT‑TGCGTGGACGACGAGTACGACGACGCGGTGACGGTTCTTCATCCAGCTc                        <  1:481639/251‑1 (MQ=255)                                                                                                                                                                                                                                     aGTAATCGGGTTGATTTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCATAATGCTTTCAGG‑GCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACGCCACTGTCAGGCGGTACCCTCACGACGTGGTGCGCGTT‑TGCGTGGCCGCCGAGTACGACGACGCGGTGACGGGTCTTCATCCAGCTCGTTGAGCATcc             >  2:486548/1‑251 (MQ=255)                                                                                                                                                                                                                                                tGATTTTGCGTCCACAGAGGTGGGGGGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATGTCGCCTAATGCTTTAAGT‑GCAGCATTGCACTGCTCCAGCGGAAAGGCCAGATTGCTGCCTATGGCAATATACGCCACTGTCATGCGGTACCCTCACGACGTGGTGCGCGTTGTG‑GTGGACGCCGAGTGCGCAGCCGAGGGGACGGTTTTTCATCCAGCTCGTTGTGAATCCCATTTTgctat  >  2:476050/1‑247 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                         CAAACAGCGGCCACAGCATAAATCCACGATTCTTCATATCGTAGTGCGGAACGGTCAGGCGTTCAGTATTTATCACTTCATTACCAAACAGCATGATGTCGAGATCCAGCGTGCGTGGTCCCCAGCGTTCAGCTTTGCGGACGCGACCTTGCTGCAATTCAATACGCTGTGTGTGATTGAGTAGCTCTTCAGGTGCAAGAGAGGTTTCCAGCGCCACGGCTGCGTTTAAGTAATCGGGTTGATCTTGCGGCCCCAGCGGTGGGGTGCGGTAAAACGAAGAAACGGTAAGAATGTGGCTTTCAGGGATATCGCCTAATGCTTTCAGG‑GCAGCATTGACCTGCTCCAGCGGAGAGGCCAGATTGCTGCCTATGGCAATATACGCCACTGTCATGCGGTACCCTCACGACGTGGTGCGCGTT‑TGCGTGGACGACGAGTACGACGACGCGGTGACGGTTCTTCATCCAGCTCGTTGAGCATCCCTTTTTGGTCT  >  NC_000913/157342‑157831

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