Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,984,090	0	T	C	34.8%	81.0 / 11.2	46	F378S (TTC→TCC)	aslB	putative anaerobic sulfatase maturation enzyme AslB
Reads supporting (aligned to +/- strand):  ref base T (26/4);  new base C (0/16);  total (26/20)
Fisher's exact test for biased strand distribution p-value = 4.89e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ACGCCTGTGATCACTATGTTTATCCGCAATATCGGCTGGGGAATATGCACCAGCAAACAATTGCAGAAATGATCGATTCCCCGCAACAGCAGGCGTTTGG‑TGAAGATAAATTTAAGCAGTTACCGGC‑GCAGTGTCGCAGTTGTAACGTGTTAAAAG‑CGTGCTGGGGAGGCTGCCCGAAACACC‑GCTTCATGCTCGATGCCAGCGGCAAACCGGGA‑CTGAATTATTT‑GTGTGCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGC  >  NC_000913/3983841‑3984339                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                          aCGCCTGTGATCACTATGTTTATCCGCAATATCGGCTGGGGAATATGCACCAGCAAACAATTGCAGAAATGATCGATTCCCCGCAACAGCAGGCGTTTGG‑TGAAGATAAATTTAAGCAGTTACCGGC‑GCAGTGTCGCAGTTGTAACGTGTTAAAAG‑CGTGCTGGGGAGGCTGCCCGAAACACC‑GCTTCATGCTCGATGCCAGCGGCAAACCGGGA‑CTGAATTATTT‑GTGTGCCGGGTATCAGCGTTAttt                                                                                                                                                                                                                                                           <  1:117656/250‑1 (MQ=255) aCGCCTGTGATCACTATGTTTATCCGCAATATCGGCTGGGGAATATGCACCAGCAAACAATTGCAGAAATGATCGATTCCCCGCAACAGCAGGCGTTTGG‑TGAAGATAAATTTAAGCAGTTACCGGC‑GCAGTGTCGCAGTTGTAACGTGTTAAAAG‑CGTGCTGGGGAGGCTGCCCGAAACACC‑GCTTCATGCTCGATGCCAGCGGCAAACCGGGA‑CTGAATTATTT‑GTGTGCCGGGTATCAGCGTTATTTc                                                                                                                                                                                                                                                          >  1:466552/1‑251 (MQ=255)            cACTATGTTTATCCGCAATATCGGCTGGGGAATATGCACCAGCAAACAATTGCAGAAATGATCGATTCCCCGCAACAGCAGGCGTTTGG‑TGAAGATAAATTTAAGCAGTTACCGGC‑GCAGTGTCGCAGTTGTAACGTGTTAAAAG‑CGTGCTGGGGAGGCTGCCCGAAACACC‑GCTTCATGCTCGATGCCAGCGGCAAACCGGGA‑CTGAATTATTT‑GTGTGCCGGGTATCAGCGTTATTTCCGCCATCTAcc                                                                                                                                                                                                                                               <  1:39030/251‑1 (MQ=255)            cACTATGTTTATCCGCAATATCGGCTGGGGAATATGCACCAGCAAACAATTGCAGAAATGATCGATTCCCCGCAACAGCAGGCGTTTGG‑TGAAGATAAATTTAAGCAGTTACCGGC‑GCAGTGTCGCAGTTGTAACGTGTTAAAAG‑CGGGCTGGGGAGGCTGCCCGAAACACC‑GCTTCATGCTCGATGCCAGCGGCAAACCGGGA‑CTGAATTATTT‑GTGTGCCGGGTATCAGCGTTATTTCCGCCATCTAcc                                                                                                                                                                                                                                               >  2:473966/1‑251 (MQ=255)              cTATGTTTATCCGCAATATCGGCTGGGGAATATGCACCAGCAAACAATTGCAGAAATGATCGATTCCCCGCAACAGCAGGCGTTTGG‑TGAAGATAAATTTAAGCAGTTACCGGC‑GCAGTGTCGCAGTTGTAACGTGTTAAAAG‑CGTGCTGGGGAGGCTGCCCGAAACACC‑GCTTCATGCTCGATGCCAGCGGCAAACCGGGA‑CTGAATTATTT‑GTGTGCCGGGTATCAGCGTTATTTCCGCCATCTACCGc                                                                                                                                                                                                                                             >  2:275167/1‑251 (MQ=255)                  gTTTATCCGCAATATCGGCTGGGGAATATGCACCAGCAAACAATTGCAGAAATGATCGATTCCCCGCAACAGCAGGCGTTTGG‑TGAAGATAAATTTAAGCAGTTACCGGC‑GCAGTGTCGCAGTTGTAACGTGTTAAAAG‑CGTGCTGGGGAGGCTGCCCGAAACACC‑GCTTCATGCTCGATGCCAGCGGCAAACCGGGA‑CTGAATTATTT‑GTGTGCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCat                                                                                                                                                                                                                                          >  1:569819/1‑250 (MQ=255)                   tttATCCGCAATATCGGCTGGGGAATATGCACCAGCAAACAATTGCAGAAATGATCGATTCCCCGCAACAGCAGGCGTTTGG‑TGAAGATAAATTTAAGCAGTTACCGGC‑GCAGTGTCGCAGTTGTAACGTGTTAAAAG‑CGTGCTGGGGAGGCTGCCCGAAACACC‑GCTTCATGCTCGATGCCAGCGGCAAACCGGGA‑CTGAATTATTT‑GTGTGCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCatat                                                                                                                                                                                                                                        >  1:244997/1‑251 (MQ=255)                                     tGGGGAATATGCACCAGCAAACAATTGCAGAAATGATCGATTCCCCGCAACAGCAGGCGTTTGG‑TGAAGATAAATTTAAGCAGTTACCGGC‑GCAGTGTCGCAGTTGTAACGTGTTAAAAG‑CGTGCTGGGGAGGCTGCCCGAAACACC‑GCTTCATGCTCGATGCCAGCGGCAAACCGGGA‑CTGAATTATTT‑GTGTGCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGAt                                                                                                                                                                                                                      >  2:283286/1‑251 (MQ=255)                                      ggggAATATGCACCAGCAAACAATTGCAGAAATGATCGATTCCCCGCAACAGCCGGCGTTTGG‑TGAAGATAAATTTAAGCAGTTACCGGC‑GCAGTGTCGCAGTTGTAACGTGTTAAAAG‑CGTGCTGGGGAGGCTGCCCGAAACACC‑GCTTCATGCTCGATGCCAGCGGCAAACCGGGA‑CTGAATTATTT‑GTGTGCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGAtt                                                                                                                                                                                                                     >  1:492928/1‑251 (MQ=255)                                      ggggAATATGCACCAGAAAAAAATTGCAGAAATGATCGATTCCCCGCAACAGC‑GGGGTTTGGGTGAAAAAAAATTAAAGCATTTACCGGC‑GCAGTGCCGCAGTGGTAACGTTTAAAAAG‑CGTGCTGGGGAGGCTGCCGAAAACCCC‑GCTTCATGCTCGATGCCAGCGGCAAACCGGGA‑CTGAATTTTTT‑GTGTGCGGGGTATCAGCGTTATTCCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGAtt                                                                                                                                                                                                                     <  1:283286/251‑1 (MQ=255)                                      ggagAATACGCACAAGCAAACAATTGCAGAAATGATCGTTCCCCCGCACCAGCAGGCTTTGGG‑TGAAGAAAAATTTAGGCATTGACGGGC‑GCGGTGTCGGATTTGTAACGTGTAAAAAG‑CGTGCTGGGGAGGCTCCCAAAAACCCC‑GCTTCATGCTCGATGCCGGCGGCAGCCCGGGA‑GTAAATTTTTT‑GTGTGCCGGGTACGAGTTTAATTCCCCCCATCTACCGCCATATCTTAAAGCAATGGCTGAtt                                                                                                                                                                                                                     <  2:373707/248‑1 (MQ=255)                                             atGCACCAGCAAACAATTGCAGAAATGATCGATTCCCCGCAACAGCAGGCGTTTGG‑TGAAGATAAATTTAAGCAGTTCCCGGC‑GCAGTGTCGCAGTTGTAACGTTTTAAAAG‑CGTGCTGGGGAGGCTGCCCGAAACCCC‑GCTTCATGCTCGATGCCAGCGGCAACCCGGGA‑CTGAATTATTT‑GTGGGCCGGGTATCAGCGTTATTCCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGgc                                                                                                                                                                                                              <  1:473966/251‑1 (MQ=255)                                                        aaaaatggCAGAAATGATCGTTTCCCCGCAACACCGGGCGTTTGG‑TAAAGATAAATTTAAGCATTTCCCGGC‑GCGGTGTCGCAGTTTTAACGTGTAAAAAG‑CGTGCTGGGGGGGCTCCCCGAACCACC‑GCTTCATGCGCGTTGCCAGCGGCAAACCGGGA‑CTGAATTATTT‑GTGTGCCGGGTATCAGCTTTATTCCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGcc                                                                                                                                                                                                   <  1:175999/244‑1 (MQ=255)                                                            aTTGCAGAAATGATCGATTCCCCGCAACAGCAGGCGTTTGG‑TGAAGATAAATTTAAGCAGTTACCGGC‑GCAGTGTCGCAGTTGTAACGTGTTAAAAG‑CGTGCTGGGGAGGCTGCCCGAAACACC‑GCTTCATGCTCGATGCCAGCGGCAAACCGGGA‑CTGAATTATTT‑GTGTGCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGc                                                                                                                                                                                               >  2:124785/1‑251 (MQ=255)                                                               gaaGAATTGATCGATTCCCCGCAACGGCAGGCGTTTGG‑TGAAGAAAATTTTAAGCATTTCCGGGC‑GCAGTGCCGCAGTTGTAACTTGTAAAAAG‑CGGGCGGGGGAGGCTCCCCGAACCCCC‑GTTTCATGCTCGATCCGAGCGGCAACCGGGGA‑TTGAATTTTTT‑TTGTGCGGGGTATAAGCTTTTTGTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCa                                                                                                                                                                                             <  1:229458/248‑1 (MQ=255)                                                                    aaTGATCGATTCCCCGCAACAGCAGGCGTTTGG‑TGAAGATAAATTTAAGCAGTTACCGGC‑GCAGTGTCGCAGTTGTAACGTGTTAAAAG‑CGTGCTGGGGAGGCTGCCCGAAACACC‑GCTTCATGCTCGATGCCAGCGGCAAACCGGGA‑CTGAATTATTT‑GTGTGCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGAc                                                                                                                                                                                        >  1:304943/1‑250 (MQ=255)                                                                     aTGATCGATCCCCCCCAACTGTGGGCGGTTGG‑AAAAGAAAAATTTAAGC‑GTTTCCGGCCGCAGTGTCGCATTTGAAAGGT‑TTAAAAGGCGTGCGGGGGGGGCGGCCCAAAACCCC‑GCTTCAGGCCCGATGCCGGCGGAAAACCGGTA‑CTAATTTATTT‑TGGGGCGGGGTATCAGCTTTATTCCCGCCTTCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACAt                                                                                                                                                                                      <  2:244997/251‑1 (MQ=255)                                                                                        agcagGCGTTTGG‑TGAAGATAAATTTAAGCAGTTACCGGC‑GCAGAGTCGCAGTTGTAACGTGTTAAAAG‑CGTGCTGGGGAGGCTGCCCGAAACACC‑GCTTCATGCTCGATGCCAGCGGCAAACCGGGA‑CTGAATTATTT‑GTGTGCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCtg                                                                                                                                                                   >  1:423617/1‑251 (MQ=255)                                                                                          cagGCGTTTGG‑TGAAGATAAATTTAAGCAGTTACCGGC‑GCAGTGTCGCAGTTGTAACGTGTTAAAAG‑CGTGCTGGGGAGGCTGCCCGAAACACC‑GCTTCATGCTCGATGCCAGCGGCAAACCGGGA‑CTGAATTATTT‑GTGTGCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCtggt                                                                                                                                                                 >  1:90953/1‑251 (MQ=255)                                                                                             gtgttttgtagaaaaTAAATTAAAGCAGTGACCGGC‑GAAGTATCGCATTTGTAACGTGTTAAAAG‑GGTGCTGGGGAGGCTGCCCGAACCCCC‑GCTTCATGCTCGATGCCAGCGGCAACCCGGGA‑CGGATTTATTG‑GTGTGCGGGTTACCGGGGTTATTCCCGCCATCTCCCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCtggtgg                                                                                                                                                               <  1:276041/237‑1 (MQ=255)                                                                                                        aaGATAAATTTAAGCAGTTACCGGC‑GCAGTGTCGCAGTTGTAACGTGTTAAAAG‑CGTGCTGGGGAGGCTGCCCGAAACACC‑GCTTCATGCTCGATGCCAGCGGCAAACCGGGA‑CTGAATTATTT‑GTGTGCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGagtaagt                                                                                                                                                      >  2:281688/1‑249 (MQ=255)                                                                                                           aTAAATTTAAGCAGTTACCGGC‑GCAGTGTCGCAGTTGTAACGTGTTAAAAG‑CGTGCTGGGGAGGCTGCCCGAAACACC‑GCTTCATGCTCGATGCCAGCGGCAAACCGGGA‑CTGAATTATTT‑GTGTGCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGaa                                                                                                                                                  >  1:148211/1‑250 (MQ=255)                                                                                                                 ttAAGCAGTTCCCGGC‑GCAGTGCCGCGTTTGAACCGTGTAAAAGG‑CGTGCGGGGGAGGCTGCCCAAAACCCC‑GCTCAATGCTCGATGCCGGCGGAAAACGGGGA‑CTAAATTATTT‑GTGTGCCGGTTATCAGCTTTATTCCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATcggc                                                                                                                                            <  2:569819/250‑1 (MQ=255)                                                                                                                      cAGTTACCGGC‑GCAGTCTCGCAGTTGTAACGTGTTAAAAG‑CGTGCTGGGGAGGCTGCCCGAAACACC‑GCTTCATGCTCGATGCCAGCGGCAAACCGGGA‑CTGAATTATTT‑GTGTGCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGAAGAAATCGGCGgccgc                                                                                                                                      >  2:465649/1‑251 (MQ=255)                                                                                                                           aCCGGG‑GCAGTGTCGCAGTTGTAACGTGTTAAAAG‑CGTGTGGGGGAGGCTGCC‑GAAACCCCCGCTTCATGCTCGGTGCCAGCGGCAACCGGGGA‑CTGAATTTTTT‑GTGTGCGGGGTACCAGCGTTTTTCCCGCCTTCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGc                                                                                                                                  <  1:108004/250‑1 (MQ=255)                                                                                                                             cGGC‑GCAGTGTCGCAGTTGTAACGTGTTAAAAG‑CGTGCTGGGGAGGCTGCCCGAAACACC‑GCTTCATGCTCGATGCCAGCGGCAAACCGGGA‑CTGAATTATTT‑GTGTGCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGt                                                                                                                               >  1:198657/1‑251 (MQ=255)                                                                                                                                c‑gcAGTGTCGCAGTTGTAACGTGTTAAAAG‑CGTGCTGGGGAGGCTGCCCGAAACACC‑GCTTCATGCTCGATGCCAGCGGCAAACCGGGA‑CTGAATTATTT‑GTGTGCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTgg                                                                                                                            >  2:231410/1‑250 (MQ=255)                                                                                                                                  gcAGTGTCGCAGTTGTAACGTGTTAAAAG‑CGTGCTGGGGAGGCTGCCCGAAACACC‑GCTTCATGCTCGATGCCAGCGGCAAACCGGGA‑CTGAATTATTT‑GTGTGCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGAt                                                                                                                           >  2:533868/1‑251 (MQ=255)                                                                                                                                          gCAGTTGTAACGTGTTAAAAG‑CGTGCTGGGGGGGCTGCCCGAAACACC‑GCTTCATGCTCGATGCCAGCGGCAACCCGGGA‑CTGAATTATTT‑GTGGGCCGGGTATCGGCGTTATTCCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGAt                                                                                                                   <  1:275167/251‑1 (MQ=255)                                                                                                                                                gTAACGTGTTAAAAG‑CGTGCTGGGGAGGCTGCCCAAAACCCC‑GCTTCATGCTCGATGCCAGCGGCAAACGGGGA‑CTGATTTATTT‑GTGTGCGGGGTATCAGCGTTATTCCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATgcggcg                                                                                                             <  1:281688/251‑1 (MQ=255)                                                                                                                                                    cGTGTTAAAAG‑CGTGCTGGGGAGGCTGCCCGAAACACC‑GCTTCATGCTCGATGCCAGCGGCAAACCGGGA‑CTGAATTATTT‑GTGTGCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGTGGCGGCGtca                                                                                                          >  2:223022/1‑248 (MQ=255)                                                                                                                                                     gtttaaaaaG‑CGTGCGGGGGAGGCTGCCCGAACCCCC‑GTTTCATGCTCGATGCCAGCGGCAAACCGGGA‑CTGAATTATTT‑TTGTGCCGGGTATCAGCGTTATTCCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAAc                                                                                                        <  1:241025/246‑1 (MQ=255)                                                                                                                                                           aaaaggggggtGGGGGGGCCGCCCGAAACCCC‑GCTTCATGCTCGATGCCAGGGGAAAACCGGGA‑C‑GAATTTTTTTGTGGGCGGGGTAACAGCGTTATTCCCCCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCtt                                                                                                   <  1:510370/241‑1 (MQ=255)                                                                                                                                                           aaaG‑CGTGCTGGGGAGGCTGCCCGAAACACC‑GCTTCATGCTCGATGCCAGCGGCAAACCGGGA‑CTGAATTATTT‑GTGTGCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCtt                                                                                                   >  1:142134/1‑250 (MQ=255)                                                                                                                                                                 gtgcgggggagccgccccGAAACACC‑GCTTCATGCTCGATGCCAGCGGCAAACCGGGA‑TTGAATTATTT‑GTGTGCGGGGTATCAGCTTTATTCCCGCCTTCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCgg                                                                                             <  1:10801/236‑1 (MQ=255)                                                                                                                                                                    gcgggggggccgcccaaaaacccctcTTCATGCGCGAGGCCGGCG‑CAAACCGGGAATTAAATTTTTT‑GTGGGCGGGGTATGGGCTTTTTTCCCGCCTTCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGcc                                                                                           <  2:394271/226‑1 (MQ=255)                                                                                                                                                                      ggggAGGCTCCCCAAAACCCC‑GTTTCATGCGCGATGCCAGCGGCAAACCGGGA‑CTGAATTATTT‑GTGTGCCGGGTATCAGCGTTATTCCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCt                                                                                          <  1:388181/249‑1 (MQ=255)                                                                                                                                                                       gggAGGCTGCCCGAAACACC‑GCTTCATGCTCGATGCCAGCGGCAAACCGGGA‑CTGAATTATTT‑GTGTGCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTAc                                                                                        >  1:333187/1‑250 (MQ=255)                                                                                                                                                                        ggAGGCGGCCAGAACCCCC‑GCTTCATTCTCGATGCCAGCGGCAATCCGGGC‑TTTATTTTTTT‑GTGTTCCGGTTATCAGCTTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAt                                                                                      <  1:508003/251‑1 (MQ=255)                                                                                                                                                                                cccGAACCCCC‑GCTTCATGCTCGATGCCAGCGGCAACCCGGGA‑CTGAATTTTTT‑GGGTGCCGGTTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGc                                                                              <  1:124785/251‑1 (MQ=255)                                                                                                                                                                                       cacC‑GCTTCATGCTCGATGCCAGCGGCAAACCGGGA‑CTGAATTATTT‑GTGTGCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGACCGTGCAAATTc                                                                        >  2:357895/1‑250 (MQ=255)                                                                                                                                                                                          ccgattCATGCTCGATGCCAGCGGCAAACCGGGA‑CGGGATTATTT‑GTGCGCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTTATTTGCTGGCGCACGGTTGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAGGTAGAATTCGGCGTCTGCCTGCGGTTGATTGCCGGATGCGGCGAAAAAGCATTATCCGGCCTACATGATAGTGCACATTctat                                                                     >  2:334678/5‑248 (MQ=255)                                                                                                                                                                                                        aTGCCAGCGGCAAACAGGGA‑ATGACTTATCT‑GTGTGACGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGAAATAACTGATTTGCTGGCGCCCGGTCGCCCGACCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGGAAGTAGAATTCGGCGGCAACCTACGGATGATTTCCGTAAGAGGCGTAAACGCCATATCCAGCCTACATGTTCGTGCAAAGTTAATAAATTGCaactgtc                                                       >  2:553203/1‑245 (MQ=255)                                                                                                                                                                                                             aGCGGCAAACCGGGA‑CTGAATTATTT‑GTGTGCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGt                                                    >  2:571977/1‑249 (MQ=255)                                                                                                                                                                                                                   aaaCCGGGA‑CTGAATTATTT‑GTGTGCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGTAGCGTTCTGTAGGCTgg                                             >  2:292188/1‑250 (MQ=255)                                                                                                                                                                                                                    aaCCGGGT‑CTGAATTATTT‑GTGTGCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATATTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGAGTTCTGTAGGCTGGAt                                           >  1:413213/1‑251 (MQ=255)                                                                                                                                                                                                                                      ttggggggctggtgttaagggTTATTCCCGCCAGCAACCGCCATATCTTAAAGCAATGTCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAgcat                            <  1:188079/230‑1 (MQ=255)                                                                                                                                                                                                                                                               ttCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGc  >  1:53733/1‑251 (MQ=255)                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                          ACGCCTGTGATCACTATGTTTATCCGCAATATCGGCTGGGGAATATGCACCAGCAAACAATTGCAGAAATGATCGATTCCCCGCAACAGCAGGCGTTTGG‑TGAAGATAAATTTAAGCAGTTACCGGC‑GCAGTGTCGCAGTTGTAACGTGTTAAAAG‑CGTGCTGGGGAGGCTGCCCGAAACACC‑GCTTCATGCTCGATGCCAGCGGCAAACCGGGA‑CTGAATTATTT‑GTGTGCCGGGTATCAGCGTTATTTCCGCCATCTACCGCCATATCTTAAAGCAATGGCTGATTTGCTGGCGCACGGTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGC  >  NC_000913/3983841‑3984339

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