Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	651,267	0	A	T	23.4%	147.3 / 12.4	64	L197* (TTA→TAA)	citC	citrate lyase synthetase
Reads supporting (aligned to +/- strand):  ref base A (5/44);  new base T (15/0);  total (20/44)
Fisher's exact test for biased strand distribution p-value = 9.72e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CCAGCCAGTAGCGCATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGT‑GCC‑TTTTAACACCAAAT‑CC‑A‑GCCGGTCTTC‑ATAGGGGAAGCGTGAAGAATCTTCTTT‑GACTAAAAACAGATG‑C‑AACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAA‑GGG‑ATTGGCGTTCATCACAATGCAGCCAAT‑CTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTG  >  NC_000913/651019‑651515                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                   ccagccagTAGCGCATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGT‑GCC‑TTTTaa                                                                                                                                                                                                                                                                   <  1:310860/250‑1 (MQ=255)   agccagTAGCGCATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGT‑GCC‑TTTTAacac                                                                                                                                                                                                                                                                <  2:223297/251‑1 (MQ=255)            cgcATATCCTGGTTGTACTGGGCAGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGTCGGT‑GCC‑TTTTAACACCAAAT‑CC‑A‑g                                                                                                                                                                                                                                                    <  1:374680/251‑1 (MQ=255)              cATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATAGGCGGT‑GCC‑TTTTTTCACCAAAA‑CC‑A‑Gcc                                                                                                                                                                                                                                                  >  2:433120/1‑251 (MQ=255)                                      aaCGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGT‑GCC‑TTTTAACACCAAAT‑CC‑A‑GCCGGTCTTC‑ATAGGGGAAGCGTgaa                                                                                                                                                                                                                          <  1:88830/250‑1 (MQ=255)                                      aaCGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGT‑GCC‑TTTTAACACCAAAT‑CC‑A‑GCCGGTCTTC‑ATAGCGGAAGCGTgaag                                                                                                                                                                                                                         <  2:448506/251‑1 (MQ=255)                                         gcgACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGT‑GCC‑TTTTAACACCAAAT‑CC‑A‑GCCGGTCTTC‑ATAGGGGAAGCg                                                                                                                                                                                                                              <  1:334648/243‑1 (MQ=255)                                         gcgACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGG‑GGC‑TTTTTACACCAAAT‑CC‑C‑GCCGGTCTTC‑ATAgggggagggggatcggg                                                                                                                                                                                                                      >  2:334648/1‑238 (MQ=255)                                                               aCAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGT‑GGC‑TTTTTACACCAAAT‑AC‑A‑ACCGGTCTTC‑CTAGGGGGAGCGGGAGGAATACTCTTT‑TGCTAAAAACaact                                                                                                                                                                                                >  2:18008/1‑247 (MQ=255)                                                                  aaGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGT‑GCC‑TTTTAACACCAAAT‑CC‑A‑GCCGGTCTTC‑ATAGGGGAAGCGTGAAGAATCTTCTTT‑GACTAAAAACAGATG‑C‑aa                                                                                                                                                                                          <  1:3722/251‑1 (MQ=255)                                                                   aGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGT‑GCC‑TTTTAACACCAAAT‑CC‑A‑GCCGGTCTTC‑ATAGGGGAAGCGTGAAGAATCTTCTTT‑GACTAAAAACAGATG‑C‑AAc                                                                                                                                                                                         <  1:136022/251‑1 (MQ=255)                                                                      ggTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGT‑GCC‑TTTTAACACCAAAT‑CC‑A‑GCCGGTCTTC‑ATAGGGGAAGCGTGAAGAATCTTCTTT‑GACTAAAAACAGATG‑C‑AACCAg                                                                                                                                                                                      <  2:378598/251‑1 (MQ=255)                                                                                         gggAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGT‑GCC‑TTTTAACACCAAAT‑CC‑A‑GCCGGTCTTC‑ATAGGGGAAGCGTGAAGAATCTTCTTT‑GACTAAAAACAGATG‑C‑AACCAGTCGCACTGTGCCGCAGCCt                                                                                                                                                                   <  1:248474/251‑1 (MQ=255)                                                                                          ggAGCGAGGTTCTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGGTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGT‑GCC‑TTTTAACACCAAAT‑CC‑A‑GCCGGTCTTC‑ATAGGGGAAGCGTGAAGAATCTTCTTT‑GACTAAAAACAGATG‑C‑AACCAGTCGCACTGTGCCGCAGCCTg                                                                                                                                                                  <  2:257296/251‑1 (MQ=255)                                                                                                    aCTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGG‑GGC‑TTTTTACACCAAAT‑CC‑C‑GCCGGTCCTC‑ATAGGGGAAGCGTGAAGAATCTTATTT‑GACTCAAAAAAGATG‑A‑AAACAGGCTGACTCTGCCGCATCCTTTTGAATc                                                                                                                                                           >  2:39225/1‑248 (MQ=255)                                                                                                       gACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGT‑GCC‑TTTTAACACCAAAT‑CC‑A‑GCCGGTCTTC‑ATAGGGGAAGCGTGAAGAATCTTCTTT‑GACTAAAAACAGATG‑C‑AACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGAt                                                                                                                                                       <  2:274759/249‑1 (MQ=255)                                                                                                          ggAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGT‑GCC‑TTTTAACACCAAAT‑CC‑A‑GCCGGTCTTC‑ATAGGGGAAGCGTGAAGAATCTTCTTT‑GACTAAAAACAGATG‑C‑AACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACgg                                                                                                                                                  <  2:323490/251‑1 (MQ=255)                                                                                                          ggAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGC‑GCC‑TTTTTACACCAAAT‑CC‑ACGCGGGT‑TTC‑ATAAGGGAAGGGCGATGAAAAATCTTT‑TCTTAAAAACAGATG‑C‑AACCAGGCGCACTGTGCCGCAGCCTGTTGAATCCGAAAACgg                                                                                                                                                  >  1:416500/1‑251 (MQ=255)                                                                                                           gAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTTGCT‑TTTTAAAACCAA‑T‑CC‑C‑CCCGGTCTTC‑CTAGGGGGAGCGTGAAGAATCTTCTTT‑TGCTAAAAACAGATG‑C‑AACCAGTCGGACAGGGGCGCAGCCTGTTGTATCAGATAAcagg                                                                                                                                                 >  2:357884/1‑248 (MQ=255)                                                                                                           gAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGT‑GCC‑TTTTTACACCAAAT‑CC‑A‑GCCCGGCTTC‑ATAGGGGGAGCGTGAAGAATCTTCTTT‑GACTAAAAACAGATG‑C‑AACCAGTCGCACTGTGCCGCAGCCCGTTGGATCAGATAACggg                                                                                                                                                 >  1:165515/1‑250 (MQ=255)                                                                                                                  ttCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGA‑GCC‑TTTTTACACCCAATACC‑G‑ACCGG‑CTTT‑AGAGGGGGAGCGTGGAGAATCTTTTTT‑TGCTAAAAACAAATGGC‑ACCCAG‑CGCACTTTGCCGCCGCCCTTTGAATCAGATAACGGGGAccctt                                                                                                                                          >  1:222486/1‑248 (MQ=255)                                                                                                                    cAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGT‑GCC‑TTTTAACACCAAAT‑CC‑A‑GCCGGTCTTC‑ATAGGGGAAGCGTGAAGAATCTTCTTT‑GACTAAAAACAGATG‑C‑AACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCg                                                                                                                                        <  1:26099/251‑1 (MQ=255)                                                                                                                    cAGATCAATTTCGGTGTAACAAGGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGAGGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGT‑GCC‑TTTTAACACCAAAT‑CC‑A‑GCCGGTCTTC‑ATAGGGGAAGCGTGAAGAATCTTCTTT‑GACTAAAAACAGAGG‑C‑AACCATTCGCACTGTGCCGCAGCCTGTTGGATCAGAGAACGGTGACCATTc                                                                                                                                         <  2:70865/250‑1 (MQ=255)                                                                                                                    cAGATCAATTTCGGTGGAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGT‑GCC‑TTTTAACACCAAAT‑CC‑A‑GCCGGTCTTC‑ATAGGGGAAGCGTGAAGAATCTTCTTT‑GACTAAAAACAGATG‑C‑AACCAGTCGCGCTGTGCCGCAGCCTGGTGAATCAGATAACGGTGACCATTc                                                                                                                                         <  1:151522/250‑1 (MQ=255)                                                                                                                     aGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGT‑GGC‑TTTTAACACCAAAT‑CC‑A‑GCCGGTCTTC‑ATAGGGGGAGCGTGAAGAATCTTCTTT‑GACTAAAAACAGATG‑C‑AACCAGTCGCACTGTGCCGCCGCCTGTTGAATCAGATAACGGGGACCATTCGt                                                                                                                                       >  1:248367/1‑251 (MQ=255)                                                                                                                         cAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGT‑GCC‑TTTTTACACCAAAA‑CC‑A‑GCCGGTCTTC‑ATAGGGGGAGCGTGAAGAATCTTCTTT‑TAATAAAAAAAGATG‑C‑AACCAATCGCACTGTGCCGCCGCCTGTTTAATCAGATatgtgggtttgttcactttt                                                                                                                                   >  1:9743/1‑232 (MQ=255)                                                                                                                           aTTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGT‑GCC‑TTTTTAAACCAAAT‑CC‑A‑GCCGGGCTTT‑CTAGGGGGAGCGTGTAGAATTTTCTTT‑TACTCAAAACAAATG‑G‑AACCAGTCTCCCTTTGCCGCAGCCTGTTGGATCAGATAACGGGGACCATTCGTTAA‑ggg                                                                                                                                >  1:449526/1‑251 (MQ=255)                                                                                                                             ttCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGT‑GCC‑TTTTTACACCAAAT‑CC‑A‑GCCGGGCTTC‑ATTGGGGGAGAGGGGAGAATCTTCTTTTGA‑TAAAAACAAAAG‑G‑AACCAGGCGCCCTGTGCCGCCGCCCGGTGAATCAGATAACGGTGACCCTTCGGAAA‑ggggt                                                                                                                              >  1:152039/1‑249 (MQ=255)                                                                                                                               cGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGT‑GCC‑TTTTAACACCAAAT‑CC‑A‑GCCGGTCTTC‑ATAGGGGAAGCGTGAAGAATCTTCTTT‑GACTAAAAACAGATG‑C‑AACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAA‑GGG‑Att                                                                                                                            <  1:18008/250‑1 (MQ=255)                                                                                                                                         aTGGTTAATGACGCTCTGTTCGTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGT‑GCC‑TTTTAACACCAAAT‑CC‑A‑GCCGGTCTTC‑ATAGGGGAAGCGTGAAGAATCTTCTTT‑GACTAAAAACAGATG‑C‑AACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAA‑GGG‑ATTGGCGTtcatca                                                                                                                 <  2:20398/251‑1 (MQ=255)                                                                                                                                         aTGGTTAATGACGCTCGGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGT‑GCC‑TTTTAACACCAAAT‑CC‑A‑GCCGGTCTTC‑ATAGGGGAAGCGTGAAGAATCTTCTTT‑GACTAAAAACAGATG‑C‑AACCAGTCGCACTGTGCCGCAGCCTGGGGAATCATATAACGGTGACCATTCGTAAA‑GGG‑ATTGGCGTtcatca                                                                                                                 <  1:196533/251‑1 (MQ=255)                                                                                                                                                aTGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGT‑GCC‑TTTTTAAACCAAAT‑TC‑A‑GCCGGGCTTCCTTAGGGGA‑GCGGGGAGAATCTTTTTT‑GGCTAAAAAAAGAAG‑C‑AACCAGTCGCCCTGTGGCGCAGCCCGTTGAATCCGATAAAGGGGGCCCATCGTAAA‑AGG‑GTTGGGGTTCATCACAATGc                                                                                                           >  1:111916/1‑250 (MQ=255)                                                                                                                                                 tGACGCTCTGTTATTTAATGAAGTAGCACGGGAACGCAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGT‑GTC‑TTTTACAACCAAAT‑CC‑C‑GCCGGGCCTC‑ATAGGGGGAGCGGGAAAAATCTTCTTT‑TGCTAAAAACAGATG‑C‑AAACAGTCGCCCTGCGCCGCAGACTGTTGAATCCGATAACGGGGAGCCTTCGTAAA‑GGG‑TGTGGGGTTCATcaccatcccg                                                                                                         >  1:379817/1‑244 (MQ=255)                                                                                                                                                        ctGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGT‑GCC‑TTTTAACACCAAAT‑CC‑A‑GCCGGTCTGC‑ATAGGGGAAGCGTGAAGAATCTTCTTT‑GACTAAAAACAGATG‑C‑AACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAA‑GGG‑ATTGGCGTTCATCACAATGCAGCCAAt                                                                                                    <  1:314960/249‑1 (MQ=255)                                                                                                                                                                tAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGT‑GCC‑TTTTAACACCAAAT‑CC‑A‑GCCGGTCTTC‑ATAGGGGAAGCGTGAAGAATCTTCTTT‑GACTAAAAACAGATG‑C‑AACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAA‑GGG‑ATTGGCGTTCATCACAATGCAGCCAAT‑CTTGTTccc                                                                                          <  2:248367/250‑1 (MQ=255)                                                                                                                                                                  aTGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGT‑GCC‑TTTTTACACCAAAT‑CC‑A‑GCCGGGCTTC‑TTAGGGGGAGCGGGAAGAATCTTCTTT‑GACTAAAAACAGATG‑C‑AACCAGCCGCACTGTGCCGCCGCCTGTTGAATCAGAt                                                                                                                                                       >  1:151978/1‑190 (MQ=255)                                                                                                                                                                  aTGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGT‑GCC‑TTTTAACACCAAAT‑CC‑A‑GCCGGTCTTC‑ATAGGGGAAGCGTGAAGAATCTTCTTT‑GACTAAAAACAGATG‑C‑AACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGAt                                                                                                                                                       <  2:151978/190‑1 (MQ=255)                                                                                                                                                                  aTGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGT‑GCC‑TTTTAACACCAAAG‑CC‑A‑GCCgg                                                                                                                                                                                                                                                <  1:380598/105‑1 (MQ=255)                                                                                                                                                                    gAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGT‑GCC‑TTTTAACACCAAAT‑CC‑A‑GCCGGTCTTC‑ATAGGGGAAGCGTGAAGAATCTTCTTT‑GACTAAAAACAGATG‑C‑AACCAGTCGCACTGTGCCGCAGCCGGTTGAATCAGATAACGGTGACCATTCGTAAA‑GGG‑ATTGGCg                                                                                                                        <  2:231378/217‑1 (MQ=255)                                                                                                                                                                    gAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGA‑GCC‑TTTTTACACCAAAT‑CC‑A‑GCCGGTCTTC‑ATAGGGGAAGCGTGAAGAATCTTCTTT‑GACTAAAAACAGATG‑C‑AAACCGTCGCCCTGTGCCGCAGCCTGTTGAATCAGATAACGGGGACCATTCGTAAA‑GGG‑ATTGGCg                                                                                                                        >  1:231378/1‑217 (MQ=255)                                                                                                                                                                     aaGTAGCAAGGGAAGGTAGCGCGGGAGAGGATGTATTCGGAGCCACGAGGCACAGTCAGGCGTGAAATATCGGCGGT‑GCC‑TTTTAACACCAAAT‑CC‑A‑GCCTGTCTTC‑ATAGGGGAAGCGTGAAGAATCTTCTTT‑GACTAAAAACAGATG‑C‑AACCAGTCGCACTGTGCCGCAGCCTTTTGAATCAGATAACGGTGACCATTCGTAAA‑GGG‑ATTGGCGTTCATCACAATGCACCCAAT‑CTTGTTCCCTGGATg                                                                                    <  2:289800/251‑1 (MQ=255)                                                                                                                                                                         aGCAAGGGAACGTAGCGCGGGAGATGAGTTATGCGGAGCCACGATGCGTAGTCAGGCGTGGAATATCGGCGGT‑GCC‑TTTTAACGCCAAAT‑CC‑A‑GCCGGGCTTC‑ATAGGGGAAGCGTGAAGAATCTTCTTT‑GACTAAAAACAGATT‑C‑AACCAGTCGTACTGTGCCGCAGCCTGGTGAATCAGATAACGGTGACCATTCGTAAA‑GTG‑ATGGGCTGTCATCACAATGCAGCCAAT‑CTTGTTCCCTGGATGACg                                                                                 <  1:317238/250‑1 (MQ=255)                                                                                                                                                                                      aGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGT‑GCC‑CTTTTATACCAAAT‑CC‑C‑CCCGGTCTTT‑CTAGGGGGAGCGGGAAGAATCTTTTTT‑CACTAAAAACACATG‑CCAACCA‑TTGCACCTTGTCGCCGCCCGTTGAATCCGATAACGGGGACCATTCGTTAA‑GGG‑GTTGGGGTTCATTCCCATGCAGCCAAA‑CTTGTTTCCTCGATGACGAAATTTTTTCTGCg                                                                   >  1:271867/1‑251 (MQ=255)                                                                                                                                                                                           gggAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGG‑GCCCTTTTAA‑AACAAAT‑CCCA‑CCCGG‑CTTT‑ATTGGGGGAGAGGGAAGAAACTTCTTT‑GACTAAAAAAAAATG‑G‑AACCAGACGCCCTGTGCCCCAGCCTGTTTAAACAGATAACGGTGACCCATCGTAAA‑GGG‑ATTGGCGTTCATTACAAATCAGCCACT‑CTTGTTCCCTGGAGGAAGAAATTTTTTCATCGATTcc                                                              >  2:393128/1‑250 (MQ=255)                                                                                                                                                                                                     gTATTCGGAGCCACGATGCCCAGTCAGGCGTGGAATATCGGCGGT‑GCC‑TTTTAACACCAAAT‑CC‑A‑GCCGGTCTTC‑ATAGGGGAAGCGTGAAGAATCTTCTTT‑GACTAAAAACAGATG‑C‑AACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAA‑GGG‑ATTGGCGTTCATCACAATGCAGCCAAT‑CTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGtt                                                    <  2:167541/251‑1 (MQ=255)                                                                                                                                                                                                     gTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGT‑GCC‑TTTTAACACCAAAT‑CC‑A‑GCCGGTCTTC‑ATAGGGGAAGCGTGAAGAATCTTCTTT‑GACTAAAAACAGATG‑C‑AACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAA‑GGG‑ATTGGCGTTCATCACAATGCAGCCAAT‑CTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGt                                                     <  1:456516/250‑1 (MQ=255)                                                                                                                                                                                                         tCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGT‑GCC‑TTTTAACACCAAAT‑CC‑A‑GCCGGTCTTC‑ATAGGGGAAGCGTGAAGAATCTTCTTT‑GACTAAAAACAGATG‑C‑AACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAA‑GGG‑ATTGGCGTTCATCACAATGCAGCCAAT‑CTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCa                                                 <  1:413388/250‑1 (MQ=255)                                                                                                                                                                                                          cGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGT‑GCC‑TTTTAACACCAAAT‑CC‑A‑GCCGGTCTTC‑ATAGGGGAATCGTGAAGAATCTTCTTT‑GACTAAAAACAGATG‑C‑AACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAA‑GGT‑AGTGGCGTTCATCACAATGCAGCCAAT‑CTTGTACCCGGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAg                                                <  1:438032/250‑1 (MQ=255)                                                                                                                                                                                                          cGGAGCCACGAGGAACAGCAAGGCGTGGAATATCGGCGGT‑GCC‑TTTTAACACCAAAT‑CC‑A‑GCCGGTCTTC‑ATAGGGGAAGCGTGAAGAATCTGCTTT‑GACTAAAAACAGATG‑C‑AACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAA‑GGG‑ATTGGCGTTCATCACAATGCAGCCAAT‑CTGGTTCCCTTGATGACGAAATGTTTGCATCGAtgctggatatcttttcag                                                <  2:431309/250‑18 (MQ=255)                                                                                                                                                                                                           ggAGCCACGATGCACAGTCAGGCGTTGAATAGCGGGGGT‑GCC‑TTTTAACGCCAGAG‑CC‑A‑GCCGGTCTTC‑ATAGGGGAGGCGTGAAGAATCTTCTTG‑GGCTAAAGACAGATG‑C‑AACCAGTCTCAGTGTGCGGCAGCCTGTTGAAGCAGATAACGGTGACCATTCGGAAT‑TGG‑GTTGGCGTTCATCAGAATGCAGCCAAT‑CTTGTCCCCTGTATTACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTc                                              <  1:410207/251‑1 (MQ=255)                                                                                                                                                                                                            gAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGT‑GCC‑TTTTAACACCAAAT‑CC‑A‑GCCGGTCTTC‑ATAGGGGAAGCGTGAAGAATCTTCTTT‑GACTAAAAACAGATG‑C‑AACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAA‑GGG‑ATTGGCGTTCATCACAATGAAGCCAAT‑CTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTcg                                             <  2:165515/251‑1 (MQ=255)                                                                                                                                                                                                              gCCACGATGCACAGTCAGGCGTGGAATATCGGCGGT‑GCC‑TTTTAACACCAAAT‑CC‑A‑GCCGGTCTTC‑ATAGGGGAAGCGTGAAGAATCTTCTTT‑GACTAAAAACAGATG‑C‑AACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAA‑GGG‑ATTGGCGTTCATCACAATGCAGCCAAT‑CTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTcgcg                                           <  1:45127/251‑1 (MQ=255)                                                                                                                                                                                                               ccACGATGCACAGTCAGGCGTGGAATATCGGCGGT‑GCC‑TTTTAACACCAAAT‑CC‑A‑GCCGGTCTTC‑ATAGGGGAAGCGTGAAGAATCTTCTTT‑GACTAAAAACAGATG‑C‑AACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAA‑GGG‑ATTGGCGTTCATCACAATGCAGCCAAT‑CTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGt                                          <  1:330720/251‑1 (MQ=255)                                                                                                                                                                                                                cACGATGCACAGTCAGGCGTGGAATATCGGCGGT‑GCC‑TTTTAACACCAAAT‑CC‑A‑GGAGGGCTTA‑ATTGGGGGAGACTGTAGAAACTTCTTT‑GTCTAAAAAAAAATG‑C‑AACCAGGAGCACTGGGCCGCCGGCTGTTGAATTAGATAACGGTGACAATTTGTAAAAGGG‑AT‑GGCGGTTATCAAAATGAAGCCAATATTTTTTCCCGGGAT‑AAGAAAATTTTTTAGCGGTTCGGCCTAGCGTTTCAgatcggc                                          >  1:188905/1‑244 (MQ=255)                                                                                                                                                                                                                cACGATGCACAGTCAGGCGTGGAATATCGGCGGT‑GCC‑TTTTAACACCAAAT‑CC‑A‑GCCGGTCTTC‑ATAGGGGAAGCGTGAAGAATCTTCTTT‑GACTAAAAACAGATG‑C‑AACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAA‑GGG‑ATGGGCGTTCATCACAATGCAGCCAAT‑CTTGTTCCCTGGATGACGAAATGTTTTCAGCGATTCGGCATAGCGTTTc                                                  <  2:188905/242‑1 (MQ=255)                                                                                                                                                                                                                  cGATGCACAGTCAGGCGTGGAATATCGGCGGT‑GCC‑TTTTAACACCAAAT‑CC‑A‑GCCGGTCTTC‑ATAGGGGAAGCGTGAAGAATCTTCTTT‑GACTAAAAACAGATG‑C‑AACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAA‑GGG‑ATTGGCGTTCATCACAATGCAGCCAAT‑CTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGgc                                       <  1:15573/251‑1 (MQ=255)                                                                                                                                                                                                                           gTCAGGCGTGGAATATCGGCGGC‑GCC‑TTTTTACACCAAAA‑CC‑C‑GCCCGGCTTC‑ATAGGGGGAGCGTGGAGAAACTTCTTT‑TACTAAAAACAGAGG‑C‑AACCAGTCGCCCTGTGCCCCAGCCTGTTGAATCAGATAACGGTGACCATTCGGAAA‑GGGGTTTGGGGTT‑ATCACCATGCAGCCCAA‑CTTGTTCCCTGGGTGACGAAATTTTTTTAGCGATTTGGCATATGGTTTCAGTCGCGGGGCGCTGttttt                              >  2:152060/1‑250 (MQ=255)                                                                                                                                                                                                                              aGGCGTGGAATATCGGCGGT‑GCC‑TTTTATCACCAAAT‑CC‑A‑GCCGGTCTGC‑ATAGGGGAAGCGTGAAGAATCTTCTTT‑GACTAAAAACAGATG‑C‑AACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGGAAA‑GGG‑ATTGGCGTTCATCACAATGCAGCCAAT‑CTTGTTCCCTGGATGACGAAATTTTTTCGGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCAt                           <  1:411867/251‑1 (MQ=255)                                                                                                                                                                                                                               ggCGTGGAATATCGGCGGT‑GCC‑TTGTAACACCAAAT‑CC‑A‑GCCGGTCTTC‑ATAGGGGAAGCGTGAAGAATCTTCTTT‑GACTAAAAACAGATG‑C‑AACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAA‑GGG‑ATTGGCGTTCATCACAATGCAGCCAAT‑CTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATc                          <  2:87882/251‑1 (MQ=255)                                                                                                                                                                                                                                     gAATATCGGCGGT‑GCC‑TTTTAACACCAAAT‑CC‑A‑GCCGGTCTTC‑ATAGGGGAAGCGTGAAGAATCTTCTTT‑GACTAAAAACAGATG‑C‑AACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAA‑GGG‑ATTGGCGTTCATCACAATGCAGCCAAT‑CTTGTTCCCTGTATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCAcc                    <  1:183774/251‑1 (MQ=255)                                                                                                                                                                                                                                       atatCGGCGGT‑GCC‑TTTTAACACCAAAT‑CC‑A‑GCCGGTCTTC‑ATAGGGGAAGCGTGAAGAATCTTCTTT‑GACTAAAAACAGAGG‑C‑AACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAA‑GGG‑AATGGCGTTCATCACAATGCAGCCAAT‑CTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCAt                  <  1:318857/251‑1 (MQ=255)                                                                                                                                                                                                                                           cggcgtT‑GCC‑GTTTAGCGCCAAAG‑CC‑A‑GCCGGTCTTC‑ATAGGGGAAGCGTGAAGAATCTTCCTT‑CACTAAAAACAGATG‑C‑ACTCAGTCTCACTGTGACTCAGCCTGTTGCATCAGATAACGGTGACCATCTGTAAA‑GGG‑ATTGTCATTCATCACAATGCAGCCAAT‑CTTGTTCCCTGGATGACGAAATTTTTTCAGCGCTTCGGCATATCGTTTCGGTCGCGTGGCGCGGTTTTCCTTCAGCACCATCACg              <  2:29821/251‑1 (MQ=255)                                                                                                                                                                                                                                           cggcggT‑GCC‑TTTTAACACCAAAT‑CC‑A‑GCCGGTCTTC‑ATAGGGGAAGCGTGAAGAATCTTCTTT‑GACTAAAAACAGATG‑C‑AACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAA‑GGG‑ATTGGCGTTCATCACAATGCAGCCAAT‑CTTGTTCCCTGGATGACGAAATTTTTTCAGCGAGTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCAc               <  1:433120/250‑1 (MQ=255)                                                                                                                                                                                                                                                         ttAACCCCAAAT‑AC‑A‑GCCGGACTTC‑ATAGGGGAAGCGTGAAGAATCTTCTTT‑GACTAAAAACAGATG‑C‑AACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAA‑GGG‑ATTGGCGTTCATCACAATGCAGCCAAT‑CTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCt   <  2:178581/250‑1 (MQ=255)                                                                                                                                                                                                                                                          tAACACCAAAT‑CC‑A‑GCCGGGCTTC‑ATAGGGGAAGCGTGAAGAATCTTCCTT‑GACTAAAAACAGATG‑C‑AACCAGTCGCACTGGGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAA‑GGG‑ATTGGCGTTCATCACAATGCAGCCAAT‑CTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCTTGGCGCTGTTGTCCATCAGCACCATCACGCCTGGTACGCTg  <  2:212119/250‑1 (MQ=255)                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                   CCAGCCAGTAGCGCATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGT‑GCC‑TTTTAACACCAAAT‑CC‑A‑GCCGGTCTTC‑ATAGGGGAAGCGTGAAGAATCTTCTTT‑GACTAAAAACAGATG‑C‑AACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAA‑GGG‑ATTGGCGTTCATCACAATGCAGCCAAT‑CTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTG  >  NC_000913/651019‑651515

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