Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,027,233	0	T	A	29.8%	81.1 / 12.2	47	intergenic (+5/+376)	fre/fadA	NAD(P)H‑flavin reductase/3‑ketoacyl‑CoA thiolase
Reads supporting (aligned to +/- strand):  ref base T (24/9);  new base A (0/14);  total (24/23)
Fisher's exact test for biased strand distribution p-value = 2.39e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GGCGCTTTCGTTGAAGCATCCTGGTCTGCAAGTGGTGCCGGTGGTTGAACAACCGGAAGCGGGCTGGCGTGGGCGTACTGGCACCGTGTTAA‑CGGCGGTATTGCAGGATCACGGTACGCTG‑GCAGAGCATGATATCTATATT‑GCCGGACG‑TTTTG‑AGATGGCGAAAATTGCCCGCGATCTGTTTT‑GCAGTGAGCGTAA‑TGCGCGGG‑AAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAAC‑‑CCGCCCCTGACAGGCGGG‑AAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAAC‑CGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAG  >  NC_000913/4026985‑4027473                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                     ggCGCTTTCGTTGAAGCATCCGGGTCTGCAAGTGGTGCCGGTGGTTGAACAACCGGAAGCGGGCTGGCGTGGGCGTACTGGCACCGTGTTAA‑CGGCGGTATTGCAGGATCACGGTACGCTG‑GCAGAGCATGATATCTATATT‑GCCGGACG‑TTTTG‑AGATGGCGAAAATTGCCCGCGATCTGTTTT‑GCAGTGAGCGTAA‑TGCGCGGG‑AAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATaa                                                                                                                                                                                                                                                    <  2:17691/251‑1 (MQ=255)          gTTGAAGCATCCTGGTCTGCAAGTGGTGCCGGTGGTTGAACAACCGGAAGCGGGCTGGCGTGGGCGTACTGGCACCGTGTTAA‑CGGCGGTATTGCAGGATCACGGTACGCTG‑GCAGAGCATGATATCTATATT‑GCCGGACG‑TTTTG‑AGATGGCGAAAATTGCCCGCGATCTGTTTT‑GCAGTGAGCGTAA‑TGCGCGGG‑AAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAAC‑‑CCg                                                                                                                                                                                                                                         <  2:178238/251‑1 (MQ=255)            tGAAGCATCCTGGTCTGCAAGTGGTGCCGGTGGTTGAACAACCGGAAGCGGGCTGGCGTGGGCGTACTGGCACCGTGTTAA‑CGGCGGTATTGCAGGATCACGGGATGCCG‑GGCGAGGCTGATATCTATATT‑TCCGGACGGTTTTG‑‑GAAGGCGAAAAATGCCCGCGATCTGTTTT‑TCAGTGAGCGCAA‑TGAGCGGG‑GAGATCGCCTGTTTGGGGGTGCGTTTGTATTTTTTTGTGATATAAAAAAAc‑‑ccccc                                                                                                                                                                                                                                       >  1:447138/1‑248 (MQ=255)                        gtgtGCAAGTGGTGCCGGTGGTTGAACAACCGGAGGCGGGCGGGCGGGGGCGTACTGGCACCGTGTTAA‑CGGCGGTATTGCAGGATCGCGGTACGCTG‑GCAGAGCATGATCTCTATATT‑GCGGGACT‑TTTTG‑AGATGGCGAAAATTGCCCGCGATCTGTTTG‑CCAGTGAGCGTAA‑TGCGCGGG‑AAGATCGCCTGTTTGGCGATGCTTTTGCATTTATCGGAGATAAAAAAAAAC‑‑CCGCCCCTGACAGGCgg                                                                                                                                                                                                                           <  1:69429/248‑1 (MQ=255)                                        ggtggtTGAACAACCGGAAGCGGGATGGCGTGGGCGTACTGGCACCGTGTTAA‑CGGCGGTATTGCAGGGACAAGGGACGCTG‑GCAGGGCATGAGATTTATATT‑GCCGGACG‑TTTTG‑AGAGGGGGGAAATTGGACGCGGACTGGTTT‑GCAGTGAGGGGAA‑AGGGGGGG‑GGGATTGCCCGTTTTGGGATGGGGTTGGATTTATCTGTGAATTAAAAAAAA‑‑CACCCCCTGACAGGCggggaaaaaaggcaata                                                                                                                                                                                                            >  1:586088/1‑236 (MQ=255)                                               aaacaccCGGAAGCGGGTTGGCGTGGGCGTACTGGC‑CCTGGTTAAACGGCGTTATGTCGGGATCCCGGTCCGCGG‑GCAGAGCATTATATCTATATT‑GCCGGCTT‑TTTTG‑AGATGGCGAAAATTGCCCGCGATCTTTTTG‑GCAGTGAGCGTAA‑TGCGCGGG‑AAGATCGCCTGTTTGGCGATGGTTTTTGTTTTATCTGAAAAAAAAAAAACC‑‑CCGCCCCTGACAGGCGGG‑AAGAACGGCAACTAAACTGtt                                                                                                                                                                                                    <  1:668755/249‑1 (MQ=255)                                                    acgggaggggggggggCGGGGGCGTACGGGCCCCGTGTTAA‑GGGCGGTATGGAAGGATCAGGGTGCGCGG‑GAAGAGCATGATATCTATATG‑GCGGGACG‑TTTTG‑AGAGGGGAAAAATGGCCCGGGATCTGTTTT‑GCGGTGAG‑GTAAAGGGGGGGG‑AAGACCGCCGGTTGGGGGAGGCGTTGGCATTTATCGGAGATAAAAAAAACC‑‑CCGCCCCTGACAGGCGGG‑AAGAACGGCAACTAAACTGTGATTCAg                                                                                                                                                                                              <  2:559994/237‑1 (MQ=255)                                                     ccggaggggggttggggggggCGTCCTGGCCCCGTGTTAA‑CGGCGTTTTTGCAGGATCACGGTACGCGG‑GCAGAGCATGATATTTATATT‑GCCGGACG‑TTTTG‑AGATGGCGAAAATTCCCCGCGATCTTTTTT‑GCAGTGAGCGTAA‑TGCGCGGG‑AAGATCGCCTTTTTGGCGATGCGTTTGCATTTATCTGAAATAAAAAAAAAC‑‑CCGCCCCTGACAGGCGGG‑AAGAACGGCAACTAAACTGTTATTCAGt                                                                                                                                                                                             <  1:136941/233‑1 (MQ=255)                                                      cGGAAGCGGGCTGGCGTGGGCGTACTGGCACCGTGTTAA‑CCGCGGACTTGACCGATTACGGGACGAAGAGC‑GACGATTATATCTAAATTTCCCGGA‑G‑CTTTG‑GGAGAGCGGAAATTGACCGGGATCTGTTTT‑GCAGTGAGCGGAC‑TGCGCGGG‑GGGATCGCCTGTTTGGCGATGGGGTTGGCTTTTTCTGAGGTATAAAAAAAA‑‑CCGCCCCCGTCAGGCGGGGaagaagggaccaaaaatggtatttag                                                                                                                                                                                               >  1:324127/1‑228 (MQ=255)                                                                               gggcccggtgtAAA‑CGGCGGTATTGAGGGATCACGGTGCGCGG‑GCAGAGCATGATATCTATATT‑GCCGGACG‑TTTTG‑AGATGGCAAAAATTGCCCGCGTTGTGTTTT‑GCAGGGAGGTTAA‑TGCGGGGA‑AAGATCGCCTGTTTGGCGATGCGTTTGCTTTTATCTGAGATAAAAAAAACC‑‑CCGCCCCTGACAGGCGGG‑AAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTgg                                                                                                                                                                   <  2:29231/244‑1 (MQ=255)                                                                                               ggcggTATTGCAGGATCACGGTACGCTG‑GCAGAGCATGATATCTATATT‑GCCGGACG‑TTTTG‑AGATGGCGAAAATTGCCCGCGATCTGTTTT‑GCAGTGAGCGTAA‑TGCGCGGG‑AAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAAC‑‑CCGCCCCTGACAGGAGGG‑AAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGAt                                                                                                                                                    >  1:480468/1‑251 (MQ=255)                                                                                                       tGCAGGATCACGGTACGCTG‑GCAGCGCATGATATCTATATT‑GCCGGACG‑TTTTG‑AGATGGCGCAACTTGCCCGCGATCTGTTTT‑GCAGTGAGAGTAA‑TGCGCGGG‑CAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAAC‑‑CCGCCCCTGACAGGCGGG‑AAGAACGGCAACTAAACTGATATTCAGAGGCATTTTGATCTTTGACGTATCTGGCAAACCCTGCCCGAGGCGAccct                                                                                                                                            >  2:695463/1‑248 (MQ=255)                                                                                                                   gTACGCGG‑GCAGAGAATGATATCTATATT‑GCGGGACG‑TTTTG‑AGAGGGCGAAAATGGCCCGCATTCTTTTTT‑GCAGGGAGCGTAA‑TGCGCGGA‑AAGACCGCCTGTTGGGCGAGGCGTTTGCTTTTATCTGAGATATAAAAAAAC‑‑CCGCCCCTGACAGGCGGG‑AAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTAt                                                                                                                                <  1:393246/251‑1 (MQ=255)                                                                                                                      ggCTG‑GCAGAGCATGATTGCTATAGT‑GCGGGACG‑TTTTG‑AGATGGCGAAAATTGCCCGCGATTTGTTTT‑GCAGTGAGCGAAA‑TGCGGGGG‑AAGATCGCCTGTTTGGCGAGGCGTTTGCATTTATCTGAAATATAAAAAAAC‑‑CCGCCCCTGACAGGCGGG‑AAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGt                                                                                                                             <  1:292677/250‑1 (MQ=255)                                                                                                                        gtG‑GCAGAGCATGATTTCTATATT‑GCCGGACG‑TTTTG‑AGAGGGCAAAAATTGCCCGCGATCTGTTTT‑GCGGTGAGCGAAA‑TGCGGGGG‑AAGATCGCCTGTTGGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAAC‑‑CCGCCCCTGACAGGCGGG‑AAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGc                                                                                                                           <  1:186537/250‑1 (MQ=255)                                                                                                                        gcgggaaaagCATGATATCTATATG‑GCCGGACG‑TTGGG‑AGAGGGGAAAAATTGCCCGCGATCTGTTTG‑GGAGTGAGCGTAA‑TGCGCGAA‑AAGACCCCCTTTTGGGCGTTGCGTTTGCATTTATTGGAGATAAAAAAAAAC‑‑CCGCCCCTGACAGGCGGG‑AAGAACGGCAACTAAACTGTTAGTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTg                                                                                                                            <  1:117889/243‑1 (MQ=255)                                                                                                                        cTG‑GCAGAGCATGATATCTATATT‑GCCGGACG‑TTTTG‑AGATGGCGAAAATTGCCCGCGATCTGTTTT‑GCAGTGAGCGTAA‑TGCGCGGG‑AAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAAC‑‑CCGCCCCTGACAGGCGGG‑AAGAACGGCAACTAAACTGTTATTCAGTGGCCTTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTATTATCGTGc                                                                                                                           >  1:97679/1‑251 (MQ=255)                                                                                                                             cAGAGCATGATATCTATATT‑GCGGCACG‑TTTTG‑AGATGGCTAAAATTCCCCGCTATCTTTTTT‑GCAGTGGGCGGAA‑TTCGGGGG‑AAGATCCCCTGTGTGGCGAGGCGTTGGCATTTATCTGAGATATAAAAAACC‑‑CCGCCCCTGACAGGCGGG‑AAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTAc                                                                                                                       <  1:356501/251‑1 (MQ=255)                                                                                                                                agCATGATATCTATTTG‑GCCGGACG‑TTTGG‑AGATGGCGAAAATGGCCCTCGTTTTGTGTT‑GAAGTGAGCGTAA‑GGCGGGGG‑AAGAGCCCCTTTTTGGCGATGCTTTGGCATTAATCTGAGATATAAAAAAAC‑‑CCGCCCCTGACAGGCGGG‑AAGAACGGCAGCTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTcgc                                                                                                                                       <  2:483414/232‑1 (MQ=255)                                                                                                                                agCATGATATCTATATT‑GCCGGACG‑TTTTG‑AGATGGCGCAAATTGCCCGCGATCTGTTTT‑GCAGTGAGCGTAA‑TGCGCGGG‑AAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAAC‑‑CCGCCCCTGACAGGCGGG‑AAGAACGGCAGCTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTcgcagatcggcagagcacacgc                                                                                                                    >  1:483414/1‑232 (MQ=255)                                                                                                                                     ataatttctaatttgcgggag‑TTTTGGAGATGGCAAAAATGCCCCGCTATCTTTTTT‑GCAGGGAGCGAAA‑TGCGGGGG‑AAGACCGCCTGTTGGGCGAGGCGTTGGAATTTATCTGAGATAAAAAAAAAC‑‑CCGCCCCTGACAGGCGGG‑AAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAg                                                                                                                 <  2:678190/231‑1 (MQ=255)                                                                                                                                      atatCTATTTT‑GCCGGACG‑TTGTG‑AGAGGGCAAAAATCCCCCG‑GATTTTTTTTTACGGGGAGGGTAA‑TGCGCGGA‑AAGATCGCTTGTTTGGGGAGGCGTTTCAATTTATCGGAGATATAAAAAAAC‑‑CCGCCCCTGACAGGCGGG‑AAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTcc                                                                                                              <  1:41952/251‑1 (MQ=255)                                                                                                                                             tttt‑GCCGGAGT‑TTTTG‑AGAGGGGAAAAATTGCCCGCGATTTGTTTT‑GCAGTGAGCGTAA‑TGCGCGGG‑AAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCGGAAATAAAAAAAACC‑‑CCGCCCCTGACAGGCGGG‑AAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAAC‑Cg                                                                                                       <  1:122080/248‑1 (MQ=255)                                                                                                                                             tatT‑GCCGGACG‑TTTTG‑AGATGGCGAAAATTGCCCGCGATCTGTTTT‑GCAGTGAGCGTAA‑TGCGCGGG‑AAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAAC‑‑CCGCCCCTGACAGGCGGG‑AAGCACGGCAAATAAACTGTTATTCAGTGGCCTTTAGATCTATGCCGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGGGCCTACAAATAATCCGAAC‑CGt                                                                                                      >  1:696750/1‑251 (MQ=255)                                                                                                                                             atggtgcgggccG‑TTTTG‑AGAGGGCGAAAATTGCCCGGGATCGGTTTT‑GCAGGGATTGTAA‑GGCGGGGG‑AAGATCGCTGGTTGGGCGATGCGTTTTCATTTATCGGATATAAAAAAAATC‑‑CCGCCCCTGACAGGCGGG‑AAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAAC‑Cg                                                                                                       <  1:622049/240‑1 (MQ=255)                                                                                                                                                t‑GCCGGACG‑TTTTG‑AGATGGCGAAAATTGCCCGCGATCTGTTTT‑GCAGTGAGCGTAA‑TGCGCGGG‑AAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAAC‑‑CCGCCCCTGACAGGCGGG‑AAGAACGGCAACTAAACTGTTATTCAGTGGCCTTTAGATCTATGACGTATCTGGCAAACCATGCCCGGTGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCCAAC‑CGCAgg                                                                                                   >  1:604627/1‑251 (MQ=255)                                                                                                                                                      gACG‑TTTTG‑AGATGGCGAAAATTGCCCGCGATCTGTTTT‑GCAGTGAGCGTAA‑TGCGCGGG‑AAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAAC‑‑CCGCCCCTGACAGGCGGG‑AAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGTTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAAC‑CGTAGGCCggg                                                                                              >  1:3975/1‑250 (MQ=255)                                                                                                                                                            tttG‑AGATGGCGAAAATTGCCCGCGATCTGTTTT‑GCAGTGAGCGTAA‑TGCGCGGG‑AAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAAC‑‑CCGCCCCTGACAGGCGGG‑AAGAACGGCAACTAAACTGTTATTCAGTGGCCTTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTa                                                                                                                        >  1:640096/1‑221 (MQ=255)                                                                                                                                                                 agaTGGCGAAAATTGCCCGCGATCTGTTTT‑GCAGTGAGCGTAA‑TGCGCGGG‑AAGATCGCCTGTTTGGCGATGAGTTTGCATTTATCTGAGATATAAAAAAAC‑‑CCGCCCCTGACAGGAGGG‑AAGAACGGCAACTAAACTGTTTTTCATTGGCATTTTAATCTATGACGTATCTGGCCAACCATGCCCGATGCGGCGATGTCGCGTCTTATCGTGCCTACAAATAGTCCGCAC‑CGTAGGCCGGTTAAGGCGtt                                                                                     >  2:448422/1‑251 (MQ=255)                                                                                                                                                                   gtggaaaaaaTGGCCCGGGATTTGTTTG‑GCAGTGAGCGTAA‑GGCGCGGA‑AAGATCGCGTGTTTGGCGATGCTTTTGCTTTTATCTGAAATAAAAAAAACC‑‑CCGCCCCTGACAGGCGGG‑AAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAAC‑CGTAGGCCGGATAAGGCGTTTa                                                                                   <  2:46350/245‑1 (MQ=37)                                                                                                                                                                                    cgATCTGTTTT‑GCAGTGAGCGTAA‑TGCGCGGG‑AAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAAC‑‑CCGCCCCTGACAGGCGGG‑AAGAAAGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACCTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAAC‑CGTAGGCCGGATAAGGCGTTTACGCCCCATCCGGCAAtt                                                                  >  2:484396/1‑251 (MQ=255)                                                                                                                                                                                        cTGTTTT‑GCAGTGAGCGTAA‑TGCGCGGG‑AAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAAC‑‑CCGCCCCTGACAGGCGGG‑AAGAAAGGCACATAAAATGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCTACGCTGTCGCGTCTTATCTTGCCTACAAATAGTCCGAAA‑CGTAGGCCGGATAAGGCGTTTATGCCGCATCCGGCAATTgggg                                                              >  2:577787/1‑249 (MQ=255)                                                                                                                                                                                          ttttt‑GCGGTGATCGAAA‑GGCGGGGG‑AAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATAAAAAAAAAC‑‑CCGCCCCTGACAGGCGGG‑AAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAAC‑CGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGt                                                               <  1:145378/247‑1 (MQ=37)                                                                                                                                                                                          gTTTT‑CCAGTGAGCG‑AA‑AGCCCGGGGAAGATCGCGTGTTTGGCTATGCGTTTGCTTTTATCTGAGATAAAAAAAACC‑‑CCGCCCCTGACAGGCGGG‑AAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAAC‑CGTAGGCCGGATAAGTCGTTTACGCCTCATCCGGCAATTGGt                                                               <  1:277467/248‑1 (MQ=37)                                                                                                                                                                                            ttt‑GCAGTGAGCGTAA‑TGCGCGGG‑AAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAAC‑‑CCGCCCCTGACAGGCGGG‑AAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCCAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGGCTACCAATATGCCGCAC‑CGGCGGCCCGATAAGGCGTTTACGCCGCATCCGGCAATTggggcctg                                                          >  2:552877/1‑245 (MQ=38)                                                                                                                                                                                                     gAGCGTAA‑TGCGCGGG‑AAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAAC‑‑CCGCCCCTGACAGGCGGG‑AAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTTCCTACAAATAGTCCGAAC‑CGTAGGCCGGATAAGGCGGTTACGCCGCATCCGGCAATTGGTGCATGATGCCTg                                                   >  1:670335/1‑250 (MQ=255)                                                                                                                                                                                                      aGCGTAA‑TGCGCGGG‑AAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAAC‑‑CCGCCCCTGACAGGCGGG‑AAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGGCTAAAAATCATCCGAAC‑CGTAGGCCGGATAAGGAGTTTACGCCGCATCCCGCAATTGGGGCCTGGTGCCTGAt                                                 >  2:573643/1‑251 (MQ=255)                                                                                                                                                                                                         gTAA‑TGCGCGGG‑AAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAAC‑‑CCGCCCCTGACAGGCGGG‑AAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCt                                                                                                                                   >  1:62099/1‑167 (MQ=255)                                                                                                                                                                                                          ggtaagggggggaaagttCGCTTGTTGGGCGAGGCGTTGACATTTATTTGAAATAAAAAAAACC‑‑CCGCCCCTGACAGGCGGG‑AAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAAC‑CGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGc                                               <  2:636621/234‑1 (MQ=21)                                                                                                                                                                                                                   ggg‑AAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAAC‑‑CCGCCCCTGACAGGCGGG‑AAGAAAGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGGATCTGGCCAACCATGCCCGATGCGAAGCTGTCCCGTCTTATCGTGCCTACAAATAGTCCGAAC‑CGGAGGCCGGATAAGGCGTTTACGCCCCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTcg                                     >  2:112220/1‑251 (MQ=37)                                                                                                                                                                                                                     g‑AAGATCGCCTGTTTGGCGCTGCGTTTGCATTTATCTGAGATATAAAAAAAC‑‑CCGCACCTGACAGGCGGG‑AAGAACGACAAATCAACCGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACAACACCCGCTGCGACGCTTGCGCGACTTACCGGGCCTCCCAATAGGCAAAACACTGAGGCCGG‑TAAAGCGTTTACGACGCAACCCGCACTTGGTGCCTGATGCCTGAGGcgccgcattccgg                                   >  2:737269/1‑244 (MQ=255)                                                                                                                                                                                                                                 gTTTGGCGATGCTTTTGCATTTATCTAAGATAAAAAAACAC‑‑CCGCCCCTGACAGGCGGG‑AAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAAC‑CGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGt                          <  1:261858/249‑1 (MQ=255)                                                                                                                                                                                                                                     ggCGATGCGTTTGCATTTATCTGAGATATAAAAAAAC‑‑CCGCCCCTGACAGGCGGG‑AAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAAC‑CGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTAc                    >  1:521726/1‑251 (MQ=255)                                                                                                                                                                                                                                            ccggtttcattctaccagaaataAAAAAACATCCGCCCCTGACAGGCGGG‑AAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGAGCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCTAAC‑CGTAGGCCGGATAAGGCGTTTTCGCCGCATCCGGCAATTGGTGCATGATGCCTCATGCGACGCTGTCGCGTCTTATCGTGCCTACaaaga               <  1:97826/231‑3 (MQ=21)                                                                                                                                                                                                                                              tttGCATTTATCTGAGATATAAAAAAAC‑‑CCGCCCCTGACAGGCGGG‑AAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGGATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAAC‑CGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTTTCGTGCCTACAAATAGTcc           >  1:583684/1‑251 (MQ=255)                                                                                                                                                                                                                                               ttGCATTTATCTGAGATATAAAAAAAC‑‑CCGCCCCTGACAGGCGGG‑AAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAAC‑CGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCg          >  1:8946/1‑251 (MQ=255)                                                                                                                                                                                                                                                         gtgggaTATAAAAAACC‑‑CCGCCCCTGACAGGCGGG‑AAGAACGGCAACTAAACTGTTATTCAGGGGCATTTAGATCTATGACGTATCGGGCAAACCATGCCCGATTCGACGGTGTCGCGTCTTGTCGTGCCTACAAAAAGTCCGAAC‑CGGAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGTGTCTTATCGTGCCTACAAATAGTCCGAACCGTAg  <  2:85571/245‑1 (MQ=25)                                                                                                                                                                                                                                                         cTGAGATATAAAAAAAC‑‑CCGCCCCTGACAGGCGGG‑AAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAAC‑CGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGGGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTAc                    >  1:286552/1‑231 (MQ=255)                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                     GGCGCTTTCGTTGAAGCATCCTGGTCTGCAAGTGGTGCCGGTGGTTGAACAACCGGAAGCGGGCTGGCGTGGGCGTACTGGCACCGTGTTAA‑CGGCGGTATTGCAGGATCACGGTACGCTG‑GCAGAGCATGATATCTATATT‑GCCGGACG‑TTTTG‑AGATGGCGAAAATTGCCCGCGATCTGTTTT‑GCAGTGAGCGTAA‑TGCGCGGG‑AAGATCGCCTGTTTGGCGATGCGTTTGCATTTATCTGAGATATAAAAAAAC‑‑CCGCCCCTGACAGGCGGG‑AAGAACGGCAACTAAACTGTTATTCAGTGGCATTTAGATCTATGACGTATCTGGCAAACCATGCCCGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAAC‑CGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGTCGCGTCTTATCGTGCCTACAAATAGTCCGAACCGTAG  >  NC_000913/4026985‑4027473

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