Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	507,782	T→G	25.0%	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	25.0%	80.9 / 2.6	51	H100P (CAC→CCC)	ybaP	TraB family protein
Reads supporting (aligned to +/- strand):  ref base T (20/16);  new base G (6/9);  total (26/25)
Fisher's exact test for biased strand distribution p-value = 3.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

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