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	96.4%	164.9 / ‑2.4	56	intergenic (‑10/‑190)	iclR/metH	DNA‑binding transcriptional repressor IclR/cobalamin‑dependent methionine synthase
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base C (15/39);  total (15/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.58e-01

AGTTCGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTG‑CGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTT‑TCGCGGG‑AATGGGTGCGACCATG‑ACAGTC‑T‑CTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTA‑ATGCAATAGTTACTGAAC‑TGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGG‑TTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTG  >  NC_000913/4223393‑4223878                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                    aGTTCGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTG‑CGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTT‑TCGCGGG‑AATGGGTGCGACCATG‑ACAGTC‑T‑Ccttttt                                                                                                                                                                                                                                                <  1:499596/251‑1 (MQ=255)     cGCCACCCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTATCCCGGCTTGTTG‑CGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGGCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTT‑TCGCGGG‑AATGGGTGCGACCATG‑ACAGTC‑T‑CCTTTTTCTGt                                                                                                                                                                                                                                            <  2:365554/251‑1 (MQ=255)      gCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTG‑CGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTT‑TCGCGGG‑AATGGGTGCGACCATG‑ACAGTC‑T‑CCTTTTTCTGTa                                                                                                                                                                                                                                           <  2:195870/251‑1 (MQ=255)              taCGCCCGAAACCCTGCTATTTCATCGTGGTTAGCGGGCGGTGGGGCGTGGAATTGGGTAACCCGTCTTGTTG‑CGCCAGTTCCGTGAGTGCCACACTGCCATTGGAGTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTT‑TCGCGGG‑AATGGGTGCGACCATG‑ACAGTC‑T‑CCTTTTTCTGTATCGTGGa                                                                                                                                                                                                                                    <  2:342704/249‑1 (MQ=255)                  cccgataCCCTGCTTTTTCACCGGGGGTAGCAGGCGGTTGGTCGTGGAATGGGGTAACCCGTGTTGTTG‑CGCCAGTTTCGTGAGGGCCACACTGCCATTGTATTCGGCAATCCACTCCAGTAAGGTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTT‑TCGCGGG‑AATGTTTGCGACCATG‑ACAGTCGT‑CTTTTT‑CTGTATCGCGGAAAtca                                                                                                                                                                                                                               <  2:548701/245‑1 (MQ=255)                    cGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTG‑CGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTT‑TCGCGGG‑AATGGGTGCGACCATG‑ACAGTC‑T‑CCTTTTTCTGTATc                                                                                                                                                                                                                                         <  1:316464/239‑1 (MQ=255)                         cccTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTG‑CGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTT‑TCGCGGG‑AATGGGTGCGACCATG‑ACAGTC‑T‑CCTTTTTCTGTATCGTGGAAAtcattttcat                                                                                                                                                                                                                        <  1:98074/251‑1 (MQ=255)                             gctgTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTG‑CGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTT‑TCGCGGG‑AATGGGTGCGACCATG‑ACAGTC‑T‑CCTTTTTCTGTATCGTGGAAAtcattttcattt                                                                                                                                                                                                                      <  1:499219/249‑1 (MQ=255)                                gtcgCATCGTGGTTAACACGCGGTGGGTCGTGGAATGGGTTATCCCTGCTTCTGG‑CGCCTGTACCGTGAGTGCAACACTGCCATTGGAGTCGGCAATCCACTCCAGTAATGGCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACTGCTGG‑TTTTCTGCCGCGTT‑TCGAGGG‑AATGGTTGCGACCATG‑ACAGTC‑T‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTAt                                                                                                                                                                                                                  <  2:208126/247‑1 (MQ=255)                                 ttGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTG‑CGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑GTTTTTGCCGCGTT‑TTGCGGG‑GATGGGGGGGACCCTG‑ACAGGC‑T‑CCCTTTTCTGTTTCGTGGGAATAatttttattttttttt                                                                                                                                                                                                                >  1:236385/1‑247 (MQ=255)                                    cATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTG‑CGCCAGTTCCGTGAGTTCCATACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTT‑TCGCGGG‑AATGGGTGCGACCATG‑TCAGTC‑T‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTa                                                                                                                                                                                                             <  1:31149/251‑1 (MQ=255)                                                                 aaTTGGGTAACCCGGCTTGTTG‑CGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTT‑TCGCGGG‑AATGGGTGCGACCATG‑ACAGTC‑T‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTA‑ATGCAATAGTTACTGAAC‑TGATCCg                                                                                                                                                                              <  1:196607/251‑1 (MQ=255)                                                                       gTAACCCGGCTTGTTG‑CGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTT‑TCGCGGG‑AATGGGTGCGACCATG‑ACAGTC‑T‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTA‑ATGCAATAGTTACTGAAC‑TGATCCGATGAGt                                                                                                                                                                        <  1:159681/251‑1 (MQ=255)                                                                         aaCCCGGCTTGTTG‑CGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGATCCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTT‑TCGCGGG‑AATGGGTGCGACCATG‑ACAGTC‑T‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTA‑ATGCAATAGTTACTGAAC‑TGATCCGATGAGTTa                                                                                                                                                                      <  1:147302/251‑1 (MQ=255)                                                                         aaCCCGGCTTGTTG‑CGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑GTTTCTGCCGCGTT‑TTGCGGG‑GATGGGGGCGGCCCTG‑ACAGAC‑T‑CCTTTTTTTGTTTCGTGGAAATCATTTTTTTTTTTTTTTTTAGCTA‑ATTCAATAGTTTATGTAC‑TGATCCGATGGGttt                                                                                                                                                                      >  1:390696/1‑250 (MQ=255)                                                                           cccGGCT‑GTTGGCGCCCGTTCTGTGAGTGCCCCACTGCCATTGGATTCGGCAATCCACGCCATTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGGGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTT‑TCGCGGG‑AATGGGTGCGACCATG‑ACAGTC‑T‑CCTTTTTCTGTATCGTGGAAATCATTTTCATGTTTATTGTTAGCTA‑ATGCAATAGTTACTGAAC‑TGATCCGATGAGTTAAt                                                                                                                                                                    <  2:233786/251‑1 (MQ=255)                                                                               gCTTGTTG‑CGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTT‑TCGCGGG‑AATGGGTGCGACCATG‑ACAGTC‑T‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTA‑ATGCAATAGTTACTGAAC‑TGATCCGATGAGTTAATGTTg                                                                                                                                                                <  1:524391/251‑1 (MQ=255)                                                                                      g‑cgcCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTT‑TCGCGGG‑AATGGGTGCGACCATG‑ACAGTC‑T‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTA‑ATGCAATAGTTACTGAAC‑TGATCCGATGAGTTAATGTTGAACaaa                                                                                                                                                          <  1:434237/250‑1 (MQ=255)                                                                                      g‑cgcCAGTGCCGTGAGTGCCACACTGCCATGGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAATGGCGTG‑TTGTCTGCCGCGTT‑TCGCGGG‑AATGGGTGCGTCTATG‑ACAGTC‑T‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTA‑ATGCGATAGTTTCTGTAC‑TGATCCGATGAGTTAATGTTGAACaaa                                                                                                                                                          <  1:479884/250‑1 (MQ=255)                                                                                         gcCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTT‑TCGCGGG‑AATGGGTGCGACCATG‑ACAGTC‑T‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTA‑ATGCAATAGTTACTGAAC‑TGATCCGATGAGTTAATGTTGAACAAAtct                                                                                                                                                       <  1:347252/251‑1 (MQ=255)                                                                                           cAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTTTGCCGCGGT‑TTGCGGG‑GATGGGGGCGGCCATGGAC‑GGA‑T‑CCTTTTTCTTTATCGTGGGAATCATTTTTCTTTTTTTTTTTTGCTA‑ATGCAATAGTTACTGTAC‑TGATCCGATGGGTTATTTTTTAACaaaaatc                                                                                                                                                      >  1:6888/1‑246 (MQ=255)                                                                                                        gCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTT‑TCGCGGG‑AATGGGTGCGACCATG‑ACAGTC‑T‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTA‑ATGCAATAGTTACTGAAC‑TGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCgtggtg                                                                                                                                         <  1:11510/250‑1 (MQ=255)                                                                                                        gCCACACTGCCACTGGATTCGGCAATGCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTT‑TCGCGGG‑AATGGGTGCGACCATG‑ACAGTC‑T‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTA‑ATGCAATAGTTACTGAAC‑TGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCgtggtg                                                                                                                                         <  2:390696/250‑1 (MQ=255)                                                                                                           acacTGCCATTGGATTCGGCAATCCACTCCAGTAATGTCAGGCCACGCGTTAAAGACTGAACCTGGCCAGTCGCTGGGGCGGTTGCAACGGCGGG‑TTTTCTGCCGCGTT‑TCGCGGG‑AATGGGTGCGACCATG‑ACAGTC‑T‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTA‑ATGCAATAGTTACTGAAC‑TGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCg                                                                                                                                     <  1:121584/251‑1 (MQ=255)                                                                                                                          tCGGCAATCCACTCCAGTAATTTCAGGCCCCGCTTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTT‑TCGCGGG‑AATGGGTGCGACCATG‑ACAGTC‑T‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTA‑ATGCAATAGTTACTGAAC‑TGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGAt                                                                                                                       <  2:504008/250‑1 (MQ=255)                                                                                                                           cggctatctctgcgaatAATTGCAGTCCACGCGTTAACTACTGAATCTGTACAGTTGGTTGGGCGGTGGCACCGGATGT‑TTTTCTGGC‑CGTTGTCGTGGG‑AATGTGTGTGACCATG‑TCATTC‑T‑TCTGTTTCTTTATCGTGTAAATCTGTTTCATTTTTATCGTTAGTTA‑ATGCAATTGTTACTGAAC‑TGATCCGATGAGTTAATGTTGAACAAATGTCATGTTGCGTGGTGTTCGTTTTTACCACATATgt                                                                                                                     <  2:200/235‑2 (MQ=255)                                                                                                                               aaTCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTT‑TCGCGGG‑AATGGGTGCGACCATG‑ACAGTC‑T‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTA‑ATGCAATAGTTACTGAAC‑TGATCCGATGAGTTAATGATGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACTACAGATGCGttt                                                                                                                 <  1:110275/251‑1 (MQ=255)                                                                                                                                      tCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGAG‑TTTTCTTCCGCGGT‑TCTCGGG‑GATAGGGGCGACCATG‑GGAGTC‑T‑CCTTTTTTTTTTTCTTGGAGAGCTCTTTTACTTTTATTGTTAGATATAT‑CACAAGTTTATGAACCTG‑TCCCATGAGGTAATGTTGGACAAATCTCATGGTGTGTGGGGGTCGTTTTTatcccaaagtggattttgtccg                                                                                                          >  1:458919/1‑230 (MQ=255)                                                                                                                                       ccAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTTTGTCGCGTT‑TCGCGGGGAAGGGGGGGGAC‑ATG‑GCAGGC‑T‑CCCTTTTTTGTATCGTGGAGATTCTTTTTTTTTTTTTTTTTAGCTA‑ATGCAATAGTTACTGAAA‑TGTTTCGGTGGGGTTATGTTGTACAAATTTCCTGTTGCGTTGTGGTCGCtttttcacccaaagcggtttttgcggg                                                                                                         >  1:111201/1‑228 (MQ=255)                                                                                                                                       ccAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTG‑TCGCGGG‑AATGGGTGCGACCATG‑ACAGTC‑T‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTA‑ATGCAATAGTTACTGAAC‑TGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGt                                                                                                         <  1:257474/251‑1 (MQ=255)                                                                                                                                               ttcatgcccagcgTTAAAGTCTGAACCTCTCCAGTCGCTGGTGCGGTGGCAACGGCGGT‑TTTTCTTCCGCGTT‑TCGAGGG‑AATGGGTTCGACCATT‑ACATTC‑T‑CCTTTTTCTGTATCGTGTAAATCATGTTCTTGTTTATTGTTAGCTA‑ATTCAACAGTTCCTGAGC‑TGATCCGATGAGTTAATGTTGAACAAAGCTCATGTTGCGTGGTGGTCACTTTTACCACAGATGCGTTTATGCCAGTATGG‑ttt                                                                                                 <  2:135914/240‑1 (MQ=255)                                                                                                                                               ttCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTT‑TCGCGGG‑AATGGGTGCGACCATG‑ACAGTC‑T‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTA‑ATGCAATAGTTACTGAAC‑TGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGc                                                                                                                     <  1:82063/231‑1 (MQ=255)                                                                                                                                               ttCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTT‑TCGCGGG‑AATGGGTGCGACCATG‑ACAGTC‑T‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTA‑ATGCAATAGTTACTGAAC‑TGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGG‑ttt                                                                                                 <  1:333538/250‑1 (MQ=255)                                                                                                                                                     ccACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTT‑TCGCGGG‑AATGGGTGCGACCATG‑ACAGTC‑T‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTA‑ATGCAATAGTTACTGAAC‑TGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTAt                                                                                                       <  1:512987/239‑1 (MQ=255)                                                                                                                                                     ccACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTT‑TCGCGGG‑AATGGGTGCGACCATG‑ACAGTC‑T‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTA‑ATGCAATAGTTACTGAAC‑TGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGG‑TTTGTTGAAt                                                                                          <  1:15768/251‑1 (MQ=255)                                                                                                                                                              aaaGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTT‑TCGCGGG‑AATGGGTGCGACCATG‑ACAGTC‑T‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTA‑ATGCAATAGTTACTGAAC‑TGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGG‑TTTGTTGAATTTttatta                                                                                  <  1:88411/250‑1 (MQ=255)                                                                                                                                                                       aCCTGTCCAGTCGCTGGTGCGGTGGCAACGTCGGG‑TTTTCTGCCGCGTT‑TCGCGGG‑AATGGGTGCGACCATG‑ACAGTC‑T‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTA‑ATGCAATAGTTACTGAAC‑TGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGAGGCGTTTATGCCAGTATGG‑TTTGTTGAATTTTTATTAAATCGGGGtt                                                                        <  1:226071/251‑1 (MQ=255)                                                                                                                                                                          tGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑GTTTCTGCCGCGTT‑TTGCGGG‑GATGGGTGCGGCCATG‑ACAGTC‑T‑CCTTTTTCTGTATCGTGGGAATCATTTTCATTTTTTTTTTTATCTT‑ATGCAATAGTTACTGAAC‑TGATCCGCTGAGTTAATGTTTTACAAATCTCATGTTTCGTGGTGGTCGCTTTTTCCACAGATGCGTTTTTGTCAGTATGGGTTTTTTGA‑TTTTTTTTAAATCTGGGTTggg                                                                     >  1:4538/1‑249 (MQ=255)                                                                                                                                                                              cAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTTGCC‑CGTT‑TCGCGGG‑GAAGGGGGCGACCATG‑ACAGTC‑T‑CCTTTTTTTGTATCGTGGGAAACATTTTTATTTTTTTTGTTAGCTA‑ATGCAATAGTTAATGAAA‑TGATCCGATGAGTTTATGTTTAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCCCAGATGCGGTTATGCCAGTATGG‑TTTTTTTAATTTTTTTTAAATCTGGGtttaggggg                                                                 >  1:293367/1‑244 (MQ=255)                                                                                                                                                                                      ggTGCGGTGGCATCGGTGGG‑TTTTCTGCCGTGTT‑TCGCGGG‑AATGGGTGCGACCAGG‑ACAGTC‑T‑CCTTTTGCTGTTTCGTGGAAATCATTTTCATTTGTATTGTTAGCGA‑ATGCAATAGGGCCGGAAC‑TGATCCGATGAGTTAATTTTGAAGAAATCTCATGTTTCGTGGTGGTCGTTTTTACCACAGATGGTTTTATGCCAGTATGG‑TTTGTTGAATTTTTATTAAAGCTGGGTTTAGCgttttgggat                                                          <  2:530005/250‑9 (MQ=255)                                                                                                                                                                                          cGGTGGCAACGGCGGG‑TTTTCTGCCGCGTT‑TCGCGGG‑AATGGGTGCGACCATG‑ACAGTC‑T‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTA‑ATGCAATAGTTACTGAAC‑TGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGG‑TTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAgt                                                     <  1:354116/251‑1 (MQ=255)                                                                                                                                                                                                 cactgcggG‑GTTTCTGCCGCCTT‑TCGCGGT‑AATTTGTGCGTCTATG‑ACTGTC‑T‑CCTTTTTCTGTTTCGTGGTAATCTGTTTCATGTTTAGTGGTAGCTA‑ATGCAATATTTACTGAAC‑TGAGCCGATGAGTTAATGGTGAACAAATCTCATGTTGCGTGGTGGTCGCTTGTTCCACAGTTGCGTTTATGCCAGTATGT‑TTTGTTGAATTTTTTTTAAATCTGGGTTGAGCGTGTCGGGTGCAAGTTTGagca                                              <  2:369389/247‑1 (MQ=255)                                                                                                                                                                                                       ggg‑TTTTCTGCCGCGTT‑TCGCGGG‑AATGGGTGCGACCATG‑ACAGTC‑T‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTA‑ATGCAATAGTTACTGAAC‑TGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGG‑TTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCaa                                          >  1:409139/1‑249 (MQ=255)                                                                                                                                                                                                        gg‑TTTTCTGCCGCGTT‑TCGCGGG‑AATGGGTGCGACCATG‑ACAGTC‑T‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTA‑ATGCAATAGTTACTGAAC‑TGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGG‑TTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGt                                       <  1:283284/251‑1 (MQ=255)                                                                                                                                                                                                         g‑TTTCCTGCCGCGTA‑TCGCGGG‑AACGGGTGCGACCATG‑ACAGTC‑T‑CCTTTTTCTGTATCGTGGCAGCCCTTTTCATTCTTATCGTTAGCTA‑ATGCAATAGTTACTGAAC‑TGTTACGCTGGGTTAATGTTGAACATATCTAATGCTACGTGGTGGTCGTCTTTACAACAGATGCGTTCTTGTCCGTATGG‑TTTGTTGAATTTTTATTAAATCTGGGTTGAGCTGGTCGTGAGCACGTGTGCGCAGCCAAGTg                                      >  2:259343/1‑251 (MQ=255)                                                                                                                                                                                                          tttttCCGCCGCGTT‑TCGCGGG‑AATTGGTGCGACCCTG‑ACCGTC‑T‑CCTTTTTCTGTATCGTGGAACTCATTTTCATTTTTATTGTTAGCTA‑ATGCAATAGTTACTGAAC‑TGATCCGATGAGTTAATGTTGAACAAATCTCATTTTGCGTGGTGGTCACTTTTACCACAGCTACGTTTATGCCCGTATGG‑TTTGTTGATTTTTTCTTAAATCTGGTTTAAGCTTGTCGGGAGCACGGGTGagctgccatgtg                                      >  2:259338/2‑242 (MQ=255)                                                                                                                                                                                                            tttCTGCCGGGTT‑TCGCGGG‑AATGGGTGCGACCATG‑CCAGTC‑T‑CCTTTTTCTGTATCGTGGAAATCATTTTCCTTTTTATTGTTAGCTA‑ATGCAATAGTTACTGAAC‑TGATCCGTTGAGTTAATGTTGAACAAATCTAATTTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCGGTATGG‑TTTGTTGAATTTTTCTTAACTCTGGTGTGGGCGTGTCGGGAGCAAGTGTGAGAAGCAAAgatg                                     >  2:339358/1‑247 (MQ=255)                                                                                                                                                                                                                 gCCGCGTT‑TCGCGGG‑AATGGTTGCGACCATG‑ACAGTC‑T‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTA‑ATGCAATAGTTACTGAAC‑TGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTAt                                                                                                       <  2:522845/180‑1 (MQ=255)                                                                                                                                                                                                                 gCCGCGTT‑TCGCGGG‑AATGGGTGCGACCATG‑ACAGTC‑T‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTA‑ATGCAATAGTTACTGAAC‑TGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTAt                                                                                                       >  1:522845/1‑180 (MQ=255)                                                                                                                                                                                                                    tcgTT‑TCGCGGG‑AATGGGTGCGACCATG‑ACAGTC‑T‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTA‑ATGCAATAGTTACTGAAC‑TGATCCGTTGATTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGTTGCGTTTATGCCAGTATGG‑TTTGTTGAATTTTTATTCAATCTGGGTTGTGCGTGTCGGGATCAAGTGTGAGCAGCAAACTGGACCAACttc                            >  2:173325/2‑249 (MQ=255)                                                                                                                                                                                                                    gcgTT‑TCGCAGG‑AATGGGAGCGACCACG‑ACAGTC‑T‑CCTTTTTCTTTATCGTTGACATCATATTCATTATTATTGTTAGCTA‑ATTCAATAGTTACTGAAC‑TGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTTCCACAGTTGCGTTTATGCCAGTATTG‑TTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTTTGAGCAGCAAAGTGGAACAACtgc                            <  1:292207/251‑1 (MQ=255)                                                                                                                                                                                                                        t‑tCGCGGG‑AATGGGTGCGACCATG‑ACAGTC‑T‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTA‑ATGCAATAGTTACTGAAC‑TGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGG‑TTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACtgcgtg                         <  1:502041/250‑1 (MQ=255)                                                                                                                                                                                                                              ggg‑AATGGGTGAGACCATG‑ACAGTC‑T‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTA‑ATGCAATAGTTACTGAAC‑TGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGTTGCGTTTATGCCTGTATGG‑TTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGAAGCAAAGTGGAACAACTGCGTGCGCAGt                   >  2:441446/1‑251 (MQ=255)                                                                                                                                                                                                                                       gTGCGACCATG‑ACAGTC‑T‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTA‑ATGCAATAGTTACTGAAC‑TGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGG‑TTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGa            >  1:473725/1‑250 (MQ=255)                                                                                                                                                                                                                                         gCGACCATG‑ACAGTC‑T‑CCGTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTA‑ATGCAATAGTTACTGAAC‑TGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGTTCGCTTTTACCACAGATGCGTTTATGCCAGTATGG‑TTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAAc          <  2:409139/250‑1 (MQ=255)                                                                                                                                                                                                                                          cgatcATG‑ACAGTT‑TTCTTTTTT‑TGTATCGTGGAAAGTATTTTCATTTTTATTGTTAGCTA‑ATGAAATAGTTACTGAAC‑TGATCCGATGAGTTAATGTTGATCAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGG‑TTTGTTGAATTTGTATTAAATCTGGGTTGAGCGTGTCGGGATCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTTCGCAGTTAAATGAACGt        <  2:236385/247‑1 (MQ=255)                                                                                                                                                                                                                                                 g‑ACAGTC‑T‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTA‑ATGCAATAGTTACTGAAC‑TGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGG‑TTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTg  <  1:279442/250‑1 (MQ=255)                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                    AGTTCGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTG‑CGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTT‑TCGCGGG‑AATGGGTGCGACCATG‑ACAGTC‑T‑CTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTA‑ATGCAATAGTTACTGAAC‑TGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGG‑TTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTG  >  NC_000913/4223393‑4223878

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: ‑