Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,354,498	T→C	18.5%	L79P (CTC→CCC)	glpB →	anaerobic sn‑glycerol‑3‑phosphate dehydrogenase membrane anchor subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,354,498	0	T	C	18.5%	99.5 / 3.5	52	L79P (CTC→CCC)	glpB	anaerobic sn‑glycerol‑3‑phosphate dehydrogenase membrane anchor subunit
Reads supporting (aligned to +/- strand):  ref base T (24/16);  new base C (3/8);  total (27/25)
Fisher's exact test for biased strand distribution p-value = 8.77e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCATTATGGGCGGCGGCCTCGCCGGATTACTCTGTGGCCTGCAACTGCAAAAACACGGCCTGC‑GCTGTGCCATTGT‑CACTCGTGGTCAAAGCGCACTGCATTTCT‑CAT‑CCGGATCGCTGG‑ATTTGCTGAGCCATCTGCCAGATGGTCAACCGGTGA‑CAGACATTCACAGTGGACTGGAATCTTTGCGTCAGCAGGCACCAGCCCATCCTTACTCCCTTCTCGAGCCACAACGCGTGCTCGATCTCGCTTGCCAGGCGCAGGCATTAATCGCTGAAAGCGGTGCGCAATTGCAGGGCAGCGTAGAACTTGCTCACCAGCGGGTTACGCCGCTCGGCACTCTGCGCTCTACCTGGCTAAGTTCGCCAGAAGTCCCCGTCTGGCCGCTGCCCGCGAAGAAAATATGTGTAGTGGGAATTAGCGGCCTGATGGATTTTCAGGCGCACC  >  NC_000913/2354279‑2354722                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                 tCATTATGGGCGGCGGCCTCGCCGGATTACTCTGTGGCCTGCAACTGCAAAAACACGGCCTGC‑GCTGTGCCATTGT‑CACTCGTGGTCAAAGCGCACTGCATTTCT‑CAT‑CCGGATCGCTGG‑ATTTGCTGAGCCATCTGCCAGATGGTCAACCGGTGA‑CAGACATTCACAGTGGACTGGAATCTTTGCGTCAGCAGGCACCAGCCCATCCTTACTCCCTTCTCGAGCcac                                                                                                                                                                                                                          >  2:412750/1‑228 (MQ=255) tCATTATGGGCGGCGGCCTCGCCGGATTACTCTGTGGCCTGCAACTGCAAAAACACGGCCTGC‑GCTGTGCCATTGT‑CACTCGTGGTCAAAGCGCACTGCATTTCT‑CAT‑CCGGATCGCTGG‑ATTTGCTGAGCCATCTGCCAGATGGTCAACCGGTGA‑CAGACATTCACAGTGGACTGGAATCTTTGACTCAGCAGGCACCACCCCAACCTTAATCCCTTCTCGAGCcaa                                                                                                                                                                                                                          >  1:19853/1‑227 (MQ=255)   aTTATGGGCGGCGGCCTCGCCGGATTACTCTGTGGCCTGCAACTGCAAAAACACGGCCTGC‑GCTGTGCCATTGT‑CACTCGTGGTCAAAGCGCACTGCATTTCC‑CAT‑CCGGATCGCTGG‑ATTGGCTGAGCCATCTGCCAGATGGTCAACCGGTGA‑CAGACATTCACATTGGACTGGAATCTTTGCGTCAGCAGCCACCAGCCCGTCCTTACTCCCTTCTCGCGCcacac                                                                                                                                                                                                                        >  2:338061/1‑227 (MQ=255)      aTGGGCGGCGGCCTCGCCGGATTACTCTGTGGCCTGCAACTGCAAAAACACGGCCTGC‑GCTGTGCCATTGT‑CACTCGTGGTCAAAGCGCACTGCATTTCT‑CAT‑CCGGATCGCTGG‑CTTTGCTGAGCCATCTGCCAGATGGGCAACCCGTGA‑CAGGCATTCACCGTGGACTGGAATCTTTGGGTACACCGGCCACAGCCCAACCTAACGTCCTCCTCGAGcccaaacc                                                                                                                                                                                                                      >  1:587063/1‑221 (MQ=255)      aTGGGCGGCGGCCTCGCCGGATTACTCTGTGGCCTGCAACTGCAAAAACACGGCCTGC‑GCTGTGCCATTGT‑CACTCGTGGTCAAAGCGCACTGCATTTCT‑CAT‑CCGGATCGCTGG‑ATTTGCTGAGCCATCTGCCAGATGGTCAACCGGTGA‑CAGCCATTCACAGTGGACTGGAATCTTTGCGTCAGCAGGCACCAGCCCAGCCTTACTCCCTTCTCGAGCCCCAAcg                                                                                                                                                                                                                      >  2:196172/1‑227 (MQ=255)         ggcggcggcCTCGCCGGATTACTCTGTGGCCTGCAACTGCAAAAACACGGCCTGC‑GCTGTGCCATTGT‑CACTCGTGGTCAAAGCGCACTGCATTTCT‑CAT‑CCGGATCGCTGG‑ATTTGCTGAGCCATCTGCCAGATGGTCAACCGGTGA‑CAGCCATTCACAGTGGACTGGAATCTTTGCGTCAGCAGGCACCAGCCCATCCTTCTTCCCTTTCCGAGCCCCAAcacg                                                                                                                                                                                                                    >  2:608505/1‑226 (MQ=255)           cggcggcCTCGCCGGATTACTCTGTGGCCTGCAACTGCAAAAACACGGCCTGC‑GCTGTGCCATTGT‑CACTCGTGGTCAAAGCGCACTGCATTTCT‑CAT‑CCGGATCGCTGT‑ATTTGCTGACCCTTCTGCCAGATGGTCACCCGGTGA‑CAGACATTCACAGTGGCCTGGAACCTTTGCGTCCGCCGGCACCAGCCCATCCTTCCTCCCTTCCCGCGCCACAACCCgcgc                                                                                                                                                                                                                 >  1:575961/1‑227 (MQ=255)           cggcggcCTCGCCGGATTACTCTGTGGCCTGCAACTGCAAAAACACGGCCTGC‑GCTGTGCCATTGT‑CACTCGTGGTCAAAGCGCACTGCATTTCT‑CAT‑CCGGATCGCTGG‑ATTTGCTGAGCCATCTGCCAGATGGTCAACCGGTGA‑CAGACATTCACAGTGGACTGGAATCTTTGCGTCAGCAGGCACCAGCCCATCCTTACTCCCTTCTCGAGCCACACAGCGTGc                                                                                                                                                                                                                 >  2:139854/1‑227 (MQ=255)           cggcggcCTCGCCGGATTACTCTGTGGCCTGCAACTGCAAAAACACGGCCTGC‑GCTGTGCCATTGT‑CACTCGTGGTCAAAGCGCACTGCATTTCT‑CAT‑CCGGATCGCTGG‑ATTTCCTGAGCCATCTGCCAGATGGTCAACCGGTGA‑CAGACATTCACAGTGGCCTGGAATCTTTGCGTCAGCAGGCACCAGCCCATCCCTACTCCCTTCCCGAGCCCCAACGCCTGc                                                                                                                                                                                                                 >  1:396293/1‑227 (MQ=255)                              cctggtgGCATGCAACTGCAAAACCACGGCCTGC‑GCTGTGCCATTGT‑CACTCGTGGTCAAAGCGCACTGCATTTCT‑CAT‑CCGGATCGCTGG‑ATTTGCTGAGCCATCTGCCAGATGGTCAACCGGTGA‑CAGACATTCACAGTGGACTGGAATCTTTGCGTCAGCAGGCACCAGCCCATCCTTACTCCCTTCTCGAGCCACAACGCGTGCTCGATCTCGCTTGCCAGgcg                                                                                                                                                                                             <  2:138150/224‑1 (MQ=255)                              cctcgggcccTGCACCCGC‑AAAACGCGGGCCGCGGCTGTGGCTGTGT‑CCCTCGGGGTCAACGCCCACGTCTTTTCT‑CAT‑CGGGATCGCTGT‑ATTGGTTGAGCCATCTGCCAGATGGTCAACC‑GTGACCATACATTCCCAGTTGACTGGAATCTTTGGGTCCGCAGGCACCAGCCCATCCTTACTCCCTCCTCGAGCCACAACGCGTGCTCGATCTCGCTTGCCAGgcg                                                                                                                                                                                             <  1:661324/220‑1 (MQ=255)                              acctGGGGCCTGGAACGGAAAAAACTCGGCCTGC‑GCGGTGACATTGT‑CACCCGGGGTCAAAGCGCCCAGCATTTC‑‑CATGCCGGACCGGTGG‑AGTTTCTGAGCCATCTGCAAGTTGGTCAACCGGTGA‑CAGACATTCCCAGTGGACTGGAATCTTTGGGCCAGCAGGCCCCAGCAAATCCTTACTCCCTTCTCGAGCCACAACGCGTGCTCGATCTCGCTTGCCAGgcg                                                                                                                                                                                             <  1:564724/226‑1 (MQ=255)                                   tgGCCTGCAACTGCAAAAACACGGCCTGC‑GCTGTGCCATTGT‑CACTCGTGGTCAAAGCGCACTGCATTTCT‑CAT‑CCGGATCGCTGG‑ATTTGCTGAGCCATCTGCCAGATGGTCAACCGGTGA‑CAGCCATTCACAGTGGCCTGGAATCTTTGCGTCAGCAGGCACCAGCCCATCCTTACTCCCTTCTCGACCCCCACCGCGTGCTCGACCCCCCTTCCCAgcccccgcc                                                                                                                                                                                        >  1:77734/1‑220 (MQ=255)                                       ctgaaaccgcaaAAACCCGGCCTGC‑GCTGTGCCATTGTCCACT‑GGGGCCAAAGCGCACTGCATTTCC‑CAT‑CCGGATCGCTGG‑ATTTGCTGAGCCATCTGCCAGATGGTCAACCGGTGA‑CAGACATTCACAGTGGACTGGAATCTTTGCGTCAGCAGGCATCAGCCCACCCTTACTCCCTTCTCGAGCCACAACGCGTGCTCGATCTCGCTTGCCAGGCGCAGGCATTa                                                                                                                                                                                    <  2:465259/228‑1 (MQ=255)                                        tgcaactgcaaAAACACGGCCTGC‑GCTGTGCCATTGT‑CACTCGTGGTCAAAGCGCACTGCATTTCT‑CAT‑CCGGATCGCTGG‑ATTTGCTGAGCCATCTGCCAGATGGTCAACCGGTGA‑CAGACATTCACAGTGGCCTGGAATCTTTGCGTCAGCAGGCACCAGCCCATCCTTACTCCCTTCTCGAGCCACAACGCGTGCTCGATCTCGCTTGCCAGGCGCAGGCACTaa                                                                                                                                                                                   >  2:42875/1‑228 (MQ=255)                                                          gCCTGC‑GCTGTGCCTTTGT‑CCCTCGGGGTAAAAGCGCACCGACTTTCT‑CAT‑CCGGATCGCTGG‑ATTTGCTGAGCCATCTGCCAGATGGTCAACCGGTGA‑CAGACATTCACAGTGGACTGGAATCTTTGCGGCAGCAGGCACCAGCCCATCCTTACTCCCTTCTCGAGCCACAACGCGTGCTCGATCTCGCTTGCCAGGCGCAGGCATTAATCGCTGAAAGCGGTgcgc                                                                                                                                                                 <  2:567842/228‑1 (MQ=255)                                                           gcccggcttgtgCCTTTGT‑CCCTCGGGGTCAAAGCGCCCGGGATTTCC‑CGT‑CGGGATGGCTGG‑ATTTGCGTGGCCATCTGCCAGATGGGCAACAGGTGA‑CAGACATTCACAGTGGACTGGAATCTTTGCGTCAGCAGGCACCAGCCCATCCTTACTCCCTTCTCGAGCCACAACGCGGGCTCGAGCTCGCTTGCCAGGCGCAGGCATTAATCGCTGAAAGCGGTgcgc                                                                                                                                                                 <  1:196172/220‑1 (MQ=255)                                                                    gtgCCATTGT‑CACTCGTGGTCAAAGCGCACTGCATTTCT‑CAT‑CCGGATCGCTGG‑ATTTGCTGAGCCATCTGCCAGATGGTCAACCGGTGA‑CAGACATTCACAGTGGACTGGAATCTTTGCGTCAGCAGGCACCAGCCCATCCTTACTCCCTTCTCGAGCCACAACGCGTGCTCGATCTCGCTTTCCAGGCGCAGGCATTAATCGCTGAAAGCGCTGCGCAATTGCAg                                                                                                                                                         >  2:359729/1‑227 (MQ=255)                                                                     tgCCATTGT‑CCCTCGTGGTCAAAGCGCACTTCATTTCT‑CAT‑CCGGATCGGTGG‑ATTTGCTGAGCCATCTGCCAGATGGTCAACCGGTGA‑CAGACATTCACAGTGGACTGGAATCTTTGCGTCAGCAGGCACCAGCCCATCCTTACTCCCTTCTCGAGCCACAACGCGTGCTCGATCTCGCTTGCCAGGCGCAGGCATTAATCGCTGAAAGCGGTGCGCAATTGCAggg                                                                                                                                                       <  2:396293/228‑1 (MQ=255)                                                                      tcccatttccccgcggtgTCAAAGCGCACTGCATTTCC‑AAT‑CCGGATCGCTGG‑ATTTGCTGAGCCATCTGCCAGATGGTCAACCGGTGA‑CAGACTTTCACAGTGGACTGGAATCTTTGCGTCAGCAGGCACCAGCCCATCCTTACTCCCTTCTCGAGCCACAACGCGTGCTCGATCTCGCTTGCCAGGCGCAGGCATTAATCGCTGAAAGCGGTGCGCAATTGCAggg                                                                                                                                                       <  2:357928/211‑1 (MQ=255)                                                                         cattgctcctctggggcAAAGCGCCCGGC‑TTTCTCCAC‑CCGGAACGGTGG‑ATTTGCTGAGCCAGCTGCCAGTGTGTCAACCGGTGA‑CAGACATTCACAGTGGCCTGGAATCTTTGCCGCAGCAGGCACCAGCCCCTCCTTCCTCCCTTCTCGAGCCACAACGCGTGCTCGATCTCGCTTGCCAGGCGCAGGCATTAATCGCTGAAAGCGGTGCGCAATTGCAGGGCAg                                                                                                                                                    <  1:388426/212‑1 (MQ=255)                                                                                aCTCGTGGTCAAAGCGCACTGCATTTCT‑CAT‑CCGGATCGCTGG‑ATTTGCTGAGCCATCTGCCAGATGGTCAACCGGTGA‑CAGACATTCACAGTGGACTGGAATCTTTGCGTCAGCAGGCACCAGCCCATCCTTACTCCCTTCTCGAGCCACAACGCGTGCTCGCTCTCGCTTGCCAGGCGCAGGCCTTAATCGCTGAAAGCGGTGCGCAATTGCAGGGCAGCGTAGaa                                                                                                                                             >  2:342389/1‑228 (MQ=255)                                                                                aCTCGTGGTCAAAGCGCACTGCATTTCT‑CAT‑CCGGATCGCTGG‑ATTTGCTGAGCCATCTGCCAGATGGTCAACCGGTGA‑CAGACATTCACAGTGGACTGGAATCTTTGCGTCAGCAGGCACCAGCCCATCCTTACTCCCTTCTCGAGCCACAACGCGTGCTCGATCTCGCTTGCCAGGCGCAGGCCTTAATCGCTGAAAGCGGTGCGCAATTGCAGGGAAGCGTaaaa                                                                                                                                             >  1:394995/1‑228 (MQ=255)                                                                                aCTCGTGGTCAAAGCGCACTGCATTTCT‑CAT‑CCGGATCGCTGG‑ATTTGCTGAGCCATCTGCCAGATGGTCAACCGGTGA‑CAGACATTCACAGTGGACTGGAATCTTTGCGTCAGCAGGCACCAGCCCATCCTTACTCCCTTCTCGAGCCACAACGCGTGCTCGATCTCCCTTGCCAGGCGCAGGCCTTAACCGCTGAACGCGGGGCGCAATTCCCGGGCAGCGCAGaa                                                                                                                                             >  1:15533/1‑228 (MQ=255)                                                                                 ctggtgtTCAAAGCGCAGTGCATTTCT‑CAA‑CCGGATCGCTGG‑ATTGGCTGATCCATCTGCCAGATGGTCACCCGGTGA‑CAGACATTCACAGTGGACTGGAATCTTTGCGTCAGCAGTCCCCCGCCCATCCTTACTCCCTTCTCGAGCCACAACGCGTGCTCGATCTCGCTTGCCAGGCGCAGGCATTAATCGCTGAAAGCGGTgcgc                                                                                                                                                                 <  2:657045/204‑1 (MQ=255)                                                                                                  cTGCATTTCT‑CAT‑CCGGATCGCTGG‑ATTTGCTGAGCCATCTGCCAGATGGTCAACCGGTGA‑CAGACATTCACAGTGGACTGGAATCTTTGCGTCAGCAGGCACCAGCCCATCCTTACTCCCTTCTCGAGCCACAACGCGTGCTCGATCTCGCTTGCCAGGCGCAGGCATTAATCGCTGAAAGCGGTGCGCAATTGCAGGGCACCGTAGAACTTGCTCACCATCGGGt                                                                                                                            >  2:62231/1‑227 (MQ=255)                                                                                                        ttgtcact‑CCGGTGCGCTGGAATTTGCAGAGCACTCGCCCAGTGGCTCACCCGGTAA‑CAGACCTTCCCTGGGCACTGCACTTTTCGCGGCAGCAGCCCCCACCCCATCCGTACTCCCTTCACGAGCCCCAACGCGTGCTCGATCTCGCTTGCCAGGCGCAGGCATTAATCGCTGAAAGCGGTGCGCAATTGCAGGGCAGCGTAGAACTTGCTCCCCAGCGGGTTAcg                                                                                                                        <  1:26961/220‑1 (MQ=255)                                                                                                              aT‑CCGGATCGCTGG‑ATTTGCTGAGCCATCTGCCAGATGGTCAACCGGTGA‑CAGACATTCACAGTGGCCTGGAATCTTTGCGTCAGCAGGCACCAGCCCATCCTTACTCCCTTCTCGAGCCACAACGCGTGCCCGATCTCGCTTGCCAGGCGCATGCCTTAATCGCTGAACGCGGTGCGCAATTGCAGGGCCGCACAGAATTTGCCCCCCACCGGGttgacccgcccg                                                                                                                 >  1:271534/1‑217 (MQ=255)                                                                                                              aT‑CCGGATCGCTGG‑ATTTGCTGAGCCATCTGCCAGATGGTCAACCGGTGA‑CAGACATTCACAGTGGACTGGAATCTTTGCGTCAGCAGGCACCAGCCCATCCTTACTCCCTTCTCGAGCCCCAACGCGTGCTCGATCTCGCTTGCCAGGCGCAGGCATTAATCGCTGAAAGCGGTGCGCACTTGCAGGGCAGCGTAGAACTTGCTCACCAGCGGGTTAggcacgccg                                                                                                                 >  2:165751/1‑221 (MQ=255)                                                                                                                          ggg‑ATTTGCTGGGCCATCTGCCAGATGGGCAACCGGTGA‑CAGACATTCACAGTGGACTGGAATCTTTGCGTCAGCAGGCACCAGCCCATCCTTACTCCCTTCTCGAGCCACAACGCGTGCTCGATCTCGCTTGCCAGGCGCAGGCATTAATCGCTGAAAGCGGTGCGCAATTGCAGGGCAGCGTAGAACTTGCTCACCAGCGGGTTACGCCGCTCGGCACTCTGCGct                                                                                                     <  2:77734/227‑1 (MQ=255)                                                                                                                                        ccATGTGCCAGGGGGTAACCCGTCGA‑CCGACCTTCCGCGTGGACGGGTATCTGTGCCTCAGTAGGCCCCAACCCCTCCTTACCCCCGGCTCGAGCCACAACGCGGGGTCGATCTCGCTTGCCAGGCGCAGGCATTAATCGCTGAAAGCGGTGCGCAATTGCAGGGCAGCGTAGAACTTGCTCACCAGCGGGTTACGCCGCTCGGCACTCTGCGCTCTACCTGGCTAAg                                                                                        <  1:608505/228‑1 (MQ=255)                                                                                                                                                cggttGGTCACCCGGGGA‑CAGCCCTTCCCGGGGGGTTGGAATCTTGGGGTCGTGCGGCCCCGGCCCTTCCTTCCCCCCTTCCCGAGCCACAACGCGTGCTCGATCTCGCTTGCCAGGCGCAGGCATTAATCGCTGAAAGCGGTGCGCAATTGCAGGGCAGCGTAGAACTTGCTCACCAGCGGGTTACGCCGCTCGGCACTCTGCGCTCTACCTGGCTAAGTTCGCCAg                                                                                <  1:42875/224‑1 (MQ=255)                                                                                                                                                                            ccaGTGGACTGGAATCTTTGCGGCAGCCGGCCCCACCCCATCCTTACGTCCTTCCCGAGCCTCAACGCGTGCTCGATCTCGCTTGCCAGGCGCAGGCATTAATCGCTGAAAGCGGTGCGCAATTGCAGGGCAGCGTAGAACTTGCTCACCAGCGGGTTACGCCGCTCGGCACTCTGCGCTCTACCTGGCTAAGTTCGCCAGAAGTCCCCGTCTGGCCGCTGCCcgcg                                                      <  2:587063/226‑1 (MQ=255)                                                                                                                                                                                tGGACTGGAATCTTTGCGTCAGCAGGCACCAGCCCATCCTTACTCCCTTCTCGAGCCACAACGCGTGCTCGATCTCGCTTGCCAGGCGCAGGCATTAATCGCTGAAAGCGGTGCGCAATTGCAGGGCAGCGTAGAACTTGCTCACCAGCGGGTTACGCCGCTCGGCACTCTGCGCTCTACCTGGCTAAGTTCGCCAGAAGTCCCCGCCTGGCCGCTGCCCCCgagcaa                                                 >  2:200362/1‑224 (MQ=255)                                                                                                                                                                                  gACTGGAATCTTTGCGTCAGCAGGCACCAGCCCATCCTTACTCCCTTCTCGAGCCACAACGCGTGCTCGATCTCGCTTGCCAGGCGCAGGCATTAATCGCTGAAAGCGGTGCGCAATTGCAGGGCAGCGTAGAACTTGCTCACCAGCGGGTTACGCCGCTCGGCACTCTGCGCTCTACCTGGCTAAGTTCGCCCGAAGTCCCCGTCTGGCCGCTGCCCGCGACGaaaa                                               >  2:409130/1‑228 (MQ=255)                                                                                                                                                                                  gACTGGAATCTTGGCGTCAGCAGGCACCAGCCCATCCTTACTCCCTTCCCGAGCCACAACGCGTGCTCGATCTCGCTTGCCAGGCGCAGGCATTAATCGCTGAAAGCGGTGCGCAATTGCAGGGCAGCGTAGAACTTGCTCACCAGCGGGTTACGCCGCTCGGCACTCTGCGCTCTACCTGGCTAAGTTCGCCAGAAGTCCCCGTCTGGCCGCTGCCCGCgaagaa                                                 <  2:394995/226‑1 (MQ=255)                                                                                                                                                                                   cttGGAATTTGTGCGGCAGGGGGCCCCCGCCCCTCCTTACCCCCCTCTCGAGCCACAACGCGTGGCCGACCTCGCTTGCCAGGCGCAGGCATTAAGCGCTGAAAGCGGTGCGCAATTGCAGGGCAGCGTAGAACTTGCTCACCAGCGGGTTACGCCGCTCGGCACTCTGCGCTCTACCTGGCTAAGTTCGCCAGAAGTCCCCGTCTGGCCGCTGCCCGCGAAGAAAat                                              <  1:342389/226‑1 (MQ=255)                                                                                                                                                                                   aCTGGAATCTTTGCGTCAGCAGGCCCCGGCCCATCCTTACTCCCATCCCGAGCCAAAACGCGTGCTCGATCTCGCTTGCCAGGCGCAGGCATTAATCGCTGAAAGCGGTGCGCAATTGCAGGGCAGCGTAGAACTTGCTCACCAGCGGGTTACGCCGCTCGGCACTCTGCGCTCTACCTGGCTAAGTTCGCCAGAAGTCCCCGTCTGGCCGCTGCCCGCGAAGAAAat                                              <  2:575961/228‑1 (MQ=255)                                                                                                                                                                                     tGGAATCTTTGCGTCAGCAGGCACCAGCCCATCCTTACTCCCTTCTCGAGCCACAACGCGTGCTCGATCTCGCTTGCCAGGCGCAGGCATTAATCGCTGAAAGCGGTGCGCAATTGCAGGGCAGCGTAGAACTTGCTCACCAGCGGGTTACGCCGCCCGGCACCCTCCGCTCTACCTGGCTCAGTCCGCCCGAGGCCCCCCTCCGGCCGCTCCCCCCCAACAAAaaat                                            >  1:40409/1‑228 (MQ=255)                                                                                                                                                                                       gAATCTTTGCGTCAGCAGCCACCAGCCCATCCTTACTCCCTTCTCGAGCCACAACGCGTGCTCGATCTCGCTTGCCAGGCCCAGGCATTAATCGCTGAAAGCGGTGCGCAATTCCAGGGCAGCGTAGAACTTGCTCACCAGCGGGGTACGCCGCTCGGCACTCTGCGCTCTACCTGCGTAAGTTCGCCAGAAGGCCCCGGCCGGCCGcgcccccagccgcaaatatgt                                          >  1:432531/1‑208 (MQ=255)                                                                                                                                                                                        aaTCTTTGCGTCAGCAGGCACCAGCCCATCCTTACCCCCTTCCCGAGCCACAACGCGTGCTCGATCTCGCTTGCCAGGCGCAGGCATTAATCGCTGAAAGCGGTGCGCAATTGCAGGGCAGCGTAGAACTTGCTCACCAGCGGGTTACGCCGCTCGGCACTCTGCGCTCTACCTGGCTAAGTTCGCCAGAAGTCCCCGTCTGGCCGCTGCCCGCGAAGAAAATAtgtg                                         <  2:20485/228‑1 (MQ=255)                                                                                                                                                                                            tttGCGTCAGCAGGCACCAGCCCATCCTTACTCCCTTCTCGAGCCACAACGCGTGCTCGATCTCGCTTGCCAGGCGCAGGCATTAATCGCTGAAAGCGGTGCGCAATTGCAGGGCAGCGTAGAACTTGCTCACCAGCGGGTTACGCCGCTCGGCACTCTGCGCTCTACCTGGCTAAGTTCGCCAGAAGTCCCCGTCCGGCCGCTGCCCGCGACGAAAATTTGGGTAg                                      >  2:217094/1‑227 (MQ=255)                                                                                                                                                                                            tttGCGTCAGCAGGCACCAGCCCATCCTTACTCCCTTCTCGAGCCACAACGCGTGCTCGATCTCGCTTGCCAGGCGCAGGCATTAATCGCTGAAAGCGGTGCGCAATTGCAGGGCAGCGTAGAACTTGCTCACCAGCGGGTTACGCCGCTCGGCACTCTGCGCTCTACCTGGCTAAGTTCGCCAGAAGTCCCCGCCTGGCCGCTGCCCGCGAAGAAAATAtgggcagg                                     >  1:56056/1‑222 (MQ=255)                                                                                                                                                                                                    atcagGCGCCAGCCCATCCTTACTCCCTTCTCGAGCCACAACGCGTGCTCGATCTCGCTTGCCAGGCGCAGGCATTAATCGCTGAAAGCGGTGCGCAATTGCAGGGCAGCGTAGAACTTGCTCACCAGCGGGTTACGCCGCTCGGCACTCTGCGCTCTACCTGGCTAAGTTCGCCAGAAGTCCCCGTCTGGCCGCTGCCCGCGAAGAAAATATGTGTAGTGGGAATTa                             <  2:271534/226‑1 (MQ=255)                                                                                                                                                                                                         gcccccgccctcctcttcccccTTGCCGAGCCCCAACGCGTGGTCGATCTCGCTTGCCAGGCGCAGGCATTAACCGCTGAAGGCGGTGCGCAGTGGCAGGGCAGCGTAGAACTTGCTCACCAGCGGGTTACGCCGCTCGGCACTCTGCGCTCTACCTGGCTAAGTTCGCCAGAAGTCCCCGTCTGGCCGCTGCCCGCGAAGAAAATATGTGTAGTGGGAATTAGCGGc                        <  1:165751/209‑1 (MQ=255)                                                                                                                                                                                                          acagaGCCCATCCTTACACCCTTCCCGAGCCACAACGCGTGCTCGATCTCGCTTGCCAGGCGCAGGCATTAATCGCTGAAAGCGGTGCGCAATTGCAGGGCAGCGTAGAACTTGCTCACCAGCGGGTTACGCCGCTCGGCACTCTGCGCTCTACCTGGCTAAGTTCGCCAGAAGTCCCCGTCTGGCCGCTGCCCGCGAAGAAAATATGTGTAGTGGGAATTAGCGGc                        <  2:116139/223‑1 (MQ=255)                                                                                                                                                                                                                 ctctctctccTCCTTTCTGGAGCCACAACGGGGGCGCGATCTCGCTTGCCAGGCGCAGGCATTAATCGCTGAAAGCGGTGCGCAATTGCAGGGCAGCGTAGAACTTGCTCACCAGCGGGTTACGCCGCTCGGCACTCTGCGCTCTACCTGGCTAAGTTCGCCAGAAGTCCCCGTCTGGCCGCTGCCCGCGAAGAAAATATGTGTAGTGGGAATTAGCGGCCTGAt                   <  1:62231/216‑1 (MQ=255)                                                                                                                                                                                                                     ccTTACTCCCTTCTCGAGCCACAACGCGTGCTCGATCTCGCTTGCCAGGCGCAGGCATTAATCGCTGAAAGCGGTGCGCAATTGCAGGGCAGCGTAGAACTTGCTCACCAGCGGGTTACGCCGCTCGGCACTCTGCGCTCTACCTGGCTAAGTTCGCCAGAAGTCCCCGTCTGGCCGCTGCCCGCGAAGCAAATATGTGTAGTGGGAATTCGCGGCCTGATGGAtttt            >  2:121356/1‑228 (MQ=255)                                                                                                                                                                                                                         aCTCCCTTCTCGAGCCACAACGCGTGCTCGATCTCGCTTGCCAGGCGCAGGCATTAATCGCTGAAAGCGGTGCGCAATTGCAGGGCAGCGTAGAACTTGCTCACCAGCGGGGTACGCCGCTCGGCACTCTGCGCTCTACCTGGCTAAGTTCCCCCGAAGTCCCCGCCTGGCCGCTGCCCGCGCAGAAAATATGTGTAGTGGGAATTAGCCGCCGGATGGGtttctgcc        >  1:116139/1‑223 (MQ=255)                                                                                                                                                                                                                           ccccTTCCTGAGCCCCAAGGCTTGCGCTATCTCGCTTGCCAGGCGCAGGCATTAAGCGCTGAAAGGGGTGCCCAATTGCAGGGCAGCGTAGAACTTGCTCACCAGCGGGTTACGCCGCTCGGCACTCTGCGCTCTTCCTGGCTAAGTTCGCCAGAAGTCCCCGTCTGGCCGCTGCCCGCGAAGAAAATATGTGTAGTGGGAATTAGCGGCCTGATGGATTTTCAGgcg      <  2:683917/227‑1 (MQ=255)                                                                                                                                                                                                                            cccTTCTCGAGCCACAACGCGTGCTCGATCTCGCTTGCCAGGCGCAGGCATTAATCGCTGAAAGCGGTGCGCAATTGCAGGGCAGCGTAGAACTTGCTCACCAGCGGGTTACGCCGCTCGGCACTCTGCGCTCTACCTGGCTAAGTTCGCCAGAAGTCCCCGTCTGGCCGCTGCCCGCGAAGAAAATATGTGTAGTGGGAATTAGCGGCCTGATGGATTTTCAGgcg      <  2:40409/227‑1 (MQ=255)                                                                                                                                                                                                                               ttCTCGAGCCACAACGCGTGCTCGATCTCGCTTGCCAGGCGCAGGCATTAATCGCTGAAAGCGGTGCGCAATTGCAGGGCAGCGTAGAACTTGCTCACCAGCGGGTTACGCCGCTCGGCACTCTGCGCTCTACCTGGCTAAGTTCGCCAGAAGTCCCCGTCTGGCCGCTGCCCGCGAAGAAAAAATGTGTAGTGGGAATTAGCGGCCTGATGGCATTTCCGGCGcccc  >  2:28425/1‑228 (MQ=255)                                                                                                                                                                                                                                tctcGAGCCACAACGCGTGCTCGATCTCGCTTGCCAGGCGCAGGCATTAATCGCTGAAAGCGGTGCGCAATTGCAGGGCAGCGTAGAACTTGCTCACCAGCGGGTTACGCCGCTCGGCACTCTGCGCTCTACCTGGCTAAGTTCGCCAGCAGTCCCCGTCTGGCCGCTGCCCGCgaaga                                                  >  2:561678/1‑179 (MQ=255)                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                 TCATTATGGGCGGCGGCCTCGCCGGATTACTCTGTGGCCTGCAACTGCAAAAACACGGCCTGC‑GCTGTGCCATTGT‑CACTCGTGGTCAAAGCGCACTGCATTTCT‑CAT‑CCGGATCGCTGG‑ATTTGCTGAGCCATCTGCCAGATGGTCAACCGGTGA‑CAGACATTCACAGTGGACTGGAATCTTTGCGTCAGCAGGCACCAGCCCATCCTTACTCCCTTCTCGAGCCACAACGCGTGCTCGATCTCGCTTGCCAGGCGCAGGCATTAATCGCTGAAAGCGGTGCGCAATTGCAGGGCAGCGTAGAACTTGCTCACCAGCGGGTTACGCCGCTCGGCACTCTGCGCTCTACCTGGCTAAGTTCGCCAGAAGTCCCCGTCTGGCCGCTGCCCGCGAAGAAAATATGTGTAGTGGGAATTAGCGGCCTGATGGATTTTCAGGCGCACC  >  NC_000913/2354279‑2354722

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