Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	77,776	C→A	G52G (GGC→GGA)	setA →	sugar exporter SetA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	77,776	0	C	A	95.5%	101.5 / ‑5.0	46	G52G (GGC→GGA)	setA	sugar exporter SetA
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base A (19/23);  total (20/26)
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.25e-01

GTGACTGAGTATTGGTGTAAAATCACCCGCCAGCAGATTATACCTGCTGGTTTTTTTTATTCTCGCCGCGCTAAAAAGGGAACGTATGATCTGGATAATGACGATGGCTCGCCGTATGAACGGTGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGT‑GGCCGGGGCGCTACAGGCTCCTACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTT‑TCTGGATCGGCCTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGTTTGATGGCTATCGGCAATGCGCTATTGTTTGCATTTAATCGTCATTATCTGACGCTTATCACCTGTGGTGTGCTTCTGGCATCTCTGGCCAATACGGCAATGCCACAGTTATTTGCTCTGGCGCGGGAATATGCGGAT  >  NC_000913/77536‑78019                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    gggACTGAGTATTGGTGTAAAATCCCCCGCCAGCAGATTATACCTGCTGGTTTTTTTTATTCTCGCCGCGCTAAAAAGGGAACGTATGATCTGGATAATGACGATGGCTCGCCGTATGAACGGTGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGT‑GGCCGGGGCGCTACAGGCTCCTACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTT‑TCTGGATCGGACTCTTTtata                                                                                                                                                                                                                                           <  1:121753/249‑1 (MQ=255)       gAGTATTGGTGTAAAATCCCCCGCCAGCAGATTATACCTGCTGGTTTTTTTTATTCTCGCCGCGCTAAAAAGGGAACGTATGATCTGGATAATGACGATGGCTCGCCGTATGAACGGTGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGT‑GGCCGGGGCGCTACAGGCTCCTACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTT‑TCTGGATCGGACTCTTTTATACGGTGa                                                                                                                                                                                                                                     <  1:322462/251‑1 (MQ=255)           aTTGGTGTAAAATCACCCGCCAGCAGATTATACCTGCTGGTTTTTTTTATTCTCGCCGCGCTAAAAAGGGAACGTATGATCTGGATAATGACGATGGCTCGCCGTATGAACGGTGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGT‑GGCCGGGGCGCTACAGGCTCCTACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTT‑TCTGGATCGGACTCTTTTATACGGTGAATGc                                                                                                                                                                                                                                 >  1:65963/1‑251 (MQ=255)               gtgtAAAATCCCCCGCCAGCAGATTATACCTGCTGGTTTTTTTTATTCACGCCGCGCTAAAAAGGGAACGTATGATCTGGCTAATGACGATGGCTCGCCGTATGAACGGTGTTTACGCGGCATTTCTGCTGGTCGCTTTTATGATGGGGGA‑GGCCGGGGCGCTATAGGCTCCTACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCACCCTT‑TCTGGATCGGACTCTTTACTCCGGTGAATGCTAtt                                                                                                                                                                                                                             >  2:225702/1‑251 (MQ=255)                     aacaccccGCCAGCAGATTATACCTGCTGGTTTTTTTTATTCTCGCCGCGCTAAAAAGGGAACGTATGATCTGGATAATGACGATGGCTCGCCGTATGAACGGTGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGT‑GGCCGGGGCGCTACAGGCTCCTACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTT‑TCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTggg                                                                                                                                                                                                                       <  1:377758/246‑1 (MQ=255)                      atcccccGCCAGCAGATTATACCTGCTGGTTTTTTTTATTCTCGCCGCGCTAAAAAGGGAACGTATGATCTGGATAATGACGATGGCTCGCCGTATGAACGGTGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGT‑GGCCGGGGCGCTACAGGCTCCTACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTT‑TCTGGATCGGACTCTGTTATACGGTGAATGCTATTGCTGGGa                                                                                                                                                                                                                      <  1:368336/247‑1 (MQ=255)                                  cagATTATACCTGCTGGTTTTTTTTATTCTCGCCGCGCTAAAAAGGGAACGTATGATCTGGATAATGACGATGGCTCGCCGTATGAACGGTGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGT‑GGCCGGGGCGCTACAGGCTCCTACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTT‑TCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGcc                                                                                                                                                                                                          <  1:124914/251‑1 (MQ=255)                                   agATTATACCTGCTGGTTTTTTTTATTCTCGCCGCGCTAAAAAGGGAACGTATGATCTGGATAATGACGATGGCTCGCCGTATGAACGGTGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGT‑GGCCGGGGCGCTACAGGCTCCTACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTT‑TCTTGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCct                                                                                                                                                                                                         >  1:250680/1‑251 (MQ=255)                                               ctgGTTTTTTTTATTCTCGCCGCGCTAAAAAGGGAACGTATGATCTGGATAATGACGATGGCTCGCCGTATGAACGGTGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGT‑GGCCGGGGCGCTACAGGCTCCTACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCCACCTT‑TCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCaa                                                                                                                                                                                              >  1:42574/1‑250 (MQ=255)                                                    tttttttATTCTCGCCGCGCTAAAAAGGGAACGTATGATCTGGATAATGACGATGGCTCGCCGTATGAACGGTGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGT‑GGCCGGGGCGCTACAGGCTCCTACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTT‑TCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGtt                                                                                                                                                                                        >  1:161603/1‑251 (MQ=255)                                                       ttttATTCTCGCCGCGCTAAAAAGGGAACGTATGATCTGGATAATGACGATGGCTCGCCGTATGAACGGTGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGT‑GGCCGGGGCGCTACAGGCTCCTACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTT‑TCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTg                                                                                                                                                                                     >  1:194120/1‑251 (MQ=255)                                                       gtttATGCTCGCGTCGCTAAAAAGGGACCGTATGTGTTGGATTAGTACGATGGCCCTCCGGTTGAACGGTGTGTACGCGGCATTTATTCTTGTCGCGTTTATGATGGGGGT‑GGCTTGGGCGCTATAGGCGCCTCCATTGAGCTTATATCTGAGTCGTGAGGTTGGCGCGCA‑CTTTGTCTGGATCGGCCTCTTTTATACGGTGAATGCTATTGCTGCGATCGGCTTAAGCTTCTTGTTGGCAAAACGTTCTg                                                                                                                                                                                     <  2:26014/250‑1 (MQ=255)                                                             tctcGCCGCGCTAAAAAGGGAACGTATGATCTGGATAATGACGATGGCTCGCCGTATGAACGGTGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGT‑GGCCGGGGCGCTACAGGCTCCTACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTT‑TCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTc                                                                                                                                                                               >  1:23180/1‑251 (MQ=255)                                                                      gcTAAAAAGGGAACGGATGAGCTGGATAATGACGATGGCTCGCCGTAGGAACGGTGTTTACGCGGCATTTATGCTTGTCGCCTTTATGATGGGGGT‑GGCCGGGGCGCTACAGGCTCCTACATTGAGCTTATTGCTGAGTCGTGAGGTTGGCGCGCAACCTT‑TCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATc                                                                                                                                                                      <  1:434204/251‑1 (MQ=255)                                                                         aaaaaGGGAACGTATGATCTGGATAATGACGATGGCTCGCCGTATGAACGGTTTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGT‑GGCCGGGGCGCTACAGGCTCCTACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTT‑TCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGc                                                                                                                                                                   <  2:121974/251‑1 (MQ=255)                                                                           aaaGGGAACGTATGATCTGGATAATGACGATGGCTCGCCGTATGAACGGTGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGT‑GGCCGGGGCGCTACAGGCTCCTACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTT‑TCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCg                                                                                                                                                                  >  2:200312/1‑250 (MQ=255)                                                                           aaaGGGAACGTATGATCTGGATAATGACGATGGCTCGCCGTATGAACGGTGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGT‑GGCCGGGGCGCTACAGGCTCCTACATTGAGCTTATTTCTGAGGCGTGAGGTTGGCGCGCAACCTT‑TCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTTGCAAAACGTTCTGACAGTCAGGGCGGTCGGCg                                                                                                                                                                  >  2:250698/1‑250 (MQ=255)                                                                              gggAACGTATGATCTGGATAATGACGATGGCTCGCCGTATGAACGGTGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGT‑GGCCGGGGCGCTACAGGCTCCTACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTT‑TCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGaaaa                                                                                                                                                              <  2:194120/251‑1 (MQ=255)                                                                                                 aaTGACGATGGCTCGCCGTATGAACGGTGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGT‑GGCCGGGGCGCTACATGCTCCTACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTT‑TCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGt                                                                                                                                            <  1:212289/250‑1 (MQ=255)                                                                                                    gACGATGGCTCGCCGTATGAACGGTGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGT‑GGCCGGGGCGCTACAGGCTCCTACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTT‑TCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGc                                                                                                                                                                   <  1:256295/224‑1 (MQ=255)                                                                                                    gACGATGGCTCGCCGTATGAACGGTGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGT‑GGCCGGGGAGCTACAGGCTCCTACATTGAGCTTATTTCTGGGTCGTGAGGTTGGCGCGCAACCTT‑TCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTACGGGCGATCGGc                                                                                                                                                                   >  2:256295/1‑224 (MQ=255)                                                                                                           gCTCGCCGTATTAACGGTGTTTACGCTGCATTTATGCTGGTCGCTTTTATGATGGGGGT‑GGTCGGGGCGCTTCAGGCTCCTACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTT‑TCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGTCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGTTTGATGGCTAt                                                                                                                                 <  2:310176/251‑1 (MQ=255)                                                                                                                           gtgtTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGT‑GGCAGGGGCGCTACAGGCTCCTACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTT‑TCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGTTTGATGGCTATCGGCAATGCGCTAttt                                                                                                                 >  1:346976/1‑250 (MQ=255)                                                                                                                             gtTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGT‑GGCCGGGGCGCTACAGGCTCCTACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTT‑TCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGTTTGATGGCTATCGGCAATGCGCTATTGtt                                                                                                               <  1:250698/251‑1 (MQ=255)                                                                                                                             ggttACGCGGCATTAATGTTGTTTGCTTTTATG‑TGGGGGTTGGCCGGGGCGCTCCAGGCTCTTACATTGAGCTTATTTCTGAGGCGTGAGGTTGGCGCGCAACCTT‑TCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAACCTGATTATATTGTGCTGTTTGATGGCTATCGGCAATGCGCTATTGtt                                                                                                               <  2:4557/249‑1 (MQ=255)                                                                                                                                           tATGCTGGTCGCTTTTATGATGGGGGT‑GGCCGGGGCGCTACAGGCTCCTACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTT‑TCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGTTTGATGGCTATCGGCAATGCGCTATTGTTTGCATTTAATCGt                                                                                                  <  1:167601/250‑1 (MQ=255)                                                                                                                                           tATGCTGGTCGCTTTTATGATGGGGGT‑GGCCGGGGCGCTACAGGCTCCTACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTT‑TCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGTTTGATGGCTATCGGCAATGCGCTATTGTTTGCATTTAATCGt                                                                                                  <  1:177887/250‑1 (MQ=255)                                                                                                                                                 ggTCGCTTTTATGATGGGGGT‑GGCCGGGGCGCTACAGGCTCCTACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTT‑TCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGTTTGATGGCTATCGGCAATGCGCTAtt                                                                                                                  >  2:406984/1‑228 (MQ=255)                                                                                                                                                 ggTCGCTTTTATGATGGGGGT‑GGCCGGGGCGCTACAGGCTCCTACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTT‑TCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGTTTGATGGCTATCGGCAATGCGCTAtt                                                                                                                  <  1:406984/228‑1 (MQ=255)                                                                                                                                                 ggTCGCTTTTATGATGGGGGT‑GGCCGGGGCGCTACAGGCTCCTACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTT‑TCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGTTTGATGGCTATCGGCAATGCGCTATTGTTTGCATTTAATCGTCATTATc                                                                                           >  2:348229/1‑251 (MQ=255)                                                                                                                                                         ttATGATGGGGGT‑GGCCGGGGCGCTACAGGCTCCTACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTT‑TCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGTTTGATGGCTATCGGCAATGCGCTATTGTTTGCATTTAATCGTCATTATCTGACGCt                                                                                    <  1:258370/250‑1 (MQ=255)                                                                                                                                                                   ggT‑GGCTGGGTCGCTACAGGCTCCTACATTGAGCTTATGTCTGAGTCGTGAGGTTGGCGCGCAACCTT‑TCTGGATCGGGCTCTTTTATACGGTGAATGCTATTGCTGGGAACGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGTTTGATGGCTATCGTCGATGCGCGATTGTTTGCATTTAATCGTCATTATCTGTCGCTTATCACCtgtg                                                                         <  1:279656/251‑1 (MQ=255)                                                                                                                                                                              gcgcTACAGGCTCCTACATTGAGCTTATTGCTGAGTCGTGAGGTTGGCGCGCAACCTT‑TCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGTTTGATGGCTATCGGCAATGCTCTATTGTTTGCATTTAATCGTCATTAGCTGACGCTTATCACCTGTGGTGTGCTTCt                                                               <  1:382674/251‑1 (MQ=255)                                                                                                                                                                                     aGGCTCCTACATTGAGCTTATTTCTGCGTCGTGAGGTTGGCGCGCAACCTT‑TCTGGACCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACGGTCAAGGCTATCGGCGAATACTGATTATATTTTGCAGTTTGCTGGCTATCGGCAATGCGCTATTGTTTGCATTTAATCGTCATCATCTGGCGCTTAGCACCTGTGGTGTGCTTCGGGCAtct                                                        >  2:94293/1‑251 (MQ=255)                                                                                                                                                                                       gCTCCTACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTT‑TCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGTTTGATGGCTATCGGCAATGCGCTATTGTTTGCATTTAATCGTCATTATCTGACGCTTATCACCTGTGGTGTGCTTCTGGCATctct                                                      <  1:30296/251‑1 (MQ=255)                                                                                                                                                                                                          ttCTGAGTCGAGAGGTTGGCGCGCAACCTT‑TCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGACAACTGCTTATATTTTGCTGTTTGATGGCTATCGGCAATGCGCTATTGTTTGCATTTAATCGTCATTATCTGACGCTTATCACCTGTGGTGTGCTTCTGGCATCTCTGGCCAATACGGCAATGccc                                   >  2:41638/1‑250 (MQ=255)                                                                                                                                                                                                            cTGAGTCGTGAGGTTGGCGCGCAACCTT‑TCTGGATCGGCCTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGTTTGATGGCTATCGGCAATGCGCTATTGTTTGCATTTAATCGTCATTATCTGACGCTTATCACCTGTGGTGTGCTTCTGGCATCTCTGGCCAATACGGCAATGCcac                                  >  2:151756/1‑250 (MQ=255)                                                                                                                                                                                                                   atGAGGGTGACGCGCACCCTT‑TCTGGATCGGACTCTTTTCTACGGTGATTGCTATTGCTGGGCTCGGCGTACGCATCTGGTTGGCAAAACGTTCTGCCAGTCTGGGCGATCGGCGCTAACTGATTATATTTTGCTGTTTGATGGCTTTCGGCAATGCGCTATTGTTTGCATTTAATCTTAATTATCTGACGCTTATCACCTGTAGTGTGATTCTTGCCTCTCTGTCCAATCCGTCCATGCCACGGTTAttt                          >  2:87815/2‑251 (MQ=255)                                                                                                                                                                                                                            gcgcgcAACCTT‑TCTGGATCGGGCTCTTTTATACGGTGGATGCTATTGCGGGGATCGGCGTAAGCCTCTGGGTGGCAAAACTTTCTGGCAGGCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGTTTGATGGCTATCGGCAATGCGCGATTGTTTGCATTGAATCGTCATTAGCTGACGCTTAGCTCCTGTGGGGTGCTTCTGGCATCTCTGGGCAATACGGCAATGCCACAGTTATTTGCTCTGgcg                 <  2:42574/251‑1 (MQ=255)                                                                                                                                                                                                                                gcAACCTT‑TCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGTTTGATGGCTATCGGCAATGCGCTATTGTTTGCATTTAATCGTCATTATCTGACGCTTATCACCTGTGGTGTGCTTCTGGCATCTCTGGCCAATACGGCAATGCCACAGTTATTTGCTCTGGCGCggg             <  1:41638/251‑1 (MQ=255)                                                                                                                                                                                                                                   aCCTT‑TCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTACGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGTTTGATGGCTATCGGCAATGCGCTATTGTTTGCATGTAATCGTCATTATCTGACGCTTATCACCTGTGGTGTGCTTCTGGCATCTCTGGCCAATACGGCAATGCCACAGTTATTTGCTCTGGCGCGGGAat          <  1:17340/251‑1 (MQ=255)                                                                                                                                                                                                                                      tt‑tCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGTTTGATGGCTATCGGCAATGCGCTATTGTTTGCATTTAATCGTCATTATCTGACGCTTATCACCTGTGGTGTGCTTCTGGCATCTCTGGCCAATACGGCAATGCCAGAGTTATTTGCTCTGGCGCGGGAATATg       >  2:313867/1‑251 (MQ=255)                                                                                                                                                                                                                                            ggttCGGACTCTTTTAGACGGTGAATGCTATTGCTGGGATCGGCGTAATCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCTGCGAAAACTGATTATATTTTGCTGTTTGATGGCTATCGGCAATGCGCTATTGTTTGCATTTAATCGTCATTATCTGACGCTTATCACCTGGGGTGTGCTTCTGGCATCTTTGGCTAATACGGCAATGCCACAGTTATTTGCGCTGGCGCGGGAATATGCg     <  2:59448/245‑1 (MQ=255)                                                                                                                                                                                                                                             gATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGTTTGATGGCTATCGGCAATGCGCTATTGTTTGCATTTAATCGTCATTATCTGACGCTTATCACCTGTGGTGTGCTTCTGGCATCTCTGGcc                                                  >  1:40980/1‑202 (MQ=255)                                                                                                                                                                                                                                             gATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGTTTGATGGCTATCGGCAATGCGCTATTGTTTGCATTTAATCGTCATTATCTGACGCTTATCACCTGTGGTGTGCTTCTGGCATCTCTGGcc                                                  <  2:40980/202‑1 (MQ=255)                                                                                                                                                                                                                                             gATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGGTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGTTTGATGGCTATCGGCCATGCGCTATTGTTTGCATTTAATCGTCATTATCTGACGCTTATCACCTGTGGTGTGCTTCTGGCATCTCTGGCCAATACGGCAATGCCACAGTTATTTGCTCTGGCGCGGGAATATGCGGAt  >  1:70458/1‑250 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    GTGACTGAGTATTGGTGTAAAATCACCCGCCAGCAGATTATACCTGCTGGTTTTTTTTATTCTCGCCGCGCTAAAAAGGGAACGTATGATCTGGATAATGACGATGGCTCGCCGTATGAACGGTGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGT‑GGCCGGGGCGCTACAGGCTCCTACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTT‑TCTGGATCGGCCTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGTTTGATGGCTATCGGCAATGCGCTATTGTTTGCATTTAATCGTCATTATCTGACGCTTATCACCTGTGGTGTGCTTCTGGCATCTCTGGCCAATACGGCAATGCCACAGTTATTTGCTCTGGCGCGGGAATATGCGGAT  >  NC_000913/77536‑78019

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