Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,665,698	0	T	C	47.5%	46.2 / 58.7	83	D38G (GAT→GGT)	gadX	acid resistance regulon transcriptional activator; autoactivator
Reads supporting (aligned to +/- strand):  ref base T (28/14);  new base C (0/38);  total (28/55)
Fisher's exact test for biased strand distribution p-value = 3.57e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACAGTGATGAACGAGTAATGAAGAAACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAG‑CGTGTCCCTTGA‑CACGAATACAAAATTTTC‑AACACAC‑TTATCTACTCGAA‑TTTGGCTTGCATCCGCAAAAA‑CC‑AGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAACATATTACTG  >  NC_000913/3665550‑3665842                                                                                                                                                          |                                                                                                                                                   acAGTGATGAACGTGTAATGAAGAAACTTTTTTATGTGCATGAAGATTTAATGCTTCCTCCTTGAGCATCGGGAGAAATAATGAAAG‑CGTGTCTCTTGA‑CACGAATACAAAATTTTT‑AACACAC‑TTATTTATTCGAA‑TTTGGCTTGGATc                                                                                                                                                   >  2:659449/1‑150 (MQ=255) acAGTGATGAACGAGTAATGAAGAAACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAG‑CGTGTCCCTTGA‑CACGAATACAAAATTTTC‑AACACAC‑TTATCTACTCGAA‑TTTGGCTTGCATc                                                                                                                                                   >  2:93168/1‑150 (MQ=255) acAGTGATGAACGAGTAATGAAGAAACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAG‑CGTGTCCCTTGA‑CACGAATACAAAATTTTC‑AACACAC‑TTATCTACTCGAA‑TTTGGCTTGCATc                                                                                                                                                   >  2:127721/1‑150 (MQ=255)    gTGATGAACGAGTAATGAAGAAAATTTTTTCTGTGAGTGAAGGTTTAATGCCTTCTCCTTGAGCATCGTGAGAAATAATGAAAG‑CGTGTCCCTTGA‑CACGAATAAAAAGTTTTC‑AACACAC‑TTATCTACTCGAA‑TTTGGCTTGCATCCGCa                                                                                                                                               >  2:227126/1‑151 (MQ=255)          aaCGAGTAATGAAGAAACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAG‑CGTGTCCCTTGA‑CACGAATACAAAATTTTC‑AACACAC‑TTATCTACTCGAA‑TTTGGCTTGCATCCGCAAAAA‑cc                                                                                                                                        <  2:15459/151‑1 (MQ=255)               gTAATGAAGAAACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAG‑CGTGTCCCTTGA‑CACGAATACAAAATTTTC‑AACACAC‑TTATCTACTCGAA‑TTTGGCTTGCATCCGCAAAAA‑CC‑AGGt                                                                                                                                   >  2:551444/1‑150 (MQ=255)                    gaagaaACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAG‑CGTGTCCCTTGA‑CACGAATACAAAATTTTC‑AACACAC‑TTATCTACTCGAA‑TTTGGCTTGCATCCGCAAAAA‑CC‑AGGTCACCg                                                                                                                              >  1:267630/1‑150 (MQ=255)                     aagaaACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAG‑CGTGTCCCTTGA‑CCCGAATACAAAATTTTC‑AACACCC‑TTATCTACTCGAA‑TTTGGCTTGCATCCGCAAAAA‑CC‑AGGTCACCTcc                                                                                                                            <  2:222699/151‑1 (MQ=255)                     aagaaACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAG‑CGTGTCCCTTGA‑CACGAATACAAAATTTTC‑AACACAC‑TTATCTACTCGAA‑TTTGGCTTGCATCCGCAAAAA‑CC‑AGGTCACCGc                                                                                                                             <  1:407618/150‑1 (MQ=255)                           cTTTTTTATGTGCATGAAGATTTACTGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAG‑CGTGTCCCTTGA‑CACGACTACAAAATTTTC‑AACACAC‑TTATCTACTCGAA‑TTTGGCTTGCATCCGCAAAAA‑CC‑AGCTCACCGCCCTTaaa                                                                                                                      >  2:949149/1‑151 (MQ=255)                           cTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAG‑CGTGTCCCTTGA‑CACGAATACAAAATTTTC‑AACACAC‑TTATCTACTCGAA‑TTTGGCTTGCATCCGCAAAAA‑CC‑AGGTCACCGCCATTaaa                                                                                                                      >  1:406690/1‑151 (MQ=255)                             tttttATGTGCATGAAGATTAAATGCCTCCTCCTTGACAATGGGGAAAAATAATAAAAG‑CGTGTCCCTTGC‑CACGAATACAAAATTTCC‑AACACCC‑TTATCTCCCAAAA‑TTTGGCTTGCATCCGCAAAAC‑CC‑AGGTCACCGCCATTAAAAt                                                                                                                    <  2:921702/151‑1 (MQ=255)                                           aaGTTTTAATGCCCCCCCTTTGGCCATCGGAAGAAAAAAAAAAAG‑CGTCCCCTTTAA‑CCCAAATAAAAAATTTTC‑ACCCCAC‑TTATCTCCCCAAT‑TTTGGTTTGCCCCCAAAAAAC‑CC‑AGGTCACCGCCATTAAAATAGCGATATTCAcc                                                                                                       <  1:1147911/150‑1 (MQ=255)                                             gATTTAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAG‑CGTGTCCCTTGA‑CACGAATACAAAATTTTC‑AACACAC‑TTATCTACTCGAA‑TTTGGCTTGCATCCGCAAAAA‑CC‑AGGTCACCGCCATTAAAATAGCGATATTCACCAtt                                                                                                    >  2:390277/1‑151 (MQ=255)                                             gATTTAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAG‑CGTGTCCCTTGA‑CACGAATACAAAATTTTC‑AACACAC‑TTATCTACTCGAA‑TTTGGCTTGCATCCGCAAAAA‑CC‑AGGTCACCGCCATTAAAATAGCGATATTCACCAtt                                                                                                    >  2:1014639/1‑151 (MQ=255)                                             gATTTAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAG‑CGTGTCCCTTGA‑CACGAATACAAAATTTTC‑AACACAC‑TTATCTACTCGAA‑TTTGGCTTGCATCCGCAAAAA‑CC‑AGGTCACCGCCATTAAAATAGCGATATTCACCAtt                                                                                                    >  1:656595/1‑151 (MQ=255)                                             gATTTAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAG‑CGTGTCCCTTGA‑CACGAATACAAAATTTTC‑AACACAC‑TTATCTACTCGAA‑TTTGGCTTGCATCCGCAAAAA‑CC‑AGGTCACCGCCATTAAAATAGCGATATTCACCAtt                                                                                                    >  2:955301/1‑151 (MQ=255)                                               tttcacgCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAG‑CGTGTCCCTTGA‑CACGAATACAAAATTTTC‑AACACAC‑TTATCTACTCGAA‑TTTGGCCTGCATCCGCAAAAA‑CC‑AGGTCACCGCCATTAAAATAGCGATATTCACCATTaa                                                                                                  >  2:1183052/7‑151 (MQ=255)                                                 tAATGCCTCCTCTTAGAACATCGGGAAAAAAAATGAAAG‑CGTGCCCCTTGA‑CACGAATACAAAATTTTC‑AACACAC‑TTATCTACTCGAA‑TTTGGCTTGCACCCGAAAAAA‑CC‑AGTTCACCGCCATTAAAATAGCGATATTCACCATTAAcc                                                                                                <  2:753836/151‑1 (MQ=255)                                                 tAATCCCTCCCCCTTGAGCATCGGGAAAAAAAAAGAAAG‑CGTGTCCTTTGA‑CACGAATAAAAAATTTTC‑ACCCCAT‑TTTTCTCCTGAAA‑TTTGCCTTCCCCCCGCAAAAC‑CC‑AGTTCCCCGCCATTAAAATAGCGATATTCACCATTAAcc                                                                                                <  1:183514/151‑1 (MQ=255)                                                 tAATCCCCCCTCTTTGACAATGGGGAAAAATAATAAAGG‑CGTTCCCCTTGA‑CACAAATAAAAAATTTTC‑AACCCCC‑TTTTCTACCCAAA‑TTTGGTTTGCACCCACAAAAC‑CC‑AGGTCCCCGCCATTAAAATAGCGATATTCACCATTAAcc                                                                                                <  1:846545/151‑1 (MQ=255)                                                 tAATCCCCCCTCCTTGAGCATGGGGAAAAAAAATGAAAG‑CGTGCCCCTTGA‑CACAAATACAAAATTTTC‑AACACCC‑TTATCTACTCAAA‑TTTGGCTTGCACCCGCAAAAA‑CC‑AGGTCACCGCCATTAAAATAGCGATATTCACCATTAAcc                                                                                                <  1:1022937/151‑1 (MQ=255)                                                 tAATCCCCCCCCCTTGAGCATGGGAAAAAATAAAGAAGG‑CGTCTCCCTTGA‑CACAAAAAAAAATTTTTC‑AACCCCT‑TTTTCTACTCGAT‑TTTGGTTTGCCCCCGAAAAAC‑CC‑TGGTCACCGCCATTAAAATAGCGATATTCACCATTAAcc                                                                                                <  1:1083769/151‑1 (MQ=255)                                                 tAATCCCCCCCCCTTGAGCATCGGGAAAAAAAAGAAAGG‑CGTGCCCCTTGA‑CACAAATAAAAAATTTTA‑AACACCC‑TTTTCTCCTCGAT‑TTGGGTTTCCACCCGAAAAAA‑CC‑GGGTCACCGCCATTAAAATAGCGATATTCACCATTAAcc                                                                                                <  1:1091064/151‑1 (MQ=255)                                                 tAATCCCCCCCCCTTGAAC‑TTGGGAAAAATAAAAAAAGCCTTTTCCCTTAA‑CACAAAAAAAAATTTTTC‑ACCCCCC‑TTATCCACTCGAA‑TTTGCCTT‑CCTCCCCAAAAACCC‑ATATCACCGCCATTAAAATAGCGATATTCACCATTAAcc                                                                                                <  1:341860/151‑1 (MQ=255)                                                 tAACGCCTCCACCTGGAGCACGGGAAGAAATAATGAAGG‑CGTG‑CCCTTGACCACGAAGAAAAAATTTTC‑AACCCAC‑TTTTTTCCCCAAA‑TTTGCCTGACA‑CCAAAAAAA‑CCAAGCTAACCGCCATTAAAATAGCGATATTCACCATTAAcc                                                                                                <  2:797159/151‑1 (MQ=255)                                                  aaTGCCCCCTCCTTGAGCATCGGGAAAAAAAATGAAAG‑CGTGTCCCTTGA‑CACAAATAAAAAATTTTC‑AACACAC‑TTATCTCCCCAAT‑TTTGGCTTGCACCCGAAAAAC‑CC‑AGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCa                                                                                               <  1:518789/151‑1 (MQ=255)                                                     gCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAG‑CGTGTCCCTTGA‑CACGAATACAAAATTTTC‑AACACAC‑TTATCTACTCGAA‑TTTGGCTTGCATCCGCAAAAA‑CC‑AGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATgg                                                                                            >  1:882080/1‑151 (MQ=255)                                                     cccccctcctTGAGCATGGGAAAAAAAAAAAAAGG‑CGTTTCCCTTGA‑CACAAAAAAAAAATTTTC‑AACCCAC‑TTTTCTACTCGAA‑TTTGGCTTGCATCCGAAAAAC‑CC‑AGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATgg                                                                                            <  1:253375/147‑1 (MQ=255)                                                       ctccccctTGAGCATGGGGAAAAAAAATAAAAG‑CGTGTCCTTTGA‑CAAAAAAAAAAAATTTTA‑AACACCT‑TTATCTACTCAAA‑TTTGGTTTGCACCCAAAAAAA‑CC‑GGGTCCCCGCCATTAAAATAGCGATATTCACCATTAACCATGGTg                                                                                          <  2:407349/151‑1 (MQ=255)                                                        tcctcctTGAGCATCGGGAGAAATAATGAAAG‑CGTGTCCCTTGA‑CACGAATACAAAATTTTC‑AACACAC‑TTATCTACTCGAA‑TTTGGCTTGCATCCGCAAAAA‑CC‑AGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTga                                                                                         >  1:438838/1‑151 (MQ=255)                                                          ctcctTGAGCATCGGGAGAAATAATGAAAG‑CGTGTCCCTTGA‑CACGAATACAAAATTTTC‑AACACAC‑TTATCTACTCGAA‑TTTGGCTTGCATCCGCAAAAA‑CC‑AGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTgaga                                                                                       >  1:464912/1‑151 (MQ=255)                                                          ctcctTGAGCATCGGGAGAAATAATGAAAG‑CGTGTCCCTTGA‑CACGAATACAAAATTTTC‑AACACAC‑TTATCTACTCGAA‑TTTGGCTTGCATCCGCAAAAA‑CC‑AGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTgaga                                                                                       >  1:123529/1‑151 (MQ=255)                                                           tcctTGAGCATCGGGAGAAATAATGAAAG‑CGTGTCCCTTGA‑CACGAATACAAAATTTTC‑AACACAC‑TTATCTACTCGAA‑TTTGGCTTGCATCCGCAAAAA‑CC‑AGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGaa                                                                                      >  1:186588/1‑151 (MQ=255)                                                            cctGGAGCATGGGGAAAAAAAATGAAAG‑CGTGTCCCTTGA‑CACAAAAAAAAAATTTTC‑AACACAC‑TTATCTACTCGAT‑TTTGGCTTGCATCCAAAAAAA‑CC‑AGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAat                                                                                     <  1:867985/151‑1 (MQ=255)                                                                 accacgggaaaaaaTAATAAAAG‑CGTGTCCCTTGA‑CACAAATAAAAAATTTTC‑AACACAC‑TTATCTCCTCAAA‑TTTGGCTTGCACCCGAAAAAC‑CC‑AGTTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATAtt                                                                                <  1:910630/140‑1 (MQ=255)                                                                 accATCGGGAAAAATAATAAAAG‑CGTGTCCCTTGA‑CACGAATACAAATTTTTC‑AACCCAC‑TTATCTCCTCAAT‑TTTGGCTTGCCTCCAAAAAAA‑CC‑AGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATAtt                                                                                <  1:158084/149‑1 (MQ=255)                                                                  gcagggggAAAAAAAAAGAAGG‑GGTGTCCCTGGA‑CACGAAAAAAAAATTTTA‑AACCCCC‑TTATCTCCTCGAT‑TTTGGCTTGCACCCGAAAAAA‑CC‑AGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATAtt                                                                                <  1:138591/145‑1 (MQ=255)                                                                   cctggggAAAAATAATAAAAG‑CGTGCCCCTTCA‑CAAAAATAAAAATTTTTC‑AACCCAC‑TTACCTCCGGGAT‑TTTGGGTTCCTCCCAAAAAAC‑CC‑AGGTCCCCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATAtt                                                                                <  1:950842/145‑1 (MQ=255)                                                                   catgggaaaaaaaaaTAAAAG‑CGTGCCCTTTAA‑CACAAAAAAAAAATTTTC‑AACCCCC‑TTTTCTCCTCAAT‑TTGGGTTTCCACCCAAAAAAC‑CC‑AGGACACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATAtt                                                                                <  1:93367/140‑1 (MQ=255)                                                                   cagcggaaaaaaTAATGAAAG‑CGTATCCCTTGA‑CACGAA‑AAAAAATTTTCAAACACCC‑TTATCTACTCAAT‑TTTGTCTTCCTCCCGCAAAAC‑CC‑GGTTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATAtt                                                                                <  2:624483/140‑1 (MQ=255)                                                                   cacggggAAAAATAAGGAAGG‑CGTGCCCCTGGA‑CACAAAAAAAAATTTTCC‑AACACAC‑TTATCTCCTCGAT‑TTTGGCTTGCACCCGAAAAAC‑CC‑AGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATAtt                                                                                <  1:913412/145‑1 (MQ=255)                                                                   cATGGGGAAAAAAAATAAAAG‑CGTTTCCTTTGA‑CACAAATAAAAAATTTTC‑AACACAC‑TTATCTACTCGAA‑TTTGGCTTGCATCCGAAAAAA‑CC‑AGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATAtt                                                                                <  1:163082/149‑1 (MQ=255)                                                                     tCGGGAAAAATAATAAAAG‑CTTGTCCTTGGA‑CACAAATA‑AAAATTTTCAAACCCCC‑TTATCTACTCAAT‑TTTGGCTTGCACCCGAAAAAC‑CC‑AGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATg                                                                            <  2:315078/151‑1 (MQ=255)                                                                      gggaaaaaaTAATGAAAG‑GGTGTCCCTTAA‑CAAGAAAAAAAAATTTCA‑AAACCAT‑TTATCTAAAGAAA‑TTGGATTTGCACCCGAAAAAA‑CC‑AGCACACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATg                                                                            <  2:832687/144‑1 (MQ=255)                                                                       gggAAAAATAAGAAAGG‑CGTGTCCTTTGA‑CACAAAAAAAAAATTTCA‑AACACAT‑TTATCTCTTAAAT‑TTTGTCTTCCACCCCAAAAAA‑CC‑AGTTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGt                                                                           <  2:928582/150‑1 (MQ=255)                                                                         gaaaaaaaaTAAAGG‑CGTGCCCCTTGA‑CACAAATAAAAAATTTCC‑AACCCAT‑TTATCTCCCCAAT‑TTTGGCTTGCACCCGAAAAAC‑CC‑GGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCt                                                                         <  1:1136199/147‑1 (MQ=255)                                                                         aaaaaaTAAGAAAGG‑CGTGTCCCTTGA‑CACAAAAAAAAAATTTTA‑AACACCT‑TTATCTACTCAAA‑TTGGGTTTGCACCCGAAAAAC‑CC‑GGTTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCt                                                                         <  2:186380/147‑1 (MQ=255)                                                                           gaAATAATGAAAG‑CGTGTCCCTTGA‑CACGAATACAAAATTTTC‑AACACAC‑TTATCTACTCGAA‑TTTGGCTTGCATCCGCAAAAA‑CC‑AGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTg                                                                       >  1:170563/1‑150 (MQ=255)                                                                               taatGAAAG‑CGTGTCCCTTGA‑CACGAATACAAAATTTTC‑AACACAC‑TTATCTACTCGAA‑TTTGGCTTGCATCCGCAAAAA‑CC‑AGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATa                                                                   >  2:1037293/1‑150 (MQ=255)                                                                                 atGAAAG‑CGTGTCCCTTGA‑CACGAATAAAAAATTTTA‑AACCCAC‑TTATCTCCTCGAT‑TTTGGTTTGCACCCGAAAAAA‑CC‑AGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGc                                                                <  1:539504/151‑1 (MQ=255)                                                                                      aG‑CGTGTCCCTTGA‑CACGAATACAAAATTTTC‑AACACAC‑TTATCTACTCGAA‑TTTGGATTGCATCCGCAAAAA‑CC‑AGGTCACCGCCATTAAAATAGCGATATTCACCATTAACAATGGTGAGAATATATTTATGTCTTGCATACGCacata                                                           >  2:512739/1‑147 (MQ=255)                                                                                        ccgttccccTTGA‑CACAAAAAAAAAATTTAC‑AACACCC‑TTATCCACCCGAA‑TTTGGCTGGCAACCGCAAAAC‑CC‑AGTTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAg                                                          <  1:1062778/145‑1 (MQ=255)                                                                                           tgtCCCTTGA‑CACGAATACAAAATTTTC‑AACACAC‑TTATCTACTCGAA‑TTTGGCTTGCATCCGCAAAAA‑CC‑AGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATAtt                                                                                >  1:1141309/1‑126 (MQ=255)                                                                                           tgtCCCTAGA‑AACAAAAAAAAATTTTAA‑AACACAT‑TTATCTA‑TCGAATTTTGCCTTGCA‑CCCCAAAAG‑CCAAGGAATCCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAGTTAGAc                                                        <  2:168232/150‑1 (MQ=255)                                                                                           tgccccTTGA‑CAAAAATAAAAAATTTTC‑AACCCCC‑TTATCTCCTCAAT‑TTTGGCTTGCACCCGAAAAAA‑CC‑AGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGAc                                                        <  1:327788/147‑1 (MQ=255)                                                                                           tcccctttGA‑CACGAAAAAAAAATTTAC‑AACACTC‑TTATCTACTAGAT‑TTTGCCTTGCACCCAAAAAAC‑CC‑AGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGAc                                                        <  2:372394/144‑1 (MQ=255)                                                                                              cccTTGA‑CACGAATACAAAATTTTC‑AACACAC‑TTATCTACTCGAA‑TTTGGCTTGCATCCGCAAAAA‑CC‑AGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAAtt                                                    >  1:489092/1‑151 (MQ=255)                                                                                                    a‑cacGAATAAAAAATTTCC‑AACACCT‑TTATCTACTCGAT‑TTTGGCTTGCACCCGAAAAAA‑CC‑AGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATg                                              <  1:314777/151‑1 (MQ=255)                                                                                                      cacaAATAAAAAATTTTC‑AACACTC‑TTATCTACTCGAA‑TTTGGCTTGCATCCGAAAAAA‑CC‑AGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGt                                             <  2:123498/151‑1 (MQ=255)                                                                                                      cacGAATACAAAATTTTC‑AACACAC‑TTATCTACTCGAA‑TTTGGCTTGCATCCGCAAAAA‑CC‑AGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGt                                             >  2:454537/1‑151 (MQ=255)                                                                                                        caaaaaaaaaaTTTTC‑AACCCCT‑TTATCTCCTCGAT‑TTTGGTTTGCACCCAAAAAAC‑CC‑AGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAg                                           <  1:938046/144‑1 (MQ=255)                                                                                                        caaaTAAAAATTTTTC‑AACACTC‑TTATCTACTCGAA‑TTTGGCTTGCCCCCGAAAAAA‑CC‑AGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAg                                           <  1:968307/149‑1 (MQ=255)                                                                                                        caaaTAAAAAATTTTA‑AACACCT‑TTATCTACTCGAT‑TTTGGCTTGCACCCGAAAAAA‑CC‑AGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAg                                           <  1:996224/149‑1 (MQ=255)                                                                                                        aaaaaaaaaaaTTTAA‑AACACTT‑TTTTCAACTCAAA‑TTTGTTTTACTCCCAAAAAAA‑CA‑AGGTCCCCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAg                                           <  1:876541/144‑1 (MQ=255)                                                                                                         caaataaaaaTTTTC‑CACCCCCTTTTTCTCCTCGAT‑TTTGTCTTGCCCCCGAAAAAC‑CC‑TGTTCCCCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAg                                           <  1:547905/145‑1 (MQ=255)                                                                                                            tACAAAATTTTC‑AACACAC‑TTATCTACTCGAA‑TTTGGCTTGCATCCGCAAAAA‑CC‑AGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGAt                                       >  1:324248/1‑151 (MQ=255)                                                                                                            aaaaaaaTTTTC‑AACACAC‑TTATCTCCTCGAA‑TTTGGCTTGCACCCGAAAAAC‑CA‑AGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGAt                                       <  1:390813/148‑1 (MQ=255)                                                                                                             aaaaaaTTTTC‑AACACAT‑TTATCTCCTCAAT‑TTTGGTTTACACCCGAAAAAA‑CC‑AGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGAtt                                      <  1:786082/149‑1 (MQ=255)                                                                                                                 aaTTTCC‑AACACAT‑TTATCTTCCCAAA‑TTTGGCTTGCCCCCGAAAAAA‑CC‑GGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCAt                                  <  1:1140210/151‑1 (MQ=255)                                                                                                                     ttacaccacac‑TTATCTCCTAAAT‑TTGGGGTAGCCCCCAAAAAAA‑CC‑GGGCCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGt                               <  2:164311/145‑1 (MQ=255)                                                                                                                     ttC‑AACACAC‑TTATCTACTCGAA‑TTTGGCTTGCATCCGCAAAAA‑CC‑AGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGtt                              >  2:573026/1‑151 (MQ=255)                                                                                                                       taaaCACCC‑TTACCTACTCGAT‑TTTGCCTTGCACCCGAAAAAA‑CC‑AGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTg                             <  2:373084/149‑1 (MQ=255)                                                                                                                          acacac‑TTATCTACTCGAT‑TTTGGCTTCCATCCGCAAAAC‑CC‑AGGTAACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACt                          <  2:196632/151‑1 (MQ=255)                                                                                                                                ctgaTCTACCCGAA‑TTTCGCGTCCCCCCGCAAAGA‑TC‑‑GCTTCCCGCAATTAAAATAGCTATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTCGTGATTGCATAGTTGACTTAatat                    <  2:831468/148‑1 (MQ=255)                                                                                                                                 ttATCTACTCGAT‑TTTGGCTTGCACCCGAAAAAC‑CC‑AGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAatat                    <  1:513830/151‑1 (MQ=255)                                                                                                                                 ttATCTACTCGAA‑TTTGGCTTGCATCCACAAAAC‑CC‑AGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAatat                    <  2:196550/151‑1 (MQ=255)                                                                                                                                  tATCTACTCGAA‑TTTGGCTTGCATCCGCAAAAA‑CC‑AGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATAtt                   >  1:666616/1‑151 (MQ=255)                                                                                                                                           gAA‑TTTGGCTTGCATCCGCAAAAA‑CC‑AGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAACa          >  2:258841/1‑151 (MQ=255)                                                                                                                                                   gCTTGCCTCCGCAAAAA‑CG‑GGGTTGCCGCCATTAAAATAGCGATATTCGCCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGGCTTAATATTACATAAGCATATTAc    <  2:169825/150‑1 (MQ=255)                                                                                                                                                   gCTTGAATCCGAAAAAA‑CC‑AGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAACATATTAc    <  2:488878/150‑1 (MQ=255)                                                                                                                                                    tttGCACCCGAAAAAA‑CG‑AGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAACATATTACTg  <  1:567349/150‑1 (MQ=255)                                                                                                                                                    cTTGCACCCGAAAAAC‑CC‑AGTTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAACATATTACTg  <  1:573309/151‑1 (MQ=255)                                                                                                                                                          |                                                                                                                                                   ACAGTGATGAACGAGTAATGAAGAAACTTTTTTATGTGCATGAAGATTTAATGCCTCCTCCTTGAGCATCGGGAGAAATAATGAAAG‑CGTGTCCCTTGA‑CACGAATACAAAATTTTC‑AACACAC‑TTATCTACTCGAA‑TTTGGCTTGCATCCGCAAAAA‑CC‑AGGTCACCGCCATTAAAATAGCGATATTCACCATTAACCATGGTGAGAATATATTTATGTCTTGCATACGCAATTAGACAATTCCCATGTAGTGATTGCATAGTTGACTTAATATTACATAAACATATTACTG  >  NC_000913/3665550‑3665842

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