Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,223,638	T→C	intergenic (‑10/‑190)	iclR ← / → metH	DNA‑binding transcriptional repressor IclR/cobalamin‑dependent methionine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,223,638	0	T	C	95.8%	140.4 / ‑3.2	48	intergenic (‑10/‑190)	iclR/metH	DNA‑binding transcriptional repressor IclR/cobalamin‑dependent methionine synthase
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base C (8/38);  total (10/38)
Fisher's exact test for biased strand distribution p-value = 3.99e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.91e-01

AGTTCGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGC‑GGGAATGG‑GTGCGACCATGACAGTCT‑C‑TTTTTTCTGTATCGTGG‑AAATC‑ATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTT‑AATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATG‑CCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAAC‑TGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTG  >  NC_000913/4223393‑4223884                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                            agtgcGCCAACCTGACGCCCGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGC‑GGGAATGG‑GTGCGACCATGACAGTCT‑C‑cttttt                                                                                                                                                                                                                                                        <  2:23848/247‑1 (MQ=255)        cAACCGGACGAACGAAACCCTGCTGTTGACTCGTGGTGAGCAGGCGGTGGGTCGGGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGC‑GGGAATGG‑GTGCGACCATGACAGTCT‑C‑CTGTTTCTGTATc                                                                                                                                                                                                                                                 <  2:304133/251‑1 (MQ=255)              gccGCACGAAACCCTTCTGTTGCATCGTGGTGAGCAGGCGGTGGGTCGGGGAATCGGGTAACCCGTCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGC‑GGGAATGG‑GTGCGACCATGACAGTCT‑C‑CTTTGTCTGTAGCGTGG‑a                                                                                                                                                                                                                                           <  2:650852/248‑1 (MQ=255)               aCGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCTCGTTTCGC‑GGGAATGG‑GTGCGACCATGACAGTCT‑C‑CTTTTTCTGTATCGTGG‑aaa                                                                                                                                                                                                                                         <  1:59771/251‑1 (MQ=255)                cGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGC‑GGGAATGG‑GTGCGACCATGACAGTCT‑C‑CTTTTTCTGTATCGTGG‑AAAt                                                                                                                                                                                                                                        <  1:438379/251‑1 (MQ=255)                  cACGAATCCCTGCTGTTGCAGCGTTGTTAGCAGGCGGTGGGTCTTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGC‑GTGAATGG‑GTGCGACCATGACAGTCT‑C‑CTTTTTCTGTATCGTGG‑AAAtc‑a                                                                                                                                                                                                                                     <  2:234586/251‑1 (MQ=255)                  cACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTTGCAACGGCGGGTTTTCTGCCGCGTTTCGC‑GGGAATGG‑GTGCGACCATGACAGTCT‑C‑CTTTTTCTTTATCGTGG‑AAAtc‑a                                                                                                                                                                                                                                     <  2:542748/251‑1 (MQ=255)                   aCGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGC‑GGGAATGG‑GTGCGACCATGACAGTCT‑C‑CTTTTTCTGTATCGTGG‑AAAtc‑at                                                                                                                                                                                                                                    <  1:500283/251‑1 (MQ=255)                      aaaCCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGC‑GGGAATGG‑GTGCGACCATGACAGTCT‑C‑CTTTTTCTGTATCGTGG‑AAAtc‑atttt                                                                                                                                                                                                                                 <  1:663700/251‑1 (MQ=255)                                    cATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGC‑GGGAATGG‑GTGCGACCATGACAGTCT‑C‑CTTTTTCTGTATCGTGG‑AAATC‑ATTTTCATTTTTATTGTTa                                                                                                                                                                                                                   <  1:366741/251‑1 (MQ=255)                                     aTCGTGGTGAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGC‑GGGAATGG‑GTGCGACCATGACAGTCT‑C‑CTTTTTCTGTATCGTGG‑AAATC‑ATTTTCATTTTTATTGTTAg                                                                                                                                                                                                                  <  1:33380/251‑1 (MQ=255)                                                    cGGTGGGTCGTGGCATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGTCCCGTTTTGC‑GGGGATGG‑GGGGGTACACGCCAGCCT‑CCTTTTTTTTGTA‑CGTGG‑TAATT‑ATTTTCATTTTTTTTGTTTTCTAATGCCATaattt                                                                                                                                                                                                   >  1:543364/1‑247 (MQ=255)                                                                          aCCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGC‑GGGAATGG‑GTGCGACCATGACAGTCT‑C‑CTTTTTCTGTATCGTGG‑AAATC‑ATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTT‑aa                                                                                                                                                                            <  1:259948/251‑1 (MQ=255)                                                                          aCCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCCCGTTTTGC‑GGGGATTG‑GTGCGGCCCTGACCGTCT‑C‑CTTTTTTTTTTTTGTGG‑GAATC‑ATTTTTTTTTTTATTTTTTACTTATTCAATTGTTTCTGAACATATCCGCTGAGttt                                                                                                                                                                              >  1:141071/1‑249 (MQ=255)                                                                              ggCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGC‑GGGAATGG‑GTGCGACCATGACAGTCT‑C‑CTTTTTCTGTATCGTGG‑AAATC‑ATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTT‑AATGtt                                                                                                                                                                        <  1:616843/251‑1 (MQ=255)                                                                                  tgttgCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGC‑GGGAATGG‑GTGCGACCATGACAGTCT‑C‑CTTTTTCTGTATCGTGG‑AAATC‑ATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTT‑AATGTTGa                                                                                                                                                                      <  1:151666/249‑1 (MQ=255)                                                                                   gtttcgctaGTTCCTTGAGTGCCAACCTGCCTGTGGATGCGGGAAGCTACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTTCGGGGGCAACGGCGGGTTTTCTGCCGCGTTTCGT‑GGGAATGG‑GTGCGACCATGACAGTCT‑C‑CTTTTTCTGTATCGTGG‑AAATC‑ATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTT‑AATGTTGAAc                                                                                                                                                                    <  2:492373/242‑1 (MQ=255)                                                                                   gttgCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGC‑GGGAATGG‑GTGCGACCATGACAGTCT‑C‑CTTTTTCTGTATCGTGG‑AAATC‑ATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTT‑AATGTTGAACa                                                                                                                                                                   <  1:408136/251‑1 (MQ=255)                                                                                        gcCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGC‑GGGAATGG‑GTGCGACCATGACAGTCT‑C‑CTTTTTCTGTATCGTGG‑AAATC‑ATTTTCATTTTTATTGTTAGCTAATGCAATAGTTTCTGAACTGATCCGATGAGTT‑AATGTTGAACAAAtct                                                                                                                                                              <  2:239041/251‑1 (MQ=255)                                                                                        gcCAGTGCCGTGAGGTCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGC‑GGGAATGG‑GTGCGACCATGACAGTCT‑C‑CTTTTTCTGTATCGTGG‑AAAGC‑ATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTT‑AATGTTGAACAAAtct                                                                                                                                                              <  1:471938/251‑1 (MQ=255)                                                                                                               gCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTCTGCCGCGTTTCGC‑GGGGAAGG‑GTGGGGGCCTGACAGGTT‑C‑CTTTTTTTGTATCGGGGGAATTC‑ATTTT‑ATTTTTATTGTTTGCTAATGCAATAGTTAATGAAATGATCCGATGAGTT‑TATGTTGAAAAAATATCATGTTGCgtggggggtgctttt                                                                                                                                       >  1:12162/1‑240 (MQ=255)                                                                                                                                 ccACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCGGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGC‑GGGAATGG‑GTGCGACCATGACAGTCT‑C‑CTTTTTCTGTATCGTGG‑AAATC‑ATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTT‑AATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATg                                                                                                                     <  2:17276/251‑1 (MQ=255)                                                                                                                                  cccTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGGCCAGACGCTGGGGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGC‑GGGAATGG‑GTGCGACCATGACAGTCT‑C‑CTTTTTCTGTATCGTGG‑AAATC‑ATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTT‑AATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTAACACAGATGCGTTTATG‑c                                                                                                                   <  2:521504/249‑1 (MQ=255)                                                                                                                                  cACTCCAGGAATCTCAGTCCACGCGTTAAAGACTGAACCTGTCCAGTCGCGGGGGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGC‑GGGAATGG‑GTGCGACCATGACAGTCT‑C‑CTTTTTCTGTATCGTGG‑AAATC‑ATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTT‑AATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATG‑c                                                                                                                   <  2:698130/251‑1 (MQ=255)                                                                                                                                   aCTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGC‑GGGAATGG‑GTGCGACCATGACAGTCT‑C‑CTTTTTCTGTATCGTGG‑AAATC‑ATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTT‑AATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATG‑cc                                                                                                                  <  1:503349/251‑1 (MQ=255)                                                                                                                                         gtatgtgcAGGCCTCGTGTTAAATAATGAACCTGTCCAGTCGCTGTGGCTGTGGTATCGGCGGGTTTTCTGCCGCGTTTCGC‑GGGTATGGAGT‑CGACCATGACATGCT‑C‑CTTTTTCTGTAACGTGG‑TCATT‑ATTTTCATTTTTATTGTTAGCTAATGCAATAGTTGCTGAACTGATCCGACGAGTT‑AATGTTGAACAAATCTCATGTTGCGTGGGTGTCGCTTGTACCACAGAttcttttaagtcagta                                                                                                               <  2:259824/242‑16 (MQ=255)                                                                                                                                           aaTTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGACTTTCGG‑GGGAATGG‑GGGGGACCCAGCCAGCCT‑C‑CCTTTTTTGTTTCGTGG‑GAATCCTTTTTCATTTT‑ATTGTTGGCTAATACAATAGTTACTGAACTGATCCGATGAGTT‑AATGTTGAAAAAAACTCCTGCTGAGTGGTGGGCCCTTTTTACACCGAGGCGTATATT‑CCCGTATGGt                                                                                                          >  1:698965/1‑251 (MQ=255)                                                                                                                                                 cggCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCTC‑GGTAATGG‑GTGCGACCATGACAGTCT‑C‑CTTTTTCTGTATCGTGG‑AAATC‑ATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTT‑AATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATG‑CCAGTATGGTttgttg                                                                                                    <  2:661804/250‑1 (MQ=255)                                                                                                                                                 aGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTTTGTCGCGGTTTGC‑GGGGATGG‑GTGCGGCCCTGACAGTCTAC‑TTTTTT‑TGTATCATGG‑AAATC‑ATTTTCATTTTTTTTTTTTGCTTATTTAATAGTTTCTGTACTGATCCGATGAGTTTAAT‑TTGTACAAATCTCATTTTTCGTGGTGGTCGCTTTTTCCACAAATGCGTTTTTT‑CCAGTTTGGTttgttt                                                                                                    >  1:252436/1‑250 (MQ=255)                                                                                                                                                            tAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTTTGCCGCGTTTCGC‑GGGGATGG‑GGGGGACCATGACAATCT‑C‑CTTTTTTTTTTTCGTGG‑GAATC‑ATTTTTCTTTTTTTTTTTTGCTAATGCAATAGTTAATTAACTGATCCGATGAGTT‑TATGTTGTACAAATTTCATGTTGGGGGGTGGGCGCTTTTTCCACAGAAGAGTTTATGTCCGGT‑TGGGTTTTTGAATTTttattt                                                                                         >  1:547545/1‑250 (MQ=255)                                                                                                                                                               aGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGC‑GGGAATGG‑GTGCGACCATGACAGTCT‑C‑CTTTTTCTGTATCGTGG‑AAATC‑ATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTT‑AATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATG‑CCAGTATGGTTTGTTGAATTTTTATTAAAt                                                                                      <  1:641647/251‑1 (MQ=255)                                                                                                                                                                    gAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGC‑GGGAATGG‑GTGCGACCATGACAGTCT‑C‑CTTTTTCTGTATCGTGG‑AAATC‑ATTTTCATTTTTATTGTTAGc                                                                                                                                                                                                                 <  2:558812/125‑1 (MQ=255)                                                                                                                                                                    gAACCTGTCCAGGCTCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGC‑GGGAATGG‑GTGCGACCATGACAGTCT‑C‑CTTTTTCTGTATCGTGG‑AAATC‑ATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTT‑AATGTTGAACAAATCTCATGTTGCGTGGTTGTCGCTTTTACCACAGATGCGTTTATG‑CCAGTATGGTTTGTTGAATTTTTATTAAATCTg                                                                                   <  1:703658/249‑1 (MQ=255)                                                                                                                                                                    gAACCTGACCGGTAGCTGTTGCGGGGTCAACTGCGGTTTTGCTGCCGCGTTCAGC‑GGGAATGG‑GTGTGTCCATGACAGGCT‑C‑CGTTGTCTTTATCGTGG‑AAATC‑AGTTTCAGTGTTATTGTTATCTAAGACTATAGTGACTGAACTGTTCCGATGAGTT‑AATGTTGATCGAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATG‑CCAATATTGTGTGTTGCATTTttattacaactg                                                                                   <  2:20819/249‑7 (MQ=255)                                                                                                                                                                          gccccgtctgtgggtctgTGGCAACGGCGTGTTTTCTGCCGCGTTCCGCTGGCAATGT‑‑TGCGAACATTACAGTCT‑T‑CTTTTTCTTTATCGGGG‑AAATC‑TTTTTCATTTTTATTGTAATGTCATGCAATTGTTACTGCACTTATCCGATGAGTT‑AATGTTGAACAAATCTCATGTTGTGTGGTGTTCGCTTTTACCACAGATGCGGTTATG‑CCAGTCTGGTTTGTTGATTTTTGATTACATCTGGTTTGAg                                                                            <  2:267873/233‑1 (MQ=255)                                                                                                                                                                            cccgTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGC‑GGGAATGG‑GTGCGACCATGACAGTCT‑C‑CTTTTTCTGTATCGTGG‑AAATC‑ATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTT‑AATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATG‑CCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCg                                                                          <  2:141071/247‑1 (MQ=255)                                                                                                                                                                                  gCTGGTGCGGTGGCAACGGCGGGTTTTTTGCCGCGGTTCGC‑GGGGATGG‑GGGCGACCATGACAGGCT‑C‑CTTTTTCTGTATCGGGG‑AAATC‑ATTTTCATTTTTTTTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATTGGTT‑AATGTTGAACAAATCTCAAGTTTCGGGGGGGTCGCTTTTACCACAGATGCGTTTATG‑CCCGTATGGGTTGTTGAATTTTTATTAAATTTTGGTTGAGGGTGGCggg                                                                   >  1:404215/1‑251 (MQ=255)                                                                                                                                                                                             ggCAACGGCGGGTTTTCTGCCGCGTTTCGC‑GGGAATGG‑GTGCGACCATGACAGTCT‑C‑CTTTTTCTGTATCGTGG‑AAATC‑ATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTT‑AATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTa                                                                                                                                      <  1:173610/174‑1 (MQ=255)                                                                                                                                                                                                 aCGGCGGGTTTTCTGCCGCGTTTCGC‑GGGAAGGG‑GTGCGACCATGACAGTCT‑C‑CTGTTTCTGTATCGTGT‑AAATC‑ATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTT‑AATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATG‑CCAGTATGGTGTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGGGagcag                                                    <  2:203881/251‑1 (MQ=255)                                                                                                                                                                                                               gCCGCGTTTCGC‑GGGAATGG‑GTGCGACCATGACAGTCT‑C‑CTTTTTCTGTATCGTGG‑AAATC‑ATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTT‑AATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATG‑CCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGaacaac                                      <  1:323126/251‑1 (MQ=255)                                                                                                                                                                                                                       tCGC‑GGGAATGG‑GTGCGACCATGACAGTCT‑C‑CTTTTTCTGGATCGTGG‑AAATC‑ATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTT‑AATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATG‑CCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAAC‑tgcgtgc                              <  2:12162/250‑1 (MQ=255)                                                                                                                                                                                                                        cgc‑gGGAATGG‑GTGCGACCATGACAGTCT‑C‑CTTTTTCTGTATCGTGG‑AAATC‑ATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGAGGAGTT‑AATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGAGGCGTTTATG‑CCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAAC‑tgcgtgcg                             <  2:407196/250‑1 (MQ=255)                                                                                                                                                                                                                              gAATGG‑GTTCGACCATGACAGTCT‑C‑CTTTTTCTGTATCGTGG‑AAATC‑ATTTTCATTTTTATTGGTAGCTAATGCAATAGTTACTGAACTTATCCGATTAGTT‑AATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATG‑CCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAAC‑TGCGTGCGCAGTTa                       <  1:178510/251‑1 (MQ=255)                                                                                                                                                                                                                                aTGG‑GTGCGACCATGACAGTCT‑C‑CTTTTTCTGTATCGTGG‑AAATC‑ATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTT‑AATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATG‑CCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAAC‑TGCGTGCGCAGTTaaa                     <  2:547545/251‑1 (MQ=255)                                                                                                                                                                                                                                   gggagcggccATGACAGTCT‑C‑CTTTTTCTGTATCGTGG‑AAATC‑ATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGCTGAGTT‑AATGTTGAACAAATCTACTGTTGCGTGGTGGTCGCATTTACCACAGATACGTTTATG‑CCAGTATGGTTTGTTGACTTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGACGAACAAAGTGGAA‑AACATGCGTGCGCAGTTAAATg                   >  2:523781/9‑251 (MQ=255)                                                                                                                                                                                                                                       gCGACCATGACAGTCT‑C‑CTTTTTCTGTATCGTGG‑AAATC‑ATTTTCATGTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTT‑AATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATG‑CCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAAC‑TGCGTGCGCAGTTAAAGGAAc                <  1:493671/250‑1 (MQ=255)                                                                                                                                                                                                                                           ccATGACAGTCT‑C‑CTTTTTCTGTATCGTGG‑AAATC‑ATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTT‑AATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATG‑CCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAAC‑TGCGTGCGCAGTTAAATGAACGTAtt           <  1:196072/251‑1 (MQ=255)                                                                                                                                                                                                                                              tGACAGTCT‑C‑CTTTTTCTGTATCGTGG‑AAATC‑ATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTT‑AATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATG‑CCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAAC‑TGCGTGCGCAGTTAAATGAACGTATTCTg        >  1:264420/1‑251 (MQ=255)                                                                                                                                                                                                                                                    tct‑c‑CTTTTTCTGTATCGTGG‑AAATC‑ATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGCTGAGTT‑AATGTTGAACAAATCTAATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATG‑CCAGAATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGCGCAAGTGTGAGCAGCAAAGTGGAACAAC‑TGCGTGCGCAGTTAAATGAACGTATTCTGGAGCTg  >  1:396594/1‑251 (MQ=255)                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                            AGTTCGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGC‑GGGAATGG‑GTGCGACCATGACAGTCT‑C‑TTTTTTCTGTATCGTGG‑AAATC‑ATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTT‑AATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATG‑CCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAAC‑TGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTG  >  NC_000913/4223393‑4223884

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