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.9%	128.7 / ‑1.8	50	G52G (GGC→GGA)	setA	sugar exporter SetA
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base A (23/24);  total (24/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 = 4.82e-01

ATGTGTGACTGAGTATTGGTGTAAAATCACCCGCCAGCAGATTATACCTGCTGGTTTTTTTTATTCTCGCCGCGCTAAAAAGGGAACGTATGATCTGGATAATGACGATGGCTCGCCGTATGAACGGTGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGTGGCCGGGGCGCTACAGGCTCCTACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGCCTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGTTTGATGGCTATCGGCAATGCGCTATTGTTTGCATTTAATCGTCATTATCTGACGCTTATCACCTGTGGTGTGCTTCTGGCATCTCTGGCCAATACGGCAATGCCACAGTTATTTGCTCTGGCGCGGGAATATGCGGA  >  NC_000913/77532‑78018                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                   attgttgacggagtttggggttAAAATCCCCCGCCAGCAGATTATTCCTGCTGGTTTTTTTTATTCTCGCCGCGCTAAAAAGGGAACGTATGATCTGGATAATGACGATGGCTCGCCGTAGGAACGGTGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGTGGCCGGGGCGCTACAGGCTCCTACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCttt                                                                                                                                                                                                                                              <  2:194017/230‑1 (MQ=255)   gtgtgACTGATTATTGGTGTAAAATCCCCCGCCAGCAGATTATACCTGCTGGTTTTTTTTATTCTCGCCGCGCTAAAAAGGGAACGTATGATCTGGATAATGACGATGGCTCGCCGTATGAACGGTGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGTGGCCGGGGCGCTACAGGCTCCTACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCtttt                                                                                                                                                                                                                                             <  1:487198/250‑1 (MQ=255)        atttaGTATTGGTTAAAAATCATCCGCCAGCATATTATACCTGCTGGTTTTTTTTATTCTCGCCGCGCTAAAAAGGGAACGTATTATCTGGATAATGACGATGGCTCGCCGTATGAACGGTGTTTATGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGTGGCCGGGGCGCTACAGGCTCCTACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACgg                                                                                                                                                                                                                                       <  1:673854/247‑1 (MQ=255)                     gtAAAATCCCCCGCCAGCATATTATACCTGCTGGTTTTTTTTATTCTCGCCGCGCTAAAAAGGGAACGTATGATCTGGATAATGACGATGGCTCGCCGTAGGAACGGTGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGTGGCCGGGGCGCTACAGGCTCCTACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGc                                                                                                                                                                                                                          <  2:326573/251‑1 (MQ=255)                                 gCCAGCAGATTATACCTGCTGGTTTTTTTTATTCTCGCCGCGCTAAAAAGGGAACGTATGATCTGGATAATGACGATGGCTCGCCGTATGAACGGTGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGTGGCCGGGGCGCTACAGGCTCCTACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGt                                                                                                                                                                                                              >  2:474806/1‑251 (MQ=255)                                    agcagATTATACCTGCTGGTTTTTTTTATTCTCGCCGCGCTAACCAGGGAACGTATGATCTGGATAATGACGATGGCTAGCCGTATGAACGGTGTTTACGCGGCATTTATGCAGGTAGCTTTTATGATGGGGGTGGCCGGGGCGCTACAGGCTCGTACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCCAACTTTCTGGATCGGATTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTa                                                                                                                                                                                                             >  1:606648/1‑249 (MQ=255)                                         attttaCCTGTTGGTTTTTTTTATTCTCGCCGCGCTAAAAAGGGAACGTATTATCTGGATAATGACGATGGCTCGCCGTATGAACGGTGTTTACGCGGCATTGATGCTGGTCGCTTTTATGATGGGGGTGGCCGGGGCGCTACAGGCTCCTACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCGTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCctct                                                                                                                                                                                                      <  2:670220/247‑1 (MQ=255)                                          ttATACCTGCTGGTTTTTTTTATTCTCGCCGCGCTAAAAAGGGAACGTATGATCTGGATAATGACGATGGCTCGCCGTCTGAACGGTGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGTGGCCGGGGCGCTACAGGCTCCTACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAAACTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTg                                                                                                                                                                                                     >  1:667129/1‑251 (MQ=255)                                               ccTGCTGGTTTTTTTTATTCTCGCCGCGCTAAAAAGGGAACGTATGATCTGGATAATGACGATGGCTCGCCGTATGAATGGTGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGTGGCCGGGGCGCTACAGGCTCCTACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTgg                                                                                                                                                                                                <  1:220751/251‑1 (MQ=255)                                                         ttttttATTCTCGCCGCGCTAAAAAGGGAACGTATGATCTGGATAATGACGATGGCTCGCCGTATGAACGGTGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGTGGCCGGGGCGCTACAGGCTCCTACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGtt                                                                                                                                                                                       <  1:112800/250‑1 (MQ=255)                                                            tttATTCTCGCCGCGCTAAAAAGGGAACGTATGATCTGGATAATGACGATGGCTCGCCGTATGAACGGTGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGTGGCCGGGGCGCTACAGGCTCCTACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTg                                                                                                                                                                                    >  1:564531/1‑250 (MQ=255)                                                               aTTCTCGCCGCGCTAAAAAGGGAACGTATGATCTGGATAATGACGATGGCTCGCCGTATGAACGGTGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGTGGCCGGGGCGCTACAGGCTCCTACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGAc                                                                                                                                                                                  >  1:708389/1‑249 (MQ=255)                                                                    cgccgcGCTAAAAAGGGAACGTATGATCTGGATAATGACGATGGCTCGCCGTATGAACGGTGTTTACGCGGCATTTATGCTGGTCGTTTTTATGATGGGGGTGGCCGGGGCGCTACAGGCTCCTACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAgg                                                                                                                                                                           <  1:383662/251‑1 (MQ=255)                                                                         cgcTAAAAAGGGAACGTATGATCTGGATAATGATGATGGCTCGCTGTATGAACGGTGTTTACGCGGCATTTATGCTGGGTGCTTTTATGATGGGGGTGGGCGGTGCGCTACAGGCTCCTACATTGAGCTTATTGCTGAGTCGTGAGGTTGGCGCTCAATCTTTCTTGATCTGTCTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCCAAACGTTCTGACATTCAGGGCgtt                                                                                                                                                                      <  2:284299/251‑3 (MQ=255)                                                                          gcTAAAAAGGGAACGTATGATCTGGATAATGACGATGGCTCGCCGTATGAACGGTGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGTGGCCGGGGCGCTACAGGCTCCTACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATc                                                                                                                                                                     <  1:461356/251‑1 (MQ=255)                                                                              aaaaGGGAACGTATGATCTGGATAATGACGATGGCTCGCCGTATGAACGGTGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGTGGCCGGGGCGCTACAGGCTCCTACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCg                                                                                                                                                                 <  1:140007/251‑1 (MQ=255)                                                                              aaaaGGGAACGTATGATATGGATCATGACGATGGCTCGCCGTATGAACGGTGTTTACGCGGCATTTAAGATGGTCGCTTTTATGATGGGGGTGGCCGGGGCGCTACAGGCTCCTACATTGTGCTTGTTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACAGTGAATGCTATTGCTGTGATCGGCGTAAGACTCTGGTTGGCACAACGTTCTGACAGTCGGGGCGAtagggg                                                                                                                                                                 >  2:304995/1‑246 (MQ=255)                                                                                         tATGATGTGGATAATGACGATGGCTCGCCGTATGAACGGTGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGTGGCCGGGGCGCTACAGGCTCCTACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGAGCAACCTTTCTGGATCGGCCTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTAAGGGCGATCGGCGAAAACTGAtt                                                                                                                                                       >  2:320683/1‑250 (MQ=255)                                                                                          aTGATCTGGATAATGACGATGGCTCGCCGTATGAACGGTGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGTGGCCGGGGCGCTACAGGCTCCTACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATtat                                                                                                                                                     >  1:702686/1‑251 (MQ=255)                                                                                          aTGATCTGGATAATGACGATGGCTCGCCGTATGAACGGAGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGTGGCCGGGGCTCTACAGGCTCCTACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGTGCATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAACACTGTTtat                                                                                                                                                     >  2:637571/1‑251 (MQ=255)                                                                                                 ggATAATGACGATGGCTCGCCGTATGAACGGTGTTTACGCGGCATTTATGCTGGTCTCTTTTATGATGGGGGTGGCCGGGGCGCTACAGGCTCCTACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTg                                                                                                                                               <  2:564531/250‑1 (MQ=255)                                                                                                      aTGACGATGGCTCGCCGTATGAACGGTGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGTGGCCGGGGCGCTACAGGCTCCTACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGCCTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTAAGGGCGATCGGCGAAAACTGATTATATTTTGCTGttt                                                                                                                                         <  1:320683/251‑1 (MQ=255)                                                                                                       tGACGATGGCTCGCCGTATGAACGGTGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGTGGCCGGGGCGCTACAGGCTCCTACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGTTTg                                                                                                                                        >  1:709012/1‑251 (MQ=255)                                                                                                           gATGGCTCGCCGTATGAACGGTGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGTGGCCGGGGCGCTACAGGCTCCTACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAg                                                                                                                                                                                <  1:259713/207‑1 (MQ=255)                                                                                                           gATGGCTCGCCGTATGAACGGTGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGTGGCCGGGGCGCTACAGGCTCCTACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAg                                                                                                                                                                                >  2:259713/1‑207 (MQ=255)                                                                                                              ggCTCGCCGTATGAACGGTGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGTGGCCGGGGCGCTACAGGCTCCTACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGTTTGATGGCt                                                                                                                                  <  1:534260/250‑1 (MQ=255)                                                                                                               gCTCGCCGTATGAACGGTGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGTGGCCGGGGCGCTACAGGCTCCTACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGCCAGTCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGTTTGATGGCTAt                                                                                                                                >  2:554371/1‑251 (MQ=255)                                                                                                                           aaCGGGGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGTGGCCGGGGCGCTACAGGCTCCTACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGTTTGATGGCTATCGGCAATGCGCt                                                                                                                    <  2:256835/251‑1 (MQ=255)                                                                                                                            ccGGTGTTTAAGCGGCATATATGCTGGGCGCTTTTATGATGGGGGTGGCCGGGGCGCTACATGCTCCTACATTGAGCTTATTGCGGGGTCGAGATGTTGGCGCGACACCTTTCTGGATAGGACTCTTTTATACGGGGAAGGCTATGTCTGGGATCGGCTTAAGCCTCTGGTTGGCACAACGTTATGAGAGTCAGAGCGTTCGGCGACAACTGGTTATATtttgatgtttgatgtatatcgcaaaggcgatt                                                                                                                   >  2:517856/2‑233 (MQ=255)                                                                                                                                      cgcgGCAGTTATGCTGGTCGCTTTTATGATGTGGGTGGCCGGGGCGCTACAGGCTCCTTCATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACTGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGTTTGATGGCTATCGGCAATGCGCTATTGTTTGCAt                                                                                                         <  2:229197/251‑1 (MQ=255)                                                                                                                                       gcgGCATTTATGCTGGTCGCTTTTATGATGGGGGTGGCCGGGGCGCTACAGGCTCCTACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGTTTGATGGCTATCGGCAATGCGCTATTGTTTGCAtt                                                                                                        <  1:593176/251‑1 (MQ=255)                                                                                                                                                    tGGTCGCTTTTATGATGGGGGTGAACGGGGCGCTACAGGCTCCTACATTGTGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGTAAACTGATTATATTTTGCTGTTTGATGGCTATCGGCAATGCGCTATTGTTTGCATTTAATCGTCATTAt                                                                                           >  2:408775/1‑251 (MQ=255)                                                                                                                                                     ggTCGCTTTTATGATGGGGGTGGCCGGGGCGCTACAGGCTCCTACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGTTTGATGGCTATc                                                                                                                               >  1:639735/1‑214 (MQ=255)                                                                                                                                                     ggTCGCTTTTATGATGGGGGTGGCCGGGGCGCTACAGGCTCCTACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGTTTGATGGCTATc                                                                                                                               <  2:639735/214‑1 (MQ=255)                                                                                                                                                     ggTCGCTTTTATGATGGGGGTGGCCGGGGCGCTACAGGCTCCTACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGTTTGATGGCTATCGGCAATGCGCTATTGTTTGCATTTAATCGTCATTATc                                                                                          <  1:142701/251‑1 (MQ=255)                                                                                                                                                                 gatgGGGGTGGCCGGGGCGCTACAGGCTCCTACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGaa                                                                                                                                                               <  2:580378/170‑1 (MQ=255)                                                                                                                                                                    gggggAGGCCGGGGCGCTACAGGCTCCTACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGTTTGATGGCTATCGGCAATGCGCTATTGTTTGCATTTAATCGTCATTATCTGACGCTTATCACCt                                                                           >  1:139407/1‑251 (MQ=255)                                                                                                                                                                        gTGGCCGGGGCGCTACAGGCTCCTACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGTTTGATGGCTATCGGCAATGCGATATTGTTTGCATTTAATCGTCATTATCTGACGCTTATCACCtgt                                                                         >  1:39849/1‑249 (MQ=255)                                                                                                                                                                         gggCCGGGGCGCTACAGGCTCCTACATTGAGCTTATTTCTGGGTCGTGAGGTTGGCGCGCAACCGTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGTTTGATGGCTATCGGCAATGCGCTATTGTTTGCATTTAATCGTCATTATCTGACGCTTATCACCTgtggt                                                                      <  2:376888/250‑1 (MQ=255)                                                                                                                                                                           gCCGGGGCGCTACAGGCTCCTACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGTTTGATGGCTATCGGCAATGCGCTATTGTTTGCATTTAATCGTCATTATCTGACGCTTATCACCTGTGgtgt                                                                    <  1:433762/251‑1 (MQ=255)                                                                                                                                                                              ggggCGCTACAGGCTCCTACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGTTTGATGGCTATCGGCAATGCGCTATTGTTTGCATTTAATCGTCATTATCTGACGCTTATCACCTGTGGTTTGCt                                                                 >  2:7055/1‑251 (MQ=255)                                                                                                                                                                               gggCGCTAAAGGCACTTACTTTGAGCGTATTTCAGATTCGTGGGGTTGGCGCGGAACTGTGCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGAAAGCCTATGGTTGGCAAAGCGGTCTGACAGTCAGTGCGATCTGCGAAAACTGATTATATTTTGCGGTTTGATGGTGATCGGCAATGCGCTATTGTTTTCATTTAATCGTAATTATCTGACGCTTATCACCTGTGGTGTGCtt                                                                <  2:212326/251‑1 (MQ=255)                                                                                                                                                                                ggCGCTACAGGCTCCTAGATTGAGCTGATTTCGGATGCGTGAGGTTGGCGCGCAACCCTTCTGGATCGGACTCTTTTATACGGCTAATGCTATTGCTGGGATCTGCGTAAGCCGCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTTATTATATTTTGCTGTGTGATGGCTATCGGCAATGCGCTATTGTTTGCATTTAATCTTCATTATCTGACGCTTATCACCTGTGGTGTGCTTc                                                               <  2:230508/251‑1 (MQ=255)                                                                                                                                                                                ggCGCTACAGGATCCTACATTGAGCTTATTTCTGAGACGTGAGGTTGGCGCGCAACATTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTAAGGGCGATCGGCGAAAACTGATTATATTTTGCTGTTTGATGGCTATCGGCAATGCGATATTGTTTGCATTTAATCGTCATTATCTGACGCTTATCACCTGTGGTGTGCttt                                                               >  2:500516/1‑250 (MQ=255)                                                                                                                                                                                            tCCTACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCCACCTTTCTGGATCGGACTCTTTTATACGGTGCATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAACCGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGTTTGATGGCTATCGGCAATGCGCTATTGTTTGCATTTAATCGTCATTATCTGACGCTTATCACCTGTGGTGTGCTTCTGGCATCTCTgg                                                   >  2:576368/1‑251 (MQ=255)                                                                                                                                                                                                                    tCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGTTTGATGGCTATCGGCAATGCGCTATTGTTTGCATTTAATCGTCATTATCTGACGCTTATCACCTGTGGTGTGCTTCTGGCATCTCTGGCCAATACGGCAATGCCACAGTTAt                           >  1:467651/1‑251 (MQ=255)                                                                                                                                                                                                                               gcgcgcAACCTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGTTTGATGGCTATCGGCAATGCGCTATTGTTTGCATTTAATCGTCATTATCTGACGCTTATCACCTGTGGTGTGCTTCTGGCATCTCTGGCCAATACGGCAATGCCACAGTTATTTGCTCTGgcg                >  2:18060/1‑251 (MQ=255)                                                                                                                                                                                                                                      aCCTTTCTGGCTCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGTTTGATGGCTATCGGCAATGCGCTATTGTTTGCATTTAATCGTCATTATCTGACGCTTATCACCTGTGGTGTGCTTCTGGCATCTCTGGCCAATACGGCAATGCCACAGTTCTTTGCTCTGGCGCGGGAat         >  2:338386/1‑251 (MQ=255)                                                                                                                                                                                                                                        cTTTCTGGATCGGACTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGTTTGATGGCTATCGGCAATGCGCTATTGTTTGCATTTAATCGTCATTATCTGACGCTTATCACCTGTGGTGTGCTTCTGGCATCTCTGGCCAATACGGCAATGCCACAGTTATTTGCTCTGGCGCGGGAatat       >  1:212051/1‑251 (MQ=255)                                                                                                                                                                                                                                              tgATCGGACTTTTTTAGACGGTGAATGATATTGCTGGGATCGGCGTAAGCCTCTGGTGGGCAAAACGTTCTGACATGCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGTTTGATGGCTATCGGCAATGCGCTATTGGTTTCATGTAATCGTCATTATCTGACGCTTAGCACCTTTGGTGTGCTGCTGGCATCTCTGGCCAATACGGCAATGCCACAGTTATTTGCTCTGGCGCGGGAATATGCGGa  <  2:680313/249‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                   ATGTGTGACTGAGTATTGGTGTAAAATCACCCGCCAGCAGATTATACCTGCTGGTTTTTTTTATTCTCGCCGCGCTAAAAAGGGAACGTATGATCTGGATAATGACGATGGCTCGCCGTATGAACGGTGTTTACGCGGCATTTATGCTGGTCGCTTTTATGATGGGGGTGGCCGGGGCGCTACAGGCTCCTACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGCCTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTCAGGGCGATCGGCGAAAACTGATTATATTTTGCTGTTTGATGGCTATCGGCAATGCGCTATTGTTTGCATTTAATCGTCATTATCTGACGCTTATCACCTGTGGTGTGCTTCTGGCATCTCTGGCCAATACGGCAATGCCACAGTTATTTGCTCTGGCGCGGGAATATGCGGA  >  NC_000913/77532‑78018

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