Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,716,988	T→G	14.8%	V148G (GTG→GGG)	yiaD →	multicopy suppressor of bamB, outer membrane lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,716,988	0	T	G	14.8%	174.1 / 2.5	68	V148G (GTG→GGG)	yiaD	multicopy suppressor of bamB, outer membrane lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (13/43);  new base G (10/2);  total (23/45)
Fisher's exact test for biased strand distribution p-value = 1.64e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGGCGCAGCTCTGGGCGGCGGCGTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCGGGGATAACATTATCCTCAATATGCCGAACAATGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAACGTGA‑TTGGTTATACCG‑ACAGCACGGGTGGTCACGACCT‑GAACATGCGTCTCTCCCAGCAAC‑GTGC‑GGA‑TTCCGTTGCCAGCGCG‑TTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTA‑CTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCAC‑CGCAGAAGGTAAGGCGCAAAACCGCCGTGTAGAAATTACCTTAAGCCCGCTGTAATCCCTTT  >  NC_000913/3716752‑3717212                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                          agGCGCAGCTCTGGGCGGCGGCGTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCGTGGATAACATTATCCTCAATATGCCGAACAATGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAACGTg                                                                                                                                                                                                                                          <  1:1101427/238‑1 (MQ=255) agGCGCAGCTCTGGGCGGCGGCGTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCGGGGATAACATTATCCTCAATATGCCGAACAATGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAACGTg                                                                                                                                                                                                                                          <  2:452382/238‑1 (MQ=255)  ggCGCAGCTCTGGGCGGCGGCGTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCGGGGATAACATTATCCTCAATATGCCGAACAATGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTCACGGGag                                                                                                                                                                                                                                        >  2:309904/1‑238 (MQ=255)    cgcAGCTCTGGGCGGCGGCGTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCGGGGATAACATTATCCTCAATATGCCGAACAATGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAACGTGA‑t                                                                                                                                                                                                                                       <  1:811065/237‑1 (MQ=255)              ggcggcggcGTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCGGGGATAACATTATCCTCAATATGCCGAACAATGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGGTAAAGTGG‑TTTGGTATAACg                                                                                                                                                                                                                            >  2:608602/1‑238 (MQ=255)               gcggcggcgTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCGGGGATAACATTATCCTCAATATGCCGAACAATGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAACGTGA‑TTGGTTATACCG‑a                                                                                                                                                                                                                          <  2:821399/238‑1 (MQ=255)                   cggcgTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCGGGGATAACATTATCCTCAATATGCCGAACAATGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAAAGTGG‑TTGGGTTTACCG‑AAAgaa                                                                                                                                                                                                                     >  2:415901/1‑237 (MQ=255)                        ttGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCGGGGATAACATTATCCTCAATATGCCGAACAATGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGGTAAAGGGG‑TTGGGTATAACG‑ACAGACCgggg                                                                                                                                                                                                                >  2:977906/1‑238 (MQ=255)                               ttaCATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCGGGGATAACATTATCCTCAATATGCCGAACAATGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAACGTGA‑TTGGTTATACCG‑ACAGCACGGGTGGTCACg                                                                                                                                                                                                         <  2:315643/239‑1 (MQ=255)                                 aCATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCGGGGATAACATTATCCTCAATATGCCGAACAATGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAACGTGA‑TTGGTTATACCG‑ACAGCACGGGTGGTCACGAc                                                                                                                                                                                                       <  2:419840/239‑1 (MQ=255)                                     ggATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCGGGGATAACATTATCCTCAATATGCCGAACAATGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAACGTGA‑TTGGTTTTAACG‑ACAGCACgggggggtcccacctga                                                                                                                                                                                                   >  2:178485/1‑228 (MQ=255)                                       aTGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCGGGGATAACATTATCCTCAATATGCCGAACAATGTGACCTTCGACAGCACCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAACGTGG‑TTGGTTTTACCC‑ACCGCACGGGGGGGCCCGACCC‑Gaa                                                                                                                                                                                                 >  2:105904/1‑238 (MQ=255)                                         gtgCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCGGGGATAACATTATCCTCAATATGCCGAACAATGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGCTTACGTGA‑TTTGTTTTAACC‑CAAGAACCGGGGGGCAACACCT‑GAACac                                                                                                                                                                                              >  2:1141026/1‑238 (MQ=255)                                            cAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCGGGGATAACATTATCCTCAATATGCCGAACAATGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTTACCGGA‑TTTGGTATAATC‑ACAGCACGGGCGGTCGCCAcccttacaagca                                                                                                                                                                                            >  2:768448/1‑229 (MQ=255)                                            cAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCGGGGATAACATTATCCTCAATATGCCGAACAATGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTGAAGGGG‑GTGGGTGTACCC‑ACAGCACCGGCGGGCACgcaccggacaaag                                                                                                                                                                                             >  1:190236/1‑226 (MQ=255)                                                 agcgaagcTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCGGGGATAACATTATCCTCAATATGCCGAACAATGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGGATCCGAAACCGGCGGTTGACGTGG‑CTGGGTTTTCCG‑CCACCCCGGGGGGGGACGACCT‑GGACATGGGcctc                                                                                                                                                                                       >  1:482780/1‑238 (MQ=255)                                                                cAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCGGGGATAACATTATCCTCAATATGCCGAACAATGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAACGTGA‑TTGGTTATACCG‑ACAGCACGGGTGGTCACGACCT‑GAACATGCGTCTCTCCCAGCAAC‑GTGC‑gg                                                                                                                                                                     <  2:814418/239‑1 (MQ=255)                                                                             aCTGGTGTTAGCGTAACCCGCAGCGGGGATAACATTATCCTCAATATGCCGAACAATGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAACGTGA‑TTGGTTATACCG‑ACAGCACGGGTGGTCACGACCT‑GAACATGCGTCTCTCCCAGCAAC‑GTGC‑GGA‑TTCCGTTGCCa                                                                                                                                                        <  1:1099226/238‑1 (MQ=255)                                                                               tGGTGTTAGCGTAACCCGCAGCGGGGATAACATTATCCTCAATATGCCGAACAATGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAACGTGA‑TTGGTTATACCG‑ACAGCACGGGTGGTCACGACCT‑GAACATGCGTCTCTCCCAGCAAC‑GTGC‑GGA‑TTCCGTTGCCAgcg                                                                                                                                                     <  2:605127/239‑1 (MQ=255)                                                                               tGGTGTTAGCGTAACCCGCAGCGGGGATAACATTATCCTCAATATGCCGAACAATGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGGTAACGGGA‑TTGGTTATAACC‑ACAACACGGGGGGGCACGAACA‑GAACAAGCGCTTCTCCCCACCAC‑GTgacgagtgccgtgacccg                                                                                                                                                        >  1:245383/1‑221 (MQ=255)                                                                                      aGCGTAACCCGCAGCGGGGATAACATTATCCTCAATATGCCGAACAATGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAACGTGA‑TTGGTTATACCG‑ACAGCACGGGTGGTCACGACCT‑GAACATGCGTCTCTCCCAGCAAC‑GTGC‑GGA‑TTCCGTTGCCAGCGCG‑TTGAt                                                                                                                                             <  2:907013/239‑1 (MQ=255)                                                                                      aGCGTAACCCGCAGCGGGGATAACATTATCCTCAATATGCCGAACAATGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAACCGGG‑TTTGGTATATCC‑ACAGACCGCGGGGGCACGACGT‑GAACAAGCGGCTCTCCCCCCAAC‑AGGC‑ggggtttgttttccgcggggtggt                                                                                                                                               >  1:916350/1‑217 (MQ=255)                                                                                                  aGCGGGGATAACATTATCCTCAATATGCCGAACAATGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAACGTGA‑TTGGTTATACCG‑ACAGCACGGGTGGTCACGACCT‑GAACATGCGTCTCTCCCAGCAAC‑GTGC‑GGA‑TTCCGTTGCCAGCGCG‑TTGATCACCCAGGGCg                                                                                                                                  <  2:853370/238‑1 (MQ=255)                                                                                                     ggggATAACATTATCCTCAATATGCCGAACAATGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAACGTGA‑TTGGTTATACCG‑ACAGCACGGGTGGTCACGACCT‑GAACATGCGTCTCTCCCAGCAAC‑GTGC‑GGA‑TTCCGTTGCCAGCGCG‑TTGATCACCCAGGGCGTg                                                                                                                                <  1:487517/237‑1 (MQ=255)                                                                                                            aCATTATCCTCAATATGCCGAACAATGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAACGTGA‑TTGGTTATACCG‑ACAGCACGGGTGGTCACGACCT‑GAACATGCGTCTCTCCCAGCAAC‑GTGC‑GGA‑TTCCGTTGCCAGCGCG‑TTGATCACCCAGGGCGTGGACGCCAGc                                                                                                                       <  2:277241/239‑1 (MQ=255)                                                                                                             cATTATCCTCAATATGCCGAACAATGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAACGTGA‑TTGGTTATACCG‑ACAGCACGGGTGGTCACGACCT‑GAACATGCGTCTCTCCCAGCAAC‑GTGC‑GGA‑TTCCGTTGCCAGCGCG‑TTGATCACCCAGGGCGTGGACGCCAGc                                                                                                                       <  2:594652/238‑1 (MQ=255)                                                                                                              aTTATCCTCAATATGCCGAACAATGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACCCCCTGACCGGCGTGGCAATGGTCCTGAAAGAGTATCCGAAAACGGCGGTTAACGTGA‑TTGGTTATACCG‑ACAGCACGGGTGGTCACGACCT‑GAACATGCGTCTCTCCCAGCAAC‑GTGC‑GGA‑TTCCGTTGCCAGCGCG‑TTGATCCCCCAGGGCGTGGACGCCAGCCg                                                                                                                     <  1:1073134/239‑1 (MQ=255)                                                                                                              aTTATCCTCAATATGCCGAACAATGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAACGTGA‑TTGGTTATACCG‑ACAGCACGGGTGGTCACGACCT‑GAACATGCGTCTCTCCCAGCAAC‑GTGC‑GGA‑TTCCGTTGCCAGCGCG‑TTGATCACCCAGGGCGTGGACGCCAGCCg                                                                                                                     <  2:272083/239‑1 (MQ=255)                                                                                                                 aTCCTCAATATGCCGAACAATGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAACGTGA‑TTGGTTATACCG‑ACAGCACGGGTGGTCACGACCT‑GAACATGCGTCTCTCCCAGCAAC‑GTGC‑GGA‑TTCCGTTGCCAGCGCG‑TTGATCGCCCAGGGCGTGGACGCCAGCCGCAt                                                                                                                  <  2:778987/239‑1 (MQ=255)                                                                                                                 aTCCTCAATATGCCGAACAATGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGGTAACGGGG‑TTTGGTTTAACG‑ACAACACGGGGGGGCACCACCTGGAACAAG‑GTCTCTCCCACCCAC‑GGGC‑GGG‑ATCCGTTTCCCGCGCG‑GTGATCACCCCGGGCGTGGACGCCAGCCGcc                                                                                                                   >  1:272083/1‑237 (MQ=255)                                                                                                                    cTCAATATGCCGAACAACGTGACCTTCGACAGCAGCAGGGGGACCCTGAAACGGGCGGGCGCGAACTCCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAACGTGA‑TTGGTTATACCG‑ACAGCACGGGTGGTCACGACCT‑GAACATGCGTCTCTCCCAGCAAC‑GTGC‑GGA‑TTCCGTTGCCAGCGCG‑TTGATCACCCAGGGCGTGGACGCCAGCCGCATCCg                                                                                                               <  1:309904/239‑1 (MQ=255)                                                                                                                     tCAATATGCCGAACAATGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTTAAGTGG‑TTGGTTATACCG‑ACAGCACGGGGGGGCACGACCC‑GAACAAGCGGCTCTCCCCACCAC‑CACC‑GGAGTTCCTTCGCC‑GCGCG‑GTGGTCGCCCCCGGGGTGGGCGCCACCCGCAcccc                                                                                                               >  1:569634/1‑237 (MQ=255)                                                                                                                       aaTATGCCGAACAATGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAACCGGG‑TTGGGTTTACCC‑ACAAGACGGGGGGGCACGAACT‑CAAACAGCGGCTCTCCCCACAAA‑ACTG‑GGG‑TTGTGCTGTCAGCGCG‑GTGATTACCCCCGGCGGGGACGCCcgccgcgaccctg                                                                                                             >  1:594652/1‑231 (MQ=255)                                                                                                                         tatGCCGAACAATGTGCCCGTCGACAGCGGCAGCGCGAACCTGAAACCGGCGGTCGCTAACACCCGGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAACGTGA‑TTGGTTATACCG‑ACAGCACGGGTGGTCACGCCCT‑GAACATGCGTCTCTCCCAGCAAC‑GTGC‑GGA‑TTCCGTTGCCAGCGCG‑TTGATCACCCAGGGCGTGGACGCCAGCCGCAGCCGTA‑CTc                                                                                                         <  1:304267/239‑1 (MQ=255)                                                                                                                          atGCCGAACAATGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAACGTGA‑TTGGTTATACCG‑ACAGCACGGGTGGTCACGACCT‑GAACATGCGTCTCTCCCAGCAAC‑GTGC‑GGA‑TTCCGTTGCCAGCGCG‑TTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTA‑CTCa                                                                                                        <  2:482780/239‑1 (MQ=255)                                                                                                                           tGCCGAACAATGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAACGTGA‑TTGGTTATACCG‑ACAGCACGGGTGGTCACGACCT‑GAACATGCGTCTCTCCCAGCAAC‑GTGC‑GGA‑TTCCGTTGCCAGCGCG‑TTGATCACCCAGGGCGTGGACGCCAGCCg                                                                                                                     <  2:1032019/226‑1 (MQ=255)                                                                                                                           tGCCGAACAATGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAACGTGA‑TTGGTTATACCG‑ACAGCACGGGTGGTCACGACCT‑GAACATGCGTCTCTCCCAGCAAC‑GTGC‑GGA‑TTCCGTTGCCAGCGCG‑TTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTA‑Ct                                                                                                          <  2:869837/236‑1 (MQ=255)                                                                                                                            gCCGAACAATGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAACGTGA‑TTGGTTATACCG‑ACAGCACGGGTGGTCa                                                                                                                                                                                                           <  1:18758/144‑1 (MQ=255)                                                                                                                            gCCGAACAATGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAAAGGGA‑TTGGTTATACCG‑ACAGCACGGGTGGTCa                                                                                                                                                                                                           >  2:18758/1‑144 (MQ=255)                                                                                                                            gCCGAACAATGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAAACTGAGTTGGTT‑TATCGCACACAACAGGGGGGTACCACCT‑‑AACACGCGTCTCTCCCCACACC‑ATGCGGGC‑TTCCGTTG‑CACCCCGCTTGAACACCAAGGGGGGGGACCCCACCCCCACCCgacccagg                                                                                                        >  1:1188547/1‑232 (MQ=255)                                                                                                                                   aaTGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAACGTGA‑TTGGTTATACCG‑ACAGCACGGGTGGTCACGACCT‑GAACATGCGTCTCTCCCAGCAAC‑GTGC‑GGA‑TTCCGTTGCCAGCGCG‑TTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTA‑CTCAGGGCCTTgg                                                                                               <  2:190236/239‑1 (MQ=255)                                                                                                                                     tgtgACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGGTAACGTGG‑TTGGGTATTACG‑ACCGCAAGGGGGGGCCAGGCCC‑TAAAAAGCGTCTCTCCCCCCAAC‑AGGC‑GGG‑GTCCGCTGCCCCCGCG‑TGGGTCCACCCGGGGGGGGGGGCCCGCCGCAccccgacccacgggcctc                                                                                                 >  1:869837/1‑221 (MQ=255)                                                                                                                                                        agcGCGACCCTGAAACCGGCGGGCGCTAACCCCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAACGTGA‑TTGGTTATACCG‑ACAGCACGGGTGGTCACGACCT‑GAACATGCGTCTCTCCCAGCAAC‑GTGC‑GGA‑TTCCGTTGCCAGCGCG‑TTGATCCCCCAGGGCGTGGACGCCAGCCGCATCCGTA‑CTCAGGGCCTTGGCCCGGCTAACCCAATCGcc                                                                            <  2:560113/237‑1 (MQ=255)                                                                                                                                                        agcGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAACGTGA‑TTGGTTATACCG‑ACAGCACGGGTGGTCACGACCT‑GAACATGCGTCTCTCCCAGCAAC‑GTGC‑GGA‑TTCCGTTGCCAGCGCG‑TTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTA‑CTCAGGGCCTTGGCCCGGCTAACCCAATCGcc                                                                            <  2:164792/237‑1 (MQ=255)                                                                                                                                                           gcgACCCTGAAACCGGCGGGCTCTAACACCCTGACCGGCGGGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAACGTGA‑TTGGTTAGACCG‑ACAGCACGGGTGGTCACGACCT‑GAACATGCGTCTCTCCCAGCAAC‑GTGC‑GGA‑TTCCGTTGCCAGCGCG‑TTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTA‑CTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCaa                                                                       <  2:634813/239‑1 (MQ=255)                                                                                                                                                                ccTGAAACCGGCGGGCGCTAACACCCTGGCCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAACGTGA‑TTGGTTATACCG‑ACAGCACGGGTGGTCACGACCT‑GAACATGCGTCTCTCCCAGCAAC‑GTGC‑GGA‑TTCCGTTGCCAGCGCG‑TTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTA‑CTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCa                                                                  <  2:919435/239‑1 (MQ=255)                                                                                                                                                                ccTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAACGTGA‑TTGGTTATACCG‑ACAGCACGGGTGGTCACGACCT‑GAACATGCGTCTCTCCCAGCAAC‑GTGC‑GGA‑TTCCGTTGCCAGCGCG‑TTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTA‑CTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCa                                                                  <  2:574482/239‑1 (MQ=255)                                                                                                                                                                 cTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAACGTGA‑TTGGTTATACCG‑ACAGCACGGGTGGTCACGACCT‑GAACATGCGTCTCTCCCAGCAAC‑GTGC‑GGA‑TTCCGTTGCCAGCGCG‑TTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTA‑CTCAGGGCCTTgg                                                                                               <  2:723547/209‑1 (MQ=255)                                                                                                                                                                                  tAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAACGTGA‑TTGGTTATACCG‑ACAGCACGGGTGGTCACGACCT‑GAACATGCGTCTCTCCCAGCAAC‑GTGC‑GGA‑TTCCGTTGCCAGCGCG‑TTGATCAcc                                                                                                                                         <  2:813871/151‑1 (MQ=255)                                                                                                                                                                                         cTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAACGTGA‑TTGGTTATACCG‑ACAGCACGGGTGGTCACGACCT‑GAACATGCGTCTCTCCCAGCAAC‑GTGC‑GGA‑TTCCGTTGCCAGCGCG‑TTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTA‑CTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCAC‑CGCAGAAGGTAAGGCGCAAAAccg                                        <  2:569634/239‑1 (MQ=255)                                                                                                                                                                                         cTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAACGTGA‑TTGGTTATACCG‑ACAGCACGGGTGGTCACGACCT‑GAACATGCGTCTCTCCCAGCAAC‑GTGC‑GGA‑TTCCGTTGCCAGCGCG‑TTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTA‑CTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCAC‑CGCAGAAGGTAAGGCGCAAAAccg                                        <  2:828676/239‑1 (MQ=255)                                                                                                                                                                                            aCCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAACGTGA‑TTGGGTTTACCG‑ACAGCACCGGGGGTCACGACCC‑GGACATTGGTCTCTCCCAGCAAC‑GTGG‑GGG‑TTCCGTTGCCCGCGCG‑GTGGGCAACCCAGGCGGGGGCGCCCGCCGCCTCCCCACCGCCGGGCCTTGGCCCGGCTACCCCAACCCCCCCCACCCCCCC‑CCCAAAAGGGAAGCGGCAAAAccccc                                      >  1:810061/1‑239 (MQ=255)                                                                                                                                                                                              cGGGGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAACGTGA‑TTGGTTATACCG‑ACAGCACGGGTGGTCACGACCT‑GAACATGCGTCTCTCCCAGCAAC‑GTGC‑GGA‑TTCCGTTGCCAGCGCG‑TTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTA‑CTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCAC‑CGCAGAAGGTAAGGCGCAAAACCGCCgtg                                   <  2:670777/239‑1 (MQ=255)                                                                                                                                                                                                gCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTTAAGGGG‑TTGGTTATACCG‑ACAGCACGGGGGGGCACGAACC‑GAAAATGCGTGTCTCCCAGCAAC‑GTGC‑GGA‑CTTCGTTGGCCGCCCG‑GTGTTTACCCACGGCGGGGGCGCCGACCGCATCCCTA‑CTCATGGCCGTGGCCCGGCTACACCAATCGCCCGCAACAACAC‑CGCCGAAGGCAAGGCGCCAAAACGCCgtgga                                 >  1:1116544/1‑237 (MQ=255)                                                                                                                                                                                                 cGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAACGTGG‑TTTGGTTTACCC‑GCCACACCGGTGGTCACCGCCC‑GTACAAGACTCCTTTCCAGCAAC‑GGGC‑GGG‑TTTCGCTGCCAACGCC‑CTGCTCCCCCAAGGCGCGGGCGGCTGGCCCCTCCCTA‑CCTATGGCCCTGGCCCGGCTAACCCCATCCCCAGCCACACCCC‑CGCCGGAGGCAAGGCACAAAACCGCCCTGCAg                                >  1:164382/1‑239 (MQ=255)                                                                                                                                                                                                  gTGGAAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAACGTGA‑TTGGTTATACCG‑ACAGCACGGGTGGTCACGACCT‑GAACATGCGTCTCTCCCAGCAAC‑GTGC‑GGA‑TTCCGTTGCCAGCGCG‑TTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTA‑CTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCAC‑CGCAGAAGGTAAGGCGCAAAACCGCCGTGTAg                                <  2:637797/238‑1 (MQ=255)                                                                                                                                                                                                      cAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAACGTGA‑TTGGTTATCCCG‑ACAGCACGGGTGGTCACGACCT‑GAACATGCTTCTCTCCCAGCAAC‑GTGC‑GGA‑TTCCGTTGCCAGCGCG‑TTGATCGCCCAGGGCGTGGACGCCAGCCGCATCCGTA‑CTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGAAACAGCAC‑CGCAGAAGGTAAGGCGCAAAACCGCCGTGTAGAAAt                            <  2:749865/238‑1 (MQ=255)                                                                                                                                                                                                      cAATGGGCCTGAAGGAGATTCCGAAAGCGCTGTTAAACGTGA‑TTGGTTCTACCG‑ACGGAACGGGTGGTCACCACCC‑GCACATGCGTCTCTCCCAGCAAC‑GTGC‑GGA‑TTCCGTTGCCACGGCG‑TGGACCACCCAGGGCGTGGACGCCAGCCGCATCCGTA‑CTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACCGCAC‑CGCAAAAGGTAAGGCGCAAAACCGCCGTGTAGAAAt                            <  1:147917/238‑1 (MQ=255)                                                                                                                                                                                                            tACTGAAAGAGTATCCGAAAACGGCGGTTAACGTGA‑TTGGTTATACCG‑ACAGCACGGGTGGTCACGACCT‑GAACATGCGTCTCTCCCAGCAAC‑GTGC‑GGA‑TTCCGTTGCCAGCGCG‑TTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTA‑CTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCAC‑CGCAGAAGGTAAGGCGCAAAAccg                                        <  2:912906/220‑1 (MQ=255)                                                                                                                                                                                                              cTGAAAGAGTATCCGAAAACGGCGGTTAACGTGA‑ATGGGTATACAG‑ACAGACCAGGGGGGTACGACCT‑TAACATTCATCTCTCCTCGGAACAGGGG‑GGA‑TT‑CGCTGCCCGCGAG‑GTGGTCACCCCGGGGGAGGAGGCCCGAAGCTTCCGTC‑ATCCAGGCCTTGGCCCCGCTAAACCCACCGCACCCAAC‑CCACACCCAGACAGTAAGGCGCAAAAACGCCAGGCACAAATTACCttca                    >  1:699403/1‑236 (MQ=255)                                                                                                                                                                                                               tGAAAGAGTATCCGAAAACGGCGGTTAACGTGA‑TTGGGTATAACG‑ACAGCACGGGGGGGCACGACCT‑GAACATGCGGCTCTTCCAACAAA‑GCGC‑GGG‑TTCCGTTGCCAGCGCG‑GTGGTCACCCAAGGCGGGGACGCCAGCCGCATCCGTA‑CCCAAGGCCCTGGCCCCGCTAAACCCATCCCCCGCAACAGCCC‑CCCAGAAGGGCAGGCGGCAAACCGCCCGGTCGAAATTAACTTAAg                   >  1:567865/1‑238 (MQ=255)                                                                                                                                                                                                                gAAAGAGTTTCCGAAAACGGCGGTGAACGTGA‑TTGGTTATCCCG‑CCAGCACGGGTGGTCCCGACCT‑GAACATGCGTCTCTCCCAGCAAC‑GTGC‑GGA‑TTCCGTTGCCAGCGCG‑TTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTC‑CTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCAC‑CGCAGAAGGTAAGGCGCAAAACCGCCGTGTAGAAATTACCTTAAGcc                 <  1:608602/239‑1 (MQ=255)                                                                                                                                                                                                                gAAAGAGTATCCGAAAACGGCGGTTAACGTGA‑TTGGTTATACCG‑ACAGCCCGGGTGGTCACGACCT‑GAGCATGCGTCTCTCCCAGCAAC‑GTGC‑GGA‑TCCCTTTGCCAGCGCG‑TTGATCCCCCAGGGCGTGGACGCCAGCCGCATCCGTA‑CTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCAC‑CGCAGAAGGTAAGGCGCAAAACCGCCGTGTAGAAATTACCTTAAGcc                 <  2:1099875/239‑1 (MQ=255)                                                                                                                                                                                                                gAAAGAGTATCCGAAAACGGCGGTTAACGTGA‑TTGGTTATACCG‑ACAGCACGGGTGGTCACGACCT‑GAACATGCGTCTCTCCCAGCAAC‑GTGC‑GGA‑TTCCGTTGCCAGCGCG‑TTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTA‑CTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCAC‑CGCAGAAGGTAAGGCGCAAAACCGCCGTGTAGAAATTACCTTAAGcc                 <  2:810061/239‑1 (MQ=255)                                                                                                                                                                                                                  aaGAGTATCCGAAAACGGCGGTTAACGTGA‑TTGGTTATACCG‑ACAGCACGGGTGGTCACGACCT‑GAACATGCGTCTCTCCCAGCAAC‑GTGC‑GGA‑TTCCGTTGCCAGCGCG‑TTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTA‑CTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCAC‑CGCAGAAGGTAAGGCGCAAAACCGCCGTGTAGAAATTACCTTAAGCCCg               <  2:207663/239‑1 (MQ=255)                                                                                                                                                                                                                         tCCGAAAACGGCGTTTAACGGGA‑TTGGTAAAGCGG‑ACGGCACGGGTGAGAGAGACCT‑GAACAGACGTCTCGCCCAGCAAC‑GGGC‑GGA‑TTCGGTTGCCAGCGCG‑TTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTA‑CTCAGGGCCTTGGCCCGGCTACCCCAATCGCCAGCAACAGCAC‑CGCAGAAGGTAAGGCGCAAAACCGCCGTGTAGAAATTACCTTAAGCCCGCTGTAAt        <  1:178485/239‑1 (MQ=255)                                                                                                                                                                                                                             aaagCGGCGTGTAACGGGA‑GTGGTTATACCA‑ACAGCACGGGTGGTCACGACCT‑GAACATGCGCCTCTCCCAGCAAC‑GTGC‑GGA‑TTCCGTTGCCGGCGCG‑GTGATCACCCAGGGCGTGGACGCCGGCCGCATCCGTA‑CTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCAC‑CGCAGAAGGTAAGGCGCAAAACCGCCGTGTAGAAATTACCTTAAGCCCGCTGTAATCCCt    <  1:1141026/239‑1 (MQ=255)                                                                                                                                                                                                                               aaCGGCGGTTAACGTGA‑TTGGTTATCCCG‑ACAGCACGGGTGGTCACGACCT‑GACCATGCGTCTCTCCCAGCAAC‑GTGC‑GGA‑TTCCGTTGCCAGCGCG‑TTGATCCCCCAGGGCGTGGACGCCAGCCGCATCCGTA‑CTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCAC‑CGCAGAAGGTAAGGCGCAAAACCGCCGTGTAGAAATTACCTTAAGCCCGCTGTAATCCCttt  <  2:887380/239‑1 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                          AGGCGCAGCTCTGGGCGGCGGCGTTGGTTATTACATGGATGTGCAGGAAGCGAAGCTGCGCGACAAAATGCGCGGCACTGGTGTTAGCGTAACCCGCAGCGGGGATAACATTATCCTCAATATGCCGAACAATGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAACGTGA‑TTGGTTATACCG‑ACAGCACGGGTGGTCACGACCT‑GAACATGCGTCTCTCCCAGCAAC‑GTGC‑GGA‑TTCCGTTGCCAGCGCG‑TTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTA‑CTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCAC‑CGCAGAAGGTAAGGCGCAAAACCGCCGTGTAGAAATTACCTTAAGCCCGCTGTAATCCCTTT  >  NC_000913/3716752‑3717212

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