Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,045,276	T→C	7.9%	D113D (GAT→GAC)	yeeJ →	putative adhesin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,045,276	0	T	C	7.9%	141.7 / 2.0	54	D113D (GAT→GAC)	yeeJ	putative adhesin
Reads supporting (aligned to +/- strand):  ref base T (28/18);  new base C (2/2);  total (30/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.06e-01

GTTCCCTATGGCTGCGGCAGCACAAGGTGTGGTAAACGCCGCAACCCAACAACCAGTTCCTGCACAAATTGCCATTGCAAATGCCAATACGGTGCCCTACACCCTTGGAGCGTTGGAATCGGCCCAAAGCGTTGCCGAACGTTTCGGTATTTCGGTGGCTGAGTTACGCAAACTCAACCAGTTTCGTACGTTTGCTCGAGGTTTTGATAATGTCCGCCAGGGTGATGAACTGGATGTCCCGGCACAAGTTAGTGAAAAAAAATTAACCCCGCCGCCGGGTAATAGCAGTGACAACCTCGAGCAACAGATAGCCAGTACTTCACAGCAAATCGGGTCTCTGCTCGCCGAAGATATGAACAGCGAGCAAGCGGCAAATATGGCGC‑GTGGATGGGCCTCTTCTCAGGCTTCAGGCGCAATGACAGACTGGTTAAGCCGCTTCGGTACCGC  >  NC_000913/2045051‑2045497                                                                                                                                                                                                                                  |                                                                                                                                                                                                                               gTTCCCTATGGCTGCGGCAGCACAAGGTGTGGTAAACGCCGCAACCCAACAACCAGTTCCTGCACAAATTGCCATTGCAAATGCCAATACGGTGCCCTACACCCTTGGAGCGTTGGAATCGGCCCAAAGCGTTGCCGAACGTTTCGGTATTTCGGTGGCTGAGTTACGCAAACTCAACCAGTTTCGTACGTTTGCTCGAGGTTTTGATAATGTCCGCCAGGGtgatga                                                                                                                                                                                                                              <  2:522251/228‑1 (MQ=255)    cccTATGGCTGCGGCAGCACAAGGTGTGGTAAACGCCGCAACCCAACAACCAGTTCCTGCACAAATTGCCATTGCAAATGCCAATACGG‑GCCCTACACCCTTGGAGCGTTGGAATCGGCCCAAAGCGTTGCCGAACGTTTCTGTATTTCGGTGGCTGCGTTACGCAAACTCAACCAGTTTCGTACCTTTGCTCGAGGTTTTGATAATGTCGGCCAGGGTGATAAACTg                                                                                                                                                                                                                          >  1:398191/1‑228 (MQ=255)       tATGGCTGCGGCAGCACAAGGTGTGGTAAACGCCGCAACCCAACAACCAGTTCCTGCACAAATGGCCATTGCAAATGCCAATACGGTGCCCTACACCCTTGGAGCGTTGGAATCGGCCCAAAGCGTTGCCGAACGTTTCGGTATTTCGGTGGCTGAGTTACGCAAACTCAACCAGTTTCGTACGTTTGCTCGAGGTTTTGATAATGTCCGCCAGGGTGATGAACTGGa                                                                                                                                                                                                                        <  1:405368/228‑1 (MQ=255)               gggCAGCACAAGGTGTGGTAAACGCCGCAACCCAACAACCAGTTCCTGCACAAATTGCCATTGCAAATGCCAATACGGTGCCCTACACCCTTGGAGCGTTGGAATCGGCCCAAAGCGTTGCCGAACGTTTCGGTATTTCGGTGGCTGAGTTACGCAAACTCAACCAGTTTCGTACGTTTGCTCGAGGTTTTGATAATGTCCGCCAGGGTGATGAACTGGATGTCCCgg                                                                                                                                                                                                                <  1:464342/227‑1 (MQ=255)               cGGCAGCACAAGGTGTGGTAAACGCCGCAACCCAACAACCAGTTCCTGCACAAATTGCCATTGCAAATGCCAATACGGTGCCCTACACCCTTGGAGCGTTGGAATCGGCCCAAAGCGTTGCCGAACGTTTCGGTATTTCGGTGGCTGAGTTACGCAAACTCAACCAGTTTCGTACGTTTGCTCGAGGTTTTGTTAATGTCCGCCCGGGTGATGAACTGGATGTCCCgg                                                                                                                                                                                                                >  1:89356/1‑228 (MQ=255)                                 tAAACGCCGCAACCCAACAACCAGTTCCTGCCCAAATTGCCATTGCAAATGCCAATCCGGTGCCCTACCCCCTTGGAGCGTTGGTATCGGCCCAAAGCGTTGCCGAACGTTTCGGTATTTCGGTGGCTGAGTTACGCAAACTCAACCAGTTTCGTACGTTTGCTCGAGGTTTTGATAATGTCCGCCAGGGTGATGAACTGGATGTCCCGGCACAAGTTAGTGaaaa                                                                                                                                                                                                <  1:93574/226‑1 (MQ=255)                                 tAAACCCCGCAACCCCACAACCAGTTCCCGCCAAAATTGCCATTGCAAATGCCAATTCGGTGCCCTACACCCTTGGAGCGTTGGAATCGGCCAAAAGCGTTGCCGAACGTTTCGGTATTTCGGTGGCTGAGTTACGCAAACTCAACCAGTTTCGTACGTTTGCTCGAGGTTTTGATAATGTCCGCCAGGGTGATGAACTGGATGTCCCGGCACAAGTTAGTGaaa                                                                                                                                                                                                 <  2:360125/225‑1 (MQ=255)                                    aCGCCGCAACCCAACAACCAGTTCCTGCACAAATTGCCATTGCAAATGCCAATACGGTGCCCTACACCCTTGGAGCGTTGGAATCGGCCCAAAGCGTTGCCGAACGTTTCGGTATTTCGGTGGCTGAGTTACGCAAACTCAACCAGTTTCGTCCGTTTGCTCGAGGTTTTGATAATGTCCGCCAGGGTGATGAACTGGATGTCCCGGCCCAAGTTAGTGaaaaaaaa                                                                                                                                                                                            >  1:243754/1‑227 (MQ=255)                                      gccgcAACCCAACAACCAGTTCCTGCACAAATTGCCATTGCAAATGCCAATACGGTGCCCTACACCCTTGGAGCGTTGGAATCGGCCCAAAGCGTTGCCGCACGTTTCGGTATTTCGGTGGCTGAGTTACGCAAACTCAACCAGTTTCGTACGTTTGCTCGAGGTTTTGATAATGTCCGCCAGGGTGATGAACTGGATTTCCCGGCACAAGTTAGTGAAAAAAAATTa                                                                                                                                                                                         >  1:253306/1‑228 (MQ=255)                                      gccgcAACCCAACAACCAGTTCCTGCACAAATTGCCATTGCAAATGCCAATACGGTGCCCTACACCCTTGGAGCGTTGGAATCGGCCCAAAGCGTTGCCGAACGTTTCGGTATTTCGGTGTCTGAGTTACGCAAACTCAACCAGTTTCGTACGTTTGCCCGAGGTTTTGATAATTACCGCCAGGGGGATGAACTGGATGTCCCGGCCCAAGATAGGTAAAAAAAATTa                                                                                                                                                                                         >  2:20609/1‑228 (MQ=255)                                                           ccTGCACAAATGGCCATTGCAAAGGCCAATACGGGGCCCTTCTCCCTGGGAGGTTTGGAATCGGCCCAAAGCGTTGCCGAACGTTTCGGTTTTTCGGTGGCTGAGTTACGCAAACTCAACCAGTTTCGTACGTTTGCTCGAGGTTTTGATAATGTCCGCCAGGGTGATGAACTGGATGTCCCGGCACAAGTTAGTGAAAAAAAATTAACCCCGCCGCCGGGTAATAGc                                                                                                                                                                    <  2:89356/228‑1 (MQ=255)                                                                            tGCAAATGCCAATACGGTGCCCTACACCCTTGGAGCGTTGGAATCGGCCCAAAGCGTTGCCGAACGTTTCGGTATTTCGGTGGCTGAGTTACGCAAACTCAACCAGTTTCGTACGTTTGCTCGGGGTTTTGATAATGTCCGCCAGGGTGATGAACTGGATGTCCCGGCCACAGTTTGTGAAAAAAAATTAACCCCGCCGCCGGGGAAAAGCAGGGACAAACCCGaaaa                                                                                                                                                   >  2:632775/1‑225 (MQ=255)                                                                                             tGCCCTACACCCTTGGAGCGTTGGAATCGGCCCAACGCGTTGCCGAACGTTTCGGTATTTCGGTGGCTGAGTTACGCAAACTCAACCAGTTTCGTACGTTTGCTCGAGGTTTTGATAATGTCCGCCAGGGCGATGAACTGGATGGCCCGGCACAAGTTAGTGAAAAAAAATTAACCCCCCCCCCGGGTAAAAGCATGGACAACCCCGAGCAACAAAAAGCCAGGACtt                                                                                                                                  >  2:337010/1‑228 (MQ=255)                                                                                             tGCCCTACACCCTTGGAGCGTTGGAATCGGCCCAAAGCGTTGCCGAACGTTTCGGTATTTCGGTGGCTGAGTTACGCAAACTCAACCAGTTTCGTACGTTTGCTCGAGGTTTTGATAATGTCCGCCAGGGTGATGAACTGGATGTCCCGGCACAAGTTAGTGAAAAAAAATTAACCCCGCCGCCGGGTACTAGCAGTGACAGACTCGAGCAAAAGATAACCAGAACtt                                                                                                                                  >  2:195795/1‑228 (MQ=255)                                                                                             tGCCCTACACCCTTGGAGCGTTGGAATCGGCCCAAAGCGTTGCCGAACGTTTCGGTATTTCGGTGGCTGAGTTACGCAAACTCAACCAGTTTCGTACGTTTGCTCGAGGCTTTGATAATGTCCGCCAGGGTGATGAACTGGATGTCCCGGCACAAGTTAGTGAAAAACAATTAACCCCTCCGCCGGGTAATAGCAGTGACAACCCCGGGCAAACAATACCCCGTACtc                                                                                                                                  >  1:593364/1‑227 (MQ=255)                                                                                              gCCCTACACCCTTGGAGCGTTGGAATCGGCCCAAAGCGTTGCCGAACGTTTCGGTATTTCGGTGGCTGAGTTACGCAAACTCAACCAGTTTCGTACGTTTGCTCGAGGTTTTGATAATGTCCGCCAGGGTGATGAACTGGATGTCCCGGCACAAGTTAGTGAAAAAAAATTAACCCCGCCGCCGGGAAATAGCAGTGACAACCTCGAGCAACAGATAGCCAGTACTTc                                                                                                                                 >  1:584742/1‑228 (MQ=255)                                                                                                         ttGGAGCGTTGGAATCGGCCCAAAGCGTTGCCGAACGTTTCGGTATTTCGGTGGCTGAGTTACGCAAACTCACCCAGTTTCGTACGTTTGCTCGAGGTTTTGATAATGTCCGCCAGGGTGATGAACTGGATGTCCCGGCACAAGTTAGTGAAAAAAAATTAACCCCGCCGCCGGGTAATAGCAGTGACAACCTCGAGCAACAGATAGCCAGTACTTCACAGCAAATc                                                                                                                       >  1:126533/1‑227 (MQ=255)                                                                                                                 ttGGAATCGGCCCAAAGCGTTGCCGAACGTTTCGGTATTTCGGTGGCTGCGTTACGCAAACTCAACCAGTTTCGTACGTTTGCTCGAGGTTTTGATAATGTCCGCCAGGGTGATGAACTGGATGTCCCGGCACAAGTTAGTGACAAAAAATTAACCCCGCCGCCGGGGAATAGCAGGGACAACCTCGAGCAACAGATAGCCAGTACTTCACATCAAATAGGGTCTCTg                                                                                                              >  1:254952/1‑228 (MQ=255)                                                                                                                 ttGGAATCGGCCCAAAGCGTTGCCGAACGTTTCGGTATTTCGGTGGCTGAGTTACGCAAACTCAACCAGTTTCGTACGTTTGCTCGAGGTTTTGATAATGTCCGCCAGGGTGATGAACTGGATGTCCCGGCACAAGTTAGTGAAAAAAAATTAACCCCGCCGCCGGGTAATAGCAGTGACAACCTCAAGCAACAGATAGCCAGCACTTCACAGCCAATCGGGtctccg                                                                                                              >  1:275142/1‑226 (MQ=255)                                                                                                                 ttGGAATCGGCCCAAAGCGTTGCCGAACGTTTCGGTATTTCGGTGGCTGAGTTACGCAAACTCAACCAGTTTCGTACGTTTGCTCGAGGTTTTGATAATGTCCGCCAGGGTGATGAACTGGATGTCCCGGCACAAGTTAGTGAAAAAAAATTAACCCCGCCGCCGGGAAATAGCAGGGACAACCTCGAGCGACAGATACCCAGTACTTCACAGGAAATAGGGTCCCTg                                                                                                              >  1:115051/1‑228 (MQ=255)                                                                                                                   ggAATCGGCCCAAAGCGTTGCCGAACGTTTCGGTATTTCGGTGGCTGAGTTACGCAAACTCAACCAGTTTCGTACGTTTGCTCGAGGTTTTGATAATGTCCGCCAGGGTGATGAACTGGATGTCCCGGCACAAGTTAGTGAAAAAAAATTAACCccgcc                                                                                                                                                                                 >  2:123197/1‑159 (MQ=255)                                                                                                                    gAATCGGCCCAAAGCGTTGCCGAACGTTTCGGTATTTCGGTGGCTGAGTTACGCAAACTCAACCAGTTTCGTACGTTTGCTCGAGGTTTTGATAATGTCCGCCAGGGTGATGAACTGGATGTCCCGGCACAAGTTAGTGAAAAAAAATTAACCccgccg                                                                                                                                                                                <  1:123197/159‑1 (MQ=255)                                                                                                                     aaTCGGCCCAAAGCGTTGCCGAACGTTTCGGTATTGCGGTGGCTGAGTTACGCAAACTCAACAAGTTTCGTACGTTTGCTCGAGGTTTTGATAATGTCCGCCAGGGTGATGAACTGGATGTCCCGGCACAAGTTAGTGAAAAAAAATTAACCCCGCCGCCGGGTAATAGCAGTGACAACCTCGAGCAACAGATAGCCAGTACTTCACAGCAAATCGGGTCTCTGCTCg                                                                                                          <  1:369248/228‑1 (MQ=255)                                                                                                                         ggCCCAAAGCGTTACCGAACGTTTCGGTATTTCGGTGGCTGAGTTACGCAAACTCAACCAGTTTCGTACGTTTGCTCGAGGTTTTGATAATGTCCGCCAGGGTGATGAACTGGATGTCCCGGCACAAGTTAGTGAAAAAAAATTAACCCCGCCGCCGGGTAATAGCAGTGACAACCTCGAGCAACAGATAGCCAGTACTTCACAGCAAATCGGGTCTCTGCTCGCCg                                                                                                       <  1:195795/227‑1 (MQ=255)                                                                                                                         ggCACAAAGCTTTGCGGAACGTTTCGGTATTTCGGTGGCTGAGTTACGCAAACTCAACCAGTTGCGTACGTTTGCTCAAGGTTTTGATAATGCCCGCCAGGGTGATGAACTGGATGGCCCGGCACAAGTTAGTGAAAAAAAATTAACCCCGCCGCCGGGTAATAGCAGTGACAACCTCGAGCAACAGATAGCCAGTACTTCACAGCAAATCGGGTCTCTGCTCGCCg                                                                                                       <  1:316474/227‑1 (MQ=255)                                                                                                                           cccAAAGCGTTGCCGAACGTTTCGGTTTTTCGGGTGCTTAGTTACGCAAACTCAACCGGTTTCGTACGTTTGCTCGAGGTTTTGATAATGTCCGCCAGGGTGATGAACTGGATGTCCCGGCACAAGTTAGTGAAAAAAAATTAACCCCGCCGCCGGGTAATAGCAGTGACAACCTCGAGCAACAGATAGCCAGTACTTCACAGCAAATCGGGTCTCTGCTCGCCGaa                                                                                                     <  2:593364/227‑1 (MQ=255)                                                                                                                                aGCGTTGCGCACCTTTTCGGTATTTCGGTGGTTGATTTACGCAAACTCAACCATTTTCGTACGTTTGCTCGAGGTTTTGATAATGTCCGCCAGGGTGATGAACTGGATGTCCCGGCACAAGTTAGTGAAAAAAAATTAACCCCGCCGCCGGGTAATAGCAGTGACAACCTCGAGCAACAGATAGCCAGTACTTCACAGCAAATCGGGTCTCTGCTCGCCGAAGATATg                                                                                               <  1:541175/228‑1 (MQ=255)                                                                                                                                     tGCCGATCGTTTCGGTATTTCGGTGGCTGAGTTACGCAAACTCAACCAGTTTCGTACGTTTGCTCGAGGTTTTGATAATGTCCGCCAGGGTGATGAACTGCATGTCCCGGCACAAGTTAGTCAAAAAAAATTAACCCCACCCCCGGGTAATAGCACTGACAACCTCGAGAACCAAAAAGCCAGTCCTTCACAGCAAATAGGGTTTATGCCGGCCGAAGAAATGaaaaa                                                                                          >  2:444699/1‑225 (MQ=255)                                                                                                                                     tGCCGAACGTTTCGGTATTTCGGTGGCTGAGTTACGCAAACTCAACCAGTTTCGTACGTTTGCTCGAGGTTTTGATAATGTCCGCCAGGGTGATGAACTGGATGTCCCGGCACAAGTTAGTGAAAAAAAATTAACCCCGCCGCCGGGTAATACCAGTGACAACCTCGAGCAACAAATAGCCAGAACTTCCCAGCAAATCGGGCCTCTCCTCCCCCGAGAAATGAACac                                                                                          >  2:381200/1‑227 (MQ=255)                                                                                                                                     tGCCGAACGTTTCGGTATTTCGGTGGCTGAGTTACGCAAACTCAACCAGTTTCGTACGTTTGCTCGAGGTTTTGATAATGTCCGCCAGGGTGATGAACTGGATGTCCCGGCACAAGTTAGTGAAAAAAAATTAACCCCGCCGCCGGGAAATAGCAGTGACAACCTCGAGCAACAGGTAGCCCGTACTTCACAGCAAATCGGGTCTCGCCCCGCCGAAGATGTGaaacg                                                                                          >  1:601750/1‑225 (MQ=255)                                                                                                                                            cGTTTCGGTATTTCGGTGGCTGAGTTACGCAAACTCAACCAGTTTCGTACGTTTGCTCGAGGTTTTGATAATGTCCGCCAGGGTGATGAACTGGATGTCCCGGCCCAAGTTAGTGAAAAAAAATTAACCCCGCCGCCGGGTAATAGCAGTGACAACCTCGAGCAACAGATAGCCAGTACTTCACACAAAATCGGGACCCTGCTCGCCGCAGATATGCACAACGac                                                                                      >  2:354636/1‑224 (MQ=255)                                                                                                                                               ttCGGTATTTCGGTGGCTGAGTTACGCAAACTCAACCAGTTTCGTACGTTTGCTCGAGGTTTTGATAATGTCCGCCAGGGTGACGAACTGGATGTCCCGGCACAAGTTAGTGAAAAAAAATTAACCCCGCCGCCGGGTAATAGCAGTGACAACCTCGAGCAACAGATAGCCAGTACTTCACAGCAAATCGGGTCTCTGCTCGCCGAAGATATGAACAGCGAGCAAGc                                                                                 >  1:572689/1‑227 (MQ=255)                                                                                                                                                tCGGTATTTCGGTGGCTGAGTTACGCAAACTCAACCAGTTTCGTACGTTTGCTCGAGGTTTTGATAATGTCCGCCAGGGTGACGAACTGGATGTCCCGGCACAAGTTAGTGAAAAAAAATTAACCCCGCCGCCGGGTAATAGCAGTGACAACCTCGAGCAACAGATAGCCAGTACTTCACAGCAAATCGGGTCTCTGCTCGCCGAAGATATGAACAGCGAGCAAGCgg                                                                               <  2:572689/228‑1 (MQ=255)                                                                                                                                                         cGGTGGCTGAGTTACGCAAACTCAACCAGTTTCGTACGTTTGCTCGAGGTTTTGATAATGTCCGCCAGGGTGATGAACTGGATGTCCCGGCACAAGTTAGTGAAAAAAAATTAACCCCGCCGCCGGGTAATAGCAGTGACAACCTCGACCAACAGAAAGCCAGTACTTCACAGCAAATCGGGTCTCTGCTCGCCGAAGATATGAACAGAGAGCAAGCGGCAAATATgg                                                                      >  1:52078/1‑228 (MQ=255)                                                                                                                                                            tggcctaGTTACGCAAACTCAACCAGTTTCGTACGTTTGCTCGAGGTTTTGATAATGTCCGCCAGGGTGATGAACTGGATGTCCCGGCACAAGTTAGTGAAAAAAAATTACCCCCGCCGCCGGGTAATAGCAGTGCCAACCTCGAGCAACAGATAGCCAGTACTTCACAGCAAATCGGGTCTCTGCTCGCCGAAGATATGAACAGCGAGCAAGCGGCAAATATGgcgc                                                                   <  1:49485/222‑1 (MQ=255)                                                                                                                                                                     tACGCAAACTCAACCCTTTTGTGTGGTTTGCTCGAGGTTTTGATAATGGCCGCGAGGTGGAGGACCTGGATTCCCCGGCACAAGTTATTGAAAAAAAATTAACCCCGCCGCCGGGTAATAGCAGTGACAACCTCGAGCAACAGATAGCCAGTACTTCACAGCAAATCGGGTCTCTGCTCGCCGAAGATATGAACAGCGAGCAAGCGGCAAATATGGCGC‑GTGGATggg                                                         <  2:243754/228‑1 (MQ=255)                                                                                                                                                                            aCTCAACCAGTTTCGTACGTTTGCTCGAGGTTTTGATAATGTCCGCCAGGGTGATGAACTGGATGTCCCGGCACAAGTTAGTGAAAAAAAATTAACCCCGCCGCCGGGTAATAGCAGTGACAACCTCGAGCAACAGATAGCCAGTACTTCACAGCAAATCGGGTCTCTGCTCGCCGAAGATATGAACACCGAGCAACCGGCAAAAAAGGCGCAGGGGAGGGCCctctc                                                   >  2:302766/1‑227 (MQ=255)                                                                                                                                                                            aCTCAACCAGTTTCGTACGTTTGCTCGAGGTTTTGATAATGTCCGCCAGGGTGACGAACTGGATGTCCCGGCACAAGTTAGTGAAAAAAAATTAACCCCGCCGCCGGGTAATAGCAGTGACAACCTCGAGCAACAGATAGCCAGTACTTCACAGCAAATCGGGTCTCTGCTCGCCGAAGATATGAACAGCGAGCAAGCGGCAAATATGGCGG‑GGGGATGGGCCtcttc                                                  >  1:462906/1‑228 (MQ=255)                                                                                                                                                                             cTCAACCAGTTTCGTACGTTTGCTCGAGGTTTTGATAATGTCCGCCAGGGTGATGAACTGGATGTCCCGGCACAAGTTAGTGAAAAAAAATTAACCCCGCCGCCGGGTAATAGCAGTGACAACCTCGAGCAACAGATAGCCAGTACTTCACAGCAAATCGGGTCTCTGCTCGCCGACGATATGACCAGCGAGCAAGCGGCAAAAATGGCGC‑GTGGATGGGCCtcttct                                                 >  2:259507/1‑228 (MQ=255)                                                                                                                                                                              accacccttttCGTTCGTGTGCTGGGGTTTTTGATAAAGTCCCCCAGGTGGATTAACGGGATGTCCCGGCACAAGTTAGTCAAAAAAAATTAACCCCGCCGCCTGGTAATAGCAGTGACAACCTCGAGCGACAGATAGCCAGTACTTCACAGCAAATCGGGTCTCTGCTCGCCGAAGATATGAACAGCGAGCAAGCGGCAAATATGGCGC‑GTGGATGGGCCtcttct                                                 <  2:115051/219‑1 (MQ=255)                                                                                                                                                                                      tttCGTCAGTTTGCTCGAGGTTTTGATAATGTCCGCCAGGGTGATGAACTGGATGTCCCGGCACAAGTTAGTGAAAAAAAATTAACCCCGCCGCCGGGTAATAGCAGTGACAACCTCGAGCAACAGATAGCCAGTACTTCACAGCAAATCGGGTCTCTGCTCGCCGAAGATATGAACAGCGAGCAAGCGGCAAATATGGCGC‑GTGGATGGGCCTCTTCTCAGGCTTCa                                        <  1:337010/228‑1 (MQ=255)                                                                                                                                                                                      tttCGTACGTTTGCTCGAGGTTTTGATAATGTCCGCCAGGGTGATGAACTGGATGTCCCGGCACAAGTTAGTGAAAAAAAATTAACCCCGCCGCCGGGTAATAGCAGTGACAACCTCGAGCAACAGATAGCCAGTACTTCACAGCAAATCGGGTCTCTGCTCGCCGAAGATATGAACATCGAACAAGCGGCAAATATGGCGC‑TTGGATAGCCCCCTTCTCAGGCTTCa                                        >  1:452956/1‑228 (MQ=255)                                                                                                                                                                                      tttCGTACGTTTGCTCGAGGTTTTGATAATGTCCGCCAGGGTGATGAACTGGATGTCCCGGCACAAGTTAGTGAAAAAAAATTAACCCCGCCGCCGGGTAATAGCAGTGACAACCTCGAGCAACAGATAGCCAGTACTTCACAGCAAATCGGGTCTCTGCTCGCCGAAGATATGAACAGCGAGCAAGCGGCAAATATGGCGC‑GTGGATGGGCCTCTTCTCAGGCTTCa                                        >  1:18182/1‑228 (MQ=255)                                                                                                                                                                                      gttgtgcggTTTGCCGGAGGTTTGGAAAATGGCCCCAAGGGGGAGGAACGGGATGCCCCGGCACAGGCTAGTGAAAAAAAATTAACCCCGCCGCCGGGTAATAGCAGGGACAACCTCGAGCAACAGATAGCCAGTACTTCACAGCAAATCGGGTCTCTGCTCGCCGAAGATATGAACAGCGAGCAAGCGGCAAATATGGCGC‑GTGGATGGGCCTCTTCTCAGGCTTCa                                        <  2:163788/220‑1 (MQ=255)                                                                                                                                                                                        acGTACGTTTGCTCGAGGTTTTGATAATGTCCGCCAGGGTGATGAACTGGATGTCCCGGCACAAGTTAGTGAAAAAAAATTAACCCCGCCGCCGGGTAATAGCAGTGACAACCTCGAGCAACAGATAGCCAGTACTTCACAGCAAATCGTGTCTCTGCTCGCCGAAGATATGAACAGCGAGCAAGCGGCAAATATGGCGC‑GTGGATGGGCCTCTTCTCAGGCGTCAgg                                      <  1:444699/227‑1 (MQ=255)                                                                                                                                                                                          tcttgTTTGTGCGAGTTTTTGAAAATTCCCGCCAGTGTGAATAACTGGATTGCCCGGCACAAATTAGTGAAAAAAAATTAACCCCGCCGCCGGGTAATAGCAGTGACAACCTCGAGCAACAGATAGCCAGTACTTCACAGCAAATCGGGTCTCTGCTCGCCGAAGATATGAACAGCGAGCAAGCGTCAAATATGGCGC‑GTGGATGGGCCGCTTCTCAGGCTTCAGgcg                                    <  2:444799/224‑1 (MQ=255)                                                                                                                                                                                                ttGCTCGAGGTTTTGATAATGTCCGCCAGGGTGATGAACTGGATGTCCCGGCACAAGTTAGTGAAAAAAAATTAACCCCGCCGCCGGGTAATAGCAGTGACAACCTCGAGCAACAGATAGCCAGTACTTCCCAGCAAATCGGGTCTCTGCTCGcc                                                                                                        >  1:583715/1‑155 (MQ=255)                                                                                                                                                                                                 tGCTCGAGGTTTTGATAATGTCCGCCAGGGTGATGAACTGGATGTCCCGGCACAAGTTAGTGAAAAAAAATTAACCCCGCCGCCGGGTAATAGCAGTGACAACCTCGAGCAACAGATAGCCAGTACTTCACAGCAAATCGGGTCTCTGCTCGCCg                                                                                                       <  2:583715/155‑1 (MQ=255)                                                                                                                                                                                                   ctgaaGGTTTTGATAATGCCAGCCAGGGTGACGAAATGGCTGTCCCGGCACAGGTTAGTGAAAAAAAATTAGCCCCGCCGCCGGGTAATAGCAGTGACAACCTCGAGCAACAGATAGCCAGTACTTCCCAGCAAATCGGGTCTCTGCTCGCCGAAGATATGAACAGCGAGCAAGCGGCAAATATGGCGC‑GTGGATGGGCCTCTTCTCAGGCTTCAGGCGCAATGACAg                           <  2:462906/224‑1 (MQ=255)                                                                                                                                                                                                                    gTCCGCCAGTGTGGTGAACTGGATGTCCCGGCACAAGTTAGGGAAAAAAAATTAACCCCGCCGCCGGGTAATAGCAGTGACAACCTCGAGCAACAGATAGCCAGTACTTCACAGCAAATCGGGTCTCTGCTCGCCGAAGATATGAACAGCGAGCAAGCGGCAAATATGGCGC‑GTGGATGGGCCTCTTCTCAGGCTTCAGGCGCAATGACAGACTGGTTAAGCCGCtt           <  2:254952/227‑1 (MQ=255)                                                                                                                                                                                                                    gTCCGCCAGGGTGATGAACTGGATGTCCCGGCACAAGTTAGTGAAAAAAAATTAACCCCGCCGCCGGGTAATAGCAGTGACAACCTCGAGCAACAGATAGCCAGTACTTCACAGCAAATCGGGTCTCTGCTCGCCGAAGATATGAACAGCGGGCAAGCGCCAAATATGGCGC‑GTGGAGTGGCCTCTTCTCAGGCTTACGGGGCAATGAAAGACCGGTTAAACAGCtcg          >  2:576790/1‑226 (MQ=255)                                                                                                                                                                                                                       gcccAGGGTGAGGAACTGGATGCCCCGGCACAAGTAAGTAAAAAAAAATGACCCCCCCCGCCGGGTAAGAGCAGTGACAACCTCGAGCAACAGATGGCCAGTACTTCACAGCAAATCGGGTCTCTGCTCGTCGAAGATATGAACAGCGAGCAAGCGGCAAATATGGCGC‑GTGGATGGGCCTCTTCTCAGGCTTCAGGCGCAATGACAGACTGGTTAAGCCGCTTCGGt       <  2:160569/226‑1 (MQ=255)                                                                                                                                                                                                                         ccAGGGTGATGAACTGGATGTCCCGGCACAAGTTAGTGAAAAAAAATTAACCCCGCCGCCGGGTAATAGCAGTGACAACCTCGAGCAACAGATAGCCAGTACTTCACAGCAAATCGGGTCTCTGCTCGCCGAAGATATGAACAGCGAGCAAGCGGCAAATATGGCGC‑GTGGATGGGCCTCTTCTCAGGCTTCAGGCGCAATGACAGACGGGTTAAGCCGCTTCgggca     >  2:330803/1‑225 (MQ=255)                                                                                                                                                                                                                         ccAGGGTGATGAACTGGATGTCCCGGCACAAGTTAGTGAAAAAAAATTAACCCCGCCGCCGGGTAATAGCAGTGACAACCTCGAGCAACAGATAGCCAGTACTTCACAGCAAATCGGGTCTCTGCTCGCCGAAGATATGAACAGCGAGCAAGCGGCAAAAATGGCGC‑GTGGATGGGCCTCTTCTCAGGCTTATGGCGCAATGACAAACTGGGTAAGCCGCTTCGGTa      >  2:540342/1‑227 (MQ=255)                                                                                                                                                                                                                            ggggGATGAACTGGATGTCCCCGCAAAAGTTAGTGAAAAAAAATTAACCCCGCCGCCGGGTAATAGCAGTGACAACCTCGAGGAACAGATAGCCAGTACTTCACAGCAAATCGGGTCTCTGCTCGCCGAAGATATGAACAGCGAGCAAGCGGCAAATATGGCGC‑GTGGATGGGCCTCTTCTCAGGCTTCAGGCGCAATGACAGACTGGTTAAGCCGCTTCGGTACCGc  <  2:275142/228‑1 (MQ=255)                                                                                                                                                                                                                                  |                                                                                                                                                                                                                               GTTCCCTATGGCTGCGGCAGCACAAGGTGTGGTAAACGCCGCAACCCAACAACCAGTTCCTGCACAAATTGCCATTGCAAATGCCAATACGGTGCCCTACACCCTTGGAGCGTTGGAATCGGCCCAAAGCGTTGCCGAACGTTTCGGTATTTCGGTGGCTGAGTTACGCAAACTCAACCAGTTTCGTACGTTTGCTCGAGGTTTTGATAATGTCCGCCAGGGTGATGAACTGGATGTCCCGGCACAAGTTAGTGAAAAAAAATTAACCCCGCCGCCGGGTAATAGCAGTGACAACCTCGAGCAACAGATAGCCAGTACTTCACAGCAAATCGGGTCTCTGCTCGCCGAAGATATGAACAGCGAGCAAGCGGCAAATATGGCGC‑GTGGATGGGCCTCTTCTCAGGCTTCAGGCGCAATGACAGACTGGTTAAGCCGCTTCGGTACCGC  >  NC_000913/2045051‑2045497

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