Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,364,402	G→A	11.9%	R533K (AGA→AAA)	yhcD →	putative outer membrane fimbrial subunit usher protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,364,402	0	G	A	11.9%	166.4 / 2.4	70	R533K (AGA→AAA)	yhcD	putative outer membrane fimbrial subunit usher protein
Reads supporting (aligned to +/- strand):  ref base G (51/8);  new base A (2/7);  total (53/17)
Fisher's exact test for biased strand distribution p-value = 1.87e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

CTTATGCGCTTGGTTTAAACTGGTCTAACTCCATTCTGGGTGGTTTCAACATCAGTGGCTATAAAACATATAGTTACGACGGTGACAATGAC‑TCAAGCAACC‑TT‑AATATTAACT‑GGAACAAAGCG‑TTCAAA‑CACGCCAC‑GGTTTCCG‑TCAA‑CTGGCAGCATCAACTTAGTGCTTCAGAAAATAA‑TGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTT‑GGCAGATCAAAC‑ACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTC  >  NC_000913/3364168‑3364636                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                            cTTATGCGTTTGGTTTAAACTGGTCTAACTCCATTCGGGGTGGTTTCAACATCAGGGGCTATAAAACATATAGTTACGACGGTGACAATGAC‑TCAAGCAACC‑TT‑AATTTTAACT‑GGAACAAAGCG‑TTCAAA‑CACGCCAC‑GGTTTCCG‑TCAA‑CTGGCAGCATCAACTTAGTGCTTCAGAAAATAA‑TGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTT‑GGCAGATCa                                                                                                                                                                                                                                         <  2:779535/239‑1 (MQ=255)            ggTTTAAACTGGTCTAACTCCATTCTGGGTGGTTTCAACATCAGTGGCTATAAAACATATAGTTACGACGGTGACAATGAC‑TCAAGCAACC‑TT‑AATATTAACT‑GGAACAAAGCG‑TTCAAA‑CACGCCAC‑GGTTTCCG‑TCAA‑CTGGCAGCATCAACTTAGTGCTTCAGAAAATAA‑TGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTT‑GGCAGATCAAAC‑ACCGCcac                                                                                                                                                                                                                             >  2:1262859/1‑239 (MQ=255)            ggTTTAAACTGGTCTAACTCCATTCTGGGTGGTTTCAACATCAGTGGCTATAAAACATATAGTTACGACGGTGACAATGAC‑TCAAGCAACC‑TT‑AATATTAACT‑GGAACAAAGCG‑TTCAAA‑CACGCCAC‑GGTTTCCG‑TCAA‑CTGGCAGCATCAACTTAGTGCTTCAGAAAATAA‑TGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTT‑GGCAGATCAAAC‑ACCGCcac                                                                                                                                                                                                                             >  2:1293694/1‑239 (MQ=255)                       gTCTAACTCCATTCTGGGTGGTTTCAACATCAGTGGCTATAAAACATATAGTTACGACGGTGACAATGCC‑TCAAGAAACC‑TT‑AATATTAACT‑GGAACAAAGCG‑TTCAAA‑CACGCCAC‑GGTTTCCG‑TCAA‑CTGGCAGCGTCAACTTAGTGCTTCAGAAACTAA‑TGAATACGCTGGTGATCTGTTCTACGTCAACATCAGTATTCCATTT‑TGCGGGTCAAAC‑AACGCTACACTGTATACt                                                                                                                                                                                                                   >  2:735582/1‑238 (MQ=255)                       gTCTAACTCCATTCTGGGTGGTTTCAACATCAGTGACTATAAAACATATAGTTAAGACGGTGACAATGAC‑TCAAGCAACC‑TT‑AATATTAACT‑GGAACAAAGCG‑TTCAAA‑CACGCCAC‑GGTTTCCG‑TCAA‑CTGGCAGCATCAACTTAGTGCTTCAGAAAATAA‑TGAGTCCGATGGGAATCTGTCCTACGTCAAAATCAGTATTCCATTT‑AGCCGCTCAAAA‑CCCGCCACACTGTATCCt                                                                                                                                                                                                                   >  2:788834/1‑238 (MQ=255)                       gTCTAACTCCATTCTGGGTGGTTTCAACATCAGCGGCTATAAAACATATAGTTACGACGGTGACAATGAC‑TCAAGCAACC‑TT‑AATATTAACT‑GGAACAAAGCG‑TTCAAA‑CACGCCAC‑GGTTTCCG‑TCAA‑CTGGCAGCATCAACTTAGTGATTCAGAAAATAA‑TGAAGACGATGGTGACATGTTCTACGTCAACCTCAGTATTCCATTT‑GGCAGATTAAAA‑ACCGCCACACTGTATAcct                                                                                                                                                                                                                  >  2:347166/1‑237 (MQ=255)                        tCTAACTCCATTCTGGGTGGTTTCAACATCAGTGGCTATAAAACATATAGTTACGACGGTGACAATGAC‑TCAAGCAACC‑TT‑AATATTAACT‑GGAACAAAGCG‑TTCAAA‑CACTCCAC‑GGTTTCCG‑TCAA‑CTGGCAGCATCAACTTATTGCTTCAGAAAATAA‑TGAAGACGAAGGTGATCTGTTCTACGTCAAAACCACTATTCCATTT‑GGCAAATCAAAC‑ACAGCCCCCCTGTAAACTc                                                                                                                                                                                                                  >  1:779618/1‑238 (MQ=255)                        tCTAACTCCATTCTGGGTGGTTTCAACATCAGTGGCTATAAAACATATAGTTACGACGGTGACAATGAC‑TCAAGCAACC‑TT‑AATATTAACT‑GGAACAAAGCG‑ATCAAA‑CACGCCAC‑GGTTTCCG‑TCAA‑CTGGCAGCATCAACTTAGTGCTTCAGAAAATAA‑TGAAGACGATGGTGATCTGTTCTACGTCAACATCAGGATTCCATTT‑GGCAGATCAAAC‑ACCGCCAAACTGTATACTc                                                                                                                                                                                                                  >  2:86753/1‑238 (MQ=255)                                 aTTCTGGGTGGTTTCAACATCAGTGGCTATAAAACATATAGTTACGACGGTGACAATGAC‑TCAAGCAACC‑TT‑AATATTAACT‑GGAACAAAGCG‑TTCAAA‑CACGCCAC‑GGTTTCCG‑TCAA‑CTGGCAGCATCAACTTAGTGCTTCAGAAAATAA‑TGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTCTTCCATTT‑GCCAGATCAAAC‑ACCGCAACACTGGCTACTCGCCATGACg                                                                                                                                                                                                         >  2:112547/1‑238 (MQ=255)                                         attttttttaaaaaaactagttaAAACACTTAGTTAAAGCGGTGAAAAAGTC‑CCAAACAACCATA‑AAAATGAACG‑GGAACAAAGCG‑TAAAAC‑GACGCCAG‑TGTTCCTG‑TCAA‑CGGGAAGCATTAACTTAGTCCTAAAGAAAATAA‑TGAAGAGGATGGGTATTTTTTCTACGCAAACATAATTATTCATTTT‑GGCATAACAAAC‑CCCGCCACACTGTATACTCGCCATGACGATCATaaa                                                                                                                                                                                                 <  1:442512/218‑1 (MQ=255)                                           gtttaaaatcagggtgtaaaaaatCATAAGGTTACACCGTTGAAAATGAC‑AAAACAAACT‑TA‑AATAGCAACT‑GGAAAAAGGCT‑TTAAAA‑CACGCCAC‑GGTTGCAC‑TCAT‑CGGGCAGCATCAACTTAGTGCTTAAGAAAAAAA‑TGAACAGGATGGTGATCTTTTCTACGTAAACTCAAGTATTCCATTT‑GCGAAATCAAAC‑ACCGCCACACTGTATACTCGCCATGACGATCATAAAAcc                                                                                                                                                                                              <  1:880575/219‑1 (MQ=255)                                            tttCAACATCAGTGGCTATAAAACATATAGTTACGCCGGCGACAATGAC‑TCAAGCAACC‑TT‑AATATAAACT‑GGAAAAAAGCG‑ATCAAA‑CACGCCAC‑GGTTTCCG‑TCAA‑CTGACAGCATCAACTTAGTGCCTCAGAAAATAA‑TGAAGACGATG‑‑GACCTGTGCTACGTCAACATCCGAATTCCATCT‑GACAGACCAAAC‑ACCCCCACACTGAAAAAACGCAAAGAAAATAATAAAACCAAc                                                                                                                                                                                           >  1:413967/1‑239 (MQ=255)                                            tttCAACATCAGTGGCTATAAAACATATAGTTACGACGGTGACAATGAC‑TCAAGCAACC‑TT‑AATATTAACT‑GGAACAAAGCG‑TTCAAA‑CACGCCAC‑GGTTTCCG‑TCAA‑CTGGCAGCATCAACTTAGTGCTTCAGAAAATAA‑TGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTT‑GGCAGATCAAAC‑ACCGCCACACTGTATACTCGCCATGACGATCATAAAAccc                                                                                                                                                                                             >  2:1212051/1‑239 (MQ=255)                                            tttCAACATCAGTGGCTATAAAACATATAGTTACGACGGTGACAATGAC‑TCAAGCAACC‑TT‑AATATTAACT‑GGAACAAAGCG‑TTCAAA‑CACGCCAC‑GGTTTCCG‑TCAA‑CTGGCAGCATCAACTTAGTGCTTCAGAAAATAA‑TGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTT‑GGCAGATCAAAC‑ACCGCCACACTGTATACTCGCCATGACGATCATAAAAcc                                                                                                                                                                                              >  2:314170/1‑238 (MQ=255)                                                     gagggcttaaaaaaCAAATAGTTCCGAGGGTGAAAATGCC‑TAAAGCAACC‑TT‑AATATT‑ACT‑GGGAAAAAAGG‑TTCAAACCCCGCCCG‑GTTTTCCG‑TCAA‑CGGGCAGCATCAACTTGGGGCTTAAAAAAATAA‑TGAAGAGGATGGTGATCTGTTCTACGCCAACATCAGTATTCCTTTT‑GGAAGATCAAAC‑ACCGCCACACTGTATACTTGCCATGACGATCATAAAACCCACTATGGTa                                                                                                                                                                                    <  2:886320/229‑1 (MQ=255)                                                      aGTGGCTATAAAACATATAGTTACGACGGTGACAATGAC‑TCAAGCAACC‑TT‑AATATTAACT‑GGAACAAAGCG‑TTCAAA‑CACGCCAC‑GGTTTCCG‑TCAA‑CTGGCAGCATCAACTTAGTGCTTCAGAAAATAA‑TGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTT‑GGCAGATCAAAC‑ACCGCCACACTGTATACTCGCCATGACGATCTTAAAACACACTATGGAAc                                                                                                                                                                                   >  2:381986/1‑239 (MQ=255)                                                        gggctaaaaacacttaggTTACGAGGGTGACAA‑GACATCAAGCACCT‑TA‑AATATAAACG‑GGAAAAAAGGT‑TTAAAA‑CACGCCAG‑GGTTCCCG‑TAAC‑CTGGAAGCATAAACTGAGTGCTTAAGAAAATAA‑GAAAGAGGAGGGTGATTTGTTCTCCGTAAACATCATTATTCCTTTT‑GCAAAATAAAAC‑ACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTAc                                                                                                                                                                                   <  2:361545/220‑1 (MQ=255)                                                         ggCTATAAAAAATAAATTTACGGTGGGAAAAATCCC‑CCAA‑CAACCTTT‑AATA‑TAACT‑GGAAAAAAGGGTTTAAAA‑CCCGCCCG‑GGTTCCCG‑TCAC‑CGGGCAAAATAAACTAAGTGCTTCAAAAAATAA‑TAAAGACTATGGTGATTTTTTCTCCGTCAACATCAGTATTCTTTTT‑GGCAAACAAAAC‑ACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACt                                                                                                                                                                                  <  2:734483/237‑1 (MQ=255)                                                               aaaaCATATAGTTACGACGGTGACAATGAC‑TCAAGCAACC‑TT‑AATATTAACT‑GGAACAAAGCG‑TTCAAA‑CACGCCAC‑GGTTTCCG‑TCAA‑CTGGCAGCATCAACTTAGTGCTTCAGAAAATAA‑TGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTT‑GGCAGATCAAAC‑ACCGCCACACTGTATACTCGCCATGACGATCATCAAACCCACTATGGTACTGGGGTCAt                                                                                                                                                                          >  2:45759/1‑239 (MQ=255)                                                                     tataGTTCCTACGGTGACAATAAC‑TCAAGCAACC‑TT‑AATATTAACT‑GGAACAACGCG‑CTCAAA‑CACGCCAC‑GGTTTCCG‑TCAC‑CTGACAGCATCAACTTAGTGCTTCAGAAACTAA‑TTCAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTT‑GGCAGACGCAAC‑AACGCCACACTGTATACGCGCCATGCCGATCTAACAACACACAATGCAACTgagataagaggaga                                                                                                                                                                    >  1:590673/1‑226 (MQ=255)                                                                     tataGTTACGACGGTGACAATGAC‑TCAAGCAACC‑TT‑AATATTAACT‑GGAACAAAGCG‑TTCAAA‑CACGCCAC‑GGTTTCCG‑TCAA‑CTGGCAGCATCAACTTAGTGCTTCAGAAAATAA‑TGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTT‑GGCAGATCAAAC‑ACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGGGCCATGGGAgt                                                                                                                                                                    >  2:510313/1‑239 (MQ=255)                                                                     tataGTTACAACTGTGACAATGAC‑TCAAGCAACC‑TT‑AATATTAACT‑GGAACAAAGCG‑TTCAAA‑CCCGCCAC‑GGTTTCCG‑TCAA‑CTGGCAGCATCAACTTAGTGCTTCAGAAAATAA‑TGAAGACGATGGTTATCTGTTCTACGTCAACATCAGTATTCTACTT‑GCCAGATCACAC‑ACCCCCACACTTCATACTCGCCATGACAGTCACAAAAACAACTaagataaacggtaaatcggaga                                                                                                                                                                    >  1:512974/1‑218 (MQ=255)                                                                      ataGTTACGACGGTGACAATGAC‑TGAAGCAACC‑TT‑AATATTAAAT‑GGAACAAAGCG‑TTCCAA‑CACGCCAC‑GGTTTCCG‑TCAA‑CTGGCAGCATCAACTTAGTGCTTCAGAAAAAAA‑TGAAGACGATCGTGATCTGTTCTACGTCAACATCAGAATTACATTC‑GGCAGATCTAAA‑AGCTCCACACTATATACCCGCAATGACGAACAAAAAACCACCCATGGAACGGGTAGCACGGGAgta                                                                                                                                                                   >  1:1066662/1‑238 (MQ=255)                                                                      ataGTTACGACGGTGACAATGAC‑TCAAGCAACC‑TT‑AATATTAACT‑GGAACAAAGCG‑TTCAAA‑CACGCCAC‑GGTTTCCG‑TCAA‑CTGGCAGCATCAACTTAGTGCTTCAGAAAATAA‑TGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTT‑GGCAGATCAAAC‑ACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATg                                                                                                                                                                                       >  2:442512/1‑219 (MQ=255)                                                                      ataGTTACGACGGTGACAATGAC‑TCAAGCAACC‑TT‑AATATTAACT‑GGAACAAAGCG‑TTCAAA‑CACGCCAC‑GGTTTCCG‑TCAA‑CTGGCAGCATCAACTTAGTGCTTCAGAAAATAA‑TGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTT‑GGCAGATCAAAC‑ACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAgtc                                                                                                                                                                   >  2:253444/1‑239 (MQ=255)                                                                      ataGTTACGACGGTGACAATGAC‑TCAAGCAACC‑TT‑AATATTAACT‑GGAACAAAGCG‑TTCAAA‑CACGCCAC‑GGTTTCCG‑TCAA‑CTGGCAGCATCAACTTAGTGCTTCAGAAAATAA‑TGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTT‑GGCAGATCAAAC‑ACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAgtc                                                                                                                                                                   >  2:788596/1‑239 (MQ=255)                                                                      ataGTTACGACGGTGACAATGAC‑TCAAGCAACC‑TT‑AATATTAACT‑GGAACAAAGCG‑TTCAAA‑CACGCCAC‑GGTTTCCG‑TCAA‑CTGGCAGCATCAACTTAGTGCTTCAGAAAATAA‑TGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTT‑GGCAGATCAAAC‑ACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGCCATGGGAgtc                                                                                                                                                                   >  2:1037485/1‑239 (MQ=255)                                                                      ataGTTACGACGGTGACAATGAC‑TCAAGCAACC‑TT‑AATATTAACT‑GGAACAAAGCG‑TTCAAA‑CACGCCAC‑GGTTTCCG‑TCAA‑CTGGCAGCATCAACTTAGTGCTTCAGAAAATAA‑TGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTT‑GGAAGATCAAAC‑AACGCCACACTGTATACTCCCCATCACGATCATAAAACACACTATGACACCGGTGGCATGagagaa                                                                                                                                                                   >  1:1194592/1‑237 (MQ=255)                                                                      ataGTTACGACGGTGACAATGAC‑TCAAGCAACC‑TT‑AATATTAACT‑GGAACAAAGCG‑TTCAAA‑CACGCCAC‑GGTTTCCG‑TCAA‑CTGGCAGCATCAACTTAGTGCTTCAGAAAATAA‑TGAAGACGATGGTGATCTGTTCTACGTCAAAATCAGTATTCCATTT‑GGCAAATCAAAC‑ACCGCCACACTGTAAAATCTCCATGACGATCATAAAAACCACAATGGTACAGGAGTCaagagaatc                                                                                                                                                                   >  1:1060492/1‑231 (MQ=255)                                                                        aGTTACGACGGTGACAATGAC‑TCAAGCAACC‑TT‑AATATTAACT‑GGAACAAAGCG‑TTCAAA‑CACGCCAC‑GGTTTCCG‑TCAA‑CTGGCAGCATCAACTTAGTGCTTCAGAAAATAA‑TGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTT‑GGCAGATCAAAC‑ACCGCCACCCTTACTACTCGCCATGCCGATCATAAAACCCCCTCTGGTACTGTGGTCATGGGAgtcc                                                                                                                                                                  >  1:1018377/1‑237 (MQ=255)                                                                        aGTTACGACGGTGACAATGAC‑TCAAGCAACC‑TT‑AATATTAACT‑GGAACAAAGCG‑TTCAAA‑CACGCCAC‑GGTTTCCG‑TCAA‑CTGGCAGCATCAACTTAGTGCTTCAGAAAATAA‑TGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTT‑GGCAGATCAAAA‑ACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGGACTGGTGTCATGGGAgtcg                                                                                                                                                                  >  2:663431/1‑238 (MQ=255)                                                                             cgacgGTGACAATGAC‑TCAAGCAACC‑TT‑AATATTAACT‑GGAACAAAGCG‑TTCAAA‑CACGCCAC‑GGTTTCCG‑TCAA‑CTGGCAGCATCAACTTAGTGCTTCAGAAAATAA‑TGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTT‑GGCAGATCAAAC‑ACCGCCACACTGTATACTCGCCATGACGTTCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCag                                                                                                                                                            >  2:1178404/1‑239 (MQ=255)                                                                             cgacgGTGACAATGAC‑TCAAGCAACC‑TT‑AATATTAACT‑GGAACAAAGCG‑TTCAAA‑CACGCCAC‑GGTTTCCG‑TCAA‑CTGGCAGCATCAACTTAGTGCTTCAGAAAATAA‑TGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTT‑GGCAGATCAAAC‑ACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAgtcgttcc                                                                                                                                                              >  2:5577/1‑234 (MQ=255)                                                                                 gggtaaaaaagcctaaacaaaCT‑TT‑AATATTAACT‑GGAACAAGGCG‑TACAAA‑CTCGCC‑T‑GGTTTCCC‑TTAACCTGGCATCATCAACTGAGTGCTTCAGAAAATAA‑TAAAGCCGATGGTGCTCTGTTTTACGTCAACATCAGTATTCCTTTT‑GGCAGACAAAAC‑ACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCagatg                                                                                                                                                         <  2:356817/220‑1 (MQ=255)                                                                                   gtaCAATCAC‑CCAAGTAACT‑TT‑AATATAAACG‑GAAACAAAGCG‑TTAAAA‑CACGCCAC‑GGTTTCCG‑TAAA‑CTGGCAGCATCAACTTAGTTCTTCAGAAAATAA‑TGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCTTTT‑GGCAGAACAAAC‑ACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCagatgag                                                                                                                                                       <  2:360783/236‑1 (MQ=255)                                                                                   aacaaaTTAC‑TAAAGAATCC‑TT‑AATATTAGCG‑GGAACAAAGCG‑TCCAAC‑CACCCCAC‑‑GTTTCCGTTCAC‑CTGGCAGCACCAATTTAGTCCTACAGAAAATAA‑TGAAGACGATGGTGATTTGTTCTCCGTCAACATCAGTATTCCTTTT‑GGCAAATCAACC‑ACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCagatgag                                                                                                                                                       <  2:1060492/234‑1 (MQ=255)                                                                                    acaatacaC‑TCAA‑CAACC‑TTAAATATAAAC‑‑GGAATAAAGCG‑TTCAAA‑AACGCCTCGGGTTTCCT‑GACA‑CTGCCTG‑ACTAAATTTGTTCATAAAAAAAAAATTGGAGGCGTTGGTGATATGTTCTCGGTACTCCTCAGTATT‑CATTT‑TGCAAGACAAACAACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCagatgagat                                                                                                                                                     <  1:1099918/232‑1 (MQ=255)                                                                                     aCAATGAC‑TCAAGCAACC‑TT‑AATATTAACT‑GGAACAAAGCG‑TTCAAA‑CACGCCAC‑GGTTTCCG‑TCAA‑CTGGCAGCATCAACTTAGTGCTTCAGAAAATAA‑TGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTT‑GGCAGATCAAAC‑ACCGCCACACTGTATAc                                                                                                                                                                                                                    >  2:788679/1‑175 (MQ=255)                                                                                     aCAATGAC‑CAAAGCATCC‑TT‑AATAATAACG‑GAAACAAAGCC‑TTAAAA‑CAGCCCAG‑GGTTTCCG‑TAAA‑CGGGCAAAATCAATTTAGTGCTTCAGAAAAAAA‑TGAAGCCGATGGTATTTTGTTCTACGTAAACATCAGTATTCCATTT‑GGCAGATCAAAC‑ACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCagatgagatg                                                                                                                                                    <  2:1018377/239‑1 (MQ=255)                                                                                          gAC‑TCAAGCAACC‑TT‑AATATTAACT‑GGAACAAAGCG‑TTCAAA‑CACGCCAC‑GGTTTCCG‑TCAA‑CTGGCAGCATCAACTTAGTGCTTCAGAAAATAA‑TGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTT‑GGCAGATCAACA‑ACCGCCACACTGTAAACCCGCCATGACGATCATAAAACCCACAATGGTACTGGTGTCAAGGGCGTCGCTCCAGATAAAATGACct                                                                                                                                                >  1:360783/1‑238 (MQ=255)                                                                                          aaaggacacaaacacatttaatttAACTGGGAACAAAG‑G‑TTAAAACCACGCAAC‑TTTTTCCG‑CAAT‑CTGGCAGCTTCAACTAAGTTCTCCAAAAAATAA‑TGAAGCGGATGGTGTTTTGTTCTCAGCTAACACCAGTATTCTTTTT‑GGAAGATAAATC‑CCCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGt                                                                                                                                                   <  2:1004684/217‑1 (MQ=255)                                                                                           aC‑TCAAGCAACC‑TT‑AATATTAACT‑GGAACAAAGCG‑TTCAAA‑CACGCCAC‑GGTTTCCG‑TCAA‑CTGGCAGCATCAACTTAGTGCTTCAGAAAATAA‑TGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTT‑GGCAGATCAAAC‑ACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCctac                                                                                                                                              >  2:518698/1‑239 (MQ=255)                                                                                                aaGCAACC‑TT‑AATATTAACT‑GGAACAAAGCG‑TTCAAA‑CACGCCAC‑GGTTTCCG‑TCAA‑CTGGCAGCATCAACTTAGTGCTTCAGAAAATAA‑TGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTT‑GGCAGATCAAAC‑ACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATTGGAGCCGTCTCAGATGGGAGTTCCTACTAtg                                                                                                                                          >  2:419141/1‑239 (MQ=255)                                                                                                      cc‑TT‑AATATTAACT‑GGAACAAAGCG‑TTCAAA‑CACGCCAC‑GGTTTCCG‑TCAA‑CTGGCAGCATCAACTTAGTGCTTCAGAAAATAA‑TGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTT‑GGCAGATCAAAC‑ACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAAt                                                                                                                                     >  2:878858/1‑238 (MQ=255)                                                                                                               tagaacgggaaaaaaaGCT‑TTCAAA‑CCCGCCCC‑GGTTCCGG‑TAAA‑CTGGAAGCATCAACTTAGTGCTTCAAAAAATAA‑GAAAGCGGTGGGTGATCTTTTCTACGTCAACATCATTATTCCATTT‑GCCAGAACAAAC‑ACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAAc                                                                                                                              <  2:868447/226‑1 (MQ=255)                                                                                                                   aCT‑GGAACAAAGCG‑TTCAAA‑CACGCCAC‑GGTTTCCG‑TCAA‑CTGGCAGCATCAACTTAGTGCTTCAGAAAATAA‑TGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTT‑GGCAGAt                                                                                                                                                                                                                                           >  2:510050/1‑126 (MQ=255)                                                                                                                   aCT‑GGAACAAAGCG‑TTCAAA‑CACGCCAC‑GGTTTCCG‑TCAA‑CTGGCAGCATCAACTTAGTGCTTCAGAAAATAA‑TGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTT‑GGCAGATCAAAC‑ACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGGGAt                                                                                                                         >  2:772537/1‑239 (MQ=255)                                                                                                                         aaCAAAGCG‑TTCAAA‑CACGCCAC‑GGTTTCCG‑TCAA‑CTGGCAGCATCAACTTAGTGCTTCAGAAAATAA‑TGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTT‑GGCAGATCAAAC‑ACCGCCACACTGTATAc                                                                                                                                                                                                                    >  2:659291/1‑143 (MQ=255)                                                                                                                         aaCAAAGCG‑TTCAAA‑CACGCCAC‑GGTTTCCG‑TCAA‑CTGGCAGCATCAACTTAGTGCTTCAGAAAATAA‑TGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTT‑GGCAGATCAAAC‑ACCGCCACACTGTATACTCGCCATGACGATCATAAAACCAACTATGGTACTGGTGTCAt                                                                                                                                                                          >  2:776001/1‑185 (MQ=255)                                                                                                                                                   ggtttcccgcaaa‑CGGGCAGTATAAACTTAGTGTTTCAAAAAATAA‑TGAAGACGATGGTGATCTTTTATAGCTCAACATCAGTATTCCTTTG‑GACATATAAAAC‑ACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAAc                                                                                              <  2:896430/228‑1 (MQ=255)                                                                                                                                                                    gcagcaTCAACTTAGTGCTTCAGAAAATAA‑TGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTT‑GGCAGATCAAAC‑ACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGCCCTCCTATGTGAATCTTGACCGCGATAACTAAAAACGTGAAACGGGCTTGAACggacgcaacatatcaa                                                                              >  1:1113357/1‑225 (MQ=255)                                                                                                                                                                       gcaTCATCTTAGGTCGTAAGAAAAAAA‑AGAAGACGATGGGTCTCTTTTCTACGGCAACAATAGTTTTCGTTTT‑GGCAGACAAAAC‑ACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAAt                                                                            <  2:612130/238‑1 (MQ=255)                                                                                                                                                                              cTTAGTGCTTCAGAAAATAA‑TGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTT‑GGCAGATCAAAC‑ACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCacg                                                                                                                     >  2:543099/1‑190 (MQ=255)                                                                                                                                                                                        cAGAAAATAA‑TGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTT‑GGCAGATCAAAC‑ACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAGACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTc                                                          >  2:438632/1‑239 (MQ=255)                                                                                                                                                                                        cAGAAAATAA‑TGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTT‑GGCAGATCAAAC‑ACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTc                                                          >  2:398911/1‑239 (MQ=255)                                                                                                                                                                                        cAGAAAATAA‑TGAAGACGAATGTGATCTGTTCTACGTCAACATCAGTATTCCATTT‑GGCAGATCAAAC‑ACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGACATGGGAGTAGTCCCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACAAGCTTGAACGGCCGAATCCGTTCCAATCTCCATTACTCCCACGTc                                                          >  2:922357/1‑239 (MQ=255)                                                                                                                                                                                                        aCGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTT‑GGCAGATCAAAC‑ACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTAACAGTTCCAATCTCCATTAAACCCacagcgcccttcccccagaa                                           >  1:1134058/1‑220 (MQ=255)                                                                                                                                                                                                           aTGGTGATCTGTTCTACGTCAACATCAGTATTCCATTT‑GGCAGATCAAAC‑ACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCa                                        >  2:276625/1‑239 (MQ=255)                                                                                                                                                                                                                 tttttttCAACGACAACATCA‑TTTTCCATTTGGGAAAAACAAAC‑CCCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCgg                                   <  2:1113357/233‑1 (MQ=255)                                                                                                                                                                                                                      ttCTACGTCAACATCAGTATTCCATTT‑GGCAGATCAAAC‑ACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGCGATCACGACGAACGAGAAACGAGCTTGAACGGCAGTATCAGTTCCAAACTCCATTACACCCAAAGCAGCCATTGCGCAGGCGCAACCGGCAggca                             >  1:756272/1‑236 (MQ=255)                                                                                                                                                                                                                        cTACGTCAACATCAGTATTCCATTT‑GGCAGATCAAAC‑ACCGCCACACTGTATACTCGCCATGACGTTCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGACCGTGAACCGAGCTTGAACGTCAGGTTCAGTTCCAATCTCCATTACACCCAATTCAGCCTTGCCCCAAGAGCAAACGGCAGGGATa                           >  1:728909/1‑239 (MQ=255)                                                                                                                                                                                                                             taaaacgaattttccctttt‑GGCAAACCAAAC‑ACCGGCACCCTGTATCCTCGGCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCATTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGt                      <  2:512974/222‑1 (MQ=255)                                                                                                                                                                                                                                aCATCAGTATTCCATTT‑GGCAGATCAAAC‑ACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACt                   >  2:251317/1‑239 (MQ=255)                                                                                                                                                                                                                                aCATCAGTATTCCATTT‑GGCAGATCAAAC‑ACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGACGTGTTAACCGGACt                   >  2:1226244/1‑239 (MQ=255)                                                                                                                                                                                                                                aCATCAGTATTCCATTT‑GGCAGATCAAAC‑ACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAACCTTGCCCCAGGAGCAAACggacggggaagccgcacc                   >  1:515911/1‑223 (MQ=255)                                                                                                                                                                                                                                     aGTATTCCATTT‑GGCAGATCAAAC‑ACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACa               >  2:50259/1‑238 (MQ=255)                                                                                                                                                                                                                                           ccATTT‑GGCAGATCAAAC‑ACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACggacc        >  2:535711/1‑236 (MQ=255)                                                                                                                                                                                                                                              ttt‑GGCAGATCAAAC‑ACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGACTTGCCGCAGCAGCAAGCGGCACGGATAGCCGAACTTACAACgggaggag     >  1:1290112/1‑233 (MQ=255)                                                                                                                                                                                                                                              ttt‑GGCAGATCAAAC‑ACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGGATCAGTTCCAATCTCCAGTACACCCAAGTCAGACTTGCCGCATGAGCAAACGGCAGTGGTAGCCGTAACTAAAACGGTAAGac     >  1:854058/1‑238 (MQ=255)                                                                                                                                                                                                                                                  ggCAGATCAAAC‑ACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCCTTACACCCAAGTCAGCCTTGCCGCAGGAGCCAACGGCAGAGATAACCGTAATTAACACAGTACGATGTc  >  1:669415/1‑239 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                            CTTATGCGCTTGGTTTAAACTGGTCTAACTCCATTCTGGGTGGTTTCAACATCAGTGGCTATAAAACATATAGTTACGACGGTGACAATGAC‑TCAAGCAACC‑TT‑AATATTAACT‑GGAACAAAGCG‑TTCAAA‑CACGCCAC‑GGTTTCCG‑TCAA‑CTGGCAGCATCAACTTAGTGCTTCAGAAAATAA‑TGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTT‑GGCAGATCAAAC‑ACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGATAGCCGTACTTACAACGGTACGATGTC  >  NC_000913/3364168‑3364636

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