Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	25,598	A→G	D131G (GAC→GGC)	lspA →	lipoprotein signal peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	25,598	0	A	G	92.9%	115.9 / 0.6	42	D131G (GAC→GGC)	lspA	lipoprotein signal peptidase
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base G (18/21);  total (20/22)
Fisher's exact test for biased strand distribution p-value = 5.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.38e-01

GTAACTATGGCGCGGCGTTTAGTTTCCTTGCCGATAGCGGCGGCTGGCAGCGTTGGTTCTTTGCCGGTATTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGT‑ATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTGCCTTCTAGAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTAC  >  NC_000913/25361‑25841                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                    gTAACTATGGCGCGGCGTTTAGTTTCCTTGCCGATAGCGGCGGCTGGCAGCGTTGGTTCTTTGCCGGTATTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGT‑ATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGGTATGATCGACTTCTACGTCGGCGGCTGGCACTTCGcc                                                                                                                                                                                                                                        >  2:161514/1‑251 (MQ=255)              ggggttttttttCTTTGCCGATAGCGGCGGCTGGCAGGGTTGGTTCTTTGCCGGTTTGGCGATTGGTATTAGCGTGTTCCTGGCAGGGATGAGTT‑ATCGCTCGAAGGCCACGAAGAAGCTAAAAAATATCGCTTACGCGCTGATTATTGGCGGCGCGCGGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTg                                                                                                                                                                                                                           <  1:9602/242‑1 (MQ=255)                cgTTTAGTTTCCTTGCCGATAGCGGCGGCTGGCAGCGTTGGTTCTTTGCCGGTATTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGT‑ATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGcc                                                                                                                                                                                                                         >  1:5529/1‑251 (MQ=255)                  tttAGTTTCCTGCCCGATAGCGGCGGTTGGCAGCGTTGGTTCTTTGCCGGTATTGCGATTGGTATTGGCGGGATCCTGGCAGTGATGATGT‑ACCGCTCGAAGGCCCCGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCTCCTTCAACCTTGCCg                                                                                                                                                                                                                        <  1:137832/250‑1 (MQ=255)                                       ggcggcTGGCAGCGTTGGTTCTTTGCCGGTATTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGT‑ATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCCCGGCTTCGTTGTCGCTATGATCCACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGCTCCTGCCCTCTGTGTCGGttg                                                                                                                                                                                                  >  2:232411/1‑249 (MQ=255)                                        gggggtgggtgggTTGGTTCTTTGCGGGTAGTGCGTTTGGTATGGGCGTGATCCGGGAAGGGATGTTGT‑ATCGCGCAAAGGCCACGAAGAAGCTAAACAATGTCGTTTCGGCTGTGTTTTTGGGGGGCGCGGGGGGCACCCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCg                                                                                                                                                                                                 <  2:219339/239‑1 (MQ=255)                                          ggcTGGCAGCGTTGGTTCTTTGCCGGTATTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGT‑ATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGc                                                                                                                                                                                               >  1:242035/1‑251 (MQ=255)                                           gtgggaggggttgttttttggCCGGTTTTGCGATGGGTATTAGCGTGATCCTGGCAGTGATGATGT‑ATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGTTTATTGGCGGCGCGCTGGGCACCCTGTCCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCa                                                                                                                                                                                              <  1:375085/231‑1 (MQ=255)                                              ggttgggtgtgtgttttGCCGGTATGGGGATGGGTATTAGCTGGATCCGGGCAGGGAAG‑TGTTATCGTTCGAAGGCACCGCAGAAGCAAAAAAATATCTTTTACGCGCTGATTAGTGGCGGCGCGCGGGGCACCCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTg                                                                                                                                                                                           <  2:1495/237‑1 (MQ=255)                                                     tttgttttttGCCTGTATTGTGATTGGTATTAGGGTGATCCGGGAAGTGATGATGT‑ACCGCTCGAAGCCCACGCAAAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCGGGGCAACCTTTTCTACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATtgtg                                                                                                                                                                                     <  1:467467/243‑1 (MQ=255)                                                     tggttttttttgctgTATTGCGATTTGTATTTGCGTTATCCTGGCAGTGATGATGT‑ATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGTTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTTTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATtgtg                                                                                                                                                                                     <  2:15710/236‑1 (MQ=255)                                                                 cGGTATTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGT‑ATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCGCCCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGttttt                                                                                                                                                                        >  1:89875/1‑251 (MQ=255)                                                                    tATTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGT‑ATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGttttttac                                                                                                                                                                     >  1:165093/1‑249 (MQ=255)                                                                     tttGCGATTGGTATTAGCGTGATCCGGGCAGTGATGGTGT‑ATCGCTGGAAGGCCCCGCAGAAGCTAAACAATATCGCTTACGCGTTGATTATTGGCGGCGCGCTGGGCACCCTTTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGttttttac                                                                                                                                                                     <  1:241834/249‑3 (MQ=255)                                                                           tttGGTATTAGCGTGATCCTGGCAGGGATGATGT‑ATCGCGCGAAGCCCACGCAGAAGCTAAACAATATCGCTTCCGCGCTGATTATTGGCGGCGCGCTGGGCATCCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTg                                                                                                                                                                                           <  1:318097/221‑1 (MQ=255)                                                                           aTTGGTATTAGCGTGATCCTGGCAGTGATGATGT‑ATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGGGCGGCACTg                                                                                                                                                                                           >  2:318097/1‑222 (MQ=255)                                                                                        tGATCCTGGCAGTGATGATGT‑ATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTCCCTTCTAAAGCGAAAAAACa                                                                                                                                                  >  2:406541/1‑250 (MQ=255)                                                                                           ccttggctgggaggatgT‑ATCGCTCGAAGGCCCCGCAGAAGCAAAACAATATCGCTTACGCGCTGTTTAGTGGGGGCGCGCTGGGCACCCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTGCCTTCTAGAGCGAAAAAACaata                                                                                                                                               <  2:242035/237‑1 (MQ=255)                                                                                                                    cGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGATGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGCAGGGTTTTTACCTTCTAAAGGGAAAAAACAATAATAAACCCTGCCGGATTCGATGCTg                                                                                                                      >  2:319913/1‑251 (MQ=255)                                                                                                                     gagggCCACGTATAACCTAATCACTATGGCTTACGCGCTTATTATTGGCGGCTCCCTGCTCACCCTTTTCGGCCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCTCTTCGCCACCTTCAACCGTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGa                                                                                                                     <  2:39271/248‑1 (MQ=255)                                                                                                                               cAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGACAAAACAATAATAAACCCTGCCGGATGCGATGCTTACGCATCTTAt                                                                                                           >  1:14212/1‑251 (MQ=255)                                                                                                                                agaagaTAAACAATATCGCTTACGCGCTGTTTATGGGGGGCGCGCGGGGCACCCTGTCCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATc                                                                                                          <  1:172820/251‑1 (MQ=255)                                                                                                                                 gAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCCCTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAACCAATAATAAACCCTGCCGGATGCGATGCTTACGCATCTTATcc                                                                                                         >  2:30149/1‑251 (MQ=255)                                                                                                                                  aaGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGATGCGCTGCTGACGCATCTTATCCg                                                                                                        >  1:314244/1‑251 (MQ=255)                                                                                                                                          cAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGTTACTGCCATCTGTGTCGGTGCGGCCCTGATTGTGCTGGAAGGTTTTTTACCTTCTTAAGCTAAACAACACTAATAAAACCTGCCGGCTGCGATGCTGAcccctcttttcctgccacc                                                                                                  >  2:113227/1‑231 (MQ=255)                                                                                                                                                                  gggcggggcgggggtaatccTGTTCGGCCGCCTGTTGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTGCCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGa                                                                        <  1:295754/233‑1 (MQ=255)                                                                                                                                                                        cgcgcTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGGGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATTCTGACGCATCTTATCCGGCCTCAAGGTTGCTTCGAAACAGTAGGCCGGATAAGGc                                                                  >  2:84051/1‑251 (MQ=255)                                                                                                                                                                                         cggcccGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGGCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCgg                                                 <  2:376832/247‑1 (MQ=255)                                                                                                                                                                                             ccGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAACAAATAATAAACCCTGACGGATGCGATGCTGACGCATCTTATACGGCCTACAGATTACTGCGAAATCATAGGCCGGCTAAAGCGTTTAAGCCTCCTCCGGcca                                              >  2:366818/1‑248 (MQ=255)                                                                                                                                                                                             ccGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCaaa                                             <  1:161514/251‑1 (MQ=255)                                                                                                                                                                                             ccGCCTGTGGCACGGCTTCGTTGGCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCaaa                                             <  2:89875/251‑1 (MQ=255)                                                                                                                                                                                               gCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGCTGCGATGCTGACGCCTCTTATCCGGCCAACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACCCCGACTCCGGCacaa                                            >  2:410767/1‑247 (MQ=255)                                                                                                                                                                                                   gtgGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCCCTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGGAAAAACAATAATAAACCCTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGTCGTTTACGCCGCATCCGGCCAAAATCCt                                       >  2:132337/1‑251 (MQ=255)                                                                                                                                                                                                   gtgGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCTTTGACGCCGCGGCCGGCAAAAATcc                                        <  1:373976/250‑1 (MQ=255)                                                                                                                                                                                                          ggCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAatat                                <  1:423488/251‑1 (MQ=255)                                                                                                                                                                                                                    gcGATATGATCGACTTCTAAGTCGTCTGCTGGCATTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAAc                      <  1:84051/250‑1 (MQ=255)                                                                                                                                                                                                                     cGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAag                             <  2:105939/243‑1 (MQ=255)                                                                                                                                                                                                                     cGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAag                             >  1:105939/1‑243 (MQ=255)                                                                                                                                                                                                                          tGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCAt                >  1:391053/1‑251 (MQ=255)                                                                                                                                                                                                                                      aCGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGGGCGGCACTGATTGTGCTGGAAGGTTTTTTCCCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCt      <  1:232411/249‑1 (MQ=255)                                                                                                                                                                                                                                        gTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCCAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGCATATGTAc  >  2:139734/1‑251 (MQ=255)                                                                                                                                                                                                                                        gTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTAc  <  2:453811/251‑1 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                    GTAACTATGGCGCGGCGTTTAGTTTCCTTGCCGATAGCGGCGGCTGGCAGCGTTGGTTCTTTGCCGGTATTGCGATTGGTATTAGCGTGATCCTGGCAGTGATGATGT‑ATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTGCCTTCTAGAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTAC  >  NC_000913/25361‑25841

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