Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,158,048	A→C	24.6%	P60P (CCA→CCC)	ptsG →	fused glucose‑specific PTS enzymes: IIB component/IIC component

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,158,048	0	A	C	24.6%	51.4 / 3.4	43	P60P (CCA→CCC)	ptsG	fused glucose‑specific PTS enzymes: IIB component/IIC component
Reads supporting (aligned to +/- strand):  ref base A (13/18);  new base C (10/2);  total (23/20)
Fisher's exact test for biased strand distribution p-value = 1.93e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGAGAACGTAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGTTCCGCGAATTTCAGCTGGCTGCCCGCCGTTGTATCGCATGTTATGGCAGAAGCAGGCGGTTCCGTCTTTGC‑AAACA‑TGCCACTGATTTTTGCGATCGGTGTCGCCCTCGGCTTTACC‑AATAACGATGGCGTATCCGCGCTGG‑CCGCAGTTGTTGCCTATGGCAT‑CATGGTTAAAAC‑CATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAATCGCCTCTAAACACCTGGCGGATACTGGCGTACTCGGAGGGATTATCTCCGGTGCGATCGCAGCGTACATGTTTAACCGTTTCTA  >  NC_000913/1157825‑1158275                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                        tgaaaaaaaaaaaaaaGCCCCCATTCTGAGGAACATTCCAAATTTTTTTTAAAAATTAATTTGTCAACGTGAAAAGGGGGGATAAATCGCTGAGGCGGCGGGTACCCGTCCTGCTTACCGCGGGTTTTCTGCGGGGGGTCGGTTCCGGAAATTTCAGCGGGTTCCCCGCCGTTGTATCGCATGTTATGGCAGAAGCAGGCGGTTCCGTCTTTGC‑AAACA‑TGCCACTGATTTTTGCGAt                                                                                                                                                                                                                           <  2:866685/229‑1 (MQ=255)     aaCGTAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGTTCCGCGAATTTCAGCTGGCTGCCCGCCGTTGTATCGCATGTTATGGCAGAAGCAGGCGGGTCCGCCTTTGC‑AAACC‑TGCCCCCGATTTTTTCGTTCGgtg                                                                                                                                                                                                                      >  1:134557/1‑239 (MQ=255)         tAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGTTCCGCGAATTTCAGCTGGCTGCCCGCCGTTGTATCGCATGTTATGGCAGAAGCAGGCGGTTCCGTCTTTTC‑AAACA‑AGCCACTGATTTTTGCgcttgggtttg                                                                                                                                                                                                                   >  1:740922/1‑228 (MQ=255)          aaaaaaaGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGTTCCGCGAATTTCAGCTGGCTGCCCGCCGTTGTATCGCATGTTATGGCAGAAGCAGGCGGTTCCCGCCTTGCAAAACT‑TGCCACCGATTTTTGCGTTCGgttggg                                                                                                                                                                                                                   >  1:866685/1‑234 (MQ=255)           aaaaaaGCACCCATACTCAGGAGCTCTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAGGTCGGTAAATCGCTGATGCGGCCGGTATCCGTACGGCCTATCGCAGGTATTCGGCTGGGCGTCGGTTCCGCGAATTTCAGCTGGCTGCCCGCCGTTGTATCGCATGTTATGGCAGAAGCAGGCGGTTCCGTCTTTGC‑AAACA‑TGCCACTGATTTTTGCGATCGGTGTCGccc                                                                                                                                                                                                                <  1:552378/239‑1 (MQ=255)                   acatattcTCAGGAGCCCTTTCAATTATTTTTAAGAATGCTTTTGCTAACCGACAAAGGTGTGGTAAATCGCGGATGCTGCCGGTATCCGTGCGGCCTGGCGCAGGTATTCTGCTGGGCGGCGGTTCCGCGAATTTCAGCGGGCTGCCCGCCGTTGTATCGCATGTTATGGCAGAAGCAGGCGGTTCCGTCTTTGC‑AAACA‑TGCCACTGATTTTTGCGATCGGTGTCGCCCTCGGCtt                                                                                                                                                                                                         <  2:213927/231‑1 (MQ=255)                               gAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGTTCCGCGAATTTCAGCTGGCTGCCCGCCGTTGTATCGCATGTTCTGCCAGAAGCCGGCGGGTTCGTCTTTGC‑AAACA‑TACCCATGCTTTTTTGGATCGGTGTGGcccccccgctttcccaaaca                                                                                                                                                                                               >  2:173989/1‑219 (MQ=255)                                  cACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGTTCCGCGAATTTCAGCTGGCTGCCCGCCGTTGTATCGCATGTTATGTCAGAAGCAGGCGGTTCCGTCTTTTC‑AAACA‑TGCACCTGATTTTTGCGGTCGGGGTCGCCCTCGGCTTTTCC‑CATAAcagaggg                                                                                                                                                                                        >  1:746552/1‑233 (MQ=255)                                  cACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGTTCCGCGAATTTCAGCTGGCTGCCCGCCGTTGTATCGCATGTTATGGCAGAAGCAGGCGGTTCCGGCTTTGC‑CAACA‑TGCCACTGATTTTTTCGATCGGTGTCGCCCCCGGCTTTACC‑AATAAccctgtg                                                                                                                                                                                        >  1:213927/1‑233 (MQ=255)                                  cACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGTTCCGCGAATTTCAGCTGGCTGCCCGCCGTTGTATCGCATGTTATGGCAGAAGCAGGCGGTTCCGCCTTTGC‑AAAAA‑TGCCCCTGATTTTTGCGATCGGGGTTGCCCTCGGGTTTACc‑aacaaactaggg                                                                                                                                                                                        >  1:799535/1‑232 (MQ=255)                                  cACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGTTCCGCGAATTTCAGCTGGCTGCCCGCCGTTGTATCGCATGTTATGGCAGAAGCAGGCGGGTCCGCCTTTTC‑AAAAC‑TGCCCCCCCTTTTTGTGGTCGGCGTTGCCCGCCGGTTTaaaaaaaaaaagggg                                                                                                                                                                                         >  2:848056/1‑225 (MQ=255)                                     tctcAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGTTCCGCGAATTTCAGCTGGCTGCCCGCCGTTGTATCGCATGTTATGGCAGAAGCAGGCGGTTCCGCCTTTGC‑AAACC‑TGCCACTGATTTTTTCGGTCGGCGTCGCCCCCGGCTTTACC‑CATACCCATTGCgg                                                                                                                                                                                      >  1:803053/1‑237 (MQ=255)                                         aaTTATGTTTAAGAATGCATTTGCTAACCTGCAAAAGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGTTCCGCGAATTTCAGCTGGCTGCCCGCCGTTGTATCGCATGTTATGGCAGAAGCAGGCGGTTCCGTCTTTTC‑CAAAA‑TGCCCCTGATTTTTTCGATCGGGGTGGCCCCCGGCCTTTCC‑CAAAACACTGGCGTATcc                                                                                                                                                                                  >  1:789545/1‑238 (MQ=255)                                         aaTTATGTTTAAGAATGCATTTGCTAACCTGCAAAAGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGTTCCGCGAATTTCAGCTGGCTGCCCGCCGTTGTATCGCATGTTATGGCAGAAGCAGGCGGTGCCGTCTTTTG‑AAACA‑TGCCACTGGTTTTTTGGACTGGTGTTGCCCTCCGCTTTCGC‑AACAACGACGGCGTAAcc                                                                                                                                                                                  >  1:560864/1‑238 (MQ=255)                                                 ttAAGAATGCATTTGCTAACCTGCAAAAGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGTTCCGCGAATTTCAGCTGGCTGCCCGCCGTTGTATCGCATGTTATGGCAGAAGCAGGCGGTTCCCCTTTTGC‑AAACA‑TGCCACTGATTTTTTGGATCCGTGTCGCCCTCCGCTTTACC‑CAAAAACAAGGCGCATCGCCGCTggg                                                                                                                                                                          >  1:10280/1‑237 (MQ=255)                                                   aaGAATGCATTTGCTAACCTGCAAAAGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGTTCCGCGAATTTCAGCTGGCTGCCCGCCGTTGTATCGCATGTTATGGCAGAAGCAGGCGGTTCCGTCTTTGC‑CAAAA‑TGCCACTGATTTTTTCGGTCGGTGTCGCCCCCGGGCTTTCC‑AATAACgcacgttccttatacacctccc                                                                                                                                                                        >  1:280439/1‑217 (MQ=255)                                                        tggcttggCTCACCTGAAAAAGGGCGGTAAATCGCGGAGCCGGCCGGTTCCGGTACGGCCTACCGCAGTTATTCTGCTGGGGGTCGGTTCCGCGAATTTCAGCGGGTTGCCCGCCGTTGTATCGCATGTTATGGCAGAAGCAGGCGGTTCCGTCTTTGC‑AAACA‑TGCCACTGATTTTTGCGATCGGTGTCGCCCTCGGCTTTACC‑AATAACGATGGCGTATCCGCGCTGG‑CCGCAgttg                                                                                                                                                                 <  1:349002/232‑1 (MQ=255)                                                                               tCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGTTCCGCGAATTTCAGCTGGCTGCCCGCCGTTGTATCGCATGTTATGGCAGACGCAGGCGGTTCCGTCTTTTG‑AACCA‑TGCCCCCTATTTTTTCGATCGGTGTCTCCCCCGGGTTTTCCAAATAAC‑ATGGCGGTTCCCCCCTGG‑CCCGCGTTGTTGCCCATGGcac‑caaaggtaaa                                                                                                                                         >  2:377695/1‑231 (MQ=255)                                                                               tCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGTTCCGCAAATTTCAGCTGGCTGCCCGCCGTTGTATCGCATGTTATGGCAGAAGCAGGCGGTTCCGTCTTTGC‑AAACC‑TGCCACTGATTTTTGCGATCGGTGGCACCCTCCGCTTTAAC‑AAAAACGATTGCGTATCCGCGCTGG‑CCGCAGTTGTTGCCTATGGccacaatggtaaa                                                                                                                                          >  2:344604/1‑227 (MQ=255)                                                                                             aTGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGTTCCGCGAATTTCAGCTGGCTGCCCGCCGTTGTATCGCATGTTATGGCAGAAGCAGGCGGTTTCCTCTTTGC‑AAAAA‑TTCCCCTGGTTTTTTCGGTTGGTGGTGCCCCCGGCTTTTCC‑CATAAACAAGGCGTTTCCCCCCCTG‑CCGCCGTTTTTGCCCTTGGCCT‑CCTGGGTaaaaaccaggcccgggg                                                                                                                           >  2:700269/1‑226 (MQ=255)                                                                                                 tgcCGGTCTCCGTTCTGCCTAGCCCAGGATTTCTTCTGGGCGCCGGTTCCGGGAATTTCAGGTGGCTCCCCTCCGTTGTACCGCATGTTATGGCAGAAGCAGGCGGTTCCGTCTTTGC‑AAACA‑TGCCACTGATTTTTGCGATCGGTGTCGCCCTCGGCTTTACC‑AATAACGATGGCGTATCCGCGCTGG‑CCGCAGTTGTTGCCTATGGCAT‑CATGGTTAAAAC‑CATGGCCGTGGTTgc                                                                                                                       <  1:464256/238‑1 (MQ=255)                                                                                                          ccGTTCTGCCTATCGCGTGTATTCTGCTGGGCGTGGGTTCCGCGAATTTCAGCTGGCTGCCCGCCGTTGTATCGCATGTTATGGCAGAAGCAGGCGGTTCCGTCTTTGC‑AAACA‑TGCCACTGATTTTTGCGATCGGTGTCGCCCTCGGCTTTACC‑AATAACGATGGCGTATCCGCGCTGG‑CCGCAGTTGTTGCCTATGGCAT‑CATGGTTAAAAC‑CATGGCCGTGGTTGCGCCACTGGTa                                                                                                             <  1:700269/239‑1 (MQ=255)                                                                                                            gTACTTCCTATCGGAGTTTTTCTGTGGGGCGTCGGTTCCGGGAATTTCAGTTGGCTGCCCGCCGTTGTATCGCATGTTATGGCAGAAGCAGGCGGTTCCGTCTTTGC‑AAACA‑TGCCACTGATTTTTGCGATCGGTGTCGCCCTCGGCTTTACC‑AATAACGATGGCGTATCCGCGCTGG‑CCGCAGTTGTTGCCTATGGCAT‑CATGGTTAAAAC‑CATGGCCGTGGTTGCGCCACTGGTACt                                                                                                           <  1:330596/239‑1 (MQ=255)                                                                                                             tACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGTTCCGCGAATTTCAGCTGGCTGCCCGCCGTTGTATCGCATGTTATGGCAGAAGCAGGCGGGTCCGTCTTTGC‑CAAAA‑TGCCACCGATTTTTTGGATCGGTGGCGCCCCCCGCTTTACC‑CATAAAGATGGCGGATCCGCTGTGG‑CCGCCGGTGTTGCCTTAGGCAA‑CAATGTTTAAAC‑AAAGGCCggggggtggcgacggggacgg                                                                                                          >  2:451864/1‑222 (MQ=255)                                                                                                                  ccTATCGCAGGTATTCTGCTGGGCGTCGGTTCCGCGAATTTCAGCTGGCTGCCCGCCGTTGTATCGCATGTTATGGCAGAAGCAGGCGGTTCCGGCTTTGC‑CAACA‑TGCCACTGATTTTTTCGATCCGCGTTGCCCTCGGCCTTACC‑CATAACGCTGTCCGTTCCCCGCTGG‑CCGCCGGTGTGGCCCAAGGCCT‑TCCGGTTTAAAC‑CCTGGCCCTGGTTGCGCCACTCGTGCCGCCttt                                                                                                     >  2:510563/1‑239 (MQ=255)                                                                                                                     aTCGCAGGTATTCTGCTGGGCGTCGGTTCCGCGAATTTCAGCTGGCTGCCCGCCGTTGTATCGCATGTTATGGCAGAAGCAGGCGGTTCCCGCTTTGC‑AAAAA‑TTCCCCTGATTTTTGGGATTGGGGTCGCCCCCGGCTTTTCC‑AATAACACAGGCGTATCCCCCCCGG‑CCGCAGTGTGTGCCCACGGCCT‑CAGCGGTAAAACACTTGGCCGGGTTGGCCCCAcggtcagtcattttc                                                                                                    >  2:781566/1‑224 (MQ=255)                                                                                                                                 ttgcgggggggctgTTCCGAGAATTGCAGTGGGTTGCCCGCCGTTGTATCGCATGTTAGGGCAGAAGGAGGCGGTTCCGTCTTTGC‑AAACA‑TGCCACTGATTTTTGCGATCGGTGTCGCCCTCGGCTTTACC‑AATAACGATGGCGTATCCGCGCTGG‑CCGCAGTTGTTGCCTATGGCAT‑CATGGTTAAAAC‑CATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAAt                                                                                      <  2:696336/226‑1 (MQ=255)                                                                                                                                         cGTCGGTTCCGCGAATTTCAGCTGGCTGCCCGCCGTTGTATCGCATGTTATGGCAGAAGCAGGCGGTTCCGTCTTTGC‑AAACA‑TGCCACTGGTTTTTTCGATCGGTGTCGCCCTCGGCTTTTCC‑AATAACGATGGCGTATCCGCGCTGG‑CCGCCGTTGTTGCCCATGGCAT‑CATGGTTAAAAA‑CATGGCCGTGGTTGCGCCCCTGGTACTGCATTTACCACCTGAAGAAATCGCctctt                                                                              >  1:14969/1‑238 (MQ=255)                                                                                                                                            cGGTTCCGCGAATTTCAGCTGGCTGCCCGCCGTTGTATCGCATGTTATGGCAGAAGCAGGCGGTTCCGCCTTTGC‑AAACA‑TACCAATGATTTTTTCCATCGGAGTCGCCCCCGCCTTTTCC‑CATAACGCTCGCGTATCCGCCCTGGCCCCCAGTGTTTGCC‑AAGGCAT‑CATGCGTAAAAA‑CACGGCCCGGCGCGCGCCACCGGTACTGCATTTcccctctacagaaactcactctaaac                                                                           >  2:380074/1‑216 (MQ=255)                                                                                                                                                         ttCAGCTGGCTGCCCGCCGTTGTATCGCATGTTATGGCAGAAGCAGGCGGTTCCGGCTTTTG‑AAAAA‑TGCCACTGATTTTTGCGATCGGGGTCGCCCTCGGCTTTACC‑AATAAAGAAGGCCGATCCGCGCCGG‑CCGCCGTTGGTGCCTACTGCATGCAGGTTTAAAAC‑C‑TGGGCGGGGGTGCCCCCCCGGGGACGCATTTTCCTTCGGAAAAAAGCACCCCCAAAAACCCCGCGGAtca                                                              >  2:94206/1‑237 (MQ=255)                                                                                                                                                               ggggtgCCCCCCGTTGTATCGCATGTTATGGCAGAAGCGGGCGGTTCCGTCTTTGC‑AAACA‑TGCCACTGATTTTTGCGAACGGTGGCGCCCTCGGCTTTACC‑AATAACGATGGCGTATCCGCGCTGG‑CCGCAGTTGTTGCCTATGGCAT‑CATGGTTAAAAC‑CATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAATCGCCTCTAAACACCTGGCGGATACTGGCg                                                         <  1:510563/234‑1 (MQ=255)                                                                                                                                                                    gCCCGCCGTTGTTCCGCATTTTATGGCAGAAGCAGGCGGTTCCGTCTTTGC‑AAACA‑TGCCACTGATTTTTGCGATCGGTGTCGCCCTCGGCTTTACC‑AATAACGATGGCGTATCCGCGCTGG‑CCGCAGTTGTTGCCTATGGCAT‑CATGGTTAAAAC‑CATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAATCGCCTCTAAACACCTGGCGGATACTGGCGTACTCg                                                   <  1:393333/239‑1 (MQ=255)                                                                                                                                                                      ccgccgTTGTATCGCATTTTATGGCAGAAGCAGGCGGTTCCGTCTTTGC‑AAACA‑TGCCACTGATTTTTGCGATCGGTGTCGCCCTCGGCTTTACC‑AATAACGATGGCGTATCCGCGCTGG‑CCGCAGTTGTTGCCTATGGCAT‑CATGGTTAAAAC‑CATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAATCGCCTCTAAACACCTGGCGGATACTGGCGTACTCgg                                                  <  1:399586/238‑1 (MQ=255)                                                                                                                                                                                    cATGTTATGGCAGAAGCAGGCGGTTCCGTCTTTGC‑AAACA‑TGCCCCTGATTTTTTCGATCGGTGTCGCCCTCGGCTTTACC‑AATAACGATGGCGGATCCCCGCTGG‑CCGCAGGTGTTGCCTATGGCAT‑CATGGGTAAAAC‑CATGGCCGTGGTTTCGCCCCCTGTACTGCATTTACCTGCTGAAGAAAACGCCTCTTAACACCTTGCGGATACTGGCGGACTCgggggggatttttcccg                                   >  1:705041/1‑228 (MQ=255)                                                                                                                                                                                       gTTATGGCAGAAGCAGGCGGTTCCGTCTTTGC‑AAACA‑TGCCACTGATTTTTGCGATCGGTGTCGCCCTCGGCTTTACC‑AATAACGATGGCGTATCCGCGCTGG‑CCGCAGTTGTTGCCTATGGCAT‑CATGGTTAAAAC‑CATGGCCGTGGTTGCGCCACTGGTACTGCATTTAcc                                                                                                  <  1:531536/173‑1 (MQ=255)                                                                                                                                                                                            ggcgaaagtgggCGGTCCCGTCTTTGC‑AAACA‑TGCCACTGATTTTTGCGATCGGTGTCGCCCTCGGCTTTACC‑AATAACGATGGCGTATCCGCGCTGG‑CCGCAGTTGTTGCCTATGGCAT‑CATGGTTAAAAC‑CATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAATCGCCTCTAAACACCTGGCGGATACTGGCGTACTCGGAGGGATTATCTCCGGTGCGATCg                           <  2:705041/229‑1 (MQ=255)                                                                                                                                                                                             tgcaagaagggggggtctcgtttttC‑AAACACTGCCACGGGTTTGTGCGACCGGTGTGGCCCTCGGTTTTACC‑AATAACGTTGGCGTATCCGCGGTGG‑ACTCAGTTGTTGCCTATGGCAT‑CATGGTTAAAAC‑CATGGCCGGGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAATCGCCTCTAAACGCCTGGCGGATACTGGCGTACTCGGAGGGATTATCTCCGGTGCGATCg                           <  2:708308/218‑1 (MQ=255)                                                                                                                                                                                               caagacgcgCGGTTCCGTCTTTGC‑AAACA‑TGCACCTTATTTTTGCGTGCGGTGGCGCCCTCGGCTTTACC‑AATAACGGTGGCGTATCCGCGCTGG‑CCGCAGTTGTTGCCTATGGCAT‑CATGGTTAAAAC‑CATGGCCGTGGTTGAGCCACTGGTACTGCATTTACCTGCTGAAGAAATCGCCTCTAAACACCTGGCGGATACTGGCGTACTCGGAGGGATTATCTCCGGTGCGATCGCa                         <  2:560864/230‑1 (MQ=255)                                                                                                                                                                                                     aGGCGGTTCCGTCTTTGC‑AAACA‑GGCCACTGATTTTTGCGATCGGTGTCGCCCTCGGCTTTACC‑AATAACGATGGCGTATCCGCGCTGG‑CCGCAGTTGTTGCCTATGGCAT‑CATGGTTAAAAC‑CATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAATCGCCTCTAAACACCTGGCGGATACTGGCGTACTCGGAGGGATTATCTCCGGTGCGATCGCAGCGTa                    <  1:872664/237‑1 (MQ=255)                                                                                                                                                                                                            tCCTTCTTTTG‑AAACA‑TTGCCCTGAGTTTTTGGATCGGGGGCGCCCTCGGGTTTCCC‑ATGAACGATGGCGTGCCCGCGCTGG‑CCGCAGGTGTTGCCTAGGGCAT‑CATGGTTAAAAC‑CAGGGCCGTGGTTGCGCCTCGGGTACTGCATTTACCGGCTGAAGAAATCGCCTCTAAACACCTGGCGGATACTGGCGTACTCGGAGGGATTATCTCCGGTGCGAACGCAGCGTACATGTTTa            <  2:563407/238‑1 (MQ=255)                                                                                                                                                                                                                    tGC‑AAACA‑TGCCATTGATTTTTGGGGTGGGTTGCCCCCTCGGCTTTACC‑AATAACGATGGCGTATCCGCGCTGG‑CCGCAGTTGTTGCCTATGGCAT‑CATGGTTAAAAC‑CATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAATCGCCTCTAAACACCTGGCGGATACTGGCGTACTCGGAGGGATTATCTCCGGTGCGATCGCAGCGTACATGTTTAACCGTTTCt   <  2:803053/239‑1 (MQ=255)                                                                                                                                                                                                                     gC‑AAGCA‑TGCCACGGATTGTTGCGATCGGGGTCGCCCTCGGCTTTACC‑AATAACGATGGCGTACCCGCGCTGG‑CCGCAGTTGTTGCCTATGGCAT‑CATGGTTAAAAC‑CATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAATCGCCTCTAAACACCTGGCGGATACTGGCGTACTCGGAGGGATTATCTCCGGTGCGATCGCAGCGTACATGTTTAACCGTTTCTa  <  1:377695/239‑1 (MQ=255)                                                                                                                                                                                                                     gC‑AAACA‑TGCCACTGATTTTTGCGATCGGTGTCGCCCTCGGCTTTACC‑AATAACGATGGCGTATCCGCGCTGG‑CCGCAGTTGTTGCCTATGGCAT‑CATGGTTAAAAC‑CATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAATCGCCTCTAAACACCTGGCGGATACTGGCGTACTCGGAGGGATTATCTCCGGTGCGATCGCAGCGTACATGTTTAACCGTTTCTa  <  1:94206/239‑1 (MQ=255)                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                        TGAGAACGTAAAAAAAGCACCCATACTCAGGAGCACTCTCAATTATGTTTAAGAATGCATTTGCTAACCTGCAAAAGGTCGGTAAATCGCTGATGCTGCCGGTATCCGTACTGCCTATCGCAGGTATTCTGCTGGGCGTCGGTTCCGCGAATTTCAGCTGGCTGCCCGCCGTTGTATCGCATGTTATGGCAGAAGCAGGCGGTTCCGTCTTTGC‑AAACA‑TGCCACTGATTTTTGCGATCGGTGTCGCCCTCGGCTTTACC‑AATAACGATGGCGTATCCGCGCTGG‑CCGCAGTTGTTGCCTATGGCAT‑CATGGTTAAAAC‑CATGGCCGTGGTTGCGCCACTGGTACTGCATTTACCTGCTGAAGAAATCGCCTCTAAACACCTGGCGGATACTGGCGTACTCGGAGGGATTATCTCCGGTGCGATCGCAGCGTACATGTTTAACCGTTTCTA  >  NC_000913/1157825‑1158275

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