Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	507,782	T→G	23.1%	H100P (CAC→CCC)	ybaP ←	TraB family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	507,782	0	T	G	23.1%	94.9 / 2.1	57	H100P (CAC→CCC)	ybaP	TraB family protein
Reads supporting (aligned to +/- strand):  ref base T (22/19);  new base G (6/10);  total (28/29)
Fisher's exact test for biased strand distribution p-value = 3.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTCAGGGAGCTGGAGCAACATGGCAATCTGGTTTTCAGCCCCTTCCAGTTCAATCACGGGTTTATGTTGTTGCTTCGCCGCCTGCAATAGCTGGTAATCGATACCGTATTCTGCCCGCAGCCCCAGTTTTTGTGCCTGCGTCGCCTGAAGAACCATCGCGATTTGCCACAGCGGTTGGGTAGAAAAGAGTGAGGGAGAAATGCCCATCTCCTGGCTAATGTGCTGCAGGTTTTGGAGTTGC‑TCCTCGCTAATGCGC‑TCTTCCAGCGCCTCGCAGGCAGGCAAATTAGCAAAAGGCGTATCGCT‑GGTGGAAACATCCGCCTCGACGATCAGCGCATCGGCGTTTTTGAGCTTTTTGAGCAAACGGGTGGGCAGGGGAGCCATATCGTGGCTACCCATATGAATACTGCCAATCAGATGAAAATGGCGATTGCCGGGGAGGGTGATATCG  >  NC_000913/507561‑508009                                                                                                                                                                                                                              |                                                                                                                                                                                                                                       tggtaGGGTGGTGGATCAACTTGGCATTCTGGTTTTCAGCCCCTTCGTGGTCAACCGCGGGTTCATTTTTTTGACGCGCCGCCTGCAATAGCGGGTAATAGACACAGGATTCTGCCCGGAGCCCCAGTTTTTGTGCCTGAGTAGTCTGAAGAACCATCGCGATTTGCCACAGCGGTTGGGTAGAAAAGAGTGAGGGAGAAATGCCCATCTCCTGGCTAATGtgctgc                                                                                                                                                                                                                                   <  1:704144/223‑1 (MQ=255) tGTCAGGGAGCTGGAGCAACATGGCAATCTGGTTTTCAGCCCCTTCCAGTTCAATCACGGGTTTATGTTGTTGCTTCGCCACCTGCAATAGCTGGTAATCGATACCGTATTCTGCCCGCAGCCCCAGTTTTTGTGCCTGCGTCGCCTGAAGAACCATCGCGATTTGCCACAGCGGTTGGGTGGAAAAGGGGGAGGGAGAAATGCCCATCCCCTGGCTAATGTGCTGCa                                                                                                                                                                                                                                  >  2:657233/1‑228 (MQ=255)            tGGAGCAACATGGCAATCTGGTTTTCAGCCCCTTCCAGTTCAATCACGGGTTTATGTTGTTGCTTCGCCGCCTGCAATAGCTGGTAATCGATACCGTATTCTGCCCGCAGCCCCAGTTTTTGTGCCTGCGTCGCCTGAAGAACCATCGCGATGTGCCACAGCGGTTGGGTAGAAAAGAGTGAGGGAGAAATGCCCATCTCCTGGCTAATGTGCTGCAGGTTTTGGAGt                                                                                                                                                                                                                       <  2:556956/228‑1 (MQ=255)             ggAGCAACATGGCAATCTGGTTTTCAGCCCCTTCCAGTTCAATCACGGGTTTATGTTGTTGCTTCGCCGCCTGCAATAGCTGGTAATCGATACCGTATTCTGCCCGCAGCCCCAGTTTTTGTGCCTGCGTCGCCTGAAGAACCATCGCGATTTGCCACAGCGGTTGGGTAGGAAAGAGGGAGGGAGAAATGCCCACAACCTGGTTAATGGGCTGCAGGTTTTGGAGtt                                                                                                                                                                                                                      >  2:670352/1‑228 (MQ=255)                   aataggcaagtgtgTTTCCAGCCCCTTACAGTGTAATCACGGGTTTATGGTTTTGCTACGCCGCCTGCAATAGCTGGTAATCGATACCGTATTCTGCCCGCAGCCCCAGTTTTTGTGCCTGCGTCGCCTGAAGAACCATCGCGATTTGCCACAGCGGTTGGGTAGAAAAGAGTGAGGGAGAAATGCCCATCTCCTGGCTAATGTGCTGCAGGTGTTGGAGTTGc‑tc                                                                                                                                                                                                                 <  2:467828/213‑1 (MQ=255)                                          ccTTCCAGTTCAATCACGGGTTTATGTTGTTGCTTCGCCGCCTGCAATAGCTGGTAATCGATACCGTATTCTGCCCGCAGCCCCAGTTTTTGTGCCTGCGTCGCCTGAAGAACCATCGCGATTTGCCACAGCGGTTGGGTAGAAAAGAGTGAGGGAGAAATGCCCATCTCCTGGCTAATGGGCTGAAGGGTTTGGAGTTGC‑TCCTCGCTAAT‑CGCGTCTTCCAgcgc                                                                                                                                                                                        >  2:528857/1‑227 (MQ=255)                                                       tCACGGGTTTATGTTGTTGCTTCGCCGCCTGCAATAGCTGGTAATCGATACCGTATTCTGCCCGCAGCCCCAGTTTTTGTGCCTGCGTCGCCTGAAGAACCATCGCGATTTGCCACAGCGGTTGGTTAGAAAAGAGTGAGGGAGAAAGGCCCATCTCCTGGCTAATGGTCTGCAGGTTTTGGAGTTGC‑TCCTCGCTAATGCGC‑GCTTCCAGCGCCTCGACggcgggca                                                                                                                                                                          >  2:328146/1‑228 (MQ=255)                                                       tCACGGGTTTATGTTGTTGCTTCGCCGCCTGCAATAGCTGGTAATCGATACCGTATTCTGCCCGCAGCCCCAGTTTTTGTGCCTGCGTCGCCTGAAGAACCATCGCGATTTGCCACAGCGGTTGGGTAGAAAAGAGTGAGGGAGAAATGCCCTTCTCCTGGCTAATGTGCTGCAGGTTTTGGAGTTGC‑TCCTCGCTAATGCGC‑TCTTCCAGCGCCTCGCAggcaggca                                                                                                                                                                          >  2:130243/1‑228 (MQ=255)                                                        cagggttttagtttgttgCTTCGCCGCGTGCAATAGCGGGTAAGCGATACCGTATTCTGCCCGCAGCCCCAGTTTTTGTGCCTGCGTCGCCTGAAGAACCATCGCGATTTGCCACAGCGGTTGGGTAGAAAAGAGTGAGGGAGAAATGCCCATCTCCTGGCTAATGTGCTGCAGGTTTTGGAGTTGC‑TCCTCGCTAATGCGC‑TCTTCCAGCGCCTCGCAGGCAGGCaa                                                                                                                                                                         <  2:276634/216‑1 (MQ=255)                                                            ggTTTATGTTGTTGCTTCGCCGCCTGCAATAGCTGGTAATCGATACCGTATTCTGCCCGCAGCCCCAGTTTTTGTGCCTGCGTCGCCTGAAGAACCATCGCGATTTGCCACAGCGGTTGGGTAGAAAAGAGTGAGGTGGAAATGCCCATCTCCTGGCTAATGGGCTGCAGGTTTTGGAGTTGC‑TCCTCGCTAATGCGC‑TCTTCCAGCCCCTCGCAGGCGGGCAAATTa                                                                                                                                                                     >  1:129151/1‑228 (MQ=255)                                                               ttttggtgttgCGTCGCCGCCTGGAATAGCTGGTAATCGATACCGTATTCTGCCCGCAGCCCCAGTTTTTGTGCCTGCGTCGCCTGAAGAACCATCGCGATTTGCCACAGCGGTTGGGTAGAAAAGAGTGAGGGAGAAATGCCCATCTCCTGGCTAATGTGCTGCAGGTTTTGGAGTTGC‑TCCTCGCTAATGCGC‑TCTTCCAGCGCCTCGCAGGCAGGCAAATTAGCa                                                                                                                                                                  <  2:671095/222‑1 (MQ=255)                                                                tttgttttTGGTGCGCCCCCGGCAAGAGCTGGTAAGCGATCCCGGATTCTGCCCGCAGCCCCAGTTTTTGTTCCTGCGTAGGCTGAAGAACCATCGCGATTTGCCACAGCGGTTGGGTAGAAAAGAGTGAGGGAGAAATGCCCATCTCCTGGCTAATGTGCTGCAGGTTTTGGAGTTGC‑TCCTCGCTAATGCGC‑TCTTCCAGCGCCTCGCAGGCAGGCAAATTAGCaa                                                                                                                                                                 <  1:541768/226‑1 (MQ=255)                                                                       ttgCTTCGCCGCCTGCAATAGCTGGTAATCGATACCGTATTCTGCCCGCAGCCCCAGTTTTTGTGCCTGCGTCGCCTGAAGAACCATCGCGATTTGCCACCGCGGTTGGGTAGAAAAGAGTGAGGGAGAAATGCCCATCTCCTGGCTAATGTGCTGCAGGTTTTGGAGTTGC‑TCCTCGCTAATGCGC‑TCTTCCAGCGCCTCGCAGGCAGGCAAATTAGCAAAAGGCGt                                                                                                                                                          >  2:121399/1‑228 (MQ=255)                                                                       ttgCTTCGCCGCCTGCAATAGCTGGTAATCGATACCGTATTCTGCCCGCAGCCCCAGTTTTTGTGCCTGCGTCGCCTGAAGAACCATCGCGATTTGCCACAGCGGTTGGGTAGCAAAGAGTGAGGGAGAAATGCCCATCTCCTGGCTAATGTGCTGCAGGTTTTGGAGTTGG‑TCCTCGCTAATGCGC‑TCTTCCAGCGCCTCGACGGCAGGCAAGTTAGCAAAAGGCGt                                                                                                                                                          >  2:677512/1‑228 (MQ=255)                                                                       ttgCTTCGCCGCCTGCAATAGCTGGTAATCGATACCGTATTATGCCCGCAGCCCCAGTTTTTGTGCCTGCGTCGACTGAAGAACAATCGCGATTTGCCACAGCGGTTGGGTAGAAAAGAGTGAGGGAGAAATGCCCATCTCCCGGCTAATGGGCTGCAGGTTTTGGAGGTGC‑GCCTCGCTAATGCGC‑TCTTCCAGCGCCTCGCAGGCAGGCCAATTAGCCAAAGACGt                                                                                                                                                          >  2:483645/1‑228 (MQ=255)                                                                            tCGCCGCCTGCAATAGCTGGTAATCGATACCGTATTCTGCCCGCAGCCCCAGTTTTTGTGCCTGCGTCGCCTGAAGAACCATCGGGATTTGCCACAGCGGTTGGGTAGAAAAGAGTAGGGGGGAAATGCCCATCCCCCGGCTAATGGGCTGCAGGGTTTGGAGTTGG‑TCCTCGCTAATGCGC‑CCTTCCACCGCCCGGCAGGCGGGGAAATTAGCACAAgagggatcg                                                                                                                                                      >  1:234475/1‑219 (MQ=255)                                                                            tCGCCGCCTGCAATAGCTGGTAATCGATACCGTATTCTGCCCGCAGCCCCAGTTTTTGTGCCTGCGTCGCCTGAAGAACCATCGCGATTTGCCACAGCGGTTGGGTAGAAAAGAGTGAGGGAGAAATGCCCATCTCCTGGCTAATGTGCTGCAGGTTTTGGAGTTGC‑TCCTCGCTAATGCGC‑TCTTCCAGCGCCTCGCAGGCAGGCAAATTAGCAAAAGGCGAATCGc                                                                                                                                                     >  2:566059/1‑228 (MQ=255)                                                                                    tGCAATAGCTGGTAATCGATACCGTATTCTGCCCGCAGCCCCAGTTTTTGTGCCTGCGTCGCCTGAAGAACCATCGCGATTTGCCACAGCGGTTGGGTAGAAAAGAGTGAGGGAGAAATGCCCATCTCCTGGCTAATGtgc                                                                                                                                                                                                                                      >  1:143421/1‑141 (MQ=255)                                                                                     gCAATAGCTGGTAATCGATACCGTATTCTGCCCGCAGCCCCAGTTTTTGTGCCTGCGTCGCCTGAAGAACCATCGCGATTTGCCACAGCGGTTGGGTAGAAAAGAGTGAGGGAGAAATGCCCATCTCCTGGCTAATGtgct                                                                                                                                                                                                                                     <  2:143421/141‑1 (MQ=255)                                                                                              ggTAATCGATACCGTATTCTGCACGCAGCCCCAGTTTTTGTGCCTGCGTCGCCTGAAGAACCATCGCGATTTGCCACAGCGGTTGGGTAGAAAGGAGTGAGGGAGAAATGCCCATCTCCTGGCTAATGTGCTGCCGGGTTTGGAGTTGC‑TCCTCCCTAATTCGC‑TCTTCCAGCTCCTGCCAGGCAGGCAAATTACCAAAAGGCGTATCGCT‑CGGGGACAAATcgcccc                                                                                                                                  >  1:240131/1‑223 (MQ=255)                                                                                                taacggATACCTTATTCTGCCCGCAGCCCCAGTTGTTGGGCCTGGGGCGCCTGAAGAACCATCGCGAATTGCCACAGCGGTTGGGTAGAAAAGAGTGAGGGAGAAATGCCCATCTCCTGGCTAATGTGCTGCAGGTTTTGGAGTTGC‑TCCTCGCTAATGCGC‑TCTTCCAGCGCCTCGCAGGCAGGCAAATTAGCAAAAGGCGTATCGCT‑GGTGGAAACATCCGCCTcg                                                                                                                                <  2:212789/223‑1 (MQ=255)                                                                                                         ccGTATTCTGCCCGCAGCCCCAGTTTTTGTGCCTGCGTCGCCTGAAGAACCATCGCGATTTGCCACAGCGGTTGGGTAGAAAAGAGTGAGGGAGAAATGCCCATCTCCTGGCTAATGTGCTGCAGGTTTTGGAGTTGC‑TCCTCGCTAATGCGC‑TCTTCCAGCGCCTCGCAGGCAGGCAAATTACCAAAAGTCGGATCGCT‑GGTGGAACACTCCGCCTCGACGATAAgc                                                                                                                       >  2:672477/1‑228 (MQ=255)                                                                                                                 tGCCCGCAGCCCCAGTTTTTGTGCCTGCGTCGCCTGAAGAACCATCGAGATTTGCCACAGCGGTTGGGGAGAAAAGAGTGAGGGAGAAATGCCCATCTCCTGGCTAATGTGCTGCAGGTTTTGGAGTTGC‑TCCTCGCTAATGCGC‑TCTTCCAGCGCCTCGCCGGCAGGCAAATTAGCAATAGGAGTATCGCT‑GGTGGAAACATCGGCCTCGACGCTCAGCGCATCGGc                                                                                                               >  2:527711/1‑228 (MQ=255)                                                                                                                   cccGCAGCCCCAGTTTTTGTGCCTGCGTCGCCTGAAGAACCATCGCGATTTGCCACAGCGGTTGGGTAGAAAAGAGTGAGGGAGAAATGCCCATCTCCTGGCTAATGTGCTGCAGGTTTTGGAGTTGC‑TCCTCGCTAATGCGC‑TCTTCCAGCGCCTCGCAGGCAGGCAAATTAGCAAAAGGCGTATCGCT‑GGTGTAAACATCCGCCTCGAGGATCAGAGCATCGGCGt                                                                                                             >  2:480314/1‑228 (MQ=255)                                                                                                                    ccGCAGCCCCAGTTTTTGTGCCTGCGTCGCCTGAAGAACCATCGCGATTTGCCACAGCGGTTGGGTAGAAAAGAGTGAGGGAGAAATGCCCATCTCCTGGCTAATGTGCTGCAGGTTTTGGGGTTGC‑TCCTCGCTAATGGGC‑CCTTCCAGCGCCTCGCAGGCCGGCAAATTAGCAAAAGGCGTATCGTC‑GGGGGAAACCTCCCCCCCGACGATCAGCGAACCAGCGtt                                                                                                            >  1:686695/1‑228 (MQ=255)                                                                                                                    ccGCAGCCCCAGTTTTTGTGCCTGCGTCGCCTGAAGAACCATCGCGATTTGCCACAGCGGTTGGGTAGAAAAGAGTGAGGGAGAAATGCCCATCTCCTGGCTAATGTGCTGCAGGCTTTGGAGTTGC‑TCCTCGCTAATGCGC‑TCTTCCAGCGCCTCGCAGGCAGGCAAATTAGCAAAAGGCGTATCGCT‑GGTGGAAACATCCGCCTCGCCGATCAGCGCATCGGCGtt                                                                                                            >  1:482964/1‑228 (MQ=255)                                                                                                                        aGCCCCAGTTTTTGTGCCTGCGTCGCCTGAAGAACCATCGCGATTTGCCACAGCGGTTGGGTAGAAAAGAGTGAGGGAGAAATGCCCATCTCCTGGCTAATGTGCTGCAGGTTTTGGAGTTGC‑TCCTCGCTAATGCGC‑TCTTCCAGCGCCTCGCAGGCAGGCAAATTAGCAAAAGGCGTATCGCT‑GGTGGAAACATCCGCCTCGCCGATCAGCGCATCGGCGtttttg                                                                                                        >  2:336881/1‑228 (MQ=255)                                                                                                                                tgtttGTGCCTGGCGCGCCCGAAGAACCATCGCGATTTGCCACAGGGGTTGGGTAGAAAAGAGTGAGGGAGAAATGCCCATTACCTGGCTAATGTGCTGCAGGTTTTGGAGTTGC‑TCCTCGCTAATGCGC‑TCTTCCAGCGCCTCGCAGGCAGGCAAATTAGCAAAAGGCGTATCGCT‑GGTGGAAACATCCGCCTCGACGATCAGCGCATCGGCGtttttgagcttttt                                                                                                <  2:522067/226‑1 (MQ=255)                                                                                                                                          tGCGTCGCCTGAAGAACCATCGCGATTTGCCACAGCGGTTGGGTAGAAAAGAGTGAGGGAGAAATGCCCATCTCCTGGCTAATGTGCTGCAGGTTTTGGAGTTGC‑TCCTCGCTAATGCGC‑TCTTCCAGCGCCTCGCAGGCAGGCAAATTAGCAAAAGGCGTATCGCT‑GGTGGAAACATCCGCCTCGACGATCAGCGCATCGGCGTTTTTGAGCTTTTTGGGCAAACgc                                                                                      >  2:160110/1‑227 (MQ=255)                                                                                                                                            gcaggCCTGAAGAACCATCGCGATTTGCCACAGCGGTTGGGTAGAAAAGAGTGAGGGAGAAATGCCCATCTCCTGGCTAATGTGCTGCAGGTTTTGGAGTTGC‑TCCTCGCTAATGCGC‑TCTTCCAGCGCCTCGCAGGCAGGCAAATTAGCAAAAGGCGTATCGCT‑GGTGGAAACATCCGCCTCGACGATCAGCGCATCGGCGTTTTTGAGCTTTTTGAGCAAACggg                                                                                     <  2:234475/223‑1 (MQ=255)                                                                                                                                            cGTCGCCTGAAGAACCATCGCGATTTGCCACAGCGGTTGGGTAGAAAAGAGTGAGGGAGAAATGCCCATCTCCTGGCTACTGTGCTGCAGGTTTTGGAGTTGC‑TCCTCGCTAATGCGC‑TCTTCCAGCGCCTCGCAGGCAGGCAAATTAGCAGAAGCCGTATCGCC‑GGTGGAAGCATCAGCCTCGACGATCAGCGCATCGGCGCTTTTGAGCTTTTTTAGCAAACggc                                                                                     >  2:245207/1‑226 (MQ=255)                                                                                                                                            cGTCGCCTGAAGAACCATCGCGATTTGCCACAGCGGTTGGGTAGAAAAGAGTGAGGGAGAAATGCCCATCTCCTGGCTAATGTGCTGCAGGTTTTGGAGTTGC‑TCCTCGCTAATGCGC‑TCTTCCAGCGCCTAGCAGGCAGGCAAATTAGCAAAAGGCGGATCGCTGGGTGGAAACATCCGCCTCGCCGATCAGCGCATCGGCGGTTTTGAGCTTTTTGAGCAAACgg                                                                                      >  2:517336/1‑227 (MQ=255)                                                                                                                                                    taaaaaCCTTCGCGATTTTCCACAGCGGTTGGGTAGAAAAGTGTGGGGGATAAAGCCCCATCTGCGGGCTAATGTGCGGCAGGTTTTGGAGTTGC‑TCCTCGCAAATGCGC‑TCTTCCAGCGCCTCGCAGGCAGGCAAATTAGCAAAAGGCGTTTCGCT‑GGTGGAAACACCCGCCTCGGCGATCAGCGCATCGGCGTTTTTGAGCTTTTTGAGCAAACGGGTGGGCAggg                                                                            <  1:421487/224‑1 (MQ=255)                                                                                                                                                    gaagtaCCATCGCGTTTTGCCGCTGGGGTGGGGTAGAAAAGAGTGAGGGAGAAATGCCCATCTCCTGGCTAAGGTGCTGCAGGGTTTGGAGTTGC‑TCCTCGCTAATGCGC‑CCTGCCAGCGCATCGCAGGCAGGCAAATTAGCAAAAGGGGTATCGTT‑GGTTGAAACATCCGCCTCGACGATCAGCGCATCGGCGTTTTTGAGCTTTTTGAGCAAACGGGTGGGCAggg                                                                            <  1:355854/228‑1 (MQ=255)                                                                                                                                                     aagcaccaacctgtttttCCACGGGGGTTGGGTGAAAAAGAGGTGGGGAGAAAGGCCCATCTCCGGGCTAAGGGGCGGCAGGGTTTGGAGTTGC‑TCCTCGATAATGGGC‑TCTTCCGGCGCCTCGCAGGCAGGGAAATTAGCAAAAGGCGTATCGCT‑GGTGGAAACATCCGCCTCGGCGATCAGCGCATCGGCGTTTTTGAGCTTTTTGAGCAAACGGGTGGGCAggg                                                                            <  1:657233/213‑1 (MQ=255)                                                                                                                                                      cgaaCCATCGCGATTTGCCACAGAGGTTGGGAAGAAAAGAGGGAGGGAGAAATGCCCATCTCCTGGCTAAGGTGCTGCAGGCTTTGGAGTTGC‑ACCTCGCTAATGCGC‑TCTTCAAGCGCCTAGCAGGCAGGCAAATTAGCAAAAGGCGTATCGCT‑GGTGGGAAAATCGGCCTCACCGATCAGAGAAACGGGGTTTTTGAGCTTTGTGATCAAACGGTGGGGGAggggg                                                                          >  1:450043/2‑227 (MQ=255)                                                                                                                                                        aaCCATCGGGATTTGCCACAGCGGTTGGGTAGAAAAGAGTGAGGGAGAAATGCCCATCTCCTGGCTAATGTGCTGCAGGTTTTGGAGTTGC‑TCCTCGCTAATGCGC‑TCTTCCAGCGCCTCGCAGGCAGGCAAATTAGCAAAAGGCGTATCGCT‑GGTGGAAACATCCGCCTCGACGATCAGCGCATCGGCGTTTTTGAGCTTTTTGAGCAAACGGGTGGGCAGGGGAg                                                                         <  2:150171/227‑1 (MQ=255)                                                                                                                                                              cgcgATTTGCCACAGCGGTTGGGTAGAAAAGAGTGAGGGAGAAATGCCCATCTCCTGGCTAATGTGCTGCAGGTTTTGGAGTTGC‑TCCTCGCTAATGCGC‑TCTTCCAGCGCCTCGCAGGCAGGCAAATTAGCAAAAGGCGTATCGCT‑GGTGGAAACATCCGCCTCGACGCTCAGCGCATAGGCGTTTTTGAGCTTTTTGTGCCAACGGGTGGGCAGGGGAGCCATATc                                                                  >  2:615558/1‑228 (MQ=255)                                                                                                                                                               gcgATTTGCCACAGCGGTTGGGTAGAAAAGAGTGAGGGAGAAATGCCCATCTCCTGGCTAATGTGCTGCAGGTTTTGGAGTTGC‑TCCTCGCTAATGCGC‑TCTTCCAGCGCCTCGCAGGCAGGCAAATTAGCAAAAGGCGTATCGCT‑GGTGGAAACATCCGCCTCGACGATCAGCGCATCGGCGTTTTTGAGCTTTTTGAGCAAACGGGTGGCCAGGGGGGCCATATCg                                                                 >  2:481996/1‑228 (MQ=255)                                                                                                                                                                    ttGCCACAGCGTTTGGGTAGAAAAGAGTGAGGGAGAAATGCCCATCTCCTGGCTAATGTGCTGCAGGTTTTGGAGTTGC‑TCCTCGCTAATGCGC‑TCTTCCAGCGCCTCGCAGGCAGGCAAATTAGCAAAAGGCGTATCGCT‑GGTGGAAACATCCGCCTCGACGATCAGCGCATCGGCGTTTTTGAGCTTTTTGAGCAAACGGGTGGGCAGGGGAGCCATATCGTgg                                                              <  2:473315/226‑1 (MQ=255)                                                                                                                                                                    tggCCACAGCTGTTGGGTAGAAAAGAGTGAGGGAGAAATGCCCGTCTCCTGGCTAATGGGCTGCAGGTTATGGAGTTGG‑TCCTCGCTAATGCGC‑TCTTCCAGCGCCGCGCAGGCAGGCAAATTAGCAAAAGGCGTATCGCT‑GGTGGAAACATCCGCCTCGACGATCAGCGCATCGGCGTTTTTGAGCTTTTTGAGCAAACGGGTGGGCAGGGGAGCCATATCGTgg                                                              <  2:30880/224‑1 (MQ=255)                                                                                                                                                                     tGCCACAGCGGTTGGGTAGAAAAGAGTGAGGGAGAAATGCCCATCTCCTGGCTAATGTGCTGCAGGTTTTGGAGTTGC‑TCCTCGCTAATGCGC‑TCTTCCAGCGCCTCGCAGGCAGGCAAATTAGCAAAAGGCGTATCGCT‑GGTGGAACCATCCGCCTCGACGATCAGCGCATCGGCGTTTTTGAGCTTTTTGAGCACACGGGTgg                                                                                  >  2:411820/1‑205 (MQ=255)                                                                                                                                                                     tGCCACAGCGGTGGGGTAGAAAAGGGTGAGGGGGAAATGCCCATCTCCTGGCTAATGTGCGGCAGGTTTTGGAGTTGC‑TCCTCGCTAATGCGC‑TCTTCCAGCGCCTCGCAGGCAGGCAAATTAGCAAAAGGCGTATCGCT‑GGTGGAAACATCCGCCTCGACGATCAGCGCATCGGCGTTTTTGAGCTTTTTGAGCAAACGGGTGGGCAGGGGAGCCATATCGTGGCTa                                                           <  2:198466/228‑1 (MQ=255)                                                                                                                                                                      tataagaggaaggggTAGAAAAGAGGGGGGGATAAATCCCCCTCGCCGGGGTAATGTGCTGCAGGTTTTGGAGTTGT‑CCCGCGCTAATGCGC‑TCTTCCAGCGCCTCGCAGGCAGGCAAATTAGCAAAAGGCGTATCGCT‑GGGGGAACCATCCGCCTCGACGAGCAGCGCATCGGCGTTTTTGAGCTTTTTGAGCAAACGGGTGGGCAGGGGAGCCATATCGTGGCTa                                                           <  1:454721/215‑1 (MQ=255)                                                                                                                                                                      gCCACGGGGGTTGGGTAAGAAAGGGGTAGGGAGAAAGCCCCATCTCCTGGCTAATGGGCTGCAGGTTTGGGAGTTGC‑TCCCCGGTAATGCGC‑TCTTCCAGCGCCTCGCAGGCAGGCAAATTAGCAAAAGGCGTATCGGT‑GGTGGAAACATCCGCCTCGACGATCAGCGCATCGGCGTTTTTGAGCTTTTTGAGCAAACGGGTggg                                                                                 <  1:411820/205‑1 (MQ=255)                                                                                                                                                                       tccgaGCGGTTGGTTGGAAAAGGGTGGGGGGGAAATGCACCTTGCCGGGTTAATAGGCTGGAGGTGTTGGAGTTGC‑TCATGGCTAATGCAG‑TCTTCCGGCGCCTCGCAGGCAGGAAAATTAGCAAAAGGGGTATCGGT‑GGTGGAACCAGCCGCCTGGACGCTCAGCGCAGCGGCGTTTTTGAGCTTTTTGAGCAAACGGGTGGGCAGGGGAGCCATATCGTGGCTAcc                                                         <  1:566059/224‑1 (MQ=255)                                                                                                                                                                       ctctgggggTTGGGTGGGAATGAGGGGGGGAGAATCGCCCATCTCGTGCTTATGGTGCTGCAGGTTGTGGAGTTGC‑TCCTGGATAAGGCGC‑TCTTCCAGCGCCTCGCAGGCTGGCAAATTAGCAAAAGGCGTATCGCT‑GGTGGAACCATCCGCCTCGACGATCAGCGCAGCGGCGTTTTTGAGCTTTTTGAGCAAACGGGTGGGCAGGGGAGCCATATCGTGGCTAcc                                                         <  1:483645/221‑1 (MQ=255)                                                                                                                                                                       ccACAGCGGTTGGGTAGAAAGGAGTGAGGGAGAAATGCCCATCTCCTGGCTAATGTGCTGCAGGTTTTGGAGTTGC‑TCCTCGCTAATGCGC‑TCTTCCAGCGCCTCGCAGGCAGGCAAATTAGCAAAAGGCGTATCGCT‑GGTGGAAACATCCGCCTCGACGATCAGCGCATCGGCGTTTTTGAGCTTTTTGAGCAAACGGGTGGGCAGGGGAGCCATATCGTGGCTAcc                                                         <  2:570641/228‑1 (MQ=255)                                                                                                                                                                          ctggggTGGGGTAGATAAGGGTGAGGGAAAAATGCCCATTCCCTGGTTAAGGGGCGGCAGGTTTGGGATGTGC‑TCCTCGCTAATGCGC‑TCTTCCGGCGCCTCGCAGGCAGGTAAATTAGCAAAAGGCGTATCGCT‑GGTGGAACCATCCGCCTCGACGATCAGCGCATCGGCGTTTTTGAGCTTTTTGAGCAAACGGGTGGGCAGGGGAGCCATATCGTGGCTACCCat                                                      <  1:458240/224‑1 (MQ=255)                                                                                                                                                                          caGCGGTGGGGTAGAAAAGAGTGAGGGAGAAATGCCCATCTCCTGGCTAATGTGCTGCAGGTTTTGGAGTTGC‑TCCTCGATAATGCGC‑TCTTCCAGCGCCTCGCAGGCAGGCAAATTAGCAAAAGGCGTATCGCT‑GGTGGAAACATCCGCCTCGACGATCAGCGCATCGGCGTTTTTGAGCTTTTTGAGCAAACGGGTGGGCAGGGGAGCCATATCGTGGCTACCCat                                                      <  2:83250/228‑1 (MQ=255)                                                                                                                                                                           aGCGGTTGGGTAGAAAAGAGTGAGGGAGAAATGCCCATCTCCTGGCTAATGTGCTGCAGGTTTTGGAGTTGC‑TCCTCGCTAATGCGC‑TCTTCCAGCGCCTCGCAGGCAGGCAAATTAGCAAAAGGCGTATCGCT‑CGGGGAACCATCCCCCTCGCCGATCACGGAATAGGCGTTTTTTAGGGTGTTGAGCAACCGGGTGGGCAGGGGGGcactatgatggggcccaaaa                                                     >  1:246626/1‑209 (MQ=255)                                                                                                                                                                               gTTGGGTAGAAAAGAGTGAGGGAGAAAGGCCCATCTCCTGGCTAATGGGCTGCAGGTTTTGGAGTTGA‑TCCTCGCTAAGGCGC‑TCTTCCAGCGCCTCGCAGGCAGGCAAATTAGCAAAAGGCGTATCGCT‑GGTGGAAACATCCGCCTCGACGATCAGCGCATCGGCGTTTTTGAGCTTTTTGAGCAAACGGGTGGGCAGGGGAGCCATATCGTGGCTACCCATATGaa                                                 <  2:380137/228‑1 (MQ=255)                                                                                                                                                                                         aaaGGGTGGGGGAAAAATCCCCATCTCCTGGCTAATGGGCGGCAGGTTTTGGGTTTGC‑TCCTCGATAATGCGC‑TCTTCCAGCGCCTCGCAGGCAGGCAAATTAGCAAAAGGCGTATCGCT‑GGTGGAAACAGCCGCCTCGACGATCAGCGCATCGGCGTTTTTGAGCTTTTTGAGCAAACGGGTGGGCAGGGGAGCCATATCGTGGCTACCCATATGAATACTGCCAAt                                       <  1:481996/228‑1 (MQ=255)                                                                                                                                                                                          aaGAGTGAGGGAGAAATGCCCATCTCTGGGCGAATGTGCTGCAGGTTTTGGAGTGGC‑TCCTAGCTAATGCGA‑TCTGCCAGCGCCTCGCAGGCAGGTAAATTAGCAAAAGGCGTATCGCT‑GGTGGAAACATCCGCCTCGACGATCAGCGCATCGGCGTTTTTGAGCTTTTTGAGCAAACGGGTGGGCAGGGGAGCCATAGCGTGGCTACCCATATGAATACTGCCAAt                                       <  2:527939/227‑1 (MQ=255)                                                                                                                                                                                             attgggggAGAAATGCCCATCTCCGGGCTAATGGGCTGCAGGTTTTGGAGTTGC‑TCCGCGCTAATGCGC‑TCTTCCAGCGCCTCGCAGGCAGGCAAATTAGCAAAAGGCGTATCGCT‑GGTGGAAACATCCGCCTCGACGATCAGCGCATCGGCGTTTTTGAGCTTTTTGAGCAAACGGGTGGGCAGGGGAGCCATATCGTGGCTACCCATATGAATACTGCCAATCAGa                                   <  2:392091/223‑1 (MQ=255)                                                                                                                                                                                             agTGAGGGAGAAATGCCCATCTCCTGGCTAATGTGCTGCAGGTTTTGGAGTTGC‑TCCTCGCTAATGCGC‑TCTTCCAGCGCCTCGCAGGCAGGCAAATTAGCAAAAGGCGTATCGCT‑GGTGGAAACATCCGCCTCGACGATCAGCGCATCGGCGTTTTTGAGCTTTTTGAGCAAACGGGTGGGCAGGGGAGCCATATCGTGGCTACCCATATGAAtctgtctcttatc                                    >  2:454721/1‑215 (MQ=255)                                                                                                                                                                                              gTGAGGGAGAAATGCCCATCTCCTGGCTAATGTGCTGCAGGTTTTGGAGTTGC‑TCCTCGCTAATGCGC‑TCTTCCAGCGCCTCGCAGGCAGGCAAATTAGCAAAAGGCGTATCGCT‑GGTGGAAACATCCGCCTCGACGATCAGAGCATAGGCGTTTTTGAGCTTTTTGAGCAACCGGGTGGGCAGGGGAGCCCGATCGGGGCTACCCATATGAATACCGcaatacacag                                  >  1:657888/1‑219 (MQ=255)                                                                                                                                                                                                         atccccATCCCCTGGCTAATGTGCTGCAGGTTTGGGAGTTGC‑CCCTCGCAAAGGCGC‑TCTTCCAGCGCCTCGCAGGCAGGCAAATTAGCAAAAGGCTTATCGCT‑GGTGGAAACATCCGCCTCGACGATCAGCGCATCGGCGTTTTTGAGCTTTTTGAGCAAACGGGTGGGCAGGGGAGCCATATCGTGGCTACCCATATGAATACTGCCAATCAGATGAAAATGGc                         <  2:540639/223‑1 (MQ=255)                                                                                                                                                                                                         acgCCCCTCTCCGTGGTAATGGGCGTCGGGTTTTGGATTTGC‑TCCGAGCTAACGTGT‑TTGCCCAGCGCCGCGCCGGCAGGGAAATTAGCAAAAGGCGTATCGCG‑GGTGGAAACATCCGCCTCGACGACCGGCGCAGCGGCGTTTTTGAGCTTTTTGAGCAAACGGGTGGGCAGGGGAGCCATATCGTGGCTACCCATATGAATACTGCCAATCAGATGAAAATGGCGa                       <  1:677512/226‑1 (MQ=255)                                                                                                                                                                                                          cccccAGCTCCTGGCTGAGTGGCTGCGGGTTGTAG‑GATGCGTCCTAGCTAATGCGC‑TCTTACAGCGCCTCGCAGGCAGGCAAATTAGCAAAGGGCGTGTCGGT‑GGTGTAAACATCCGCCTCGACGATCAGCGCATCGGAGTTTTTGAGCTTTTTGAGCAAACGGGTGGGCAGGGGAGCCATATCGTGGCTACCCATATGAATACTGCCAATCAGATGAAAATGGCGa                       <  1:506594/225‑1 (MQ=255)                                                                                                                                                                                                                            tgtgCTGCAGGTTTTGGAGTTGC‑TCCTCGCTAATGCGC‑TCTTCCAGCGCCTCGCAGGCAGGCAAATTAGCAAAAGGCGTATCGCT‑GGTGGAAACATCCGCCTCGACGATCAGCGCATCGGCGTTTTTGAGCTTTTTGAGCCAACGGGTGGGCAGGGGAGCCATATCGTGGCTACCCATATGAATACTGCCAATCAGAGGAAAATGGCGAATGCAGGGGGGGGGGaaaa    >  1:12232/1‑225 (MQ=255)                                                                                                                                                                                                                              tgctgcAGGTTTTGGAGTTGC‑TCCTCGCTAATGCGC‑TCTTCCAGCGCCTCGCAGGCAGGCAAATTAGCAAAAGGCGTATCGCT‑GGTGGAAACATCCGCCTCGACGATCAGCGCATCGGCGTTTTTGAGCTTTTTGAGCAAACGGGTGGGCAGGGGAGCCCTATCGTGGCTACCCATATGAATACTGCCAATCAGATGAAAATGGCGATTGCCGGGGCGGGGGataccc  >  2:487854/1‑225 (MQ=255)                                                                                                                                                                                                                              |                                                                                                                                                                                                                                       TGTCAGGGAGCTGGAGCAACATGGCAATCTGGTTTTCAGCCCCTTCCAGTTCAATCACGGGTTTATGTTGTTGCTTCGCCGCCTGCAATAGCTGGTAATCGATACCGTATTCTGCCCGCAGCCCCAGTTTTTGTGCCTGCGTCGCCTGAAGAACCATCGCGATTTGCCACAGCGGTTGGGTAGAAAAGAGTGAGGGAGAAATGCCCATCTCCTGGCTAATGTGCTGCAGGTTTTGGAGTTGC‑TCCTCGCTAATGCGC‑TCTTCCAGCGCCTCGCAGGCAGGCAAATTAGCAAAAGGCGTATCGCT‑GGTGGAAACATCCGCCTCGACGATCAGCGCATCGGCGTTTTTGAGCTTTTTGAGCAAACGGGTGGGCAGGGGAGCCATATCGTGGCTACCCATATGAATACTGCCAATCAGATGAAAATGGCGATTGCCGGGGAGGGTGATATCG  >  NC_000913/507561‑508009

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