Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	25,677	G→A	L157L (TTG→TTA)	lspA →	lipoprotein signal peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	25,677	0	G	A	95.2%	116.0 / NA	44	L157L (TTG→TTA)	lspA	lipoprotein signal peptidase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  major base A (24/16);  minor base C (0/2);  total (24/20)
Fisher's exact test for biased strand distribution p-value = 1.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATA‑TCGCTTACGCGCTGATT‑ATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTGCCTTCTAGAGCGAAAAAACAATAATAAACCCTGCCGGA‑TGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTAAAACTCGACGATGGCACCACCGCCGAGTCTACCCGCAACAACGGT  >  NC_000913/25436‑25921                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                      ttGGTATTAGCGTGATCCTGGCAGTGCTGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATA‑TCGCTTACGCGCTGATT‑ATTGGCGGCGAGCTAGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCAGGCACTTCGCCACCTTCAACCTTTCCGATACTGCCATCGTTGTCGGTGCGGCAATGATTGTGCTTGTAGGTTTTTTACCttttaa                                                                                                                                                                                                                                               >  2:90449/1‑246 (MQ=255)     tATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATA‑TCGCTTACGCGCTGATT‑ATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGGTATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGGTTTTTACCTTCTAAAGCgg                                                                                                                                                                                                                                          >  2:235443/1‑250 (MQ=255)                 ccTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATA‑TCGCTTACGCGCTGATT‑ATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAGGCGAAAAAACaataat                                                                                                                                                                                                                              >  1:277015/1‑251 (MQ=255)                   tGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATA‑TCGCTTACGCGCTGATT‑ATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAataata                                                                                                                                                                                                                             >  1:41282/1‑250 (MQ=255)                     gCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATA‑TCGCTTACGCGCTGATT‑ATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAAc                                                                                                                                                                                                                          >  1:464882/1‑251 (MQ=255)                        gggttgtgggatctcgcagaGGCCACGCAGAAGC‑AAACAAAAGTCGCTT‑CGGGGTGATTTATGGGCGGCGCGCGGGACAACCTGTCCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACGTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTTCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCt                                                                                                                                                                                                                       <  2:464882/232‑1 (MQ=255)                        gTGAGGATGTTTCGCTCGAAGGCCCGGAAGAAGTTAAACAATA‑GCGTTTACGCGCTAATT‑ATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCt                                                                                                                                                                                                                       <  1:396011/251‑1 (MQ=255)                                gtaccGCTCGAGGGCACCGCAGAAGCAAAACAATA‑TCTCTTACGCGCTGATT‑ATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGga‑tg                                                                                                                                                                                                              <  1:155693/247‑1 (MQ=255)                                        cGAAGGCCACGCAGAAGCTAAACAATA‑TAGCTTACGCGCTGATT‑ATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACATTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACATTAATAAACCCTGCCggctgagctgctg                                                                                                                                                                                                       >  1:313479/1‑240 (MQ=255)                                                     gAAGCGAAAAAAAA‑TCGCTGACGCGGTGATT‑TTCGGCGGCGGGTTGGGCAGCCTGTAGGACCGCCTGTGGCACGGCTTCGTTGTCGATATGAGCGACTTCTACGTCGGGGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGA‑TGCGATGCTGACGCATCTTATcc                                                                                                                                                                                         <  1:367503/251‑1 (MQ=255)                                                         cTAAACAATA‑TCGCTTACGCGCTGATT‑ATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGA‑TGCGATGCTGACGCATCTTATCCGGcc                                                                                                                                                                                     >  1:356349/1‑251 (MQ=255)                                                         cTAAACAATA‑TCGCTTACGCGCTGATT‑ATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGAACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGGTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGA‑TGCGATGCTGCAGCATCTTATCCGGcc                                                                                                                                                                                     >  1:70507/1‑251 (MQ=255)                                                                 atagcgtttACGCGTTGATT‑ATGGGCTGGGCGCGTGGCAACCTTTTCAGCCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTGCCTTCTAGAGCGAAAAAACAATAATAAACCCTGCCGGA‑TGCGATGCTGACGCATCTTATCCGGCCTACAGAt                                                                                                                                                                              <  1:94457/244‑1 (MQ=255)                                                                                                      gCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGA‑TGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGCAATCGTAGGCCGGATAAGGCGTTTAcgc                                                                                                                                          >  1:346049/1‑251 (MQ=255)                                                                                                        atcgtgctcgcccGCCTGTGGCTCGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGA‑TGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTAcgccg                                                                                                                                        <  2:363819/240‑1 (MQ=255)                                                                                                           cTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGA‑TGCGCTGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCAt                                                                                                                                     >  2:500353/1‑251 (MQ=255)                                                                                                             gTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGA‑TGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATcc                                                                                                                                   <  1:205020/251‑1 (MQ=255)                                                                                                                      cctggtgccccggttCGTTGGCGATTTGATTGTCTTCTGCGTCGGTGGCTGGCTATGCGCCACCTTCAACCTTGCCGTTACTGCCATCTGTGTCGGTGCGGGAATGATTGTGCTGGAAGGTTTTTTTCCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGA‑TGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAAGCGTAGTCCGTATAAGGCGTTTACGCCGCATCCGGCAAAAAt                                                                                                                          <  2:490973/238‑1 (MQ=255)                                                                                                                       cTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGA‑TGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCAGGATAAGGCGTTTACACCGCATCCGGCAAAAATc                                                                                                                         >  1:237425/1‑251 (MQ=255)                                                                                                                         gtgGCACGGCTTCGTTGTCGATATGAGCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGA‑TGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATc                                                                                                                         <  1:406612/249‑1 (MQ=255)                                                                                                                                     cGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGA‑TGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAagag                                                                                                           <  2:437068/251‑1 (MQ=255)                                                                                                                                         gTCGATATTATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGA‑TGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCaaa                                                                                                       <  1:378299/251‑1 (MQ=255)                                                                                                                                            gttaggagcGACTTCTACGGCGGCGGCTGGCACTTCGCCACCTTCAACCTGGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGA‑TGCGATGCTGACGCATCTTAGCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAACCCt                                                                                                    <  2:99220/243‑1 (MQ=255)                                                                                                                                              tatGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGA‑TGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGCAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCATTAAATATAAGAGAAAACCTGc                                                                                                  >  2:429861/1‑251 (MQ=255)                                                                                                                                               atGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGA‑TGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGc                                                                                                  <  1:234391/250‑1 (MQ=255)                                                                                                                                               atGATCGAATTCTAAGACGGCGGCTGGAACTTCGCCACCTTCAACCTTGCCGGTACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAACAACAATAATAAACCCTGCCGGA‑AGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCAGTCACACATCCTTAAATATAAGAGCAAACCTGcc                                                                                                 >  2:336976/1‑250 (MQ=255)                                                                                                                                               aggatcgctttCTACGTCGGCGGCTTGCACTTCGCCCCCTTCACCCTTGCCGATACTGCCATCTGTGTCGGTGCTGCACTGATTGTGCTGGAAGGTGTTTTCCCTTCTAAAGCGAAAAAACAATAATAAACCCGGCCGGA‑TGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGc                                                                                                  <  2:277015/241‑1 (MQ=255)                                                                                                                                                                cggcggCTGGCACTTCGCCACCTTCAACCTTGCCGATAATGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAGACAATAATAAACCCTGCCGGC‑TTCGATGCTGACGCATCTTATCAGGCCTACCGATTGCTGCGAACTCGTAGGCCAGATAAGGCGTTTACGCCGCCTCCTGCACACAGCCTTCAATATAACCGCCAACCTGCATGTCTGAATCAGTAc                                                                                  >  2:178241/1‑249 (MQ=255)                                                                                                                                                                        ggCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGA‑TGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATa                                                                         <  1:121967/250‑1 (MQ=255)                                                                                                                                                                          cacgtcgccccatttatcctggtcGATACTGCCCTCTGTGTCGGTGCGGCACGGATATTGTGGGAAGTTTTTTTACCTTCTAAAGCGAAAAAACAATAATAACCCCTGCCGGA‑TGCGATGCTGACGCATCTTATCCTGCCTATAGATTGCTGCGATATCGTAGGCCTGATAAGGCGTTTACGCTTCATCCGGCAAAAATCCATAAATATAAGAGCAAACCTGCATGTCTGATTCTGTACAGAGCAATAgcg                                                                      <  2:434684/228‑1 (MQ=255)                                                                                                                                                                            cTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGA‑TGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTATGAATCTGTACAGAGCACTAGCGcc                                                                    >  2:22733/1‑251 (MQ=255)                                                                                                                                                                                gCCAACTTCAACCTTGCCGATACTGCCATCTGTGTCGGGGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGA‑TGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAAAGCAAACCTGCATGTCTTAATCTGTACAGAGCAATCGCGCCGTcc                                                                >  2:264450/1‑251 (MQ=255)                                                                                                                                                                                          aCCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAACACAATAATAAACCCTGCCGGA‑TGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGACTCTGTACAGAGCAATAGCGCCGTCCTGGTGCACtt                                                      >  2:399330/1‑251 (MQ=255)                                                                                                                                                                                                   aTACTGCCATCTGTTTCGGTGCGGCACTGATTGTTCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGA‑TGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTaaa                                             <  2:266671/251‑1 (MQ=255)                                                                                                                                                                                                         ccATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGA‑TGCGATGCTGACGCCTCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTAAACCTcg                                        >  2:157466/1‑250 (MQ=255)                                                                                                                                                                                                           ctctgtgtcggtaagtaaCTGATTGTGCTGTAAGGCTTTTTACCTTACAAAGCGAACAAACCATAATAAACCCTGCCGGA‑TGCGATGCTGAAGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGGTAAGGCGTTTAAGCCGCATCCGGCAAAAATCCTTAGATATAAGTGAAAACCTGCCTGTCTAAATCTGTACAAAGCCCTATCGCCGTACTGGGGCTTTTAACGCTCAAACTcgact                                     >  2:467314/18‑250 (MQ=255)                                                                                                                                                                                                                   cggagCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGA‑TGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGCACAGAGCACAATCGCCGTCCTGGTGAACTGCAAACTACAACTCGACGATGGaacc                             >  2:263186/5‑251 (MQ=255)                                                                                                                                                                                                                       gCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGA‑TGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTGGACGCTaaaaagagaagatggaaccacc                          <  1:222826/250‑19 (MQ=255)                                                                                                                                                                                                                       gCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGA‑TGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAAGCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGGCCTGGTGCACTTCACTCTAAAACTCGACGATGGCACCAccg                         <  2:295408/251‑1 (MQ=255)                                                                                                                                                                                                                          gCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGA‑TGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGAAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTAAAACTCGACGATGGCACCAccgccg                      >  2:224434/1‑251 (MQ=255)                                                                                                                                                                                                                            aCTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGA‑TGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTAAAACTCGACGATGGCACCAccgccg                      >  1:58867/1‑249 (MQ=255)                                                                                                                                                                                                                                aTTGTGTTGGAAGGTTTTTGCCCTTCTAAAGCAAAAAAACAATAATAACCCCGGCCGGAGTG‑GATGCTGACGCATCTTACCCGGCCTACAGATTGCTGCGAAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCAAGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCAAAAACTCGACGAGGGCACCACCGCCGAGTCTa                <  1:277569/251‑1 (MQ=255)                                                                                                                                                                                                                                  tgtgCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGA‑TGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTAAAACTCGACGATGACACCACCGCCGAGTCTAcc              >  1:440145/1‑251 (MQ=255)                                                                                                                                                                                                                                              ttttttACCTTCTAAAGCGAAAAAAAAATAATAAACCCTGCCGGA‑TGCGATGCTGACGCATCTTATACGGCATACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTCAAACTCGACGATTGCACAACCGCCGAGGCTACCAGCAACAACGGt  >  2:342154/1‑251 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                      TTGGTATTAGCGTGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATA‑TCGCTTACGCGCTGATT‑ATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTGCCTTCTAGAGCGAAAAAACAATAATAAACCCTGCCGGA‑TGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTAAAACTCGACGATGGCACCACCGCCGAGTCTACCCGCAACAACGGT  >  NC_000913/25436‑25921

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