Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	82,633	T→G	I299I (ATA→ATC)	leuA ←	2‑isopropylmalate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	82,633	0	T	G	95.6%	137.8 / NA	45	I299I (ATA→ATC)	leuA	2‑isopropylmalate synthase
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base G (24/19);  total (25/20)
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.67e-01

ATGAAGGCCAGCGCCTCCAGATCGTAATCAAACACCTGACCTTTTTTGTCCGCCAGCTTCAGGAAAGCATCGTACAAATTGTCTAAATTATATTCACTTTCTTTATACCCCATCTCATCCATGCGATGTTTCACCGCCGCACGCCCCGAACGAGAGGTCAGATTCAGCTGGATTTGGTTCAGACCAATAGATTCTGGTGTCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGTATACCGGAGGAGTGTGCGAATGCGCCGCTGCCAACAATGGCTTTGTTTGCCGGGATCGGCATATTACAAATCTGGCTAACTAACTGGCTGGTGCGCCATATCTCCTGGTGATTAATGGCGGTGTGGACGTTGAGAATATCCTTACGAAC‑TTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTGCGCCTTCCACCTGGCGTGCACCGGCAT  >  NC_000913/82388‑82874                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   aTGAAGGCCAGCGCCTCAAGATCGTAATCAAACACCTGACCTTTTTTGTCCGCCAGCTTCAGGAAAGCATCGTACAAATTGTCTAAATTATATTCCCTTTCTTTATACCCCATCTCATCCATGCGATGTTTCACCGCCGCACGCCCCGAACGAGAGGTCAGATTCAGCTGGATTTGGTTCAGACCAATAGATTCTGGTGTCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGGATAcc                                                                                                                                                                                                                                               <  1:460949/251‑1 (MQ=255)   gAAGGCCAGCGCCTCCAGATCGTAATCAAACACCTGACCTTTTTTGTCCGCCAGCTTCAGGAAAGCATCGTACAAATTGTCTAAATTATATTCACTTTCTTTATACCCCATCTCATCCATGCGATGTTTCACCGCCGCACGCCCCGAACGAGAGGTCAGATTCAGCTGGATTTGGTTCAGACCAATAGATTCTGGTGTCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGGATACCgg                                                                                                                                                                                                                                             >  1:165264/1‑251 (MQ=255)   gAAGGCCAGCGCCTCCAGATCGTAATCAAACACCTGACCTTTTTTGTCCGCCAGCTTCAGGAAAGCATCGTACAAATTGTCTAAATTATATTCACTTTCTTTATACCCCATCTCATCCATGCGATGTTTCACCGCCGCACGCCCCGAACGAGAGGTCAGATTCAGCTGGATTTGGTTCAGACCAATAGATTCTGGTGTCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGGATACCg                                                                                                                                                                                                                                              <  1:273657/250‑1 (MQ=255)   gAAGGCCAGCGCCTCCAGATCGTAATCAAACACCTGACCTTTTTTGTCCGCCAGCTTCAGGAAAGCATCGTACAAATTGTCTAAATTATATTCACTTTCTTTATACCCCATCTCATCCATACGATGTTTCACCGCCGCACGCCCCGAACGAGAGGTCAGATTCAGCTGGATTTGGTTCAGACCAATAGATTCTGGTGTCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGGATACCgg                                                                                                                                                                                                                                             >  1:38636/1‑251 (MQ=255)         cAGCGCCTCCATATCTAAATCAAATACTTGACTTTTTTTGTCCGCCAGCTTCAGGAAAGCTTCGGCCAAATTGTTTAAATTTGAGACACTGTCTTTATACCCCACCTCATCCAGTCGATGTTGCACCGCCGCACGCCCCGAGTGAGTGGTCAGATTCAGCTGGATTTGGTTCAGACCAATAGATTCTGGTGTCATGATTTCGTAGTTTTCGCGGGTTTACAGCACGCCTTCCTGGTGTATACCGCAGGAgt                                                                                                                                                                                                                                       <  2:315289/251‑1 (MQ=255)             gcCTCCAGATCGTAATCAAACACCTGACCTTTTTTGTCCGCCAGCTTCAGGAAAGCATCGTACAAATTGTCTAAATTATATTCACTTTCTTTATACCCCATCTCATCCATGCGATGTTTCACCGCCGCACGCCCCGAACGAGAGTTCAGTTCCAGTTGGATTTGTTTCTGACCAATAGATTCTGGGGTCATGATTTCGAAGTTTTCGCGTTTTTTAAGCACGCAATCCTGGTGGTTACCGGAGGAGggtgg                                                                                                                                                                                                                                   >  2:291402/1‑250 (MQ=255)                 caaGATCGTAATCAAACACCTGACCTTTTTTGTCCGCCAGCTTCAGGAAAGCATCGTACAAATTGTCTAAATTATATTCCCTTTCTTTATACCCCAGCTCATCCAGGCGATGTTTCACCGCCGCACGCCCCGAACGAGAGGTCAGATTCAGCTGGATTTGGTTCAGACCAATAGATTCTGGTGTCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGGATACCGGAGGAGTGTGCGa                                                                                                                                                                                                                                 <  1:448511/247‑1 (MQ=255)                  cagctctaaatcaatcCCCTTACCTTTTTTGTCCGCCTGCTTCAGTTAAGCATCGTACAAATTGTCTAAATTTTATTCTCTTTCTTTATACCCCATCTCATCCATTCGATGTTTCACCGCCGCACGCCCCGAACGAGATGTCTGATTCAGCTGGATTTGGTTCAGTCCAATAGATTCTGGTGTCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGGATACCGGAGGAGTGTGCGAATg                                                                                                                                                                                                                              <  2:83026/243‑1 (MQ=255)                      tCGTAATCAAACACCTGTCCTTTTTTGTCCGCCAGCTTCAGGAAAGCATCGTACAAATTGTCTAAATTATATTCACTTTCTTTATACCCCATCTCATCCATGCGATGTTTCACCGCCGCACGCCCCGAACGAGAGGTCAGATTCAGCTGGATTTGGTTCAGACCCATAGATTCTGGTGTCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGGATACCGGTGGAGTGTGCGAGTGcgcc                                                                                                                                                                                                                          >  2:396792/1‑251 (MQ=255)                         tAATCAAACACCTGACCTTTTTTGTCCGCCAGCTTCAGGAAAGCATAGTACAAATTGTCTAAATTATATTCACTTTCTTTATACCCCATCTCATCCATGCGATGTTTCACCGCCGCACGCCCCGAACGAGAGGTCAGATTCAGCTGGATTTGGTTCAGACCAATAGATTCTGGTGTCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGACATCCTGGTGGATACCGGAGGAGTGTGCGAATGcgccgcc                                                                                                                                                                                                                       >  2:336200/1‑250 (MQ=255)                                  acCTGACCTTTTTTGTCCGCCAGCTTCAGGAAAGCATCGTACAAATTGTCTAAATTATATTCACTTTCTTTATACCCCATCTCATCCATGCGATGTTTCACCGCCGCACGCCCCGAACGAGAGGTCAGATTCAGCTGGATTTGGTTCAGACCAATAGATTCTGGTGTCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGGATACCGGAGGAGTGTGCGAATGCGCCGCTGCCAACAAt                                                                                                                                                                                                              >  1:102018/1‑251 (MQ=255)                                     tGACCTTTTTTGTCCGCCAGCTTCAGGAAAGCATCGTACAAATTGTCTAAATTATATTCACTTTCTTTATACCCCATCTCATCCATGCGATGTTTCACCGCCGCACGCCCCGAACGAGAGGTCAGATTCAGCTGGATTTGGTTCAGACCAATAGATTCTGGTGTCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGGATACCGGAGGAGTGTGCGAATGAGCCGCTGCCAACAATGGc                                                                                                                                                                                                           >  1:130651/1‑251 (MQ=255)                                      gACCTTTTTTGTCCGCCAGCTTCAGGAAAGCATCGTACAAATTGTCTAAATTATATTCACTTTCTTTATACCCCATCTCATCCATGCGATGTTTCACCGCCGCACGCCCCGAACGAGAGGTCAGATTCAGCTGGATTTGGTTCAGACCAATAGATTCTGGAGTCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGGATACCGGAGGAGTGTGCGAATGCGCCGCTGCCAACAATGGCt                                                                                                                                                                                                          >  1:81975/1‑251 (MQ=255)                                                                ttccccacccgccccttttccgttctataTTCCCTTTCTTTATCCCCCATCCCATCCAGGCGATGTTTCACCGCCGCACGCCCCGAACGAGAGGTCAGATTCAGCTGGATTTGGTTCAGACCAATAGATTCTGGTGTCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGGATACCGGAGGAGTGTGCGAATGCGCCGCTGCCAACAATGGCTTTGTTTGCCGGGATCGGCATATTAc                                                                                                                                                                                 <  2:100731/225‑1 (MQ=255)                                                                                          atatTCACTTTCTTTATACCCCATCTCATCCATGCGATGTTTCACCGCCGCACGCCCCGAACGAGAGGTCAGATTCAGCTGGATTTGGTTCAGACCAATAGATTCTGGTGTCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGGATACCGGAGGAGTGTGCGAATGCGCCGCTGCCAACAATGGCTTTGTTTGCCGGGATCGGCATATTAc                                                                                                                                                                                 >  1:100731/1‑224 (MQ=255)                                                                                                actttcggtatcccacATCTCATCCATGCGATGTTTCACAGCCGCACGCACCGAACGAGAGGTCAGATTCAGCTGGATTTGGATCAGACCAATAGATTCTGGGGTCATGATTTCGTAGTTTTCGAGGTGTTTCAGCACGACATCCTGGTGGAAAACGGAGGAGTGTGCGAATGAGCCGCTGCCAAAAATGGCGTTGTTTGCCGGGCTCGGCATATTACACATCTGGCTAACTCACTGGCTGGAGCACCat                                                                                                                                                 >  2:50133/16‑250 (MQ=255)                                                                                                        tataCCACAACTCATCCCTGCGATGTTTCACCGCCGCACCCCCCGACCGACAGGTCAGTTTAAGCTGGATTTGGTTCAGACCACTCGATTATGGTGCCATGCTTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCTGTTGTATACCGGTGGTGTGTGCAAATGCGCCGCTGCCAACAATGGCTTTGTTTGCCGGGCTCGGCTTACTACAACTCAGGCCACCTAACGGGCTGGTGAGCCaccccacctgg                                                                                                                                        >  2:309118/1‑241 (MQ=255)                                                                                                              cccaccccaccATGCGATTTTTTCCCGCCGCACGCCCCGAACGAGAGGTCAGATTCAGCTGGATTTGGTTCAGACCAATAGATTCTGGTGTCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGGATACCGGAGGAGTGTGCGAATGCGCCGCTGCCAACAATGGCTTTGTTTGCCGGGATCGGCATATTACAAATCTGGCTAACTAACTGGCTGGTGCGCCATATCTCCTGGTGATTa                                                                                                                                  <  1:291402/242‑1 (MQ=255)                                                                                                              ccATCTCATCCATGCGATGTTTCACCGCCGCACGCCCCGAACGAGAGGTCAGATTCAGCTGGATTTGGTTCAGACCAATAGATTCTGGTGTCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGGATACCGGAGGAGTGTGCGAATGCGCCGCTGCCAACAATGGCTTTGTTTGCCGGGATCGGCATATTACAAATCTGGCTAACTAACTGGCTGGTGCGCCATATCTCCTGGTGAtt                                                                                                                                   >  1:17658/1‑250 (MQ=255)                                                                                                              acATCTCATCCATGCGATGATTCACCGCCGCACGCCACGAACGAGAGGTCAGATACAGCTGGATTTGGTTCCGACCCATAGATTCTGGTGTCATGATTTCGTAGTTTTCGCGGTTTTTACGCACGCCATCATGGTGGATACCGGAGGAGTGTGCGAATGCGCCGCTGCCAACAATGGCTTTGTTTGCCGGGATCGGCATATTACAAATCTGGCTAACTAACTGGCTGGGGCGCCatcttctctggtgattt                                                                                                                                  >  2:249328/2‑236 (MQ=255)                                                                                                                        cATGCGATGTTTCACCGCCGCACGCCCCGAACGAGAGGTCAGATTCAGCTGGATTTGGGTCAGACCAATAGATTCTGGTGTCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGGATACCGGAGGAGTGTGCGAATGCGCCGCTGCCAACAATGGCTTTGTTTGCCGGGATCGGCATATTACAAATCTGGCTAACTAACTGGCTGGTGCGCCATATCTCCTGGTGATTAATGGCgggtg                                                                                                                        >  2:181450/1‑248 (MQ=255)                                                                                                                             gATGTTTAGCCGCCGCACGCCCCGAACGAGAGGTCAGATTCAGCTGGATTTGGTTCAGACCAATAGATTCTGGTGGCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGGATACCGGAGGAGTGTGCGAATGCGCCGCTGCCAACAATGGCTTTGTTTGCCGGGATCGGCATATTACAAATCTGGCTAACTAACTGGCTGGTGCGCCATATCTCCTGGTGATTAATGGCGGTGTGGACg                                                                                                                   <  2:73458/251‑1 (MQ=255)                                                                                                                                      ccgccgCACGCCCCGAACGAGAGGTCAGATTCAGCTGGATTTGGTTCAGACCACTAGATTCTGGTGTCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGGATACCGGAGGAGTGTGCGAATGCGCCGCTGCCAACAATGGCTTTGTTTGCCGGGATCGGCATATTACAAATCTGGCTAACTAACTGGCTGGTGCGCCATATCTCCTGGTGCTTAATGGCGGGGTGGACGtgtggaaa                                                                                                           >  2:439844/1‑243 (MQ=255)                                                                                                                                           gcACGCCCCGAACGAGAGGTCAGATTCAGCTGGATTTGGTTCAGACCAATAGATTATGGTGTCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGGATACCGGAGGAGTGTGCGAATGCGCCGCTGCCAACAATGGCTTTGTTTGCCGGGATCGGCATATTACAAATCTGGCTAACTAACTGGCTGGTGCGCCATATCTCCTGGTGATTAATGGCGGTGTGGACGGTGAGACAATCCtt                                                                                                     >  1:310333/1‑251 (MQ=255)                                                                                                                                                                  aTTCAGCTGGATTTGGTTCAGACCAATAGATTCTGGTGTCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGGATACCGGAGGAGTGTGCGAATGCGCCGCTGCCAACAATGGCTTTGTTTGCCGGGATCGGCATATTACAAATCTGGCTAACTAACTGGCTGGTGCGCCATATCTCCTGGTGATTAATGGCGGTTTGGAAGTTTGGGATATCCTTAAAAAC‑ttttgtagcccttgtgg                                                                             >  1:221521/1‑237 (MQ=255)                                                                                                                                                                   ttCAGCTGGATTTGGTTCAGACCAATAGATTCTGGTGTCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGGATACCGGAGGAGTGTGCGAATGCGCCGCTGCCAACAATGGCTTTGTTTGCCGGGATCGGCATATTACAAATCTGGCTAACTAACTGGCTGGTGCGCCATATCTCCTGGTGATTAATGGCGGTGGGGAGGGTGGGAATATTCTTACGAACATTTG‑TCGCCAAGATGAc                                                                            >  1:74390/1‑251 (MQ=255)                                                                                                                                                                     cAGCTGGATTTGGGGCAGACCAATAGATTCTGGTGTCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGGATACCGGAGGAGTGTGCGAATGCGCCGCTGCCAACAATGGCTTTGTTTGCCGGGATCGGCATATTACAAATCTGGCTAACTAACTGGCTGGTGCGCCATATCTCCTGGTTATTAATGGCGGTGTGAACGGTGAGAATATCCTTACGAAC‑TTTGATCGCCATGATGAct                                                                           <  2:310333/250‑1 (MQ=255)                                                                                                                                                                       gCTGGATTTGGTTCAGACCAATAGATTCTGGTGTCATGATGTTGTAGTTTTCGCGGTTTTTCAGAACGGCATCCTGGTGGATACCGGAGGAGTGTGCGAATGCGCCGATGCCAACAAAGGATTTGTTGGCCGGGATCGGCATATTACAAATCTGGCTAACTAACTGGCTGGTGCGCCATATCTCCTGGTGATTAATGGCGGTGTGGACGTTGAGAATATCCTTACGAAC‑TTTGATCGCCATGATGActttt                                                                        <  2:135960/251‑3 (MQ=255)                                                                                                                                                                           gATTTGGTTCAGACCAATAGATTCTGGTGTCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGGATACCGGAGGAGTGTGCGAATGCGCCGCTGCCAACAATGGCTTTGTTTGCCGGGATCGGCATATTACAAATCTGGCTAACTAACTGGCTGGTGCGCCATATCTCCTGGTGATTAATGGCGGTGTGGACGTTGAGAATATCCTTACGAAC‑TTTGATCGCCATGATGACTTCTTCCa                                                                    <  1:439844/251‑1 (MQ=255)                                                                                                                                                                                   tCAGACCAATAGATTCTGGTGTCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGGATACCGGAGGAGTGTGCGAATGCGCCGCTGCCAACAATGGCTTTGTTTGCCGGGATCGGCATATTACAAATCTGGCTAACTAACTGGCTGGTGCGCCATATCTCCTGGTGATTAATGGCGGTGTGGACGTTGAGAATATCCTTACGAAC‑TTTGATCGCCATGATGACTTCTTCCAGGGAACa                                                             <  1:355970/250‑1 (MQ=255)                                                                                                                                                                                    cAGACCAATAGATTCTGGTGTCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGGATACCGGAGGAGTGTGCGAATGCGCCGCTGCCAACAATGGCTTTGTTTGCCGGGATCGGCATATTACAAATCTGGCTAACTAACTGGCTGGTGCGCCATATCTCCTGGTGATTAATGGCGGTGTGGACGTTGAGAATATCCTTACGAAC‑ttt                                                                                           <  1:191377/219‑1 (MQ=255)                                                                                                                                                                                      gACCAATAGATTCTGGTGTCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGGATACCGGAGGAGTGTGCGAATGCGCCGCTGCCAACAATGGCTTTGTTTGCCGGGATCGGCATATTACAAATCTGGCTAACTAACTGGCTGGTGCGCCATATCTCCTGGTGATTAATGGCGGTGTGGGCGTTGAGAATATCCTTTCGAAC‑TTTTATCGCCATGATGACTTTTTTCCGGGAACAGttt                                                         >  1:121077/1‑251 (MQ=255)                                                                                                                                                                                      gACCAATAGATTCTGGTGTCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGGATACCGGAGGAGTGTGCGAATGCGCCGCTGCCAACAATGGCTTTGTTTGCCGGGATCGGCATATTACAAATCTGGCTAACTAACTGGCTGGTGCGCCATATCTCCTGGTGATTAATGGCGGTGTGGACGTTGAGAATATCCTTACGAAC‑TTTGATCGCCATGATGACTTCTTCCAGGGAACAGttt                                                         <  1:320015/251‑1 (MQ=255)                                                                                                                                                                                        ccAATAGATTCTGGTGTCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGGATACCGGAGGAGTGTGCGAATGCGCCGCTGCCAACAATGGCTTTGTTTGCCGGGATCGGCATATTACAAATCTGGCTAACTAACTGGCTGGTGCGCCATATCTCCTGGTGATTAATGGCGGTGTGGACGTTGAGAATATCCTTACGAAC‑TTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTcc                                                       <  1:232135/251‑1 (MQ=255)                                                                                                                                                                                            tAGATTCTGGTGTCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGGATACCGGAGGAGTGTGCGAATGCGCCGCTGCCAACAATGGCTTTGTTTGCCGGGATCGGCATATTACAAATCTGGCTAACTAACTGGCTGGTGCGCCATATCTCCTGGTGATTAATGGCGGTGGGGGAGGTGAGAATATCCTTACGGAC‑TTTGATCGCCAATATGACTTCTTcccaggaaccgattcccgcc                                                   >  1:65107/1‑233 (MQ=255)                                                                                                                                                                                                    gttgtCATGATTTCGTAGTTTTCGCGGTTTTTACGCCCGCCCGCCTGTTGGATACCGGAGGAGGGTGCGAATGCGCCTCTGCCAACAATGGCTTTGGTTGCCGGGATCGGCATATTACAAATCTGGCTAACTAACTGGCTGGTGCGCCATATCTCCTGGTGATTAATGGCGGTGTGGACGTTGAGAATATCCTTACGAAC‑GTTGATCGCCATGATGACTTCGTCCAGGGTACAGTTTCCGGCACGCTc                                              <  2:58446/246‑1 (MQ=255)                                                                                                                                                                                                       gtCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGGATACCGGAGGAGTGTGCGAATGCGCCGCTGCCAACAATGGCTTTGTTTGCCGGGATCGGCATATTACAAATCTGGCTAACTAACTGGCTGGTGCGCCATATCTCCTGGTGATTAATGGCGGTGTGGGCGTTGAGAATATCCTTACGGAC‑TTTGATCGCCCTGATAACTTTTTTCAGGGGACAGTTTCCGGCCCGCTCGcccat                                        >  1:134682/1‑248 (MQ=255)                                                                                                                                                                                                                  gTAGGTTTCGTGGTTTTTCAGCTCGCCATCCTGGTGGATACCGGAGTAGTGTGCGAATGCGCCGCTGCCAACAATGGCTTTGTTTGCCGGGATCGGCATATTACAAAACTGGCTAACTAACTGGCTGGTGCGCCATATCTCCTGGTGATTAATGGCGGTGTGGACGTTGAGAATATCCTTACGAAC‑TTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCGCGCCGATCCCGTTCATTg                             <  2:133145/251‑1 (MQ=255)                                                                                                                                                                                                                    aGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGGATACCGGAGGAGTGTGCGAATGCGCCGCTGCCAACAATGGCTTTGTTTGCCGGGATCGGCATATTACAAATCTGGCTAACTAACTGGCTGGTGCGCCATATCTCCTGGTGATTAATGGCGGTGTGGACGTTGAGAATATCCTTACGGAC‑TTTGATCGCCATGATGACTTCTTCCCGGGAACAGTTTCCGGCACGCTCGCCGATCCCCTTCCTTgcg                           >  1:26266/1‑251 (MQ=255)                                                                                                                                                                                                                    aGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGGAGTCCGGAGGAGTGGGCGAATGCGCCGCTGCCAACAATGGCTTTGTGTGCCGGGATCGGCATAGTACAAATCTGGCTAACTAACTGGCTGGTGCGCCATATCTCCTGGTGATTAATGGCGGTGTGGACGTTGAGAATATCCTTACGAAC‑TTTGATCGCCATGATGACTTCTTCCAGGGAACAGGTTCCGGCACGCTCGCCGATCCCTGTCATTgc                            <  2:42780/250‑1 (MQ=255)                                                                                                                                                                                                                     gTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGGATACTGGAGGAGTGTGCGAATGCGCCGCTGCCAACAAGGGCTTTGTTTGCCGGGATCGGCAGATTACAAATCTGGCTAACTAACTGGCTGGTGCGCCTTATCTCCTGGTGATTAATGGCGGTGTGGACGTTGAGAATATCCTTACGAAC‑TTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCTCGTTCATTgcg                           <  2:134682/250‑1 (MQ=255)                                                                                                                                                                                                                                          gCCATCCTGGTGGATACCGGAGGAGGGTGCGAAGGCGCCTCTGCCAACAATGGCTTTGTTTGCCGGGATCGGCATATTACAAATCTGGCGAACTAACTGGCTGGTGCGCCATATCTCCTGGTGATTAATGGCGGTGTGGACGTTGAGAATATCCTTACGAAC‑TTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTGCGCCTTCCACCTGGCGTGCACCgg     <  2:257985/251‑1 (MQ=255)                                                                                                                                                                                                                                           ccATCCTGGTGGATACCGGAGGCGTGTGCGAATGCGCCGCTGCCAACAATGGCTTTGTTTGCCGGGATCGGCATATTACAAATCTGGCTAACTAACTGGCTGGTGCGCCATATCTCCTGGTGATTAATGGCGGGGTGGGCGGTTAGAAAATCCTTTAGCAC‑TTTTATTGCCATGGGGGCTTCTTCCCTGGTACAGATTTCAGCCCGCCCGCCGATCCCGTTCATTTTGCCCTCCAcccggggtgacccggc    >  2:391869/1‑237 (MQ=255)                                                                                                                                                                                                                                           ccATCCTGGTGGATACCGGAGGAGTGTGCGAATGCGCCGCTGCCCACAATGGCTTTGTTTGCCGGGATCGGCATATTACAAATCTGGCTAACTAACTGGCTGGTGCGCCATATCTCCTGGTGATTAATGGCGGTGTGGGCGTTGAGAATATCCTTACGAAC‑TTTTATAGCCATGATGACTTTTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCCTTCATTGCGCCTTCCACCTGGCGTGCACCGGc    >  1:414747/1‑251 (MQ=255)                                                                                                                                                                                                                                             aTCCTGGTGGATACCGGAGGAGTGTGCGAATGCGCCGCTGCCAACAATGGCTTTGTTTGCCGGGATCGGCATATTACAAATCTGGCTAACTAACTGGCTGGTGCGCCATATCTCCTGGTGATTAATGGCGGTGTGGACGTTGAGAATATCCTTACGAAC‑TTTGATCGCCATGATGACTTCTTCCAGGGAAAAGTTTCCGGCACGCCCGCCGATCCCGTTCATTGCGCCTTCCCCCTGGCGGGCACCGGCAt  >  1:278203/1‑251 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   ATGAAGGCCAGCGCCTCCAGATCGTAATCAAACACCTGACCTTTTTTGTCCGCCAGCTTCAGGAAAGCATCGTACAAATTGTCTAAATTATATTCACTTTCTTTATACCCCATCTCATCCATGCGATGTTTCACCGCCGCACGCCCCGAACGAGAGGTCAGATTCAGCTGGATTTGGTTCAGACCAATAGATTCTGGTGTCATGATTTCGTAGTTTTCGCGGTTTTTCAGCACGCCATCCTGGTGTATACCGGAGGAGTGTGCGAATGCGCCGCTGCCAACAATGGCTTTGTTTGCCGGGATCGGCATATTACAAATCTGGCTAACTAACTGGCTGGTGCGCCATATCTCCTGGTGATTAATGGCGGTGTGGACGTTGAGAATATCCTTACGAAC‑TTTGATCGCCATGATGACTTCTTCCAGGGAACAGTTTCCGGCACGCTCGCCGATCCCGTTCATTGCGCCTTCCACCTGGCGTGCACCGGCAT  >  NC_000913/82388‑82874

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