Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	65,785	A→T	intergenic (‑5/+70)	polB ← / ← araD	DNA polymerase II/L‑ribulose‑5‑phosphate 4‑epimerase AraD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	65,785	0	A	T	94.4%	163.0 / ‑0.2	56	intergenic (‑5/+70)	polB/araD	DNA polymerase II/L‑ribulose‑5‑phosphate 4‑epimerase AraD
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base T (20/31);  total (24/32)
Fisher's exact test for biased strand distribution p-value = 5.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.47e-01

TACACCGGCTGGCGGTGAAAA‑TCCTTTAACGCCAGCGGTGTCAGGCGAAAGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGGAATAAACGCCACGGACTCTTGCGGTG‑CAAGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCT‑TGCGGGGTGTCCCGCCAGTGTCGGGTTAAGA‑TAAAACCTGCCTGCGCCACGCTG‑AAAATCC‑ATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGC‑AGATAGTGTTTATC‑CAGCAGCGTTTGCTGC‑ATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCCATGCAA  >  NC_000913/65545‑66030                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                          aacacCGGCTGGCG‑TGAAAAACCCTTTAACGCCAGCTGTGTCAGGCGATCGCCTTGTTCACCCTGAAAAATAGCCTGAGCGCGGGGAACCTGATCGGCGGGAATAAACGCCACGGACTCTTGCGGTG‑CAAGCGTAACCGTAAACGTCCCGGTGGCCGTCGCCAGCCATAAGGAGACTGCTGTCCCT‑TGCGGTGTGTCCCGCCAGTGTCGGGTTAAGA‑TAAAGCCTGCCTGCGCCACGCTG‑TAATTCC‑ATCa                                                                                                                                                                                                                                                <  2:403798/250‑1 (MQ=255)  cccccgggggggggTGAAAA‑TCCTTAAACGCCAGCGGGTTCAGGCGAACCCCTTGTTCCCCCTGAAAAATATGCTGGGCGGGGGGAACCTGACCGGGGGGAATAACCGCCAGGGACTCTTGCGGTG‑AAAGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCT‑TGCGGGGTGCCCCGCCAGTGTCGGGTTAAGA‑TAAAACCTGCCTGCGCCACGCTG‑TAATTCC‑ATCaa                                                                                                                                                                                                                                               <  2:138616/239‑1 (MQ=255)  acacCGGCTGGCGGTGAAAA‑TCCTTTAACGCCAGCGGTGTCAGGCGAACGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGGAATAAACGCCACGGACTCTTGCGGTG‑CAAGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCT‑TGCGGGGTGTCCCGCCAGTGTCGGGTTAAGA‑TAAAACCTGCCTGCGCCACGCTG‑TAATTCC‑ATCa                                                                                                                                                                                                                                                <  1:533494/250‑1 (MQ=255)  acacCGGCTGGCGGTGAAAA‑TCCTTTAACGCCAGCGGTGTCAGGCGAACGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGGAATAAACGCCACGGACTCTTGCGGTG‑CAAGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCT‑TGCGGGGTGTCCCGCCAGTGTCGGGTTAAGA‑TAAAACCTGCCTGCGCCACGCTG‑TAATTCC‑ATCa                                                                                                                                                                                                                                                <  1:484242/250‑1 (MQ=255)              ggTGAAAA‑TCCTTTAACGCCAGCGGTGTCAGGCGAACGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGGAATAAACGCCACGGACTCTTGCGGTG‑CAAGCGTAACCTGCAACAGCCCGTTGTCCGTCGCCAGCAAGAAGGAGACTTCTGTCCCT‑TGCGGGGTGTCCCGCCAGTGTCGGGGTAAGA‑TAAAAACCGCCTGCGCCCCGCTG‑TAATTCCCTTCAAAAAAC‑AGGCtt                                                                                                                                                                                                                                   >  1:380146/1‑251 (MQ=255)                            tAACGCCAGCGGTGTCAGGCGAACGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGGAATAAACGCCACGGACTCTTGCGGTG‑AAAACGTAACCCGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCT‑TGCGGGGTGTCCCGCCCGTGTCGGGGTAAAA‑TAAAACCTGCCTGCGCCACGCGGTAATTTCC‑A‑CAAAAAACCAGGCTTTATTATAGcccgg                                                                                                                                                                                                                      >  1:326304/1‑248 (MQ=255)                            tAACGCCAGCGGTGTCAGGCGAACGCCTTGTTAACCCTGCAAAATATGCGGAGCGCGGGGAACCTGACCGGCGGGAATAAACGCCACGGACTCTTGCGGTG‑CAAGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCT‑TGCGGGGTGTCCCGCCAGTGTCGGGTTAAGA‑TAAAACCTGCCTGCGCCACGCTG‑TAATTCC‑ATCAAAAAACCAGGCTTGATTATAGCCTg                                                                                                                                                                                                                       <  1:483051/250‑1 (MQ=255)                                  cAGCGGTGTCAGGCGAACGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGACCGGCGGGAATAAACGCCACGGACTCTGGCGGTG‑CAAGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCT‑TGCGGGGTGTCCCGCCAGTGTCGGGTTAAGA‑TAAAACCTGCCTGCGCCACGCTG‑TAATTCC‑ATCAAAAAACCAGGCTTGATTATAGCCTGGttt                                                                                                                                                                                                                   <  1:655918/248‑1 (MQ=255)                                                gatcgCCTTGTGCACCCTGAAAAATTGGCTGGGCGCGGGGAACCGGATCGGCGGGAATAAACGCCACGGCCGCGTGTGGGG‑AAAGGTTAACCTGCAACGGCCCGTTGGCCGGCGCCAGCAAGAAGGAGATTTCTGGCCCT‑TGCGGGGTGTCCCGCCAGTGTCGGGTTAAGA‑TAAATCCGGCCTGCGCCACGCTG‑TAATTCC‑ATCAAAAAACCAGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTgg                                                                                                                                                                                                  <  2:688177/247‑1 (MQ=255)                                                gaacgCCTTGTTCCCCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGGAATAAACGCCACGGACTCTTGCGGTG‑CAAGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCT‑TGCGGGGTGTCCCGCCAGTGTCGGGTTAAGA‑TAAAACCTGCCTGCGCCACGCTG‑TAATTCC‑ATCAAAAAACCAGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTgg                                                                                                                                                                                                  <  1:601165/247‑1 (MQ=255)                                                gaacgCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGGAATAAACGCCACGGACTCTTGCGGTG‑CAAGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCT‑TGCGGGGTGTCCCGCCAGTGTCGGGTTAAGA‑TAAAACCTGCCTGCGCCACGCTG‑TAATTCC‑ATCAAAAAACCAGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTgg                                                                                                                                                                                                  <  1:20548/247‑1 (MQ=255)                                                          gtCACCCTGCAAAATATGCGGAGCGCGGGTAGCCTGATCGGCGGAAATAAACGCCACGGACTCTTGCGGTG‑CAAGCGTAACCTGCAACGGCCCTTTGTCCGTCGCCAGCCAGAAGGAGACTGCTGTCCCT‑GGCGGGGTGTCCCGCCAGTGTCGGGTTAAGA‑TAAAACCTGCCTGCGCCACGCTG‑TAATTCC‑ATCAAAAAACCAGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAg                                                                                                                                                                                        <  1:109550/250‑1 (MQ=255)                                                             ccccggcaaaatatgtggggggcgGGGAACCTGAGCGGGGGGAATAAACGCCACGGACTCTGGCGTGG‑CAAGCGTAACCGGAAACGGCCCGTTGTCCGTCGCCAGCCAGAGGGAGACTTCTGTCCCT‑TGCGGGGTGTCCCGCCAGTGTCGGGTTAAGA‑TAAAACCTGCCTGCGCCACGCTG‑TAATTCC‑ATCAAAAAACCAGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCa                                                                                                                                                                                     <  2:677838/230‑1 (MQ=255)                                                               ccggaaaaaTATGCGGGGCGCGGGGAACGTGTCCGGGGGAAATAACCGCCCTGGCCTCTTCGGGGG‑CAAGCGTAACTTGCACCGGCCGGTGGTCCGCTGCCAGCCAGAAGGAGACTTCTGTCCCT‑TGCGGGTTGTCCCGCCAGTGTCGGGTTAAGA‑TAAAACCTGCCTGCGCCACGCTG‑TAATTCC‑ATCAAAAAACCAGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCAtt                                                                                                                                                                                   <  2:75618/246‑1 (MQ=255)                                                               agccgtgagctcttccgatcgggggAGCCTGATCGGCGGGAATAAACGCCACGGACGCTTG‑GGTGGCAAGCGAAGCATGCAACGGGCGGTTGTCCGTCGCCAGCAAGAAGGAGACTTCTGTCCCT‑TGCGGGGTGCCCCGCCAGTGTCGGGTTAAGA‑TAAAACCTGCCTGCGCCACGCTG‑TAATTCC‑ATCAAAAAACCAGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCAtt                                                                                                                                                                                   <  1:74293/230‑1 (MQ=255)                                                                          tGCTGAGCGCGGGGAACCTGATCGGCGGGAATAAACGCCACGGACTCTTGCGGTG‑CAAGCGTAACCTGCAACGGCCCGGTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCT‑TGCGGGGTGTCCCGCCCGTGTCGGGTTTAGA‑TAAAAACTGCCTGCGCCCCCCTG‑GAATTTC‑ATCAAAAAACCAAGCTTGATTTTATCCTGGGTTTGTTTGTTTTTCTGTGGGTTTTTACAGttcttattggccgta                                                                                                                                                                         >  1:521394/1‑236 (MQ=255)                                                                           gCTGAGCGCGGGGAACCTGATCGGCGGGAATAAACGCCACGGACTCTTGCGGTG‑CAAGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCT‑TGCGGGGTGTCCCGCCAGTGTCGGGTTAAGA‑TAAAACCTGCCTGCGCCACGCTG‑TAATTCC‑ATCAAAAAACCAGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTaa                                                                                                                                                                        <  1:64940/250‑1 (MQ=255)                                                                            cTGAGCGCGGGGAACCTGATCGGCGGGAATAAACGCCACGGACTCTTGCGGGG‑CAAACGTAACCTGCAACGGCCCGGTGGCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCC‑TGCGGGGTGGCCCGCCCGTGTCGGGTTAAGAGAAAAACCT‑CCTTCGCCACCCAG‑TAATTTC‑ATCAAAAAAACAGGCTTTATTATAACCTGGGTTCGTTTGATTGGGTGGGGTTTTATACAggtattaattccccttat                                                                                                                                                                       >  1:501655/1‑233 (MQ=255)                                                                                                gggcggGAATAAACGCCACGGACTCTGGCGGTG‑CAAGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCT‑TGCGGGGTGTCCCGCCAGTGTCGGGTTAAGA‑TAAAACCTGCCTGCGCCACGCTG‑TAATTCC‑ATCAAAAAACCAGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTAc                                                                                                                                                  <  1:277487/250‑1 (MQ=255)                                                                                                 ggcggAAATAAACGCCACGGACTCTGGCGGTG‑CAAGCGTAACCTGCAACGGCCCGTGGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCT‑TGCGGGGTGTCCCGCCAGTGTCGGGTTAAGA‑TAAAACCTGCCTGCGCCACGCTG‑TAATTCC‑ATCAAAAAACCAGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACg                                                                                                                                                 <  1:323794/251‑1 (MQ=255)                                                                                                          aaaCGCCACGGACTCTTGCGGTG‑CACGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCT‑TGCGGGGCGTCCCGCCAGTGACGGGGTAAGA‑TAAAACCTGCCTGCGCCCCGCTG‑TAATTCC‑ATCAACAAACCAGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGGGGTTTTATACATTCCTTACTGCCCGTAATAACACTTCGCGCCATGCTTACgacacaaggg                                                                                                                                        >  2:368701/1‑242 (MQ=255)                                                                                                          aaaCGCCACGGACTCTGGCGGTG‑CAAGCGTAACCTGCAACGGCCCGTTGTCCGGCGCCAGCCAGAAGGAGACTTCTGTCCCT‑TGCGGGGTGTCCCGCCAGTGTCGGGTTAAGA‑TAAAACCTGCCTGCGCCACGCTG‑TAATTCC‑ATCAAAAAACCAGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTta                                                                                                                                                                                             <  1:705372/198‑1 (MQ=255)                                                                                                            ccGCCACGTACTCTTGCGGTG‑CAAGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCTT‑TGCGGGATGTCCCGCCAGTGTCGGGTTAAGA‑TAAAACCTGCCTGCGCCACGCTG‑TAATTCC‑ATCAAAAAACCAGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGC‑AGATAGTGtt                                                                                                                                     <  1:121297/250‑1 (MQ=255)                                                                                                              gCCACGGACTCTTGCGGTG‑CAAGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCT‑TGCGGGGTGTCCCGCCAGTGTCGGGTTAAGA‑TAAAACCTGCCTGCGCCACGCTG‑TAATTCC‑ATCAAAAAACCAGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGC‑AGATAGTGtttt                                                                                                                                   >  1:134433/1‑250 (MQ=255)                                                                                                                 aCGGACTCTTGCGGTG‑CAAGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCT‑TGCGGGGTGTCCCGCCAGTGTCGGGTTAAGA‑TAAAACCTGCCTGCGCCACGCTG‑TAATTCC‑ATCAAAAAACCAGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCtt                                                                                                                                                                >  1:631855/1‑220 (MQ=255)                                                                                                                 aCGGACTCTTGCGGTG‑CAAGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCT‑TGCGGGGTGTCCCGCCAGTGTCGGGTTAAGA‑TAAAACCTGCCTGCGCCACGCTG‑TAATTCC‑ATCAAAAAACCAGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGC‑AGATAGTGTTTAttc                                                                                                                                >  1:473885/1‑249 (MQ=255)                                                                                                                                  caaacGTAACCTGCAACGGCCCGTTGTCCGCCGCCAGCCAGAAGGAGACTTCTGTCCCT‑TGCGGGGTGTCCCGCCAGTGTCGGGTTAAGA‑TAAAACCTGCCTGCGCCACGCTG‑TAATTCC‑ATCAAAAAACCAGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGC‑AGATAGTGTTTATT‑CAGCAGCGTTTGCTGC‑a                                                                                                              <  1:587924/247‑1 (MQ=255)                                                                                                                                  cAAGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCT‑TGCGGGGTGTCCCGCCAGTGTCGGGTTAAGA‑TAAAACCTGCCTGCGCCACGCTG‑TAATTCC‑ATCAAAAAACCAGGCTTGGTTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCtt                                                                                                                                                                >  1:368622/1‑204 (MQ=255)                                                                                                                                         aaCCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCT‑TGCGGGGTGTCCCGCCAGTGTCGGGTTAAGA‑TAAAACCTGCCTGCGCCACGCTG‑TAATTCC‑ATCAAAAAACCAGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGC‑AGATAGTGTTTATT‑CAGCAGCGTTTGCTGC‑ATATCCgg                                                                                                       <  1:89367/251‑1 (MQ=255)                                                                                                                                              gCAACGGCCCTTTGTCCTGCGCCAGCCAGAAGGAGACTTCTGTCCCT‑TGCGGGGTTTCCCGCCAGTGTCGGGTTAAGA‑TAAAACCTGCCTGCGCCACGCTG‑TAATTCC‑ATCAAAAAACCAGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGC‑AGATAGTGTTTATT‑CAGCAGCGTTTGCTGC‑ATATCCGGTaa                                                                                                    <  1:303213/249‑1 (MQ=255)                                                                                                                                                tcccgatcttggcccGTCCCCAGCCAGAAGGAGACTTC‑GGCCCTGTGCGGGTTGTCCCGCCATTGTCGGGTTAAGA‑TAAAACCTGCCTGCGCCACGCTG‑TAATTCC‑ATCAAAAAACCAGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGC‑AGATAGTGTTTATT‑CAGCAGCGTTTGCTGC‑ATATCCGGTAACTg                                                                                                 <  2:444271/237‑1 (MQ=255)                                                                                                                                                     cccGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCT‑TGCGGGGTGTCCCGCCAGTGTCGGGTTAAGA‑TAAAACCTGCCTGCGCCCCGCTG‑TAATTCC‑ATCAAAAAACCAGGCTTGATTATACCCTGGTTTCGTGTGATTGGCTGCGGTTTTATACAGCCATTACTGCCCCTAATATGCCTTTGCTCCCTGCTTCAGC‑AGATAGACTTTATT‑CAGCAGCGATCGCTGCCAT‑TCCGGCAACTGCGgcgc                                                                                           >  1:42895/1‑251 (MQ=255)                                                                                                                                                         ttGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCT‑TGCGGGGTGTCCCGCCAGTGTCGGGTTAAGA‑TAAAACCTGCCTGCGCCACGCTG‑TAATTCC‑ATCAAAAAACCAGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGC‑AGATAGTGTTTATT‑CAGCAGCGTTTGCTGC‑ATATCCGGTAACTg                                                                                                 >  1:444271/1‑241 (MQ=255)                                                                                                                                                           gTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCT‑TGCGGGGTGTCCCGCCAGTGTCGGGTTAAGA‑TAAAACCTGCCTGCGCCACGCTG‑TAATTCC‑ATCAAAAAACCAGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGC‑AGATAGTGTTTATT‑CAGCAGCGTTTGCTGC‑ATATCCGGTAACTGCGGCGCTAACTg                                                                                     >  1:376790/1‑251 (MQ=255)                                                                                                                                                            cccGTCGCCAGCCAGAAGGAGACTTCTGTCCCT‑TGCGGGGTGTCCCGCCAGTGTCGGGTTAAGA‑TAAAACCTGCCTGCGCCACGCTG‑TAATTCC‑ATCAAAAAACCAGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGC‑AGATAGTGTTTATT‑CAGCAGCGTTTGCTGC‑ATATCCGGTAACTGCGGCGCTAACTGa                                                                                    <  1:533309/250‑1 (MQ=255)                                                                                                                                                                            aGGAGACTTCTGTCCCT‑TGCGGGGTGTCCCGCCAGTGTCGGGGTAAGA‑TAAAACCTGCCTGCGCCCCGCTG‑TAATTCC‑ATCAAAAAACCAGGCTTGATTATAGCCCAGTTTCGTTTTATTTGCTGTGGTTTTATAAAGGCATTACTGACCGTAATATGCCTTCGCGCCATTCTTTCGCCAG‑TAGTTTTTTTTTCA‑CAGCTGTTGCTGC‑ATATCCGGGTACTTCGGCGCTAACTGACGGCAGaaattccccc                                                                    >  2:458397/1‑243 (MQ=255)                                                                                                                                                                            aGGAGACTTCTCTCCCT‑TGCGCGGTGTCCCCACAGCGTCGGGTTAAGC‑TAAAACCTGCCCTCGCCACGCAG‑CAATTCC‑CTCACCAAACCAGGCTTGATTATAGCCTGGTTTCGCTTGAATGGCTGTGGTTTTATAAAATAATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGC‑ATATAGAGTTTATT‑CAGAAGCGTTTGGTGC‑ATCTCCGGTAACTGCGGCGCCAACAGAAGGCAGAatttcaccc                                                                    >  2:653594/1‑244 (MQ=255)                                                                                                                                                                                       gTCCCT‑TGGGGGGTGACCCGCCAGTGTCGGGTTAAGA‑TAAAACCTGCCTGCGCCACGCTG‑TAATTCC‑ATCAAAAAACCAGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGC‑AGATAGGGTTTATT‑CAGCAGCGTTTGCTGC‑ATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGa                                                           <  1:400670/249‑1 (MQ=255)                                                                                                                                                                                       gTCCCT‑GGGGGGGTGTCCCGCCAGTGGCGGGTTAAGA‑TAAAACCGGCCTGCGCCACGCTG‑TAATTCC‑ATCAAAAAACCAGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGC‑AGATAGTGTTTATT‑CAGCAGCGTTTGCTGC‑ATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGAcc                                                         <  2:114550/251‑1 (MQ=255)                                                                                                                                                                                                gggggTGCCCCGACAGTGTCGGGTTAAGA‑TAAAACCTGCCTGCGCCACGCTG‑TAATTCC‑ATCAAAAAACCAGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGC‑AGATAGTGTTTATT‑CAGCAGCGTTTGCTGC‑ATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAg                                                 <  2:276447/250‑1 (MQ=255)                                                                                                                                                                                                  gggTGTCCCGCCAGTGTCGGGTTAAGA‑TAAAACCTGCCTGCGCCACGCTG‑TAATTCC‑ATCAAAAAACCAGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGC‑AGATAGTGTTTTTT‑CAGCAGCGTTTGCTGC‑ATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGc                                                >  1:476957/1‑250 (MQ=255)                                                                                                                                                                                                     tgtCCCGCCAGTGTCGGGTTAAGA‑TAAAACCTGCCTGCGCCACGCTG‑TAATTCC‑ATCAAAAAACCAGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGC‑AGATAGTGTTTATT‑CAGCAGCGTTTGCTGC‑ATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACg                                             >  1:35960/1‑250 (MQ=255)                                                                                                                                                                                                        cccGCCAGTGTCGGGTTAAGA‑TAAAACCTGCCTGCGCCACGCTG‑AAAATCC‑ATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGC‑AGATAGTGTTTATC‑CAGCAGCGTTTGCTGC‑ATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATgg                                         <  1:39847/251‑1 (MQ=255)                                                                                                                                                                                                        cccGCAAGTGTCGGGTTAAGA‑TAAAACCTGCCTGCGCCACGCTG‑TAATTCC‑ATCAAAAAACCAGGCTTTATTATAGCCTGGTTGCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGC‑AGATAGTGTTTATT‑CAGCAGCGTTTGCTGC‑ATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTGCCAGCACGATgg                                         <  1:572596/251‑1 (MQ=255)                                                                                                                                                                                                         ccGCCAGTGTCGGGTTAAGA‑TAAAACCTGCCTGCGCCACGCTG‑TAATTCC‑ATCAAAAAACCAGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAatat                                                                                                                                                                     >  2:309577/1‑129 (MQ=255)                                                                                                                                                                                                         ccGCCAGTGTCGGGTTAAGA‑TAAAACCTGCCTGCGCCACGCTG‑TAATTCC‑ATCAAAAAACCAGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAatat                                                                                                                                                                     <  1:309577/129‑1 (MQ=255)                                                                                                                                                                                                            cacgtggagggTTACGA‑TAACACCTGCCTGCGCCACGCTG‑TAATTCC‑ATCAAAAAACCAGGCGTGATTATAGACTGTGTTCTTTTGATTGGCTGTGGTTTTATACAGGCATTACTGCCCGTAATATGCCGTCGCGCCATGCTTACGC‑AGATAGTGTTTATT‑CAGCAGCGTTTGCTGC‑ATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCAGATAAGCGTCCTCTTCCAGCGCGATGGCGtt                                     <  2:119869/243‑1 (MQ=255)                                                                                                                                                                                                                tgtCGGGTTAAGA‑TAAAACCTGCCTGCGCCACGCTG‑AAAATCC‑ATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGC‑AGATAGTGTTTATC‑CAGCAGCGTTTGCTGC‑ATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGc                                 >  1:214703/1‑251 (MQ=255)                                                                                                                                                                                                                      gTTAAGA‑TAAAACCTGCCTGCGCCACGCTG‑TAATTCC‑ATCAAAAAACCAGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGC‑AGATAGTGTTTATT‑CAGCAGCGTTTGCTGC‑ATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGc                            >  1:707532/1‑250 (MQ=255)                                                                                                                                                                                                                       ttAAGA‑TAAAACCTGCCTGCGCCACGCTG‑TAATTCC‑ATCAAAAAACCAGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGC‑AGATAGTGTTTATT‑CAGCAGCGTTTGCTGC‑ATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCAt                          >  2:143621/1‑251 (MQ=255)                                                                                                                                                                                                                         aaGA‑TAAACCCTTCCGGCGCCACGCTG‑TAATTCC‑ATAAAAAAACCAGGCGTGATTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGC‑AGATAGTGTTTATT‑CAGCAGCGTTTGCTGC‑ATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCCCCGCATcg                        <  2:642351/251‑2 (MQ=255)                                                                                                                                                                                                                           gA‑TAAAACCTGCCTGCGCCACGCTG‑TAATTCC‑AGCAAAAAACCAGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGC‑AGATAGTGTTTATT‑CAGCAGCGTTTGCTGC‑ATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTc                      <  1:466218/251‑1 (MQ=255)                                                                                                                                                                                                                                  aCCTGCCTGCGCCACGCTG‑TAATTCC‑ATCAAAAAACCAGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGC‑AGATAGTGTTTATT‑CAGCAGCGTTTGCTGC‑ATATCCGGTAACCGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCAtt                >  1:265468/1‑251 (MQ=255)                                                                                                                                                                                                                                    ctgcctgcGCCACGCTG‑TAATTCC‑ATCAAAAAACCAGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGC‑AGATAGTGTTTATT‑CAGCAGCGTTTGCTGC‑ATATCCGGTAACTGCGGCGCTAACTGACGGCAGCATATCCCCATATAAGCGACCTCTTCCAGCACACTGGCGTTATGCACCGCATCTTCGGCTtttt              >  2:147062/1‑251 (MQ=255)                                                                                                                                                                                                                                     tgcctgcGCCACGCTG‑TAATTCC‑ATCAAAAAACCAGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGC‑AGATAGTGTTTATT‑CAGCAGCGTTTGCTGC‑ATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATCTAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCAttttt             >  2:291300/1‑251 (MQ=255)                                                                                                                                                                                                                                                cGCTG‑TAATTCC‑ATCAAAAAACCAGGCTTGATTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGC‑AGATAGTGTTTATT‑CAGCAGCGTTTGCTGC‑ATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCC‑ATGCaa  >  1:342363/1‑250 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                          TACACCGGCTGGCGGTGAAAA‑TCCTTTAACGCCAGCGGTGTCAGGCGAAAGCCTTGTTCACCCTGCAAAATATGCTGAGCGCGGGGAACCTGATCGGCGGGAATAAACGCCACGGACTCTTGCGGTG‑CAAGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCT‑TGCGGGGTGTCCCGCCAGTGTCGGGTTAAGA‑TAAAACCTGCCTGCGCCACGCTG‑AAAATCC‑ATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGC‑AGATAGTGTTTATC‑CAGCAGCGTTTGCTGC‑ATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCCCATATAAGCGACCTCTTCCAGCACGATGGCGTTATGCACCGCATCTTCGGCATTTTTGCCCCATGCAA  >  NC_000913/65545‑66030

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