Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,520,199	A→G	G158G (GGT→GGC)	gltX ←	glutamate‑‑tRNA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,520,199	0	A	G	80.3%	131.4 / 20.2	62	G158G (GGT→GGC)	gltX	glutamate‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base A (7/5);  new base G (25/24);  total (32/30)
Fisher's exact test for biased strand distribution p-value = 7.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.18e-01

TTGGACAGTTTTTTACCGTCATCGCCATTGATCATAGAAACGTGCGCGTAAACCGGCACCGGCGCTTTCAGGGCCTTAAGAATGTTGATCTGGCGTGGCGTGTTGTTGATATGGTCTTCGCCACGGATAACGTGGGTGATTTCCATATCCCAGTCATCGACAACCACACAGAAGTTATAGGTTGGGGAACCATCGGTACGGCGGATAATAAGATCGTCCAGTTCCTGGTTGCTGAACTCGATCGGACCACGGATCTGATCGTCAAAAACAACAGAACCTTCCTGCGGGTTAGCAAAACGTACAACACACGGTTCATCATCAGCATGATGCTCATGGCTGTGGCGGCAGCGACCGTCATAACGCGGCTTCTCACCTTTCGCCATTTGCTCTTCGCGCAGCGCTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTG  >  NC_000913/2519954‑2520445                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       ttGGACAGTTTTTTACCGTCATCGCCATTGATCATAGAAACGTGCGCGTAAACCGGCACCGGCGCTTTCAGGGCCTTAAGAATGTTGATCTGGCGTGGCGTGTTGTTGATATGGTCTTCGCCACGGATAACGTGGGTGATTTCCATATCCCAGTCATCGACAACCACACAGAAGTTATAGGTTGGGGAACCATCGGTACGGCGGATAATAAGATCGTCCAGTTCCTGGTTGCTGAACTCGATCGGGCCAc                                                                                                                                                                                                                                                    <  1:403540/250‑1 (MQ=255)      cAGTTTTTTACCGTCATCGCCATTGATCATAGAAACGTGCGCGTAAACCGGCACCGGCGCTTTCAGGGCCTTAAGAATGTTGATCTGGCGTGGCGTGTTGTTGATATGGTCTTCGCCACGGATAACGTGGGTGATTTCCATATCCCAGTCATCGACAACCACACAGAAGTTATAGGTTGGGGAACCATCGGTACGGCGGATAATAAGATCGTCCAGTTCCTGGTTGCTGAACTCGATCGGGCCACGGATCt                                                                                                                                                                                                                                              >  1:262809/1‑251 (MQ=255)       aGTTTTTTACCGTCATCGCCATTGATCATAGAAACGTGCGCGTAAACCGGCACCGGCGCTTTCAGGGCCTTAAGAATGTTGATCTGGCGTGGCGTGTTGTTGATATGGTCTTCGCCACGGATAACGTGGGTGATTTCCATATCCCAGTCATCGACAACCACACAGAAGTTATAGGTTGGGGAACCATCGGTACGGCGGATAATAAGATCGTCCAGTTCCTGGTTGCTGAACTCGATCGGGCCACGGATCTg                                                                                                                                                                                                                                             <  1:360837/251‑1 (MQ=255)          tttttACCGTCATCGCCATTGATCATAGAAACGTGCGCGTAAACCGGCACCGGCGCTTTCAGGGCCTTAAGAATGTTGATCTGGCGTGGCGTGTTGTTGATATGGTCTTCGCCACGGATAACGTGGGTGATTTCCATATCCCAGTCATCGACAACCACACAGAAGTTATAGGTTGGGGAACCATCGGTACGGCGGATAATAAGATCGTCCAGTTCCTGGTTGCTGAACTCGATCGGGCCACGGATCTGAt                                                                                                                                                                                                                                           <  1:349458/250‑1 (MQ=255)          tttctacctgcATCACCATTGAGCATAGAAACGTGCGCGTCAACCGGCCCCGTCGCTTTGAGCGCCTTAAGAATGTTGATCTGGCGTGGCGTGTTGTTGATATGGTCTGCGCCACGGATAACGGGGGTGATTTCCATATCCCAGTCATCGACAACCACACAGAAGTTATAGTTTGGGAAACCATCGGTACGGCGGATAATAAGATCGTCCAGTTCCTGGTTGCTGAACTCGATCGGACCACGGATCTGAt                                                                                                                                                                                                                                           <  2:284618/240‑1 (MQ=255)              tACCGTCATCGCCATTGATCATAGAAACGTGCGCGTAAACCGGCACCGGCGCTTTCAGGGCCTTAAGAATGTTGATCTGGCGTGGCGTGTTGTTGATATGGTCTTCGCCACGGATAACGTGGGTGATTTCCATATCCCAGTCATCGACAACCACACAGAAGTTATAGGTTGGGGAACCATCGGTACGGCGGATAATAAGATCGTCCAGTTCCTGGTTGCTGAACTCGATCGGGCCACGGATCTGATCGTc                                                                                                                                                                                                                                       >  1:85957/1‑250 (MQ=255)               aCCGTCATCGCCATTGATCATAGAAACGTGCGCGTAAACCGGCACCGGCGCTTTCAGGGCCTTAAGAATGTTGATCTGGCGTGGCGTGTTGTTGATATGGTCTTCGCCACGGATAACGTGGGTGATTTCCATATCCCAGTCATCGACAACCACACAGAAGTTATAGGTTGGGGAACCATCGGTACGGCGGATAATAAGATCGTCCAGTTCCTGGTTGCTGAACTCGATCGGACCACGGATCTGATCGTCa                                                                                                                                                                                                                                      <  1:340940/250‑1 (MQ=255)                   tCATCGCCATTGATCATAGAAACGTGCGCGTAAACCGGCACCGGCGCTTTCAGGGCCTTAAGAATGTTGATCTGGCGTGGCGTGTTGTTGATATGGTCTTCGCCACGGATAACGTGGGTGATTTCCATATCCCAGTCATCGACAACCACACAGAAGTTATAGGTTGGGGAACCATCGGTACGGCGGATAATAAGATCGTCCAGTTCCTGGTTGCTGAACTCGATCGGGCCACGGATCTGATCGTCaaaaa                                                                                                                                                                                                                                  >  1:226409/1‑250 (MQ=255)                              gATCATAGAAACGTGCGCGTAAACCGGCACCGGCGCTTTCAGGGCCTTAAGAATGTTGATCTGGCGTGGCGTGTTGTTGATATGGTCTTCGCCACGGATAACGTGGGTGATTTCCATATCCCAGTCATCGACAACCACACAGAAGTTATAGGTTGGGGAACCATCGGTACGGCGGATAATAAGATCGTCCAGTTCCTGGTTGCTGAACCCGATCGGGCCACGGATCTGATCGTCAAAAACAACAGAACCtt                                                                                                                                                                                                                      >  2:480044/1‑251 (MQ=255)                                 cATAGAAACGTGCGCGTAAACCGGCACCGGCGCTTTCAGGGCCTTAAGAATGTTGATCTGGCGTGGCGTGTTGTTGATATGGTCTTCGCCACGGATAACGTGGGTGATTTCCATATCCCAGTCATCGTCAACCACACAGAAGTTATAGGTTGGGGAACCATCGGTACGGCGGATAATAAGATCGTCCAGTTCCTGGTTGCTGAACTCGATCGGGCCACGGATCTGATCGTCAAAAAAAACAGAACCTTCCt                                                                                                                                                                                                                   >  1:456941/1‑251 (MQ=255)                                   tAGAAACGTGCGCGTAAACCGGCACCGGCGCTTTCAGGGCCTTAAGAATGTTGATCTGGCGTGGCGTGTTGTTGATATGGTCTTCGCCACGGATAACGTGGGTGATTTCCATATCCCAGTCATCGACAACCACACAGAAGTTATAGGTTGGGGAACCATCGGTACGGCGGATAATAAAATCGTCCAGTTCCTGGTTGCTGAACTCGATCGGGCCAAGGATCTGATCGTCCAAAAAAACAGAACCTTCCTGc                                                                                                                                                                                                                 >  2:56816/1‑251 (MQ=255)                                         cGTGCGCGTAAACCGGCACCGGCGCTTTCTATGCCTTAAAAATGTTGATCTGGCGTGGCGTGTTGTTGATATGGTCTTCGCCACGGATAACGTGGGTGATTTCCATATCCCAGTCATCGGCAACAACAGAGAAGTTATAGGCTGGGGAACAATCGGTACGTCGGAGAATAAGCGCGTCCAGTTCCTGGTAGCTGACCTCGATCGGTCCAAGGATCTGATCGTCAAAGACAACATAGCCTGACTGCGGGTTa                                                                                                                                                                                                           <  2:318250/251‑1 (MQ=255)                                                       ggcaccggcGCTTTCAGGGCATTAAGAATGTTGATCTGGCGTGGCGTGTTGTTGATATGGTCTTCGCCACGGATAACGTGGGTGATTTCCATATCCCAGTCATCGACAACCACACAGAAGTTATAGGTTGGGGAACCATCGGTACGGCGGATAATAAGATCGTCCAGTTCCTGGTTGCTGAACTCGATCGGGCCACGGATCTGATCGTCAAAAACAACAGAACCTTCCTGCGGGTTAGCAAAACGTacaa                                                                                                                                                                                              >  1:151841/1‑250 (MQ=255)                                                          accggcGCTTTCAGGGCCTTAAGAATGTTGATCTGGCGTGGCGTGTTGTTGATATGGTCTTCGCCACGGATAACGTGGGTGATTTCCATATCCCAGTCATCGACAACCACACAGAAGTTATAGGTTGGGGAACCATCGGTACGGCGGATAATAAGATCGTCCAGTTCCTGGTTGCTGAACTCGATCGGACCACGGATCTGATCGTCAAAAACAACAGAACCTTCCTGCGGGTTAGCAAAACGTacaaaac                                                                                                                                                                                           >  2:276666/1‑247 (MQ=255)                                                           caggcGCTTTCAGGGCCTTAAGAATGTTGATCTGGCGTGGCGTGTTGTTGATATGGTCTTCGCCACGGATAACGTGGGTGATTTCCATATCCCAGTCATCGACAACCACACAGAAGTTATAGGTTGGGGAACCATCGGTACGGCGGATAATAAGATCGTCCAGTTCCTGGTTGCTGAACTCGATCGGGCCACGGATCTGAGCGTCAAAAACAACAGAACCTTCCTGCGGGTTAGCAAAACGTACAacacac                                                                                                                                                                                         <  1:174855/249‑1 (MQ=255)                                                                  tttCAGGGCCTTAAGAATGTTGATCTGGCGTGGCGTGTTGTTGATATGGTCTTCGCCACGGATAACGTGGGTGATTTCCATATCCCAGTCATCGACAACCACACAGAAGTTATAGGTTGGGGAACCATCGGTACGGCGGATAATAAGATCGTCCAGTTCCTGGTTGCTGAACTCGATCGGACCACGGATCTGATCGTCCAAAACAACAGAACCTTCCTTCGGGTTAGCAAAACGTACAACAACCGGTtca                                                                                                                                                                                   >  2:86637/1‑250 (MQ=255)                                                                                  aTGTTGATCTGGCGTGGCGTGTTGTTGATATGGTCTTCGCCACGGATAACGTGGGTGATTTCCATATCCCAGTCATCGACAACCACACAGAAGTTATAGGTTGGGGAACCATCGGTACGGCGGATAATAAGATCGTCCAGTTCCTGGTTGCTGAACTCGATCGGGCCACGGATCTGATCGTCAAAAACAACAGAACCTTCCTGCGGGGTAGCAAAACGTACAACACACGGGTCATCATCAGCATGATGCTc                                                                                                                                                                  >  1:352011/1‑251 (MQ=255)                                                                                    gTTGATCTGGCGTGGCGTGTTGTTGATATGGTCTTCGCCACGGATAACGTGGGTGATTTCCATATCCCAGTCATCGACAACCACACAGAAGTTATAGGTTGGGGAACCATCGGTACGGCGGATAATAAGATCGTCCAGTTCCTGGTTGCTGAACTCGATCGGGCCACGGATCTGATCGTCAAAAACAACAGAACCTTCCTGCGGGTTAGCAAAACGTACAACACACGGTTCATCATCAGCATGATGCTCa                                                                                                                                                                 <  1:113330/250‑1 (MQ=255)                                                                                     ttGATCTGGCGTGGCGTGTTGTTGATATGGTCTTCGCCACGGATAACGTGGGTGATTTCCATATCCCAGTCATCGACAACCACACAGAAGTTATAGGTTGGGGAACCATCGGTACGGCGGATAATAAGATCGTCCAGTTCCTGGTTGCTGAACTCGATCGGGCCACGGATCTGATCGTCAAAAACAACAGAACCTTCCTGCGGGTTAGCAAAACGTACAACACACGGTTCATCATCAGCCTGCTGCTCATg                                                                                                                                                               >  2:177246/1‑251 (MQ=255)                                                                                     ttGATCTGGCGTGGCGTGTTGTTGATATGGTCCTCGCCACGGATAACGTGGGTGATTTCCATATCCCAGTCATCGACAACCACACAGAAGTTATAGGTTGGGGAACCATCGGTAAGGCGGATAATAAGATCGTCCAGTTCCTGGTTGCTGAACTCGATCGGACCACGGATCTGATCGTCAAAAACAACAAAACCTTCCTGCGGGTTACCAAAAAGCACAATACACGGCTCATAATCAGCATGATCCTCATg                                                                                                                                                               >  1:113708/1‑251 (MQ=255)                                                                                             gcgtggcgtgTTGTTGATATGGTCTTCGCCACGGATAACGTGGGTGATTTCCATATCCCAGTCATCGACAACCACACAGAAGTTATAGGTTGGGGAACCATCGGTACGGCGGATAATAAGATCGTCCAGTTCCTGGTTGCTGAACTCGATCGGGCCACGGATCTGATCGTCAAAAACAACAGAACCTTCCTGCGGGTTAGCAAAACGTACAACACACGGTTCATCATCAGCATGATGCTCATGGCTGTggc                                                                                                                                                       <  1:474573/251‑1 (MQ=255)                                                                                               gtggcgtgTTGTTGATATGGTCTTCGCCACGGATAACGTGGGTGATTTCCATATCCCAGTCATCGACAACCACACAGAAGTTATAGGTTGGGGAACCATCGGTACGGCGGATAATAAGATCGTCCAGTTCCTGGGTGCTGAACTCGATCGGGCCACGGATCTGATCGTCAAAAACAACAGAACCTTCCTGCGGGTTAGCAAAACGTACAACACACGGTTCATCATCAGCATGATGCTCATGGCTGTggcgg                                                                                                                                                     >  1:44685/1‑251 (MQ=255)                                                                                                 ggcgtgTTGTTGATATGGTCTTCGCCACGGATAACGTGGGTGATTTCCATATCCCAGTCATCGACAACCACACAGAAGTTATAGGTTGGGGAACCATCGGTACGGCGGATAATAAGATCGTCCAGTTCCTGGTTGCTGAACTCGATCGGGCCACGGATCTGATCGTCAAAAACAACAGAACCTTCCTGCGGGTTAGCAAAACGTACAACACACGGTTCATCATCAGCATGATGCTCATGGCTGTggcggc                                                                                                                                                    >  1:356804/1‑250 (MQ=255)                                                                                                       ttgttgATATGGTCTTCGCCACGGATAACGTGGGTGATTTCCATATCCCAGTCATCGACAACCACACAGAAGTTATAGGTTGGGGAACCATCGGTACGGCGGATAATAAGATCGTCCAGTTCCTGGTTGCTGAACTCGATCGGGCCACGGATCTGATCGTCAAAAACAACAGAACCTTCCTGCGGGTTAGCACAACGTACAACACACGGTTCATCATCAGCATGATGCTCATGGCTGTGGGGGCAGCGAcc                                                                                                                                             >  1:85385/1‑251 (MQ=255)                                                                                                        tgttgATATGGTCTTCGCCACGGATAACGTGGGTGATTTCCATATCCCAGTCATCGACAACCACACAGAAGTTATAGGTTGGGGAACCATCGGTACGGCGGATAATAAGATCGTCCAGTTCCTGGTTGCTGAACTCGATCGGGCCACGGATCTGATCGTCAAAAACAACAGAACCTTCCTGCGGGTTAGCAAAACGTACAACACACGGTTCATCATCAGCATGATGCTCATGGCTGTGGCGGCAGCGAcc                                                                                                                                             <  2:240803/250‑1 (MQ=255)                                                                                                               attgTCTTCGCCACGGATAACGTGGGTGATTTCCATATCCCAGTCATCGACAACCACACAGAAGTTATAGGTTGGGGAACCATCGGTACGGCGGATAATAAGATCGTCCAGTTCCTGGTTGCTGAACTCGATCGGACCACGGATCTGATCGTCAAAAACAACAGAACCTTCCTGCGGGTTAGCAAATCGTACAACACACGGTTCATCATCAGCATGATGCTCATGGCTGTGGCGGCAGCGACCGTCATaa                                                                                                                                      <  2:422212/247‑1 (MQ=255)                                                                                                                      tCGCCACGGATAACGTGGGTGATTTCCATATCCCAGTCATCGACAACCACACAGAAGTTATAGGTTGGGGAACCATCGGTACGGCGGATAATAAGATCGTCCAGTTCCTGGTTGCTGAACTCGATCGGGCCACGGATCTGATCGTCAAAAACAACAGAACCTTCCTGCGGGTTAGCAAAACGTACAACACACGGTTCATCATCAGCATGATGCTCATGGCTGTGGCGGCAGCGACCGTCATAACGCGGCt                                                                                                                               <  2:352011/250‑1 (MQ=255)                                                                                                                         ccACGGATAACGTGGGTGATTTCCATATCCCAGTCATCGACAACCACACAGAAGTTATAGGTTGGGGAACCATCGGTACGGCGGATAATAAGATCGTCCAGTTCCTGGTTGCTGAACTCGATCGGGCCACGGATCTGATCGTCAAAAACAACAGAACCTTCCTGCGGGTTAGCAAAACGTACAACACACGGTTCATCATCAGCATGATGCTCATGGCTGTGGCGGCAGCGACCGTCATAACGCGGCTtctc                                                                                                                           <  1:232886/251‑1 (MQ=255)                                                                                                                            cGGATAACGTGGGTGATTTCCATATCCCAGTCATCGACAACCACACAGAAGTTATAGGTTGGGGAACCATCGGTACGGCGGATAATAAGATCGTCCAGTTCCTGGTTGCTGAACTCGATCGGGCCACGGATCTGATCGTCAAAAACAACAGAACCTTCCTGCGGGTTAGCAAAACGTACAACACACGGTTCATCATCAGCATGATGCTCATGGCTGTGGCGGCAGCGACCGTCATAACGCGGCTTCTCAcc                                                                                                                        <  1:328795/251‑1 (MQ=255)                                                                                                                               aTAACGTGGGTGATTTCCATATCCCAGTCATCGACAACCACACAGAAGTTATAGGTTGGGGAACCATCGGTACGGCGGATAATAAGATCGTCCAGTTCCTGGTTGCTGAACTCGATCGGGCCACGGATCTGATCGTCAAAAACAACAGAACCTTCCTGCGGGTTAGCAAAACGTACAACACACGGTTCATCATCAGCATGATGCTCATGGCTGTGGCGGCAGCGACCGTCATAACGCGGCTTCTCACCttt                                                                                                                     <  2:262809/251‑1 (MQ=255)                                                                                                                                 aaCGTGGGTGATTTCCATATCCCAGTCATCGACAACCACACAGAAGTTATAGGTTGGGGAACCATCGGTACGGCGGATAATAAGATCGTCCAGTTCCTGGTTGCTGAACTCGATCGGGCCACGGATCTGATCGTCAAAAACAACAGAACCTTCCTGCGGGTTAGCAAAACGTACAACACACGGTTCATCATCAGCATGATGCTCATGGCTGTGGCGGCAGCGACCGTCATAACGCGGCTTCTCACCTTTc                                                                                                                    <  1:480044/250‑1 (MQ=255)                                                                                                                                       ggTGATTTCCATATCCCAGTCATCGACAACCACACAGAAGTTATAGGTTGGGGAACCAACGATACGGCGGATAATAAGATCGTGCAGTTCCTGGTTGCTGAACTCGATCGGGACACGGATCTGATCGTCAAAAACAAGAGAACATTAATACAAGTTAAAAAAACGTACACCACACGGTTCATCATCAGCATGATGCTCATGGCTGTGGCTGCAGCGACCGTCATAACGCGGCTTCTCACCTTTCGCCCttt                                                                                                             <  2:387451/251‑1 (MQ=255)                                                                                                                                             ttCCATATCCCAGTCATCGACAACCACACAGAAGTTATAGGTTGGGGAACCATCGGTACGGCGGATAATAAGATCGTCCAGTTCCTGGTTGCTGAACTCGATCGGGCCACGGATCTGATCGTCAAAAACAACAGAACCTTCCTGCGGGTTAGCAAAACGGACAACACACGGTTCATCATCAGCATGATGCTCATGGCTGTGGCGGCAGCGACCGTCATAACGCGGCTTCTCACCTTTCGACATTTGCTCtt                                                                                                       >  1:94666/1‑251 (MQ=255)                                                                                                                                             ttCCATATCCCAGTCATCGACAACCACACAGAAGTTATAGGTTGGGGAACCATCGGTACGGCGGATAATAAGATCGTCCAGTTCCTGGTTGCTGAACTCGATCGGACCACGGATCTGATCGTCAAAAACAACAGAACCTTCCTGCGGGTTAGCAAAACGTACAACACACGGTTCATCATCAGCATGATGCTCATGGCTGTGGCGGCCGCGCACGTCATAACGCGGCCTCTCCCCTTTCCCCCTTTGCNCtt                                                                                                       >  2:439724/1‑251 (MQ=255)                                                                                                                                                     cccAGTCATCGACAACCACACAGAAGTTATAGGTTGGGGAACCATCGGTACGGCGGATAATAAGATCGTCCAGTTCCTGGTTGCTGAACTCGATCGGGCCACGGATCTGATCGTCAAAAACAACAGAACCTTCCTGCGGGTTAGCAAAACGTACAACACACGGTTCATCATCAGCATGATGCTCATGGCTGTGGCGGCAGCGACCGTCATAACGCGGCTTCTCACCTTTCGCCATTTGCTCTTCGCGCAgc                                                                                               <  2:326582/251‑1 (MQ=255)                                                                                                                                                     cccAGTCATCGACAACCACACAGAAGTTATAGGTTGGGGAACCATCGGTACGGCGGATAATAAGATCGTCCAGTTCCTGGTTGCTGAACTCGATCGGGCCACGGATCTGATCGTCAAAAACAACAGAACCTTCCTGCGGGTTAGCAAAACGTACAACACACGGTTCATCATCAGCATGATGCTCATGGCTGTGGCGGCAGCGACCGTCATAACGCGCCTTCTCACCTTTCGCCATTTGCTCTTCGCGCAgc                                                                                               >  1:125388/1‑251 (MQ=255)                                                                                                                                                            aTCGACAACCACACAGAAGTTATAGGTTGGGGAACCATCGGTACGGCGGATAATAAGATCGTCCAGTTCCTGGTTGCTGAACTCGATCGGGCCACGGATCTGATCGTCAAAAACAACAGAACCTTCCTGCGGGTTAGCAAAACGTACAACACACGGTTCATCATCAGCATGATGCTCATGGCTGTGGCGGCAGCGACCGTCATAACGCGGCTTCTCACCTTTCGCCATTTGCTCTTCGCGCAGCGCTTCCa                                                                                        <  2:226409/251‑1 (MQ=255)                                                                                                                                                            aTCGACAACCACACAGAAGTTATAGGTTGGGGAACCATCGGTACGGCGGATAATAAGATCGTCCAGTTCCTGGTTGCTGAACTCGATCGGACCACGGATCTGATCGTCAAAAACAACAGAACCTTCCTGCGGGTTAGCAAAACGTACAACACACGGTTCATCATCAGCATGATGCTCATGGCTGTGGCGGCAGCGACCGTCATAACGCGGCTTCTCACCTTTCCCCCTTTGCTCTTCGCGCAGCGCTTccc                                                                                        >  2:118979/1‑250 (MQ=255)                                                                                                                                                                      acacaGAAGTTATAGGTTGGGGAACCATCGGTACGGCGGATAATAAGATCGTCCAGTTCCTGGTTGCTGAACTCGATCGGGCCACGGATCTGATCGTCAAAAACAACAGAACCTTCCTGCGGGTTAGCAAAACGTACAACACACGGTTCATCATCAGCATGATGCTCATGGCTGTGGCGGCAGCGACCGTCATAACGCGGCTTCTCACCTTTCGCCATTTGCTCTTCGCGCAGCGCTTCCAGGCGCTCttt                                                                              <  1:56816/251‑1 (MQ=255)                                                                                                                                                                           gAAGTTATAGGTTGGGGAACCATCGGTACGGCGGATAATAAGATCGTCCAGTTCCTGGTTGCTGAACTCGATCGGGCCACGGATCTGATCGTCAAAAACAACAGAACCTTCCTGCGGGTTAGCAAAACGTACAACACACGGTTCATCATCAGCATGATGCTCATGGCTGTGGCGGCAGCGACCGTCATAACGCGGCTTCTCACCTTTCGCCATTTGCTCTTCGCGCAGCGCTTCCAGGCGCTCTTTAGAGc                                                                         <  1:275381/251‑1 (MQ=255)                                                                                                                                                                            aaGTTATAGGTTGGGGAACCATCGGTACGGCGGATAATAAGATCGGCCAGTTCCTGGTTGCTGAACTCGATCGGGCCACGGATCTGATCGTCAAAAACAACAGAACCTTCCTGCGGGTTAGCAAAACGTACAACACACGGTTCATCATCAGCATGATGCTCATGGCTGGGGCGGCAGCGACCGTCATAACGCGGCTTCTCACCTTTCGCCATTTGCTCTTCGCGCAGCGCTTCCAGGCGCTCTTTAGAGCa                                                                        <  2:77157/251‑1 (MQ=255)                                                                                                                                                                                tataGGTTGGGGAACCATCGGTACGGCGGATAATAAGATCGTCCAGTTCCTGGTTGCTGAACTCGATCGGGCCACGGATCTGATCGTCAAAAACAACAGAACCTTCCTGCGGGTTAGCAAAACGTacaa                                                                                                                                                                                              <  2:53036/129‑1 (MQ=255)                                                                                                                                                                                tataGGTTGGGGAACCATCGGTACGGCGGATAATAAGATCGTCCAGTTCCTGGTTGCTGAACTCGATCGGGCCACGGATCTGATCGTCAAAAACAACAGAACCTTCCTGCGGGTTAGCAAAACGTacaa                                                                                                                                                                                              >  1:53036/1‑129 (MQ=255)                                                                                                                                                                                 ataGGTTGGGGAACCATCGGTACGGCGGATAATAAGATCGTCCAGTTCCTGGTTGCTTAACTCGATCGGACCCCGGATCTGATCGTCAAAAACAACAGAACCTTCCTGCGGGTTAGCAAAACGTACAACACACGGTTCATCATCAGCATGATGCTCATGGCTGTGGCGGCAGCGACCGTCATAACGCGGCTTCTCACCTTTCGCCATTTGCGCTTCGCGCAGCGCTTCCAGGCGCTCTTTAGAGCAATAg                                                                    <  2:471431/250‑1 (MQ=255)                                                                                                                                                                                  taGGTTGGGGAACCATCGGTACGGCGGATAATAAGATCGTCCAGTTCCTGGTTGCTGAACTCGATCGGGCCACGGATCTGATCGTCAAAAACAACAGACCCTTCCTGCGGGTTAGCAAAACGTACAACACACGGTtc                                                                                                                                                                                    >  1:86435/1‑137 (MQ=255)                                                                                                                                                                                  taGGTTGGGGAACCATCGGTACGGCGGATAATAAGATCGTCCAGTTCCTGGTTGCTGAACTCGATCGGGCCACGGATCTGATCGTCAAAAACAACAGAACCTTCCTGCGGGTTAGCAAAACGTACAACACACGGTtc                                                                                                                                                                                    <  2:86435/137‑1 (MQ=255)                                                                                                                                                                                     gTTGGGGAACCATCGGTACGGCGGATAATAAGATCGTCCAGTTCCTGGTTGCTGAACTCGATCGGGCCACGGATCTGATCGTCAAAAACAACAGAACCTTCCTGCGGGTTAGCAAAACGTACAACACACGGTTCATCATCAGCATGATGCTCATGGCTGTTGCGGCAGCGACCGTCATAACGCGGCTTCTCACCTTTCGCCATTTGCTCTTCGCGCAGCGCTTCCAGGCGCTCTTTAGAGCAATAGCAttt                                                               >  1:281913/1‑251 (MQ=255)                                                                                                                                                                                             aCCATCGGTACGGCGGATAATAAGATCGTCCAGTTCCTGGTTGCTGAACTCGATCGGGCCACGGATCTGATCGTCAAAAACAACAGAACCTTCCTGCGGGTTAGCAAAACGTACAACACACGGTTCATCATCAGCATGATGCTCATGGCTGTGGCGGCAGCGACCGTCATAACGCGGCTTCTCACCTTTCCCCATTTGCTCTTCGCGCAGCGCTTCCAGGCGCTCTTTAGAGCAATAGCATTTATTAGCa                                                        >  2:273643/1‑250 (MQ=255)                                                                                                                                                                                                       cggcggATAATAAGATCGTCCAGTTCCTGGTTGCTGAACTCGATCGGGCCACGGATCTGATCGTCAAAAACAACAGAACCTTCCTGCGGGTTAGCAAAACGTACAACACACGGTTCATCATCAGCATGATGCTCATGGCTGTGGCGGCAGCGACCGTCATAACGCGGCTTCTCACCTTTCGCCATTTGCTCTTCGCGCAGCGCTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTc                                             <  2:85957/251‑1 (MQ=255)                                                                                                                                                                                                          cggATAATAAGATCGTCCAGTTCCTGGTTGCTGAACTCGATCGGGCCACGGATCTGATCGTCAAAAACAACAGAACCTTCCTGCGGGTTAGCAAAACGTACAACACACGGTTCATCATCAGCATGATGCTCATGGCTGTGGCGGCAGCGACCGTCATAACGCGGCTTCTCACCTTTCGCCATTTGCTCTTCGCGCAGCGCTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTcc                                            >  1:235518/1‑249 (MQ=255)                                                                                                                                                                                                               aataaGATCGTCCAGTTCCTGGTTGCTGAACTCGATCGGGCCACGGATCTGATCGTCAAAAACAACAGAACCTTCCTGCGGGTTAGCAAAACGTACAACACACGGTTCATCATCAGCATGATGCTCATGGCTGTGGCGGCAGCGACCGTCATAACGCGGCTTCTCACCTTTCGCCATTTGCTCTTCGCGCAGCGCTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCt                                     >  2:459731/1‑251 (MQ=255)                                                                                                                                                                                                                 taaGATCGTCCAGTTCCTGGTTGCTGAACTCGATCGGGCCACGGATCTGATCGTCAAAAACAACAGAACCTTCCTGAGGGTTAGCAAAACGTACAACACACGGTTCATCATCACCATGATGCTCATGGCTGTGGCGGCAGCGACCGTCATAACGCGGCTTCTCACCTTTCGCCATTTGCTCTTCGCGCAGCGCTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTg                                    >  1:488586/1‑250 (MQ=255)                                                                                                                                                                                                                   aGATCGTCCAGTTCCTGGTTGCTGAACTCGATCGGGCCACGGATCTGATCGTCAAAAACAACAGAACCTTCCTGCGGGTTAGCAAAACGTACAACACACGGTTCATCATCAGCATGATGCTCATGGCTGTTGCGGCAGCGACCGTCATAACGCGGCTTCTCACCTTTCGCCATTTGCTCTTCGCGCAGCGCTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTgatc                                 <  2:281913/251‑1 (MQ=255)                                                                                                                                                                                                                   aGATCGTCCAGTTCCTGGTTGCTGAACTCGATCGGGCCACGGATCTGATCGTCAAAAACAACAGAACCTTCCTGCGGGTTAGCAAAACGTACAACACACGGTTCATCATCAGCATGATGCTCATGGCTGTGGCGGCAGCGACCGTCATAACGCGGCTTCTCACCTTTCGCCATTTGCTCTTCGCGCAGCGCTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTgatc                                 <  2:28101/251‑1 (MQ=255)                                                                                                                                                                                                                      tCGTCCAGTTCCTGGTTGCTGAACTCGATCGGGCCACGGATCTGATCGTCAAAAACAACAGAACCTTCCTGCGGGTTAGCAAAACGTACAACACACGGTTCATCATCAGCATGATGCTCATGGCTGTGGCGGCAGCGACCGTCATAACGCGGCTTCTCACCTTTCGCCATTTGCTCTTCGCGCAGCGCTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTgatcgat                              >  1:261945/1‑251 (MQ=255)                                                                                                                                                                                                                       cGTCCAGTTCCTGGTTGCTGAACTCGATCGGGCCACGGATCTGATCGTCAAAAACAACAGAACCTTCCTGCGGGTTAGCAAAACGTACAACACACGGTTCATCATCAGCATGATGCTCATGGCTGTGGCGGCAGCGACCGTCATAACGCGGCTTCTCACCTTTCGCCATTTGCTCTTCGCGCAGCGCTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCCTCCAGCATCTgatcgatc                             >  2:417778/1‑251 (MQ=255)                                                                                                                                                                                                                        gTCCAGTTCCTGGTTGCTGAACTCGATCGGACCACGGATCTGATCGTCAAAAACAACAGAACCTTCCTGCGGGTTAGCAAAACGTACAACACACGGTTCATCATCAGCATGATGCTCATGGCTGTGGCGGCAGCGACCGTCATAACGCGGCTTCTCACCTTTCGCCCTTTGCTCTTCGCCCAGCGCTTCCAGGCGCTCTTTATAGCAATAGCATTTATAAGCAGTGCCCTCTTCCaccctctgttccatcc                            >  2:264924/1‑236 (MQ=255)                                                                                                                                                                                                                                         tGAACTCGATCGGGCCACGGATCTGATCGTCAAAAACAACAGAACCTTCCTGCGGGTTAGCAAAACGTACAACACACGGTTCATCATCAGCATGATGCTCATGGCTGTGGCGGCAGCGACCGTCATAACGCGGCTTCTCACCTTTCGCCATTTGCTCTTCGCGCAGCGCTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATc            >  1:115036/1‑250 (MQ=255)                                                                                                                                                                                                                                               cGATCGGGCCACGGATCTGATCGTCAAAAACAACAGAACCTTCCTGCGGGTTAGCAAAACGTACAACACACGGTTCATCATCAGCATGATGCTCATGGCTGTGGCGGCAGCGACCGTCATAACGCGGCTTCTCACCTTTCGCCATTTGCTCTTCGCGCAGCGCTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGGTGTAGCGATCAAAACGt     >  1:23928/1‑251 (MQ=255)                                                                                                                                                                                                                                               cGATCGGGCCACGGATCTGATCGTCAAAAACAACAGAACCTTCCTGCGGGTTAGCAAAACGTACAACACACGGTTCATCATCAGCATGATGCTCATGGCTGTGGCGGCAGCGACCGTCATAACGCGGCTTCTCACCTTTCGCCATTTGCTCTTCGCGCAGCGCTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCGTCCAGCATCTGAGCGATCACCGCGTTGTAGCGAGCAAAAc       <  2:44685/249‑1 (MQ=255)                                                                                                                                                                                                                                                gATCGGACCACGAATCTGATCGTCAAAAACAACAGAACCGTCCTGCTGGGTAGCAAAACGTACAACACACAGTTCATCATCAGCATGATGCTCATGGCTGTGGCGGCAGCGACCGGCATAACGCGGCTTCTCATCTTTCGCCATTTGCTCTTCGCGCAGCGCTTCCAGGCGCTCTTTAGAGCAGTAGCATTTGTAAGCCGGGCCCTCTTCCAGCATCTTCTCGATCACCGCGTTGTAGCTATCAAAActtt    >  2:177821/1‑248 (MQ=255)                                                                                                                                                                                                                                                   cGGACCACGGATCTGATCGTCAAAAACAACAGAACCTTCCTGCGGGTTAGCAAAACGTACAACACACGGTTCATCATCAGCATGATGCTCATGGCTGTGGCGGCAGCGACCGTCATAACGCGGCTTCTCACCTTTCGCCATTTGCTCTTCGCGCAGCGCTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTg  <  1:15529/250‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       TTGGACAGTTTTTTACCGTCATCGCCATTGATCATAGAAACGTGCGCGTAAACCGGCACCGGCGCTTTCAGGGCCTTAAGAATGTTGATCTGGCGTGGCGTGTTGTTGATATGGTCTTCGCCACGGATAACGTGGGTGATTTCCATATCCCAGTCATCGACAACCACACAGAAGTTATAGGTTGGGGAACCATCGGTACGGCGGATAATAAGATCGTCCAGTTCCTGGTTGCTGAACTCGATCGGACCACGGATCTGATCGTCAAAAACAACAGAACCTTCCTGCGGGTTAGCAAAACGTACAACACACGGTTCATCATCAGCATGATGCTCATGGCTGTGGCGGCAGCGACCGTCATAACGCGGCTTCTCACCTTTCGCCATTTGCTCTTCGCGCAGCGCTTCCAGGCGCTCTTTAGAGCAATAGCATTTATAAGCAGTGCCCTCTTCCAGCATCTGATCGATCACCGCGTTGTAGCGATCAAAACGTTTG  >  NC_000913/2519954‑2520445

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