Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	39,025	0	G	T	29.8%	71.0 / 10.8	47	Q31K (CAA→AAA)	caiB	gamma‑butyrobetainyl‑CoA:carnitine CoA transferase
Reads supporting (aligned to +/- strand):  ref base G (13/20);  new base T (13/1);  total (26/21)
Fisher's exact test for biased strand distribution p-value = 9.63e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CACAGTACTTCATCGGTAATGCCACGACGGGCAAAGGCCGGACCTTTACTGGCTTCGATGAAGATATCGGTGGTTTCCATTAATTTCAGAAACGCTTCGCGGCCTTCATCTTTGAAAATATTTAACGACAGCGCGTGCAAATTGCGGCGGGAGAGTTGCGGGTAG‑TTCGG‑TTGAACGCGAATGGTGTC‑GGCCCAGGCGAC‑GTTCTC‑GATCCAGATAACTTCCGCGCCCCATTCTGCGAACATTTGCCCGGCAAACGGTCCGGCGATTTCG‑ATACCGGAGAAGACAACGCG‑CAATCCGGCCAACGGCCCGAATTTCGGCATGGGTAGATGATCCATTATTTGCTCCTGAAAAATTTATGTAGCGCATGACTGCCGGATGCGGCGTAAACGCTTTATCCGGCCTACATTCGTGCTCCCGTAGGCCTGATAAGACGCATCAGCGTCGCATCAGGCAGCGCACGGACTTAGCGGTATTGCTTCAGCACCGCACGACCC  >  NC_000913/38781‑39274                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                            cacaGTACTTCATCGGTAATGCCACGACGGGCAAAGGCCGGACCTTTACTGGCTTCGATGAAGATATCGGTGGTTTCCATTAATTTCAGAAACGCTTCGCGGCCTTCATCTTTGAAAATATTTAACGACAGCGCGTGCAAATTGCGGCGGGGGAGTTGCGGGGAG‑TTCGG‑TTGAACGCGAAAGGGGTC‑GGCCCCGGCGAC‑GGTCTC‑GGTCCCGATAAATTTCGCGCCCCCTTCTTCCAACAATTTCCCGGc                                                                                                                                                                                                                                                       >  1:468154/1‑251 (MQ=255) cacaGTACTTCATCGGTAATGCCACGACGGGCAAAGGCCGGACCTTTACTGGCTTCGATGAAGATATCGGTGGTTTCCATTAATTTCAGAAACGCTTCGCGGCCTTCATCTTTGAAAATATTTAACGACAGCGCGTGCAAATTGCGGCGGGGGAGGTGCGGGTAG‑ATTGG‑GTGGACGCGGAAGGTGTC‑GGCCCAGGCGAC‑GTTCTC‑GATCCAGCTAAACTCCGCCCCCCCTTCTTCGAACCTTTTCCCGGc                                                                                                                                                                                                                                                       >  1:468246/1‑251 (MQ=255)            aTCGGTAATGCCACGACGGGCAAAGGCCGGACCTTTACTGGCTTCGATGAAGATATCGGTGGTTTCCATTAATTTCAGAAACGCTTCGCGGCCTTCATCTTTGAAAATATTTAACGACAGCGCGTGCAAATTGCGGCGGGAGAGTTGCGGGTAG‑TTCGG‑TTGAACGCGAATGGTGTC‑GGCCCAGGCGAC‑GTTCTC‑GATCCAGATAACTTCCGCGCCCCATTCTGCGAACATTTGCCCGGCAAACGGTCCgg                                                                                                                                                                                                                                            <  2:326659/251‑1 (MQ=255)              cGGTAATGCCACGACGGGCAAAGGCCGGACCTTTACTGGCTTCGATGAAGATATCGGTGGTTTCCATTAATTTCAGAAACGCTTCGCGGCCTTCATCTTTGAAAATATTTAACGACAGCGCGTGCAAATTGCGGCGGGAGAGTTGCGGGTAG‑TTCGG‑TTGAACGCGAATGGTGTC‑GGCCCAGGCGAC‑GTTCTC‑GATCCAGATAACTTCCGCGCCCCATTCTGCGAACATTTGCCCGGCAAACGGTCCGGCg                                                                                                                                                                                                                                          <  1:322884/251‑1 (MQ=255)                gTAATGCCACGACGGGCAAAGGCCGGACCTTTACTGGCGTCGATGAAGATATCGGTGGTTTCCATTAATTTCAGAAACGCTTCGCGGCCTTCATCTTTGAAAATATTTAACGACAGCGCGTGCAAATTGCGGCGGGAGAGTTGCGGGTAG‑TTCGG‑TTGAACGCGAATGGTGTC‑GGCCCAGGCGAC‑GTTCTC‑GATCCAGATAACTTCCGCGCCCCATTCTGCGAACATTTGCCCGGCAAACGGTCCGGCGAt                                                                                                                                                                                                                                        <  1:478461/251‑1 (MQ=255)                  aaTGCCACGACGGGCAAAGGCCGGACCTTTACTGGCTTCGATGAAGATATCGGTGGTTTCCATTAATTTCAGAAACGCTTCGCGGCCTTCATCTTTGAAAATATTTAACGACAGCGCGTGCAAATTGCGGCGGGAGAGGTGGGGGGAGGTTGGG‑TTAACCGGGAAGGGGGT‑‑GGCCCCGGCGGC‑GTTCTC‑GGTCCCGGTAACTTCCGCCCCCCCCTCTTCCAACATTTTCCCGGCCAACGGTCCGGCGCttt                                                                                                                                                                                                                                      >  2:161033/1‑251 (MQ=255)                       cACGACGGGCAAAGGCCGGACCTTTACTGGCTTCGATGAAGATATCGGTGGTTTCCATTAATTTCAGAAACGCTTCGCGGCCTTCATCTTTGAAAATATTTAACGACAGCGCGTGCAAATTGCGGCGGGAGAGGTGCGGGGAG‑GTCGG‑GTGAAAGCGGATGGGGTC‑GGGCCAGGCGGC‑GGTCTC‑GAAACAAATAACTTCCGCGCCCCCATCTTCGGACATTTTCCCGGCCAACGGTACGGCGGttttggt                                                                                                                                                                                                                                  >  1:198867/1‑246 (MQ=255)                           acgGGCAAAGGCCGGACCTTTACTGGCTTCGATGAAGATATCGGTGGTTTCCATTAATTTCAGAAACGCTTCGCGGCCTTCATCTTTGAAAATATTTAACGACAGCGCGTGCAAATTGCGGCGGGAGAATTGCGGGTAG‑GTCGG‑GTGAAAACGAATGGGGTCCGCCCCGGGGGAA‑GT‑CTC‑TATCCCGATAAACTCCGCGCCCCCCTCTCCGAAAAATTGCCCCGCAAACGGTcccgccgtttccgtcacgg                                                                                                                                                                                                                             >  1:398686/1‑234 (MQ=255)                             gggCAAAGGCCGGACCTTTACTGGCTTCGATGAAGATATCGGTGGTTTCCATTAATTTCAGAAACGCTTCGCGGCCTTCATCTTTGAAAATATTTAACGACAGCGCGTGCAAATTGCGGCGGGAGGGGTGCGGGGAG‑GTCGG‑TTGAACGCGGATGGGGTC‑GGCCCCGGCGGC‑GTTCTC‑GATCCAGATAAATTTCGCGCCCCCTTCTGCGAACATTTTCCCGGCAAACGGGCCGGCGGTTTTG‑ATAcccg                                                                                                                                                                                                                            >  1:310453/1‑247 (MQ=255)                                     gCCGGACCTTTACTGGCTTCGATGAAGATATCGGTGGTTTCCATTAATTTCAGAAACGCTTCGCGGCCTTCATCTTTGAAAATATTTAACGACAGCGCGTGCAAATTGCGGCGGGAGAGTTACGGGGAG‑GTCGGGTTGAA‑GCGGATGGGGTC‑GGCCCCGACGAC‑GGTCTCTGACCCA‑AAAACCTCCGCGCCCCATTCTTCGGACATTTGTCCGGCAAACGGTGCGGCGAATTCGCATACCG‑CGAAGACaa                                                                                                                                                                                                                   >  1:79769/1‑250 (MQ=255)                                         gACCTTTACTGGCTTCGATGAAGATATCGGTGGTTTCCATTAATTTCAGAAACGCTTCGCGGCCTTCATCTTTGAAAATATTTAACGACAGCGCGTGCAAATTGCGGCGGGAGAGTTTCGGGGAG‑TTCGGGTTGAA‑GCGAAAGGGGTC‑GGCCCCGGCGAC‑GTTCTC‑GATCCAGAAAACTTCCGCGCCCCCTTCTGCGAACATTTTCCCGGCCAACGGTCCGGCGGTTTCG‑GTACCGGGGAAGACAAAgcg‑c                                                                                                                                                                                                             >  1:304639/1‑251 (MQ=255)                                           ccTTTACTGGCTTCGATGAAGATATCGGTGGTTTCCATTAATTTCAGAAACGCTTCGCGGCCTTCATCTTTGAAAATATTTAACGACAGCGCGTGCAAATTGCGGCGGGAGAGGTGCGGGGAG‑TTCGG‑TTGGACGCGAAAGGTGGC‑GGCCCACGCGAC‑GTTCTC‑GATTCAGATAAATTCCGCGCCCCATTCTGCGGACATTTTCCCGGCAAAAGGTCCGGCGATTTTG‑ATACCGGGGAAGAAAAAGCG‑Ca                                                                                                                                                                                                            >  1:270914/1‑250 (MQ=255)                                                 cTGGCTTCGATGAAGATATCGGTGGTTTCCATTAATTTCAGAAACGCTTCGCGGCCTTCATCTTTGAAAATATTTAACGACAGCGCGTGCAAATTGCGGCGGGAGAGTTACGGGTAG‑GTCGG‑TTGAAAGCGCATGGTGTC‑CGCCCCGGGGAC‑GCTCTC‑TGTCCCGATCAATTCCACCCCCCCTTCTCCGGACATTTTCCCCGCAAACGGGCCCGCGGTTTTG‑GTACCGGAGGAGACAACacg‑ccaacccgc                                                                                                                                                                                                     >  1:242506/1‑243 (MQ=255)                                                         gATGAAGATATCGGTGGTTTCCATTAATTTCAGAAACGCTTCGCGGCCTTCATCTTTGAAAATATTTAACGACAGCGCGTGCAAATTGCGGCGGGAGAGTTGCGGGGAG‑TTCGG‑GTGGACGCGAATGGGGGC‑GGCCCAAGCGGC‑GTTCTC‑GATCCAGATAACTTCCCCGCCCCCTTCTTCGAACCTTTTCCCGGCAAACCGGCCGGCGGTTTCG‑GTACCCGGGAAGAAAACCCC‑CCATCCCGCCCACGGcc                                                                                                                                                                                             >  1:23202/1‑251 (MQ=255)                                                                    cGGTGGGTTCCATTAATTTCAGAAACGCGTCGCGGCCTGCATCTTTGAAAATATTTAACGACAGCGCGTGCAAATTGCGGCGGGAGAGTTGCGGGTAG‑TTCGG‑TTGAACGCGAATGGTGGC‑GGCCCAGGCGAC‑GTTCTC‑GATCCAGATAACTTCCGCGCCCCATTCTGCGAACATTTGCCCGGCAAACGGTCCGGCGATTTCG‑ATACCGGAGAAGACAACGCG‑CAATCCGGCCAACGGCCCGAATTTCgg                                                                                                                                                                                   <  2:89510/250‑1 (MQ=255)                                                                            tCCATTAATTTCAGAAACGCTTCGCGGCCTTCATCTTTGAAAATATTTAACGACAGCGCGTGCAAATTGCGGCGGGAGAGGTGCGGGGAG‑TTTGG‑TTGAACGCGAATGGGGTC‑GGCCCCAGCGAC‑GCTTTC‑GATCCAGCTAAATTCCCCGCCCCCTTCTGCGAACATTTGCCCGGCAAAAGGTTCGGCGAATTCG‑CTACCGGGGAAGACAACGCGGCAA‑CCGGCCAAAGGccccgattttggcctggggcg                                                                                                                                                                          >  1:445447/1‑233 (MQ=255)                                                                                       cAGAAACCCTTCGCGGCCTTCATGTTTCAAAATATTTACCGACAGCGCGTGCAAATTGCGGCGGGAGAGTTGCGGGTAA‑TTCGG‑TTGAACGCGCATGGTTTC‑GGCCCAGGCGAC‑GTTCGC‑GATCCAGATAACTTCCGCGCCCCATTCTGCGAACATTTGCCCGGCAAACGGTCCGGCGATTTCG‑ATACCGGAGAAGACAACGTG‑CAATCCGGCCAACGGCCCGAATTTCGGCATGGGTAGATGATCCatt                                                                                                                                                                <  2:317745/250‑1 (MQ=255)                                                                                          aaaCGCTTCGCGGCCTTCATCTTTGAAAATATTTAACGACAGCGCGTGCAAATTGCGGCGGGAGAGTTGCGGGTAG‑TTCGG‑TTGAACGCGAATGGTGTC‑GGCCCAGGCGAC‑GTTCTC‑GATCCAGATAACTTCCGCGCCCCATTCTGCGAACATTTGCCCGGCAAACGGTCCGGCGATTTCG‑ATACCGGAGAAGACAACGCG‑CAATCCGGCCAACGGCCCGAATTTCGGCATGGGTAGATGATCCATTAttt                                                                                                                                                            <  2:370354/251‑1 (MQ=255)                                                                                             cGCTTCGCGGCCTTCATCTTTGAAAATATTTNACGACAGCGCGTGCAAATTGCGGCGGGAGAGTTGCGGGTAG‑TTCGG‑TTGAACGCGAATGGTGTC‑GGCCCAGGCGAC‑GTTCTC‑GATCCAGATAACTTCCGCGCCCCATTCTGCGAACATTTGCCCGGCAAACGGTCCGGCGATTTCG‑ATACCGGAGAAGACAACGCG‑CAATCCGGCCAACGGCCCGAATTTCGGCATGGGTAGATGATCCATTATTTGCt                                                                                                                                                         <  2:218749/251‑1 (MQ=255)                                                                                                 tCGCGGCCTTCATCTTTGAAAATATTTAACGACAGCGCGTGCAAATTGCGGCGGGAGAGTTGCGGGTAG‑TTCGG‑TTGAACGCGAATGGTGTC‑GGCCCAGGCGAC‑GTTCTC‑GATCCAGATAACTTCCGCGCCCCATTCTGCGAACATTTGCCCGGCAAACGGTCCGGCGATTTCG‑ATACCGGAGAAGACAACGCG‑CAATCCGGCCAACGGCCCGAATTTCGGCATGGGTAGATGATCCATTATTTGCTCCt                                                                                                                                                      <  1:152155/250‑1 (MQ=255)                                                                                                   gcgGCCTTCATCTTTGAAAATATTTAACGACAGCGCGTGCAAATTGCGGCGGGAGAGGTTCGGGGAG‑TTCGG‑TTGAACGCGAATGGGGTC‑GAACCAGGCGAC‑GTTCTC‑GATCCAGATAAATTCCCCGCCCCCTTTTCCGAACATTTTTCCGGCAAAAGGTCCGGCGCTTTCG‑ATACCGGAGAAAACAACCCG‑CAATCCCGCCAACAGCCCGGAATTCTGCATGGGGCGATGATCCCTTATTTGCTCCtta                                                                                                                                                    >  1:301626/1‑248 (MQ=255)                                                                                                   gcgGCCTTCATCTTTGAAAATATTTAACGACAGCGCGTGCAAATTGCGGCGGGAGAGGTGCGGGGAG‑TTCGG‑GTGAACCCGAATGGGGTC‑GGCCCCGGCGGC‑GGTCTC‑GATCCAGATAAATTCCGCGCCCCAATCTGCGAAAATTTTCCCCGCAAACGGTCCCGCGGTTTCG‑ATACCGGGGAAGACAACCCC‑CCATCCGGCCCAAGGCCCGCATTTTGGCCAGGGGAgatgaaccctttttttcccccga                                                                                                                                                    >  1:409485/1‑232 (MQ=255)                                                                                                          tCATCTTTGAAAATATTTAACGACAGCGCGTGCAAATTGCGGCGGGAGAGTTGCGGGGAG‑TTCGG‑GTGAACGCGAAAGGTGTC‑GGGCCAGGCGAC‑GTTCTC‑GATCCCGATAACTTCCGCGCCCCCTTCTGCGGACATTTTCCCGGCAAACGGTCCGGCGGTTTCG‑GTACCGGAGAAGAAAAAGCG‑CAAACCCGCCCACGGCCCGAATTTCGGCCTGGGGAGATGAACCATTTTTTGCTCCCGAAAAAtttt                                                                                                                                            >  1:91996/1‑250 (MQ=255)                                                                                                                    aaaTATTTAACGACAGCGCGTGCAAATTGCGNGGGGAGAGTTGCGGGTAG‑TTCGG‑TTTAACGCGAATGGTGTC‑GGCCCAGGCGAC‑GTTCTC‑GATCCAGATAACTTCCGCGCCCCATTCTGCGAACATTTGCCCGGCAAACGGTCCGGCGATTTCG‑ATACCGGAGAAGACAACGCG‑CAATCCGGCCAACGGCCCGAATTTCGGCATGGGTAGATGATCCATTATTTGCTCCTGAAAAATTTATGTAGCGCAt                                                                                                                                  <  2:242506/251‑1 (MQ=255)                                                                                                                          ttgaCGACAGCTCGTGCAAATGGCGGCGGGAGAGTTGCGGTGAG‑TTCGG‑TTGAACGGGAATTGTGTC‑GGCCCAGGCGAC‑GTTCTC‑GATCCAGATAACTTCCGCGCCCCATTCTGCGAACATTTGCCCGGCAAACGGTCCGGCGATTTCG‑ATACCGGAGAAGACAACGCG‑CAATCCGGCCAACGGCCCGAATTTCGGCATGGGTAGATGATCCATTATTTGCTCCTGAAAAATTGATGTAGCGCATGACTg                                                                                                                             <  2:110464/247‑1 (MQ=255)                                                                                                                           tAACGACAGCGCGGGCAAAGGGCTGCGGGANAGTTGCGGGTAG‑TTCGG‑TTGAACGCGAATGGTGTC‑GGCCCAGGCGAC‑GTTCCC‑GATCCAGATAACTGCCGCGACCCATTCTGCGAACATTTGCGCGGCAAACGGTCCGGCGATTTCG‑ATCCCGGAGAAGACAACGCG‑CACTCCGGCCAACGGCCCGAATTTCGGCATGGGTAGATGATCCATTATTTGCTCCTGAAAAATTTATGTAGCGCATGACTGc                                                                                                                            <  2:211722/250‑1 (MQ=255)                                                                                                                            aaCGACAGCGCGTGCAAATTGCGGCGGGAGAGTTGCGGGTAG‑TTCGG‑TTGAACGCGAATGGTGTC‑GGCCCAGGCGAC‑GTTCTC‑GATCCAGATAACTTCCGCGCCCCATTCTGCGAACATTTGCCCGGCAAACGGTCCGGCGATTTCG‑ATACCGGAGAAGACAACGCG‑CAATCCGGCCAACGGCCCGAATTTCGGCATGGGTAGATGATCCATTATTTGCTCCTGAAAAATTTATGTAGCGCATGACTGCCg                                                                                                                          <  1:136333/251‑1 (MQ=255)                                                                                                                                                      ggAGAGTTGCGGGTAG‑TTCGG‑TTGAACGCGAATGGTGTC‑GGCCCAGGCGAC‑GTTCTC‑GATCCAGATAACTTCCGCGCCCCATTCTGCGAACATTTGCCCGGCAAACGGTCCGGCGATTTCG‑ATACCGGAGAAGACAACGCG‑CAATCCGGCCAACGGCCCGAATTTCGGCATGGGTAGATGATCCATTATTTGCTCCTGAAAAATTTATGTAGCGCATGACTGCCGGATGCGGCGTAAACGCTTTATCCGGc                                                                                                <  2:295167/251‑1 (MQ=255)                                                                                                                                                       gctggtgacGGGTAG‑TTCGG‑TTGAACGCGAATGGTGCC‑GGCCCAGGCGGCTGTGC‑C‑GTGCCCGATATCTTCCGCGCCCCAGTCTGCGAACATTTGTCCGGCAAACGGTCCGGCTATTTCG‑ATACCGTAGAAGGCAACGCG‑CAATCCGGCCAACGGCCTGGATTTCGGCATGGGTGGATGATCCATTATTTGCTCCTGAAAAATTTATGTCGCTCATGGCTTCAGGATGCGGCGTAAACGCTTTATCCGGca                                                                                               <  2:263066/243‑2 (MQ=255)                                                                                                                                                                    aG‑TTCGG‑TTGAACGCGAATGGTGTC‑GGCCCAGGCGGC‑GTTCTC‑GATCCAGATAACTTCCGCGCCCCATTCTGCGAACATTTGCCCGGCAAACGGTCCGGCGATTTCG‑ATACCGGAGAAGACAACGCG‑CAATCCGGCCAACGGCCCGAATTTCGGCATGGGTAGATGATCCATTATTTGCTCCTGAAAAATTTATGTAGCGCATGACTGCCGGATGCGGCGTAAACGCTTTATCCGGCCTACATTCGTGCt                                                                                   <  1:482477/250‑1 (MQ=255)                                                                                                                                                                                 aCGCGAATGGTGTC‑GGCCCAGGCGAC‑GTTCTC‑GATCCAGATAACTTCCGCGCCCCATTCTGCGAACATTTGCCCGGCAAACGGTCCGGCGATTTCG‑ATACCGGAGAAGACAACGCG‑CAATCCGGCCAACGGCCCGAATTTCGGCATGGGTAGATGATCCATTATTTGCTCCTGAAAAATTTATGTAGCGCATGACTGCCGGATGCGGCGTAAACGCTTTATCCGGCCTACATTCGTGCTCCCGTAGGCCTg                                                                       >  1:455427/1‑251 (MQ=255)                                                                                                                                                                                     gAATGGTGTC‑GGCCCAGGCGAC‑GTTCTC‑GATCCAGATAACTTCCGCGCCCCCTTCTTCGAACATTTTCCCTTCAAACGGTCCTTCTATTTCG‑ATCCCTTCGAAGACACCTCT‑CCATCCGGCCCCCTTCCCGAATTTCTTCCTGGGTAGATTCTCCCTTATTTGCTCCTTCAAAATTTCCTTCGCGCATGACCGCACGATGCGGCGTAAACGCTTTATCCGGCCTACTTTCGTGCTCCCGTAGGCCTGATa                                                                    >  1:251291/1‑250 (MQ=255)                                                                                                                                                                                      aaTGGTGTC‑GGCCCATGCGAC‑GTTCTC‑GATCCAGATAACTTCCGCGCCCCATTCTGCGAACATTTGCCCGGCAAACGGTCCGGCGATTTCG‑ATACCGGAGAAGACAACGCG‑CAATCCGGCCAACGGCCCGAATTTCGGCATGGGTAGATGATCCATTATTTGCTCCTGAAAAATTTATGTAGCGCATGACTGCCGGATGCGGCGTAAACGCTTTATCCGGCCTACATTCGTGCTCCCGTAGGCCTGATAAg                                                                  <  2:12345/251‑1 (MQ=255)                                                                                                                                                                                            gtC‑GGCCCAGGCGAC‑GTTCTC‑GATCCAGATAACTTCCGCGCCCCATTCTGCGAACATTTGCCCGGCAAACGGTCCGGCGATTTCG‑ATACCGGAGAAGACAACGCG‑CAATCCGGCCAACGGCCCGAATTTCGGCATGGGTAGATGATCCATTATTTGCTCCTGAAAAATTTATGTAGCGCATGACTGCCGGATGCGGCGTAAACGCTTTATCCGGCCTACATTCGTGCTCCCGTAGGCCTGATAAGACGCAt                                                            <  2:440926/251‑1 (MQ=255)                                                                                                                                                                                                  cccAGGCGAC‑GTTCTC‑GATCCAGATAACTTCCGCGCCCCATTCTGCGAACATTTGCCCGGCAAACGGTCCGGCGATTTCG‑ATACCGGAGAAGACAACGCG‑CAATCCGGCCAACGGCCCGAATTTCGGCATGGGTAGATGATCCATTATTTGCTCCTGAAAAATTTATGTAGCGCATGACTGCCGGATGCGGCGTAAACGCTTTATCCGGCCTACATTCGTGCTCCCGTAGGCCTGATAAGACGCATCAGCg                                                       >  1:57492/1‑251 (MQ=255)                                                                                                                                                                                                     aGGCGAC‑GTTCTC‑GATCCAGATAACTTCCGCGCCCCATTCTGCGAACATTTGCCCGGCAAACGGTCCGGCGATTTCG‑ATACCGGAGAAGACAACGCG‑CAATCCGGCCAACGGCCCGAATTTCGGCATGGGTAGATGATCCATTATTTGCTCCTGAAAAATTTATGTAGCGCATGACTGCCGGATGCGGCGTAAACGCTTTATCCGGCCTACATTAGTGCTCCCGTAGGCCTGATAAGACGCATAAGCGTCg                                                    >  2:435709/1‑251 (MQ=255)                                                                                                                                                                                                      ggCGAC‑GTTCTC‑GATCCAGATAACTTCCGCGCCCCATTCTGCGAACATTTGCCCGGCAAACGGTCCGCCGATTTCG‑ATACCGGAGAAGACAACGCG‑CAATCCGGCAAACGGCCCGGATTTCGGCATGGGTAGATGATCCCTTATTTGCTCCTGAAAAATTTATGTACCGCATGACTGCCGGATGCGGCGTAAACCCTTTATCCGGCCTACATTAGTGCTCCAGTAGGCCTGATAAGACGCATCAGCgccac                                                   >  1:443971/1‑247 (MQ=255)                                                                                                                                                                                                         gAC‑GTTCTC‑TATCCAGATAACTTCCGCGCCCCATTCTGCGAACATTTGCCCGGCAAACGGTCCGGCGATTTCG‑ATACCGGAGAAGACAACGCG‑CAATCCGGCCAACGGCCCGAATTTCGGCATGGGTAGATGATCCATTATTTGCTCCTGAAAAATTTATGTAGCGCATGACTGCCGGATGCGGCGTAAACGCTTTATCCGGCCTACATTCGTGCTCCCGTAGGCCTGATAAGACGCATCAGCGTCGCATc                                                <  1:207181/251‑1 (MQ=255)                                                                                                                                                                                                             gTTCTC‑GATCCAGATAGCTTCCGCGCCCCATTCTGCGAACATTTTCCCGGCAAACGGTCCTGCGATTTCG‑ATACCGGAGAAGACAACTCG‑CAATCCGGCCAACGGCCCGAATTTCGGCATGGGTAGATGATCCATTATTTGCTCCTGATAAATTTATGTAGCGCATGACTGCCGGATGCGGCGTAAACGCTTGATCCGTCCTACATTCGTGCTCCCGTAGGCCTGATAAGACGCATCAGCGTCGCATCAgg                                             <  1:172853/251‑1 (MQ=255)                                                                                                                                                                                                                     aTCCAGATAACTTCCGCGCCCCATTCTGCGAACATTTGCCCGGCAAACGGTCCGGCGATTTCG‑ATACCGGAGAAGACAACGCG‑CAATCCGGCCAACGGCCCGAATTTCGGCATGGGTAGATGATCCATTATTTGCTCCTGAAAAATTTATGTAGCGCATGACTGCCGGATGCGGCGTAAACGCTTTATCCGGCCTACATTCGTGCTCCCGTAGGCCTGATAAGACGCATCAGCGTCGCATCAGGCAGCGCa                                      <  1:210369/251‑1 (MQ=255)                                                                                                                                                                                                                       ccAGATAACTTCCGCGCCCCATTCTGCGAACATTTGCCCGGCAAACGGTCCGGCGATTTCG‑ATACCGGAGAAGACAACGCG‑CAATCCGGCCAACGGCCCGAATTTCGGCATGGGTAGATGATCCATTATTTGCTCCTGAAAAATTTATGTAGCGCATGACTGCCGGATGCGGCGTAAACGCTTTATCCGGCCTACATTCGTGCTCCCGTAGGCCTGATAAGACGCATCAGCGTCGCATCAGGCAGCGCACg                                    >  1:400945/1‑251 (MQ=255)                                                                                                                                                                                                                           aTAACTTCCGCGCCCCATTCTGCGAACATTTGCCCGGCAAACGGTCCGGCGATTTCG‑ATACCGGAGAAGACAACGCG‑CAATCCGGCCAACGGCCCGAATTTCGGCATGGGTAGATGATCCATTATTTGCTCCTGAAAAATTTATGTAGCGCATGACTGCCGGATGCGGCGTAAACGCTTTATCCGGCCTACATTCGTGCTCCCGTAGGCCTGATAAGACGCATCAGCGTCGCATCAGGCAGCGCACGGACt                                >  1:100563/1‑251 (MQ=255)                                                                                                                                                                                                                            tAACTTCCGCGCCCCATTCTGAGAACATTTGCCCGGCAAACGGTCCGGCGATTTCG‑ATACCGGAGAAGACAACGCG‑CAATCCGGCCAACGGCCCGAATTTCGGCATGGGTAGATGATCCATTATTTGCTCCTGAAAAATTTATGTAGCGCATGACTGCCGGATGCGGCGTAAACGCTTTATCCGGCCTACATTCGTGCTCCCGTAGGCCTGATAAGACGCATCAGCGTCGCATCAGGCAGCGCACGGACtt                               >  1:347995/1‑251 (MQ=255)                                                                                                                                                                                                                                  ccGCGCCCCATTCTGCGAACATTTGCCCGGCAAACGGTCCGGCGATTTCG‑ATACCGGAGAAGACAACGCG‑CAATCCGGCCAACGGCCCGAATTTCGGCATGGGTAGATGATCCATTATTTGCTCCTGAAAAATTTATGTAGCGCATGACTGCCGGATGCGGCGTAAACGCTTTATCCGGCCAACATTCGTGCTCCCGTAGGCCTGATAATACGCATCAGCGTCGCCTCAGGCAGCGCACGGCCTTAGctgt                         >  1:229379/1‑248 (MQ=255)                                                                                                                                                                                                                                   cgcgCCCCATTCTGCGAACATTTGCCCGGCAAACGGTCCGGCGATTTCG‑ATACCGGAGAAGACAACGCG‑CAATCCGGCCAACGGCCCGAATTTCGGCATGGGTAGATGATCCATTATTTGCTCCTGAAAAATTTATGTAGCGCATGACTGCCGGATGCGGCGTAAACGCTTTATCCGGCCTACATTCGTGCTCCCGTAGGCCTGATAAGACGCATCAGCGTCGCATCAGGCAGCGCACGGACTTAGCGGt                         >  1:367107/1‑250 (MQ=255)                                                                                                                                                                                                                                                gCGAACATTTGCCCGGCAAACGGTCCGGCGATTTCG‑ATACCGGAGAAGACAACGCG‑CAATCCGGCCAACGGCCCGAATTTCGGCATGGGTAGATGATCCATTATTTGCTCCTGAAAAATTTATGTAGCGCATGACTGCCGGATGCGGCGTAAACGCTTTATCCGGCCTACATTCGTGCTCCCGTAGGCCTGATAAGACGCATCAGCGTCGCATCAGGCAGCGCACGGACTTAGCGGTATTGCTTCAGCAc            <  2:270914/250‑1 (MQ=255)                                                                                                                                                                                                                                                      aTTTGCCCGGCAAACGATCCGGCGATTACG‑ATACCGGAGAAGACAACGCG‑CAATCCAGCCAACGGCCCGAATTTCGGCATGGGCAGATGATCCATTATTTGATCCTGAAAAATTTATGTAGCGCATGACTGCCGGATGCGGCGTAAACGCTTTATCCGGACTACATTCGAGCTCCCGCAGGACTGATAAGCCGCCTCAGCGTCGCATCAGGCAGCGCACGGAATTCGCGGAATTGCTTCAGAACCGCACg      >  2:34642/1‑250 (MQ=255)                                                                                                                                                                                                                                                          gCCCGGCAAACGGTCCTGCGATTTCG‑ATACCGGAGAAGACAACGCG‑CAATCCGGCCAACGGCCCGAATTTCGGCATGGGTAGATGATCCATTATTTGCTCCTGAAAAATTTATGTAGCGCATGACTGCCGGATGCGGCGTAAACGCTTTATCCGGCCTACATTCGTGCTCCCGTAGGCCTGATAAGACGCATCAGCGTCGCATCAGGCAGCGCACGGACTTAGCGGTATTGCTTCAGCACCGCACGAccc  <  2:88542/250‑1 (MQ=255)                                                                                                                                                                                                                                                          gCCCGGCAAACGGTCCGGCGATTTCG‑ATACCGGAGAAGACAACGCG‑CAATCCGGCCAACGGCCCGAATTTCGGCATGGGTAGATGATCCATTATTTGCTCCTGAAAAATTTATGTAGCGCATGACTGCCGGATGCGGCGTAAACGCTTTATCCGGCCTACATTCGTGCTCCCGTAGGCCTGATAAGACGCATCAGCGTCGCATCAGGCAGCGCACGGACTTAGCGGTATTGCTTCAGCACCGCACGAccc  >  1:164189/1‑250 (MQ=255)                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                            CACAGTACTTCATCGGTAATGCCACGACGGGCAAAGGCCGGACCTTTACTGGCTTCGATGAAGATATCGGTGGTTTCCATTAATTTCAGAAACGCTTCGCGGCCTTCATCTTTGAAAATATTTAACGACAGCGCGTGCAAATTGCGGCGGGAGAGTTGCGGGTAG‑TTCGG‑TTGAACGCGAATGGTGTC‑GGCCCAGGCGAC‑GTTCTC‑GATCCAGATAACTTCCGCGCCCCATTCTGCGAACATTTGCCCGGCAAACGGTCCGGCGATTTCG‑ATACCGGAGAAGACAACGCG‑CAATCCGGCCAACGGCCCGAATTTCGGCATGGGTAGATGATCCATTATTTGCTCCTGAAAAATTTATGTAGCGCATGACTGCCGGATGCGGCGTAAACGCTTTATCCGGCCTACATTCGTGCTCCCGTAGGCCTGATAAGACGCATCAGCGTCGCATCAGGCAGCGCACGGACTTAGCGGTATTGCTTCAGCACCGCACGACCC  >  NC_000913/38781‑39274

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