Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,354,500	T→C	18.6%	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,500	0	T	C	18.6%	91.3 / 3.4	48	L79P (CTC→CCC)	glpB	anaerobic sn‑glycerol‑3‑phosphate dehydrogenase membrane anchor subunit
Reads supporting (aligned to +/- strand):  ref base T (23/14);  new base C (3/7);  total (26/22)
Fisher's exact test for biased strand distribution p-value = 8.64e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCATTATGGGCGGCGGCCTCGCCGGATTACTCTGTGGCCTGCAACTGCAAAAACACGGCCTGC‑GCTGTGCCATTGT‑CACTCGTGGTCAAAGCGCACTGCATTTCTCAT‑CCGGATCGCTGG‑ATTTGCTGAGCCATCTGCCAGATGGTCAACCGGTGA‑CAGACATTCACAGTGGACTGGAATCTTTGCGTCAGCAGGCACCAGCCCATCCTTACTCCCTTCTCGAGCCACAACGCGTGCTCGATCTCGCTTGCCAGGCGCAGGCATTAATCGCTGAAAGCGGTGCGCAATTGCAGGGCAGCGTAGAACTTGCTCACCAGCGGGTTACGCCGCTCGGCACTCTGCGCTCTACCTGGCTAAGTTCGCCAGAAGTCCCCGTCTGGCCGCTGCCCGCGAAGAAAATATGTGTAGTGGGAATTAGCGGCCTGATGGATTTTCAGGCGCACC  >  NC_000913/2354281‑2354724                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                 tCATTATGGGCGGCGGCCTCGCCGGATTACTCTGTGGCCTGCAACTGCAAAAACACGGCCTGC‑GCTGTGCCATTGT‑CACTCGTGGTCAAAGCGCACTGCATTTCTCAT‑CCGGATCGCTGG‑ATTTGCTGAGCCATCTGCCAGATGGTCAACCGGTGA‑CAGACATTCACAGTGGACTGGAATCTTTGCGTCAGCAGGCACCAGCCCATCCTTACTCCCTTCTCGAGCcac                                                                                                                                                                                                                          >  1:412750/1‑228 (MQ=255) tCATTATGGGCGGCGGCCTCGCCGGATTACTCTGTGGCCTGCAACTGCAAAAACACGGCCTGC‑GCTGTGCCATTGT‑CACTCGTGGTCAAAGCGCACTGCATTTCTCAT‑CCGGATCGCTGG‑ATTTGCTGAGCCATCTGCCAGATGGTCAACCGGTGA‑CAGACATTCACAGTGGACTGGAATCTTTGACTCAGCAGGCACCACCCCAACCTTAATCCCTTCTCGAGCcaa                                                                                                                                                                                                                          >  2:19853/1‑227 (MQ=255)   aTTATGGGCGGCGGCCTCGCCGGATTACTCTGTGGCCTGCAACTGCAAAAACACGGCCTGC‑GCTGTGCCATTGT‑CACTCGTGGTCAAAGCGCACTGCATTTCCCAT‑CCGGATCGCTGG‑ATTGGCTGAGCCATCTGCCAGATGGTCAACCGGTGA‑CAGACATTCACATTGGACTGGAATCTTTGCGTCAGCAGCCACCAGCCCGTCCTTACTCCCTTCTCGCGCcacac                                                                                                                                                                                                                        >  1:338061/1‑227 (MQ=255)      aTGGGCGGCGGCCTCGCCGGATTACTCTGTGGCCTGCAACTGCAAAAACACGGCCTGC‑GCTGTGCCATTGT‑CACTCGTGGTCAAAGCGCACTGCATTTCTCAT‑CCGGATCGCTGG‑CTTTGCTGAGCCATCTGCCAGATGGGCAACCCGTGA‑CAGGCATTCACCGTGGACTGGAATCTTTGGGTACACCGGCCACAGCCCAACCTAACGTCCTCCTCGAGcccaaacc                                                                                                                                                                                                                      >  2:587063/1‑221 (MQ=255)      aTGGGCGGCGGCCTCGCCGGATTACTCTGTGGCCTGCAACTGCAAAAACACGGCCTGC‑GCTGTGCCATTGT‑CACTCGTGGTCAAAGCGCACTGCATTTCTCAT‑CCGGATCGCTGG‑ATTTGCTGAGCCATCTGCCAGATGGTCAACCGGTGA‑CAGCCATTCACAGTGGACTGGAATCTTTGCGTCAGCAGGCACCAGCCCAGCCTTACTCCCTTCTCGAGCCCCAAcg                                                                                                                                                                                                                      >  1:196172/1‑227 (MQ=255)         ggcggcggcCTCGCCGGATTACTCTGTGGCCTGCAACTGCAAAAACACGGCCTGC‑GCTGTGCCATTGT‑CACTCGTGGTCAAAGCGCACTGCATTTCTCAT‑CCGGATCGCTGG‑ATTTGCTGAGCCATCTGCCAGATGGTCAACCGGTGA‑CAGCCATTCACAGTGGACTGGAATCTTTGCGTCAGCAGGCACCAGCCCATCCTTCTTCCCTTTCCGAGCCCCAAcacg                                                                                                                                                                                                                    >  1:608505/1‑226 (MQ=255)           cggcggcCTCGCCGGATTACTCTGTGGCCTGCAACTGCAAAAACACGGCCTGC‑GCTGTGCCATTGT‑CACTCGTGGTCAAAGCGCACTGCATTTCTCAT‑CCGGATCGCTGT‑ATTTGCTGACCCTTCTGCCAGATGGTCACCCGGTGA‑CAGACATTCACAGTGGCCTGGAACCTTTGCGTCCGCCGGCACCAGCCCATCCTTCCTCCCTTCCCGCGCCACAACCCgcgc                                                                                                                                                                                                                 >  2:575961/1‑227 (MQ=255)           cggcggcCTCGCCGGATTACTCTGTGGCCTGCAACTGCAAAAACACGGCCTGC‑GCTGTGCCATTGT‑CACTCGTGGTCAAAGCGCACTGCATTTCTCAT‑CCGGATCGCTGG‑ATTTGCTGAGCCATCTGCCAGATGGTCAACCGGTGA‑CAGACATTCACAGTGGACTGGAATCTTTGCGTCAGCAGGCACCAGCCCATCCTTACTCCCTTCTCGAGCCACACAGCGTGc                                                                                                                                                                                                                 >  1:139854/1‑227 (MQ=255)           cggcggcCTCGCCGGATTACTCTGTGGCCTGCAACTGCAAAAACACGGCCTGC‑GCTGTGCCATTGT‑CACTCGTGGTCAAAGCGCACTGCATTTCTCAT‑CCGGATCGCTGG‑ATTTCCTGAGCCATCTGCCAGATGGTCAACCGGTGA‑CAGACATTCACAGTGGCCTGGAATCTTTGCGTCAGCAGGCACCAGCCCATCCCTACTCCCTTCCCGAGCCCCAACGCCTGc                                                                                                                                                                                                                 >  2:396293/1‑227 (MQ=255)                              cctggtgGCATGCAACTGCAAAACCACGGCCTGC‑GCTGTGCCATTGT‑CACTCGTGGTCAAAGCGCACTGCATTTCTCAT‑CCGGATCGCTGG‑ATTTGCTGAGCCATCTGCCAGATGGTCAACCGGTGA‑CAGACATTCACAGTGGACTGGAATCTTTGCGTCAGCAGGCACCAGCCCATCCTTACTCCCTTCTCGAGCCACAACGCGTGCTCGATCTCGCTTGCCAGgcg                                                                                                                                                                                             <  1:138150/224‑1 (MQ=255)                              cctcgggcccTGCACCCGC‑AAAACGCGGGCCGCGGCTGTGGCTGTGT‑CCCTCGGGGTCAACGCCCACGTCTTTTCTCAT‑CGGGATCGCTGT‑ATTGGTTGAGCCATCTGCCAGATGGTCAACC‑GTGACCATACATTCCCAGTTGACTGGAATCTTTGGGTCCGCAGGCACCAGCCCATCCTTACTCCCTCCTCGAGCCACAACGCGTGCTCGATCTCGCTTGCCAGgcg                                                                                                                                                                                             <  2:661324/220‑1 (MQ=255)                              acctGGGGCCTGGAACGGAAAAAACTCGGCCTGC‑GCGGTGACATTGT‑CACCCGGGGTCAAAGCGCCCAGCATTTC‑CATGCCGGACCGGTGG‑AGTTTCTGAGCCATCTGCAAGTTGGTCAACCGGTGA‑CAGACATTCCCAGTGGACTGGAATCTTTGGGCCAGCAGGCCCCAGCAAATCCTTACTCCCTTCTCGAGCCACAACGCGTGCTCGATCTCGCTTGCCAGgcg                                                                                                                                                                                             <  2:564724/226‑1 (MQ=255)                                   tgGCCTGCAACTGCAAAAACACGGCCTGC‑GCTGTGCCATTGT‑CACTCGTGGTCAAAGCGCACTGCATTTCTCAT‑CCGGATCGCTGG‑ATTTGCTGAGCCATCTGCCAGATGGTCAACCGGTGA‑CAGCCATTCACAGTGGCCTGGAATCTTTGCGTCAGCAGGCACCAGCCCATCCTTACTCCCTTCTCGACCCCCACCGCGTGCTCGACCCCCCTTCCCAgcccccgcc                                                                                                                                                                                        >  2:77734/1‑220 (MQ=255)                                       ctgaaaccgcaaAAACCCGGCCTGC‑GCTGTGCCATTGTCCACT‑GGGGCCAAAGCGCACTGCATTTCCCAT‑CCGGATCGCTGG‑ATTTGCTGAGCCATCTGCCAGATGGTCAACCGGTGA‑CAGACATTCACAGTGGACTGGAATCTTTGCGTCAGCAGGCATCAGCCCACCCTTACTCCCTTCTCGAGCCACAACGCGTGCTCGATCTCGCTTGCCAGGCGCAGGCATTa                                                                                                                                                                                    <  1:465259/228‑1 (MQ=255)                                        tgcaactgcaaAAACACGGCCTGC‑GCTGTGCCATTGT‑CACTCGTGGTCAAAGCGCACTGCATTTCTCAT‑CCGGATCGCTGG‑ATTTGCTGAGCCATCTGCCAGATGGTCAACCGGTGA‑CAGACATTCACAGTGGCCTGGAATCTTTGCGTCAGCAGGCACCAGCCCATCCTTACTCCCTTCTCGAGCCACAACGCGTGCTCGATCTCGCTTGCCAGGCGCAGGCACTaa                                                                                                                                                                                   >  1:42875/1‑228 (MQ=255)                                                          gCCTGC‑GCTGTGCCTTTGT‑CCCTCGGGGTAAAAGCGCACCGACTTTCTCAT‑CCGGATCGCTGG‑ATTTGCTGAGCCATCTGCCAGATGGTCAACCGGTGA‑CAGACATTCACAGTGGACTGGAATCTTTGCGGCAGCAGGCACCAGCCCATCCTTACTCCCTTCTCGAGCCACAACGCGTGCTCGATCTCGCTTGCCAGGCGCAGGCATTAATCGCTGAAAGCGGTgcgc                                                                                                                                                                 <  1:567842/228‑1 (MQ=255)                                                           gcccggcttgtgCCTTTGT‑CCCTCGGGGTCAAAGCGCCCGGGATTTCCCGT‑CGGGATGGCTGG‑ATTTGCGTGGCCATCTGCCAGATGGGCAACAGGTGA‑CAGACATTCACAGTGGACTGGAATCTTTGCGTCAGCAGGCACCAGCCCATCCTTACTCCCTTCTCGAGCCACAACGCGGGCTCGAGCTCGCTTGCCAGGCGCAGGCATTAATCGCTGAAAGCGGTgcgc                                                                                                                                                                 <  2:196172/220‑1 (MQ=255)                                                                    gtgCCATTGT‑CACTCGTGGTCAAAGCGCACTGCATTTCTCAT‑CCGGATCGCTGG‑ATTTGCTGAGCCATCTGCCAGATGGTCAACCGGTGA‑CAGACATTCACAGTGGACTGGAATCTTTGCGTCAGCAGGCACCAGCCCATCCTTACTCCCTTCTCGAGCCACAACGCGTGCTCGATCTCGCTTTCCAGGCGCAGGCATTAATCGCTGAAAGCGCTGCGCAATTGCAg                                                                                                                                                         >  1:359729/1‑227 (MQ=255)                                                                     tgCCATTGT‑CCCTCGTGGTCAAAGCGCACTTCATTTCTCAT‑CCGGATCGGTGG‑ATTTGCTGAGCCATCTGCCAGATGGTCAACCGGTGA‑CAGACATTCACAGTGGACTGGAATCTTTGCGTCAGCAGGCACCAGCCCATCCTTACTCCCTTCTCGAGCCACAACGCGTGCTCGATCTCGCTTGCCAGGCGCAGGCATTAATCGCTGAAAGCGGTGCGCAATTGCAggg                                                                                                                                                       <  1:396293/228‑1 (MQ=255)                                                                                aCTCGTGGTCAAAGCGCACTGCATTTCTCAT‑CCGGATCGCTGG‑ATTTGCTGAGCCATCTGCCAGATGGTCAACCGGTGA‑CAGACATTCACAGTGGACTGGAATCTTTGCGTCAGCAGGCACCAGCCCATCCTTACTCCCTTCTCGAGCCACAACGCGTGCTCGCTCTCGCTTGCCAGGCGCAGGCCTTAATCGCTGAAAGCGGTGCGCAATTGCAGGGCAGCGTAGaa                                                                                                                                             >  1:342389/1‑228 (MQ=255)                                                                                aCTCGTGGTCAAAGCGCACTGCATTTCTCAT‑CCGGATCGCTGG‑ATTTGCTGAGCCATCTGCCAGATGGTCAACCGGTGA‑CAGACATTCACAGTGGACTGGAATCTTTGCGTCAGCAGGCACCAGCCCATCCTTACTCCCTTCTCGAGCCACAACGCGTGCTCGATCTCGCTTGCCAGGCGCAGGCCTTAATCGCTGAAAGCGGTGCGCAATTGCAGGGAAGCGTaaaa                                                                                                                                             >  2:394995/1‑228 (MQ=255)                                                                                aCTCGTGGTCAAAGCGCACTGCATTTCTCAT‑CCGGATCGCTGG‑ATTTGCTGAGCCATCTGCCAGATGGTCAACCGGTGA‑CAGACATTCACAGTGGACTGGAATCTTTGCGTCAGCAGGCACCAGCCCATCCTTACTCCCTTCTCGAGCCACAACGCGTGCTCGATCTCCCTTGCCAGGCGCAGGCCTTAACCGCTGAACGCGGGGCGCAATTCCCGGGCAGCGCAGaa                                                                                                                                             >  2:15533/1‑228 (MQ=255)                                                                                 ctggtgtTCAAAGCGCAGTGCATTTCTCAA‑CCGGATCGCTGG‑ATTGGCTGATCCATCTGCCAGATGGTCACCCGGTGA‑CAGACATTCACAGTGGACTGGAATCTTTGCGTCAGCAGTCCCCCGCCCATCCTTACTCCCTTCTCGAGCCACAACGCGTGCTCGATCTCGCTTGCCAGGCGCAGGCATTAATCGCTGAAAGCGGTgcgc                                                                                                                                                                 <  1:657045/204‑1 (MQ=255)                                                                                                  cTGCATTTCTCAT‑CCGGATCGCTGG‑ATTTGCTGAGCCATCTGCCAGATGGTCAACCGGTGA‑CAGACATTCACAGTGGACTGGAATCTTTGCGTCAGCAGGCACCAGCCCATCCTTACTCCCTTCTCGAGCCACAACGCGTGCTCGATCTCGCTTGCCAGGCGCAGGCATTAATCGCTGAAAGCGGTGCGCAATTGCAGGGCACCGTAGAACTTGCTCACCATCGGGt                                                                                                                            >  1:62231/1‑227 (MQ=255)                                                                                                       ttgtcact‑CCGGTGCGCTGGAATTTGCAGAGCACTCGCCCAGTGGCTCACCCGGTAA‑CAGACCTTCCCTGGGCACTGCACTTTTCGCGGCAGCAGCCCCCACCCCATCCGTACTCCCTTCACGAGCCCCAACGCGTGCTCGATCTCGCTTGCCAGGCGCAGGCATTAATCGCTGAAAGCGGTGCGCAATTGCAGGGCAGCGTAGAACTTGCTCCCCAGCGGGTTAcg                                                                                                                        <  2:26961/220‑1 (MQ=255)                                                                                                             aT‑CCGGATCGCTGG‑ATTTGCTGAGCCATCTGCCAGATGGTCAACCGGTGA‑CAGACATTCACAGTGGCCTGGAATCTTTGCGTCAGCAGGCACCAGCCCATCCTTACTCCCTTCTCGAGCCACAACGCGTGCCCGATCTCGCTTGCCAGGCGCATGCCTTAATCGCTGAACGCGGTGCGCAATTGCAGGGCCGCACAGAATTTGCCCCCCACCGGGttgacccgcccg                                                                                                                 >  2:271534/1‑217 (MQ=255)                                                                                                             aT‑CCGGATCGCTGG‑ATTTGCTGAGCCATCTGCCAGATGGTCAACCGGTGA‑CAGACATTCACAGTGGACTGGAATCTTTGCGTCAGCAGGCACCAGCCCATCCTTACTCCCTTCTCGAGCCCCAACGCGTGCTCGATCTCGCTTGCCAGGCGCAGGCATTAATCGCTGAAAGCGGTGCGCACTTGCAGGGCAGCGTAGAACTTGCTCACCAGCGGGTTAggcacgccg                                                                                                                 >  1:165751/1‑221 (MQ=255)                                                                                                                         ggg‑ATTTGCTGGGCCATCTGCCAGATGGGCAACCGGTGA‑CAGACATTCACAGTGGACTGGAATCTTTGCGTCAGCAGGCACCAGCCCATCCTTACTCCCTTCTCGAGCCACAACGCGTGCTCGATCTCGCTTGCCAGGCGCAGGCATTAATCGCTGAAAGCGGTGCGCAATTGCAGGGCAGCGTAGAACTTGCTCACCAGCGGGTTACGCCGCTCGGCACTCTGCGct                                                                                                     <  1:77734/227‑1 (MQ=255)                                                                                                                                       ccATGTGCCAGGGGGTAACCCGTCGA‑CCGACCTTCCGCGTGGACGGGTATCTGTGCCTCAGTAGGCCCCAACCCCTCCTTACCCCCGGCTCGAGCCACAACGCGGGGTCGATCTCGCTTGCCAGGCGCAGGCATTAATCGCTGAAAGCGGTGCGCAATTGCAGGGCAGCGTAGAACTTGCTCACCAGCGGGTTACGCCGCTCGGCACTCTGCGCTCTACCTGGCTAAg                                                                                        <  2:608505/228‑1 (MQ=255)                                                                                                                                               cggttGGTCACCCGGGGA‑CAGCCCTTCCCGGGGGGTTGGAATCTTGGGGTCGTGCGGCCCCGGCCCTTCCTTCCCCCCTTCCCGAGCCACAACGCGTGCTCGATCTCGCTTGCCAGGCGCAGGCATTAATCGCTGAAAGCGGTGCGCAATTGCAGGGCAGCGTAGAACTTGCTCACCAGCGGGTTACGCCGCTCGGCACTCTGCGCTCTACCTGGCTAAGTTCGCCAg                                                                                <  2:42875/224‑1 (MQ=255)                                                                                                                                                                           ccaGTGGACTGGAATCTTTGCGGCAGCCGGCCCCACCCCATCCTTACGTCCTTCCCGAGCCTCAACGCGTGCTCGATCTCGCTTGCCAGGCGCAGGCATTAATCGCTGAAAGCGGTGCGCAATTGCAGGGCAGCGTAGAACTTGCTCACCAGCGGGTTACGCCGCTCGGCACTCTGCGCTCTACCTGGCTAAGTTCGCCAGAAGTCCCCGTCTGGCCGCTGCCcgcg                                                      <  1:587063/226‑1 (MQ=255)                                                                                                                                                                               tGGACTGGAATCTTTGCGTCAGCAGGCACCAGCCCATCCTTACTCCCTTCTCGAGCCACAACGCGTGCTCGATCTCGCTTGCCAGGCGCAGGCATTAATCGCTGAAAGCGGTGCGCAATTGCAGGGCAGCGTAGAACTTGCTCACCAGCGGGTTACGCCGCTCGGCACTCTGCGCTCTACCTGGCTAAGTTCGCCAGAAGTCCCCGCCTGGCCGCTGCCCCCgagcaa                                                 >  1:200362/1‑224 (MQ=255)                                                                                                                                                                                 gACTGGAATCTTTGCGTCAGCAGGCACCAGCCCATCCTTACTCCCTTCTCGAGCCACAACGCGTGCTCGATCTCGCTTGCCAGGCGCAGGCATTAATCGCTGAAAGCGGTGCGCAATTGCAGGGCAGCGTAGAACTTGCTCACCAGCGGGTTACGCCGCTCGGCACTCTGCGCTCTACCTGGCTAAGTTCGCCCGAAGTCCCCGTCTGGCCGCTGCCCGCGACGaaaa                                               >  1:409130/1‑228 (MQ=255)                                                                                                                                                                                 gACTGGAATCTTGGCGTCAGCAGGCACCAGCCCATCCTTACTCCCTTCCCGAGCCACAACGCGTGCTCGATCTCGCTTGCCAGGCGCAGGCATTAATCGCTGAAAGCGGTGCGCAATTGCAGGGCAGCGTAGAACTTGCTCACCAGCGGGTTACGCCGCTCGGCACTCTGCGCTCTACCTGGCTAAGTTCGCCAGAAGTCCCCGTCTGGCCGCTGCCCGCgaagaa                                                 <  1:394995/226‑1 (MQ=255)                                                                                                                                                                                  cttGGAATTTGTGCGGCAGGGGGCCCCCGCCCCTCCTTACCCCCCTCTCGAGCCACAACGCGTGGCCGACCTCGCTTGCCAGGCGCAGGCATTAAGCGCTGAAAGCGGTGCGCAATTGCAGGGCAGCGTAGAACTTGCTCACCAGCGGGTTACGCCGCTCGGCACTCTGCGCTCTACCTGGCTAAGTTCGCCAGAAGTCCCCGTCTGGCCGCTGCCCGCGAAGAAAat                                              <  2:342389/226‑1 (MQ=255)                                                                                                                                                                                  aCTGGAATCTTTGCGTCAGCAGGCCCCGGCCCATCCTTACTCCCATCCCGAGCCAAAACGCGTGCTCGATCTCGCTTGCCAGGCGCAGGCATTAATCGCTGAAAGCGGTGCGCAATTGCAGGGCAGCGTAGAACTTGCTCACCAGCGGGTTACGCCGCTCGGCACTCTGCGCTCTACCTGGCTAAGTTCGCCAGAAGTCCCCGTCTGGCCGCTGCCCGCGAAGAAAat                                              <  1:575961/228‑1 (MQ=255)                                                                                                                                                                                    tGGAATCTTTGCGTCAGCAGGCACCAGCCCATCCTTACTCCCTTCTCGAGCCACAACGCGTGCTCGATCTCGCTTGCCAGGCGCAGGCATTAATCGCTGAAAGCGGTGCGCAATTGCAGGGCAGCGTAGAACTTGCTCACCAGCGGGTTACGCCGCCCGGCACCCTCCGCTCTACCTGGCTCAGTCCGCCCGAGGCCCCCCTCCGGCCGCTCCCCCCCAACAAAaaat                                            >  2:40409/1‑228 (MQ=255)                                                                                                                                                                                       aaTCTTTGCGTCAGCAGGCACCAGCCCATCCTTACCCCCTTCCCGAGCCACAACGCGTGCTCGATCTCGCTTGCCAGGCGCAGGCATTAATCGCTGAAAGCGGTGCGCAATTGCAGGGCAGCGTAGAACTTGCTCACCAGCGGGTTACGCCGCTCGGCACTCTGCGCTCTACCTGGCTAAGTTCGCCAGAAGTCCCCGTCTGGCCGCTGCCCGCGAAGAAAATAtgtg                                         <  1:20485/228‑1 (MQ=255)                                                                                                                                                                                           tttGCGTCAGCAGGCACCAGCCCATCCTTACTCCCTTCTCGAGCCACAACGCGTGCTCGATCTCGCTTGCCAGGCGCAGGCATTAATCGCTGAAAGCGGTGCGCAATTGCAGGGCAGCGTAGAACTTGCTCACCAGCGGGTTACGCCGCTCGGCACTCTGCGCTCTACCTGGCTAAGTTCGCCAGAAGTCCCCGTCCGGCCGCTGCCCGCGACGAAAATTTGGGTAg                                      >  1:217094/1‑227 (MQ=255)                                                                                                                                                                                           tttGCGTCAGCAGGCACCAGCCCATCCTTACTCCCTTCTCGAGCCACAACGCGTGCTCGATCTCGCTTGCCAGGCGCAGGCATTAATCGCTGAAAGCGGTGCGCAATTGCAGGGCAGCGTAGAACTTGCTCACCAGCGGGTTACGCCGCTCGGCACTCTGCGCTCTACCTGGCTAAGTTCGCCAGAAGTCCCCGCCTGGCCGCTGCCCGCGAAGAAAATAtgggcagg                                     >  2:56056/1‑222 (MQ=255)                                                                                                                                                                                                   atcagGCGCCAGCCCATCCTTACTCCCTTCTCGAGCCACAACGCGTGCTCGATCTCGCTTGCCAGGCGCAGGCATTAATCGCTGAAAGCGGTGCGCAATTGCAGGGCAGCGTAGAACTTGCTCACCAGCGGGTTACGCCGCTCGGCACTCTGCGCTCTACCTGGCTAAGTTCGCCAGAAGTCCCCGTCTGGCCGCTGCCCGCGAAGAAAATATGTGTAGTGGGAATTa                             <  1:271534/226‑1 (MQ=255)                                                                                                                                                                                                         acagaGCCCATCCTTACACCCTTCCCGAGCCACAACGCGTGCTCGATCTCGCTTGCCAGGCGCAGGCATTAATCGCTGAAAGCGGTGCGCAATTGCAGGGCAGCGTAGAACTTGCTCACCAGCGGGTTACGCCGCTCGGCACTCTGCGCTCTACCTGGCTAAGTTCGCCAGAAGTCCCCGTCTGGCCGCTGCCCGCGAAGAAAATATGTGTAGTGGGAATTAGCGGc                        <  1:116139/223‑1 (MQ=255)                                                                                                                                                                                                                ctctctctccTCCTTTCTGGAGCCACAACGGGGGCGCGATCTCGCTTGCCAGGCGCAGGCATTAATCGCTGAAAGCGGTGCGCAATTGCAGGGCAGCGTAGAACTTGCTCACCAGCGGGTTACGCCGCTCGGCACTCTGCGCTCTACCTGGCTAAGTTCGCCAGAAGTCCCCGTCTGGCCGCTGCCCGCGAAGAAAATATGTGTAGTGGGAATTAGCGGCCTGAt                   <  2:62231/216‑1 (MQ=255)                                                                                                                                                                                                                    ccTTACTCCCTTCTCGAGCCACAACGCGTGCTCGATCTCGCTTGCCAGGCGCAGGCATTAATCGCTGAAAGCGGTGCGCAATTGCAGGGCAGCGTAGAACTTGCTCACCAGCGGGTTACGCCGCTCGGCACTCTGCGCTCTACCTGGCTAAGTTCGCCAGAAGTCCCCGTCTGGCCGCTGCCCGCGAAGCAAATATGTGTAGTGGGAATTCGCGGCCTGATGGAtttt            >  1:121356/1‑228 (MQ=255)                                                                                                                                                                                                                        aCTCCCTTCTCGAGCCACAACGCGTGCTCGATCTCGCTTGCCAGGCGCAGGCATTAATCGCTGAAAGCGGTGCGCAATTGCAGGGCAGCGTAGAACTTGCTCACCAGCGGGGTACGCCGCTCGGCACTCTGCGCTCTACCTGGCTAAGTTCCCCCGAAGTCCCCGCCTGGCCGCTGCCCGCGCAGAAAATATGTGTAGTGGGAATTAGCCGCCGGATGGGtttctgcc        >  2:116139/1‑223 (MQ=255)                                                                                                                                                                                                                          ccccTTCCTGAGCCCCAAGGCTTGCGCTATCTCGCTTGCCAGGCGCAGGCATTAAGCGCTGAAAGGGGTGCCCAATTGCAGGGCAGCGTAGAACTTGCTCACCAGCGGGTTACGCCGCTCGGCACTCTGCGCTCTTCCTGGCTAAGTTCGCCAGAAGTCCCCGTCTGGCCGCTGCCCGCGAAGAAAATATGTGTAGTGGGAATTAGCGGCCTGATGGATTTTCAGgcg      <  1:683917/227‑1 (MQ=255)                                                                                                                                                                                                                           cccTTCTCGAGCCACAACGCGTGCTCGATCTCGCTTGCCAGGCGCAGGCATTAATCGCTGAAAGCGGTGCGCAATTGCAGGGCAGCGTAGAACTTGCTCACCAGCGGGTTACGCCGCTCGGCACTCTGCGCTCTACCTGGCTAAGTTCGCCAGAAGTCCCCGTCTGGCCGCTGCCCGCGAAGAAAATATGTGTAGTGGGAATTAGCGGCCTGATGGATTTTCAGgcg      <  1:40409/227‑1 (MQ=255)                                                                                                                                                                                                                              ttCTCGAGCCACAACGCGTGCTCGATCTCGCTTGCCAGGCGCAGGCATTAATCGCTGAAAGCGGTGCGCAATTGCAGGGCAGCGTAGAACTTGCTCACCAGCGGGTTACGCCGCTCGGCACTCTGCGCTCTACCTGGCTAAGTTCGCCAGAAGTCCCCGTCTGGCCGCTGCCCGCGAAGAAAAAATGTGTAGTGGGAATTAGCGGCCTGATGGCATTTCCGGCGcccc  >  1:28425/1‑228 (MQ=255)                                                                                                                                                                                                                               tctcGAGCCACAACGCGTGCTCGATCTCGCTTGCCAGGCGCAGGCATTAATCGCTGAAAGCGGTGCGCAATTGCAGGGCAGCGTAGAACTTGCTCACCAGCGGGTTACGCCGCTCGGCACTCTGCGCTCTACCTGGCTAAGTTCGCCAGCAGTCCCCGTCTGGCCGCTGCCCGCgaaga                                                  >  1:561678/1‑179 (MQ=255)                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                 TCATTATGGGCGGCGGCCTCGCCGGATTACTCTGTGGCCTGCAACTGCAAAAACACGGCCTGC‑GCTGTGCCATTGT‑CACTCGTGGTCAAAGCGCACTGCATTTCTCAT‑CCGGATCGCTGG‑ATTTGCTGAGCCATCTGCCAGATGGTCAACCGGTGA‑CAGACATTCACAGTGGACTGGAATCTTTGCGTCAGCAGGCACCAGCCCATCCTTACTCCCTTCTCGAGCCACAACGCGTGCTCGATCTCGCTTGCCAGGCGCAGGCATTAATCGCTGAAAGCGGTGCGCAATTGCAGGGCAGCGTAGAACTTGCTCACCAGCGGGTTACGCCGCTCGGCACTCTGCGCTCTACCTGGCTAAGTTCGCCAGAAGTCCCCGTCTGGCCGCTGCCCGCGAAGAAAATATGTGTAGTGGGAATTAGCGGCCTGATGGATTTTCAGGCGCACC  >  NC_000913/2354281‑2354724

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