Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,665,701	C→A	36.1%	A37S (GCG→TCG)	gadX ←	acid resistance regulon transcriptional activator, autoactivator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,665,701	0	C	A	36.1%	40.9 / 19.5	57	A37S (GCG→TCG)	gadX	acid resistance regulon transcriptional activator, autoactivator
Reads supporting (aligned to +/- strand):  ref base C (33/3);  new base A (2/19);  total (35/22)
Fisher's exact test for biased strand distribution p-value = 4.34e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.45e-01

AGCATTTCTGCATAGCGAAGATTCTTATTCTGCGATAGTTGCGCAACTTCCTGAAAAACAGG‑AATATCTCTACTACAGTGATGAACGAGTAATGAAGAAACTTTTTT‑ATGTG‑CATGAAGATTTAATGCCTCCTCCT‑TG‑AG‑CATCGGGAGAAATAATGAAAGCGTGTCCC‑TTGACACG‑AATACAAAATTT‑TC‑AAC‑A‑CACT‑TATCTACTCGAATTTGGCTTGCAT‑CCGCAAAAACCAGG‑TCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAACATATTA  >  NC_000913/3665475‑3665838                                                                                                                                                                                                                                                 |                                                                                                                                           aGCATTTCTGCATAGCGAAGATTCTTATTCTGCGATAGTTGCGCAACTTCCTGAAAAACAGG‑AATATCTCTACTACAGTGATGAACGAGTAATGAAGAAACTTTTTT‑ATGTG‑CATGAAGATTTAATGCCTCCTCCT‑TG‑AG‑CATCGGGAGAAATAATGAAAGCGTGTCCC‑TTGACACG‑AATACAAAATTT‑TC‑AAC‑A‑CACT‑AACATACTCGAATTTGGCTTGAAT‑CCGCa                                                                                                                                           >  2:710204/1‑228 (MQ=255) aGCATTTCTGCATAGCGAAGATTCTTATTCTGCGATAGTTGCGCAACTTCCTGAAAAACAGG‑AATATCTCTACTACAGTGATGAACGAGTAATGAAGAAACTTTTTT‑ATGTG‑CATGAAGATTTAATGCCTCCTCCT‑TG‑AG‑CATCGGGAGAAATAATGAAAGCGTGTCCC‑TTGACACG‑AATACAAAATTT‑TC‑AAA‑A‑CACT‑TATCTACTCGAATTTGGCTTGCAT‑CCGCac                                                                                                                                          >  1:464856/1‑228 (MQ=255)   cATTTCTGCATAGCGAAGATTCTTATTCTGCGATAGTTGCGCAACTTCCTGAAAAACAGG‑AATATCTCTACTACAGTGATGAACGAGTAATGAAGAAACTTTTTT‑ATGTG‑CATGAAGATTTAATGCCTCCTCCT‑TG‑AG‑CATCGGGAGAAATAATGAAAGCTTGTCCC‑TTGCCACG‑AATACAAAATTT‑TC‑AAA‑A‑AACT‑ACTCACTTCGAATTTGGCTTGCAT‑CCCCaaa                                                                                                                                         >  2:23026/1‑228 (MQ=255)   cATTTCTGCATAGCGAAGATTCTTATTCTGCGATAGTTGCGCAACTTCCTGAAAAACAGG‑AATATCTCTACTACAGTGATGAACGAGTAATGAAGAAACTTTTTT‑ATGTG‑CATGAAGATTTAATGCCTCCTCCT‑TG‑AG‑CATCGGGAGAAATAATGAAAGCGTTTCCC‑TTGACACG‑AATACAAAATTT‑TC‑AAC‑A‑AACT‑TAACAACTCAAATTTGGCTTGCAC‑CCGAaaaa                                                                                                                                        >  2:350564/1‑229 (MQ=255)   cATTTCTGCATAGCGAAGATTCTTATTCTGCGATAGTTGCGCAACTTCCTGAAAAACAGG‑AATATCTCTACTACAGTGATGAACGAGTAATGAAGAAACTTTTTT‑ATGTG‑CATGAAGATTTAATGCCTCCTCCT‑TG‑AG‑CATCGGGAGAAATAATGAAAGCGTGTCCC‑TTGACACG‑AATACAAAATTT‑TC‑AAC‑A‑CACT‑TATCTACTCGAATTTGGCTTGCAT‑CCGCaaa                                                                                                                                         >  1:714788/1‑228 (MQ=255)   cATTTCTGCATAGCGAAGATTCTTATTCTGCGATAGTTGCGCAACTTCCTGAAAAACAGG‑AATATCTCTACTACAGTGATGAACGAGTAATGAAGAAACTTTTTT‑ATGTG‑CATGAAGATTTAATGCCTCCTCCT‑TG‑AG‑CATCGGGAGAAATAATGAAAGCGTGTCCC‑TTGACACG‑AATACAAAATTT‑TC‑AAC‑A‑CACT‑TATCTACTCGAATTTGGCTTGCAT‑CCGCaaa                                                                                                                                         >  1:100113/1‑228 (MQ=255)   cATTTCTGCATAGCGAAGATTCTTATTCTGCGATAGTTGCGCAACTTCCTGAAAAACAGG‑AATATCTCTACTACAGTGATGAACGAGTAATGAAGAAACTTTTTT‑ATGTG‑CATGAAGATTTAATGCCTCCTCCT‑TG‑AG‑CATCGGGAGAAATAATGAAAGCGTGTCCC‑TTGACACG‑AATACAAAATTT‑TA‑AAC‑A‑CACT‑TATCTACTCGAATTTGGCTTGCAT‑CCGCaaa                                                                                                                                         >  2:141627/1‑228 (MQ=255)   cATTTCTGCATAGCGAAGATTCTTATTCTGCGATAGTTGCGCAACTTCCTGAAAAACAGG‑AATATCTCTACTACAGTGATGAACGAGTAATGAAGAAACTTTTTT‑ATGTG‑CATGAAGATTTAATGCCTCCTACT‑TG‑AG‑CATCGGGAGACATACTTAAAGCTGGTCCC‑TTGACACG‑AATACAAAATTT‑TC‑AAA‑A‑CACT‑TATCTACTCGAATTTGGCTTGAAT‑CCGCaaa                                                                                                                                         >  2:476190/1‑228 (MQ=255)                cGAAGATTCTTATTCTGCGATAGTTGCGCAACTTCCTGAAAAACAGG‑AATATCTCTACTACAGTGATGAACGAGTAATGAAGAAACTTTTTT‑ATGTG‑CATGACGATTTAATGCATCCTCCT‑TG‑AG‑CATCGGGAGAAATAATGAAAGCGTGTCCC‑TTGACACG‑AATACAAAATTT‑TC‑AAC‑A‑CACT‑TATCTACTCGAATTTGGCTTGCAC‑ACGCAAAAACCAGG‑TCACCGc                                                                                                                          >  1:62175/1‑229 (MQ=255)                cGAAGATTCTTATTCTGCGATAGTTGCGCAACTTCCTGAAAAACAGG‑AATATCTCTACTACAGTGATGAACGAGTAATGAAGAAACTTTTTT‑ATGTG‑CATGAAGATTTAATGCCTCCTCCT‑TG‑AG‑CATCGGGAGAAATAATGAAAGCGTGTCCC‑TTGACCCG‑CATACAAAATTT‑TC‑AAA‑A‑AACT‑TATCTACTCGAATTTGGCTTGCAT‑CCCAAAAAAcaaaataaacac                                                                                                                           >  2:84597/1‑218 (MQ=255)                cGAAGATTCTTATTCTGCGATAGTTGCGCAACTTCCTGAAAAACAGG‑AATATCTCTACTACAGTGATGAACGAGTAATGAAGAAACTTTTTT‑ATGTG‑CATGAAGATTTAATGCCTCCTCCT‑TG‑AG‑CATCGGGAGAAATAATGAAAGCGTGTCCC‑TTGACACG‑AATCCAAAATTT‑TC‑AAC‑A‑CACT‑TATCTACTCGAATTTGACTTGCAT‑CCGCAAAAACAAGG‑TCACCGc                                                                                                                          >  1:706218/1‑229 (MQ=255)                cGAAGATTCTTATTCTGCGATAGTTGCGCAACTTCCTGAAAAACAGG‑AATATCTCTACTACAGTGATGAACGAGTAATGAAGAAACTTTTTT‑ATGTG‑CATGAAGATTTAATGCCTCCTCCT‑TG‑AC‑CATCGGGAGAAATAATGAAAGCGTGTCCC‑TTGACACG‑AATACAAAATTT‑TC‑AAC‑A‑CACT‑TATCTACTCCATTTTGGCTTGCAT‑CCTCAAAAACCAGG‑ACACCGc                                                                                                                          >  1:348372/1‑229 (MQ=255)                     aTTCTTATTCTGCGATAGTTGCGCAACTTCCTGAAAAACAGG‑AATATCTCTACTACAGTGATGAACGAGTAATGAAGAAACTTTTTT‑ATGTG‑CATGAAGATTTAATGCCTCCTCCT‑TG‑AG‑CATCGGGAGAAATAATGAAAGCGTGTCCC‑TTGACACG‑AATACAAAATTT‑TC‑AAC‑A‑AACT‑TATCTACTCGAATTTGGCTTGCAT‑CCGCAAAAACCAGG‑TCACCCCCATTa                                                                                                                     >  2:721307/1‑229 (MQ=255)                           aTTCTGCGATAGTTGCGCAACTTCCTGAAAAACAGG‑AATATCTCTACTACAGTGATGAACGAGTAATGAAGAAACTTTTTT‑ATGTG‑CATGAAGATTTAATGCCTCCTCCT‑TG‑AG‑CATCGGGAGAAATAATGAAAGCGTGTCCC‑TTGACACG‑AATACAAAATTT‑TC‑AAC‑A‑CACT‑TATCTACTCGAATTTGGCTTGCAT‑CCGCAAAAACCAGG‑TCACCCCAATTAAAAAAg                                                                                                               >  2:299879/1‑229 (MQ=255)                                 cGATCGTTGCGCAACTTCCTGAAAAACAGG‑AATATCTCTACTACAGTGATGAACGAGTAATGAAGAAACTTTTTT‑ATGTG‑CATGAAGATTTAATGCCTCCTCCT‑TG‑AG‑CATCGGTAGAAATAATGAAAGCGTGTCCC‑TTGACACG‑AATACAAAATTT‑TC‑AAC‑A‑CACT‑TATCTACTCGAATTTGGCTTGCAT‑CCGCAAAAACCAGG‑TCACCGCCATTAAAATACCGatat                                                                                                         >  1:242491/1‑229 (MQ=255)                                 cGATAGTTGCGCAACTTCCTGATAAACAGG‑AATATCTCTACTACAGTGATGAACGAGTAATGAAGAAACTTTTTT‑ATGTG‑CATGAAGATTTAATGCCTCCTCCT‑TG‑AG‑CATCGGGAGAAATAATGAAAGCGTGTCCC‑TTGACCCG‑AATACAAAATTT‑TC‑AAC‑A‑CACT‑TATATACCCGAATTTGGTCTCATT‑CCGCAAAAACCAGG‑TCACCACCATTAAAAAAGcaaaa                                                                                                          >  2:340144/1‑224 (MQ=255)                                 cGATAGTTGCGCAACTTCCTGAAAAACAGG‑AATATCTCTACTACAGTGATGAACGAGTAATGAAGAAACTTTTTT‑ATGTG‑CATGAAGATTTAATGCCTCCTCCT‑TG‑AG‑CATCGGGAGAAATAATGAAAGCGTGTCCC‑TTGACACG‑AATACAAAATTT‑TC‑AAC‑A‑CACT‑TATCTACTCGAATTTGGCTTGCAT‑CCGCAAAAACCAGG‑TCACCGCCATTAAAATAGCGata                                                                                                          >  1:282403/1‑228 (MQ=255)                                 cGATAGTTGCGCAACTTCCTGAAAAACAGG‑AATATCTCTACTACAGTGATGAACGAGTAATGAAGAAACTTTTTT‑ATGTG‑CATGAAGATTTAATGCCTCCTCCT‑TG‑AG‑CATCGGGAGAAATAATGAAAGCGTGTCCC‑TTGACACG‑AATACAAAATTT‑TC‑AAC‑A‑CACT‑TATCTACTCGAATTTGGCTTGCAT‑CCGCAAAAACCAGG‑TACCCCCCCTTAAAAAAGCggaat                                                                                                         >  2:307521/1‑225 (MQ=255)                                 cGATAGTTGCGCAACTTCCTGAAAAACAGG‑AATATCTCTACTACAGTGATGAACGAGTAATGAAGAAACTTTTTT‑ATGTG‑CATGAAGATTTAATGCCTCCTCCT‑TG‑AG‑CATCGGGAGAAATAATGAAAGCGTGTCCC‑TTGACACG‑AATACAAAATTT‑TC‑AAC‑A‑CACT‑TATCTACTCGAATTTGGCTTGCAT‑CCGCAAAAACCAGG‑TAACCGCCATTAAAATAACGata                                                                                                          >  2:474987/1‑228 (MQ=255)                                  aaTAGTGGCGAAATTTCCTGAAAAACAGA‑TATTTCTCTACTAAAGGGAAGAAGGGGTAATGAAGAAACTTTTTT‑ATGGG‑CATGATGTTTAAATGCTTCCCCCT‑TA‑AG‑CTTAGGGAAAAATAATGAAAGCGTGCCCC‑TTGACTCA‑AATACAAAATTA‑AC‑ACC‑A‑CCCT‑TCTCTCCTCAAATTTGGCTTGCCC‑CCGAAAAAACCAGG‑TCACCGCCATTAAAATAGCGATAtt                                                                                                        <  1:466076/228‑1 (MQ=255)                                                ctcacaaaaaaCAGGAAATCTCTTTACTACATGTATGAACGAGTAAGGAAAAACTTTTTTT‑AGGGG‑AATAAAGATTTAATCCCTCCCCCT‑TG‑AC‑AATGGGGAAAAATAAGAAAAGCTTGCCCC‑TTGACACA‑AATAAAAATTTT‑TC‑AAC‑C‑CTTT‑TATCTACTCAAATTTGGCTGGCAC‑CCAAAAAAACCAGT‑TCACCGCCATTAAAATAGCGATATTCACCATTAACCAt                                                                                           <  1:721308/223‑1 (MQ=255)                                                  cataaaaaaaAGG‑AATACCTCTAATACAGTGAGGAACAAAAAAGGAAAAAATTTTTTA‑ATGTG‑CATGAAGATTGAATCCCTCCCCCT‑TA‑AG‑CACCGGGAAAAAAAATGAAGGCGTGTCCC‑TTGACAC‑‑AAAAAAAAATTT‑TCAAAC‑A‑CATT‑TATCTACTCAATTTGGGCTTCCAC‑CCAAAAAAACAAGT‑TCACCGCCATTAAAATAGCGATATTCACCATTAACCATgg                                                                                         <  1:368925/224‑1 (MQ=255)                                                              g‑AATATCTCTACTACAGTGATGAACGAGTAATGAAGAAACTTTTTT‑ATGTG‑CATGAAGATTTAATGCCTCCTCCT‑TG‑AG‑CATCGGGAGAAATAATGAAAGCGTGTCCC‑TTGACACG‑AATACAAAATTT‑TC‑AAC‑A‑CACT‑TATCTACTCGAATTTGGCTTGAAT‑CCGCAAAAACCAGG‑TCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGGGAGAATAAAttt                                                                            >  1:541373/1‑229 (MQ=255)                                                                   ttctaTACTACAGGGATAAACGAGTAATAAAAAATTTTTTTT‑ATGGG‑AAGAAAGTTTGAATGCCCCCCCCT‑TG‑GC‑CATCGAGAAAAAAAATAAAGGGTTGTCCC‑TTGAACCA‑AATAAAAAATTT‑TC‑ACC‑A‑CTCT‑TATCTCCTCGATTTTGGTTTGCAC‑CCGAAAAAACCAGG‑TCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGt                                                                        <  1:299879/228‑1 (MQ=255)                                                                          ctacAGTGATGAACGAGTAATGAAGAAACTTTTTT‑ATGTG‑CATGAAGATTTAATGCCTCCTCCT‑TG‑AG‑CATCGGGAGAAATAATGAAAGCGTGTCCC‑TTGACACG‑AATACAAAATTT‑TC‑AAC‑A‑CACT‑TATCTACTCGAATTTGGCTTGCAT‑CCGCAAAAACCAGG‑TCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCAt                                                                 >  1:379602/1‑229 (MQ=255)                                                                           taaattgaagaactataaaagaagaaCTTTTTTT‑AAGTG‑AATAAAGATTAAACGCCTCATCTT‑TA‑CG‑CACCGGGAAAAATAAGAAAAGCGTGTCCC‑TTGACAGA‑AATAAAAAATTG‑TC‑AAC‑A‑CTTT‑TATCTCCGCAATTTTGGCTTCCAT‑CCACAAAACCCAGC‑TCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATa                                                                <  1:23026/209‑1 (MQ=255)                                                                           taaaggagaaaaggAGTAAAGAAAAAGTTTTTAT‑AAGTG‑CATGAAGATTTAATGCCTCCTCGT‑TG‑AG‑CATCGGAAAAAAAAATGAAGGCGT‑TTCC‑TGAAAACGAAATAAAAATTTT‑TC‑AAC‑A‑CTCT‑TATCTACTCGATTTTGGCTTGCAT‑CCAAAAAAACGAGG‑TCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATa                                                                <  1:506942/216‑1 (MQ=255)                                                                             cAGTGATGAACGAGTAATGAAGAAACTTTTTT‑ATGTG‑CATGAAGATTTAATGCCTCCTCCT‑TG‑AG‑CATCGGGAGAAATAATGAAAGCGTGTCCC‑TTGACACG‑AATACAAAATTT‑TC‑AAC‑A‑CACT‑TATCTACTCGAATTTGGCTTGCAT‑CCGCAAAAACCAGG‑TCACCGCCATTAAAATAGCGATATTCCCCATTAACCATGGGGAGAAAATATATAAGTCTTGCaaaac                                                              >  2:715340/1‑225 (MQ=255)                                                                                tgatgaACGAGTAATGAAGAAACTTTTTT‑ATGTG‑CATGAAGATTTAATGCCTCCTCCT‑TG‑AG‑CATCGGGAGAAATAATGAAAGCGTGTCCC‑TTGACACG‑AATACAAAATTT‑TC‑AAC‑A‑CACT‑TATCTACTCGAATTTGGCTTGCAT‑CCGCAAAAAc                                                                                                                                      >  1:317847/1‑155 (MQ=255)                                                                                tgatgaACGAGTAATGAAGAAACTTTTTT‑ATGTG‑CATGAAGATTTAATGCCTCCTCCT‑TG‑AG‑CATCGGGAGAAATAATGAAAGCGTGTCCC‑TTGACACG‑AATACAAAATTT‑TC‑AAC‑A‑CACT‑TATCTACTCGAATTTGGCTTGCAT‑CCGCAAAAACCAGG‑TCACCGCCATTAAAATAGCGATAtt                                                                                                        >  1:288907/1‑184 (MQ=255)                                                                                tgatgaACGAGTAATGAAGAAACTTTTTT‑ATGTG‑CATGAAGATTTAATGCCTCCTCCT‑TG‑AG‑CATCGGGAGAAATAATGAAAGCGTGTCCC‑TTGACACG‑AATACAAAATTT‑TC‑AAC‑A‑CACT‑TATCTACTCGAATTTGGCTTGCAT‑CCGCAAAAACCAGG‑TCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCa                                                            >  1:214380/1‑228 (MQ=255)                                                                                tgatgaACGAGTAATGAAGAAACTTTTTT‑ATGTG‑CATGAAGATTTAATGCCTCCTCCT‑TG‑AG‑CATCGGGAGAAATAATGAAAGCGTGTCCC‑TTGACACG‑AATACAAAATTT‑TC‑AAC‑A‑CACT‑TATCTACTCGAATTTGGCTTGCAT‑CCGCAAAAACCAGG‑TCACCGCCATTAAAAAAGCGATATTCACCATTAACCATGGTGAAAATAAATTTATTTTTTGCAAACGCa                                                            >  2:709772/1‑228 (MQ=255)                                                                                  atgaACGAGTAATGAAGAAACTTTTTT‑ATGTG‑CATGAAGATTTAATGCCTCCTCCT‑TG‑AG‑CATCGGGAGAAATAATGAAAGCGTGTCCC‑TTGACACG‑AATACAAAATTT‑TC‑AAC‑A‑CACT‑TATCTACTCGAATTTGGCTTGCAT‑CCGCAAAAACCAGG‑TCACCGCCATTAAAATAg                                                                                                               >  1:429384/1‑175 (MQ=255)                                                                                  atgaACGAGTAATGAAGAAACTTTTTT‑ATGTG‑CATGAAGATTTAATGCCTCCTCCT‑TG‑AG‑CATCGGGAGAAATAATGAAAGCGTGTCCC‑TTGACACG‑AATACAAAATTT‑TC‑AAC‑A‑CACT‑TATCTACTCGAATTTGGCTTGCAT‑CCGCAAAAACCAGG‑TCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTAGGAATATATTAATGTCTTGCATACGCAAt                                                          >  1:102645/1‑228 (MQ=255)                                                                                  atgaACGAGTAATGAAGAAACTTTTTT‑ATGTG‑CATGAAGATTTAATGCCTCCTCCT‑TG‑AG‑CATCGGGAGAAATAATGAAAGCGAGTCCC‑TTGACACG‑AATACAAAATTT‑TC‑AAC‑A‑CACT‑TATCTACTCGAATTTGGCTTGCAT‑CCGCAAAAACCAGG‑TACCCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTATTGCAAACGCACtt                                                         >  2:379579/1‑229 (MQ=255)                                                                                     aaagaaaaataaacaatttttttt‑ATGGT‑CATAAATTTTTAATGCC‑CCTCCG‑TG‑AGCCAGTGGGAAAAAAAATAAAAGCTCGTCTC‑TTGACACA‑AAAACAAATTTT‑CC‑AAC‑‑‑CTCTGTACCTACAAGATTTTGGCTGACAC‑CCAAAAAACCAAGC‑TAACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGa                                                      <  2:430663/211‑1 (MQ=255)                                                                                       cgggTAAAGAAAAAACTTTTTT‑ATTTGACAGAAAGTTTTATGGCCTCCCCAT‑TA‑AA‑CATCGGAAAAAATAATGAAAGCGTGTCCC‑TTGACACA‑AAAACAAAATTT‑TC‑ACC‑C‑CATT‑TATCTACTCAAATTGGGCTTGCAT‑CACAAAAACCCGGG‑TACCCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGAc                                                     <  1:709773/226‑1 (MQ=255)                                                                                        gttatattgaaattttttttt‑ATGGG‑AAAAAATTTTTAACGCCTCCCTCT‑TG‑AG‑ACAGGGAAAAAATAAAAAAGGCGTGTCCC‑TTGAAACA‑AATA‑AAAATTT‑TCAAAC‑A‑CATT‑TCTTTACTTGATTTTGGTTTGCAT‑CCACAAAAACGAGG‑TCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACaa                                                   <  2:299309/214‑1 (MQ=255)                                                                                        gaaaaaagaaaaaATTTTTTT‑ATGTG‑CAAAAAGATTTAACGCCCCCTCTT‑GG‑AG‑ACGCGGAAAAAATAATGAAGGGGTGCCCC‑TTGACACA‑AATACAAAATTT‑AC‑AAC‑A‑CTTT‑TATCTACTAGATTTTGGCTTCCAC‑CCGAAAAAACCAGG‑TCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACaa                                                   <  1:487594/222‑1 (MQ=255)                                                                                        aattaataaaaaatttgtttt‑AGGTG‑AAAGAAAATTTACGCCCTCCTCCT‑GG‑AC‑CACGGGGAAAAAAAATGAATCGTTGTCCA‑TGGAAAGA‑AAAACAAAATTT‑CA‑AAC‑C‑CTCT‑TATCTACTCAAATTTGTCTGGCAC‑CCAAAAAAACCAGG‑TAACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGCCaa                                                   <  1:476190/212‑1 (MQ=255)                                                                                         aattaTGAAAAAAATTTTTTCATGGG‑AATGAAGTTTGAATCCCCCCTGCT‑GG‑AA‑CATCGGGAAAAAAAATGAAAGCTTCTCCT‑TAAACACG‑AATACAAAATTT‑AC‑AAC‑A‑CTAT‑TATTACCTGGAATTGGGTTTGCAC‑CCAAAAAACCCAGA‑TATCCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACaa                                                   <  1:340144/225‑1 (MQ=255)                                                                                          gTAATGAAGAAACTTTTTT‑ATGTG‑CATGAAGATTTAATGCCTCCTCCT‑TG‑AG‑CATCGGGAGAAATAATGAAAGCGTGTCCC‑TTGACACG‑AATACAAAATTT‑TC‑AAC‑A‑CACT‑TATCTACTCGAATTTGGCTTGCAT‑CCGCAAAAACCACG‑TCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATCTATTTATGTTTTACATAACCAATTAAAAAAtt                                                 >  2:76146/1‑229 (MQ=255)                                                                                           taataaAAAATCTTTTTT‑ATGTG‑AATAAATATTTAACCCCCCCCTCT‑AG‑AC‑CATCGGAAAAAATAATGAAAGGGTTCCCC‑TAGACACG‑AATAAAAAATTT‑CC‑AAC‑C‑CACT‑TATCTACTCAAATTTGGTTTGCAT‑CCAAAAAACCCAGT‑TCTCCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACaa                                                   <  1:350564/226‑1 (MQ=255)                                                                                                 aaaacttttttt‑AGGTG‑CATGAAGTTTTAATCCCCCCTCAT‑TG‑AA‑CATAGGGAAAAAAAATGAAGGCGTGTCCC‑TTAAAAAG‑AATACAAAATTT‑TC‑AAC‑A‑CACT‑TATCACCGAGAATTTGG‑TTTCATCCCGAAAAAACCAGG‑TCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGt                                          <  1:636742/223‑1 (MQ=255)                                                                                                        tttttttaattgaatagaaaTTTAAAGCCCCCCCCTGTG‑AG‑CATGGGAAAAAATAAGGAAAGCGTTCCCC‑TTAACACA‑AATAAAAAATTT‑AA‑ACC‑C‑CTCT‑TATCTACTCGATTTTGGCTTGACT‑CCAAAAAAACCAGT‑TAACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGt                                       <  1:307521/209‑1 (MQ=255)                                                                                                          ttt‑ATGTG‑CATGAAGATTTAATGCCTCCTCCT‑TG‑AG‑CATCGGGAGAAATAATGAAAGCGTGTCCC‑TTGACACG‑AATACAAAATTT‑TC‑AAC‑A‑CACT‑TATCTACTCGAATTTGGCTTGCAT‑CAGCAAAAACCAGG‑TCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATACATTTATGTCTTGCATACTCAATTAGACAATTCCCATGTAGAGATTga                                 >  2:98701/1‑228 (MQ=255)                                                                                                                gtg‑AATGAACTTTTAATGCCCCCCCTT‑GG‑AG‑CAGGGGAAAAAATAATGAAGGGGTGCCCG‑TTAACACG‑AATACAAATTTT‑CC‑AAC‑C‑CATT‑TATCTACTAAAATTTGGTTTGCAC‑CCAAAAAACCCAGC‑TCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGt                            <  1:98701/229‑1 (MQ=255)                                                                                                                 tgaaagaaatctttatcgCCTCCTCAT‑GG‑AC‑CATTGGAAAAAATAAT‑AAAACGTTCCCCGTTAACACA‑AATAAAAAATTC‑TA‑ACC‑A‑CCTT‑TATCTACTCGATTTTGGCTTCCAC‑CCAAAAAAACCAGT‑ACACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGt                            <  1:710205/212‑1 (MQ=255)                                                                                                                  tgtaaaaataTTTTATGCCCCCACCC‑TGAAG‑AATCGGGAGAAAAAAAGAACGGGTG‑CCC‑CTGACAAG‑AAAAAAAAATTTCAC‑AAC‑C‑ATCT‑TCTAACCATGAGTGTGGTTGCCAC‑CCACAAAAACGTGT‑ACACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGt                            <  2:167328/220‑1 (MQ=255)                                                                                                                  g‑CATGAAGATTTAATGCCTCCTCCT‑TG‑AG‑CATCGGGAGAAATAATGGAAGCGTGTCCC‑TTGACACG‑AATACAAAATTT‑TC‑AAC‑A‑TTCT‑TATCTTCTCGAATTTGGCTTGCAT‑CCGCAAAAACCAGG‑TCACCGCCATCAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTAAGTCTTGCTTACCAACTTAGACAACTCCCATggtatatttccatagttg                          >  2:74296/1‑212 (MQ=255)                                                                                                                  g‑CATGAAGATTTAATGCCTCCTCCT‑TG‑AG‑CATCGGGAGAAATAATGAAAGCGTGTCCC‑TTGACACG‑AATACAAAATTT‑TC‑AAC‑A‑CACT‑TATCTACTCGAATTTGGCTTGCAT‑CCGCAAAAACCAGG‑TCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTg                          >  1:555811/1‑229 (MQ=255)                                                                                                                                 tGCCTCCTCGT‑TG‑AG‑CAGGGGGAAAAAAAAGGAGAGGTTGCCCG‑TTCACAAA‑AATA‑AAATTTT‑TC‑CACAA‑CAAT‑TATTAACACAAATGGGGGTTGAAC‑CCAAAAAAACAAGG‑TCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACAt            <  1:74296/229‑1 (MQ=255)                                                                                                                                 tGCCTCCTCCT‑TG‑AG‑AGTGGGAAAAAATAAAGAAAGCTGCTCCC‑GTGACA‑G‑AAAAAAAAATTT‑TC‑AAA‑ACCTCT‑TACCTCCTAGATTTTGGTTTCAAT‑CCGAAAAAACAAGG‑TCGCCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACAt            <  1:76146/229‑1 (MQ=255)                                                                                                                                 tGCCTCCCCCT‑AG‑AA‑CATCGGGAAAAATAAAGAAGCCTTGCCCT‑GTAAAAAC‑AATAAAAAATTT‑AC‑ACC‑A‑CTCT‑TATCAACTAAATTTTGGCTTGCAC‑CCGAAAAACCGGGG‑ACACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACAt            <  1:224154/229‑1 (MQ=255)                                                                                                                                     tcctcct‑TG‑AG‑CATTGGGAGAAATAATGAAAGAGTGTCCC‑TTGACACG‑AATACAAAATTT‑TC‑AAC‑A‑CACT‑TATCTACTCGAATTTGGCTTGCAT‑CCGCAAAAACCAGG‑TCACCGCCATTAAAATAGCGTTATTCACCATTAACCATGGTGAGAATATATTTATGGCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCCTAGTTGACTTAAAATTACATaaca        >  1:160036/1‑227 (MQ=255)                                                                                                                                              tgaagatttgtcatatataatGATCGCGTTTCTT‑TTGACACA‑AATAAAAAATTT‑CA‑AAC‑A‑CATT‑TATCTATTAG‑TTTTGCTGTGACA‑CCCAAAAAACCGGGCAAACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAACATATTa  <  2:283497/213‑1 (MQ=255)                                                                                                                                                      cGGGAGAAATAATGAAAGCGTGTCCC‑TTGACACG‑AATACAAAATTT‑TC‑AAC‑A‑CACT‑TATCTACTCGAATTTGGCTTGCAT‑CCGCAAAAACCAGG‑TCACCGCCATTAAAATAGCGATATTCACCATTAAcc                                                                                             >  2:408626/1‑130 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                           AGCATTTCTGCATAGCGAAGATTCTTATTCTGCGATAGTTGCGCAACTTCCTGAAAAACAGG‑AATATCTCTACTACAGTGATGAACGAGTAATGAAGAAACTTTTTT‑ATGTG‑CATGAAGATTTAATGCCTCCTCCT‑TG‑AG‑CATCGGGAGAAATAATGAAAGCGTGTCCC‑TTGACACG‑AATACAAAATTT‑TC‑AAC‑A‑CACT‑TATCTACTCGAATTTGGCTTGCAT‑CCGCAAAAACCAGG‑TCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAACATATTA  >  NC_000913/3665475‑3665838

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