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_0009134,223,6380TC98.5% 214.1 / ‑4.1 68intergenic (‑10/‑190)iclR/metHDNA‑binding transcriptional repressor IclR/cobalamin‑dependent methionine synthase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (21/46);  total (22/46)
Fisher's exact test for biased strand distribution p-value = 3.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.31e-01

AGTTCGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGG‑AATGGGTGCGACC‑ATGACAGTCT‑CTTTTTTCTGTATCGTGG‑AAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTT‑ACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGT  >  NC_000913/4223393‑4223880
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    
agtgcGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGG‑AATGGGTGCGACC‑ATGACAGTCT‑Ccttttt                                                                                                                                                                                                                                                 <  2:85369/247‑1 (MQ=255)
aGTTCGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGGTTCGCGGG‑GATGGGGGGGGCC‑ATGACAGTCT‑Ccttttt                                                                                                                                                                                                                                                 >  1:383012/1‑251 (MQ=255)
    cGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGG‑AATGGGTGCGACC‑ATGACAGTCT‑CCTTTTTCTg                                                                                                                                                                                                                                              <  1:422901/250‑1 (MQ=255)
    cGCCAAACTGACGCACGAAACCCTGCTGTTGCATCGCGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCAGTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGCTTTCTGCCGCGTTTCGCGGG‑AATGGGTGCGACC‑ATGACAGTCT‑CCTTTTTCTGt                                                                                                                                                                                                                                             <  2:168135/251‑1 (MQ=255)
      ccAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACCCTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGCTTTGCGGG‑GAAGGGGGGGACC‑AAGAGAGTCT‑CCTTTTTTTGTAt                                                                                                                                                                                                                                           >  2:678593/1‑251 (MQ=255)
      ccAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGTTTCGCCGG‑AATGGGGGCGACC‑ATGACAGTCA‑CCTTTTTTTGTAt                                                                                                                                                                                                                                           >  1:76513/1‑251 (MQ=255)
          ccggACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGGGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACTCTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGTCAACGGCGGGTTTTCTGCCGCGTTTCGCGGG‑AATGGGTGCGACC‑ATGACAGTCT‑CCTTTTTCTGGATCGTg                                                                                                                                                                                                                                       <  1:482704/248‑1 (MQ=255)
          ccTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGG‑AATGGGTGCGACC‑ATGACAGTCT‑CCTTTTTCTGTATCGTg                                                                                                                                                                                                                                       <  2:257949/251‑1 (MQ=255)
           cTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTCTGCCGCGTTTCGCGGG‑GATGGGTGCGGCC‑ATGACAGTCT‑CCTTTTTCTGTATCGTgg                                                                                                                                                                                                                                      >  1:411194/1‑251 (MQ=255)
             gACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATCCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGG‑AATGGGTGCGACC‑ATGACAGTCT‑CCTTTTTCTGTATCGTgg                                                                                                                                                                                                                                      <  1:453508/249‑1 (MQ=255)
             gACGCACGAAACCCTGCTGGTGCAGCGTGGTGAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGGTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCTCTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAGCCTGTCCAGTCGCTGTTGCGGTGGCATCGGCGTGTTTTCTGCCGCGGGTCGCGGT‑AATGGGTGCGACC‑ATTACATTCT‑CCTTTCTCTGGAGCGTGG‑aa                                                                                                                                                                                                                                   <  2:483339/251‑1 (MQ=255)
                        cccTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGG‑AATGGGTGCGACC‑ATGACAGTCT‑CCTTTTTCTGTATCGGGG‑AAAtcagtttcat                                                                                                                                                                                                                        <  2:67032/251‑1 (MQ=255)
                            gcttttcaatagTGGTTAACAGGCGGTGGGTCGCGGAATTAGGTATCCCGGCTTGTTGCTCCAGTCCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCGGGCCACTCGTGAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCATCGGCGGTTTTTCTGCCGCGTTTGGCGGG‑AATTGGTGCGATC‑GTGACAGTCT‑CCTTTTTCTGTATCGTGG‑AAATCATTGTCAttttt                                                                                                                                                                                                                    <  2:631754/240‑1 (MQ=255)
                            gctgTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGG‑AATGGGTGCGACC‑ATGACAGTCT‑CCTTTTTCTGTATCGTGG‑AAAtcattttcattt                                                                                                                                                                                                                      <  2:207351/249‑1 (MQ=255)
                               gTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGG‑AATGGGTGCGACC‑ATGACAGTCT‑CCTTTTTCTGTATCGTGG‑AAATCATTTTCATTTTTAtt                                                                                                                                                                                                                 <  1:276263/251‑1 (MQ=255)
                                      cGTGGTTAGCAGGCGGTGGGTCGTGTAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGG‑AATGGGTGCGACC‑ATGACAGTCT‑CCTTTTTCTGTATCGTGG‑AAATCATTTTCATTTTTATTGTTAGCt                                                                                                                                                                                                          <  1:442270/251‑1 (MQ=255)
                                         ggTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGTTTCGCGGG‑GATGGGGGCGGCC‑ATGACAGGCT‑CCTTTTTTTGTATCGTGG‑AAATCATTTTTATTTTTTTTTTTTGCTAAt                                                                                                                                                                                                       >  1:402425/1‑251 (MQ=255)
                                         ggTTAGCAGGCGGTGGGTCGTGGAATGGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGG‑AATGGGTGCGACC‑ATGACAGTCT‑CCTTTTTCTGTATCGTGG‑AAATCATTTTCATTTTTATTGTTAGCTAAt                                                                                                                                                                                                       <  2:383633/251‑1 (MQ=255)
                                            tAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGATTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGG‑AATGGGTGCGACC‑ATGACAGTCT‑CCTTTTTCTGTATCGTGG‑AAATCATTTTCATTTTTATTGTTAGCTAATGc                                                                                                                                                                                                     <  2:413511/250‑1 (MQ=255)
                                             aGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGTCGCGTTTTGCGGG‑AGTGGGGGCGGCC‑AAGACAGTCT‑CCTTTTTTTGTATCGTGG‑GAAtcatttttctttttttttttttcttattcct                                                                                                                                                                                                   >  2:468457/1‑227 (MQ=255)
                                                       gggTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTTTGCCGCGTTTCGCGGG‑GAAGGGGGCGGCC‑ATGACCGGCT‑CCTTTTTTTTTTTTGTGGGAAACCTTTTTC‑TTTTTTTTTTTAGCTAATGCTATAGtttttt                                                                                                                                                                                          >  2:337176/1‑246 (MQ=255)
                                                        ggTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGG‑AATGGGTGCGACC‑ATGACAGTCT‑CCTTTTTCTGTATCGTGG‑AAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTactgaa                                                                                                                                                                                        <  1:270008/251‑1 (MQ=255)
                                                              ggAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCAAGTACTTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTTTGCCCCGTTTCGCGGG‑AAAGGGGGCGACC‑CTGCCAGTCT‑CCTTTTTCTGTATCGTTG‑AAATCCTTTTCATTTTTTTTGTTAGCTACAGCAATAGTTACTGCACTGATc                                                                                                                                                                                  >  1:725039/1‑251 (MQ=255)
                                                                aaTTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTTTGCCGCGGTTCGCGGG‑GATGGGGGCGACC‑ACGACAGACT‑CCTTTTTTTTTTTCGTGG‑GAATCATTTTTCTTTTTTTTTTTAGCTAATGCAATAGTTTCTGATATGATccc                                                                                                                                                                                >  2:609771/1‑250 (MQ=255)
                                                                     ggTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCCGTTCGCGGG‑GATGGGGGCGGCC‑CTGGCAGTCT‑CCCTTTTCTGTTTAGTGG‑GAATCATTTTTTCTTTTATTTTTTTCTTATTCAATAGTTacttaacttgtccgcatgg                                                                                                                                                                           >  1:452283/1‑240 (MQ=255)
                                                                       tAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGCCTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGATTTCTGCCCCGTTTCGCGGG‑GATGGGGGCGACC‑CCGACAGGCT‑CCTTTTTCTGTATTGTGG‑AGATCATTTTCCTTTTTATTGTTTTCTAATGCCATAATTTCTGTAATGATCCGATGGGtt                                                                                                                                                                         >  2:330363/1‑251 (MQ=255)
                                                                       tAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTCTTCCGCGTTTCGCGGG‑AAAGGGTGCGGCC‑CTGACAGGCT‑CCTTTTTTTTTATTGTGG‑AAATCATTTTCATTTTTTTTTTTAGCTTATTCAATAGTTAATGAACTGATCCGATGGGtt                                                                                                                                                                         >  1:579557/1‑251 (MQ=255)
                                                                         aCCCGTCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGG‑AATGGGTGCGACC‑ATGACAGTCT‑CCTTTTTCTGTATCGTGG‑AAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTaa                                                                                                                                                                       <  2:860/251‑1 (MQ=255)
                                                                           ccGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGG‑AATGGGTGCGACC‑ATGACAGTCT‑CCTTTTTCTGTATCGTGG‑AAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATg                                                                                                                                                                     <  1:408481/251‑1 (MQ=255)
                                                                           ccGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGG‑AATGGGTGCGACC‑ATGACAGTCT‑CCTTTTTCTGTATCGTGG‑AAATCAGTTTCATTTTTATTGTTAGCTAATGCAATATTTACTGAACTGATCCGAGGAGTTAATg                                                                                                                                                                     <  2:477488/251‑1 (MQ=255)
                                                                               cTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGG‑AATGGGTGCGACC‑ATGACAGTCT‑CCTTTT‑CTGTATCGTGG‑AAATCATTTTCATTTTTATTGTTAGCTAATGc                                                                                                                                                                                                     <  1:58778/214‑1 (MQ=255)
                                                                                 tgttgCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGG‑AATGGGTGCGACC‑ATGACAGTCT‑CCTTTTTCTGTATCGTGG‑AAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTCGaa                                                                                                                                                                <  1:141872/250‑1 (MQ=255)
                                                                                      cgcCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGG‑AATGGGTGCGACC‑ATGACAGTCT‑CCTTTTTCTGTATCGTGG‑AAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTGGAACAAAtc                                                                                                                                                          <  2:116867/251‑1 (MQ=255)
                                                                                         cAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGG‑AATGGGTGCGACC‑ATGACAGTCT‑CCTTTTTCTGTATCGTGG‑AAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGAt                                                                                                                                                                              <  1:25273/228‑1 (MQ=255)
                                                                                         cAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTCTGCCGCGTTTTGCGGG‑GATGGGGGGGGCC‑AAGACAGTCT‑CCTTTTTTTGTATCGTGG‑GAATCATTTTTCTTTTTTTTGTTTGCTAATGCAATAGTTTCTGAACTTATCCggtagatcggaagagcgtttggtagg                                                                                                                                                       >  2:25273/1‑226 (MQ=255)
                                                                                              ccGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTTGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGG‑AATGGGTGCGACC‑ATGACAGTCT‑CCTTTTTCTGTATCGTGG‑AAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGtt                                                                                                                                                   <  1:678593/250‑1 (MQ=255)
                                                                                              ccGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTCTGGCGCGGTTTGCGGG‑GATGGGGGGGAGC‑ATGACAGTCT‑CCTTTTTCTGTATCGGGG‑GAATCATTTTCATTTTTTTTTTTAGATAATGCAATAGTTACTGAACTGATTCGATTAGGTAATGTTTGACAAATCTCATGTTg                                                                                                                                                  >  2:206528/1‑251 (MQ=255)
                                                                                              ccGTGAGTGCCACACTGACATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGG‑AATGGGTGCGACC‑ATGACAGTCT‑CCTTTTTCTGTATCGTGG‑AAATCATTTTCATTTGTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTg                                                                                                                                                  <  2:119353/251‑1 (MQ=255)
                                                                                                 tGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTCTGGCGCGTTTCGCGGG‑GATGGGGGGGGCC‑AAGACAGTCT‑CCATTTTCTGTTTCGTGG‑GAATAATTTTCATTTTTTTTGTTTGCTTATTCGATGATTTCTGAAATCATCCGCTGAGTTTATGTTTAAAAAATCTCATTTTgtgg                                                                                                                                               >  2:125078/1‑248 (MQ=255)
                                                                                                  gAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGGCCTGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGT‑AATGGGTGCGACC‑ATGTCAGTCT‑CCTTTTGCTGTATCGTGG‑AAATCATTGTCATTTTTATTGTTAGCTAATGCAATAGTGACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCgtg                                                                                                                                              <  1:635706/251‑1 (MQ=255)
                                                                                                      gCCAGACTGCCAGTGGATTCGGCAAGCCACTCCAGTAATTTCAGGCCACGCGTGAAAGACTGAACCTGTCCAGTCGCTGGTGTGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGG‑AATGGGTGCGACC‑ATGACAGTCT‑CCTTTTCCTGGATCGTGG‑AGATGATTTTCAGTTTTATTGTTATCTAATGCAATATGTACGGAACTGATCCGAGGAGTTAATTTTGAACAAATCGCATTTTGGTtggtgg                                                                                                                                          <  2:294941/251‑1 (MQ=255)
                                                                                                        cacaCTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTTCCGCGTTTCGCGGG‑AATGGGTGCGACC‑ATGACAGTCT‑CCTTTTTCTGTATCGTGG‑AAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTc                                                                                                                                        <  1:705480/251‑1 (MQ=255)
                                                                                                              gCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTTTGCCGCGTTTCGCGGG‑GATGGGTGCGGCC‑ATGACAGTCT‑CCTTTTTTTGTATCGTGG‑GAATCATTTTTATTTTTTTTGTTAGCTAATGCCATAGTTACTGAACTGGTCCGATGAGTTAATGTTGAAAAAATCTCATGTTGCGTGGTGGGCGCcttt                                                                                                                                  >  1:359558/1‑251 (MQ=255)
                                                                                                                                     ccAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGG‑AATGGGTGCGACC‑ATGACAGTCT‑CCTTTTTCTGTATCGTGG‑AAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTT‑ACCACAGATGCGTTTATGCCa                                                                                                            <  1:468457/249‑1 (MQ=255)
                                                                                                                                     ccAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGG‑AATGGGTGCGACC‑ATGACAGTCT‑CCTTTTTCTGTATCGTGG‑AAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTT‑ACCACAGATGCGTTTATGCCAGt                                                                                                          <  1:441358/251‑1 (MQ=255)
                                                                                                                                      cAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCTGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGG‑AATGTGTGCGACC‑ATGACAGTCT‑CCTTTTTCTGTATCGTGG‑AAATCATTTTCATTTTTATTGGTAGCTAAGGCAATAGTTACTGAACTGAGCCGATGAGTTAATGTTGAACAAAGCTCATGTTGCGTGGTGGTCGCTTTT‑ACCACAGATGCGTTTATGCCAGTa                                                                                                         <  1:98471/251‑1 (MQ=255)
                                                                                                                                       aGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCGGGTGCGGGGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGG‑AATGGGTGCGACC‑ATGACAGTCT‑CCTTTTTCTGTATCGTGG‑AAATCATTTTCATTTTTATTGTGAGCTAATGCAATAGTTACTGAACTGATCTGATGAGTGAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTT‑ACCACAGATGCGTTTATGCCAGTAt                                                                                                        <  2:238310/251‑1 (MQ=255)
                                                                                                                                            tttCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTTCCGCGTTTCGCGGG‑GATGGGTGCGACC‑ATGACAGGCT‑CCTTTTTTTGTTTCGTGG‑GAATCATTTTCATTTTTTTTGTTAGCTAATGCAATAGTTTCTGAACTGATCCGCTGAGTTAATTTTGGACAAATCTCATGTTTCGTGGGGGTTGCTTTT‑TCCACAGATGCGTTTTTTCCCGTATGGttt                                                                                                   >  2:362030/1‑251 (MQ=255)
                                                                                                                                               cAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCAGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGG‑AATGGGTGCGACC‑ATGACAGTCT‑CCTTTTTCTGTATCGTGG‑AAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTT‑ACCACAGATGCGTTTATGCCAGTATGGTttgt                                                                                                 <  1:647630/250‑1 (MQ=255)
                                                                                                                                               cAGGCCACGCGTTAAAGACTGAACCTGTCCAGGCGCTGGTGCGGGGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGT‑AATGGGTGCGACC‑ATGACAGTCT‑CCTTTTTCTGTATCGTGG‑AAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAATAAATCTCATGTTGCGTGGTGGTCGTTTTT‑ACCACAGATGCGTTTATGCCAGTATGgtgtgt                                                                                                 <  2:616836/250‑1 (MQ=255)
                                                                                                                                                  gCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGACCCGGTTTGCGGG‑GATCGGTGGGGCC‑ATGGCCATCT‑CCTTTTTCTTTTTCGTGG‑AAATTCTTTTCATTTTTTTTGTTTGCTAATGCAATTGTTTCTTAAATAATCCGATGAGTTAATGTTGTACAAATCTCATGTTTCGTGGTGGTCGCTTTT‑ACCCCAGATGAGTTTATGCCAGTATGGTttgttta                                                                                              >  2:314127/1‑248 (MQ=255)
                                                                                                                                                            aaaGACTTAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTTGCGGG‑AATGGGTGCGACC‑ATGACAGTCT‑CCTTTTTCTGTATCGTGG‑AAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTT‑ACCACAGAGGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTaa                                                                                   <  1:604736/251‑1 (MQ=255)
                                                                                                                                                               gACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGG‑AATGGGTGCGACC‑ATGACAGTCT‑CCTTTTTCTGTATCGTGG‑AAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTT‑ACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATc                                                                                <  1:314127/251‑1 (MQ=255)
                                                                                                                                                                      ccTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGG‑AATGGGTGCGACC‑ATGACAGTCT‑CCTTTTTCTGTATCGTGG‑AAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTT‑ACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTg                                                                         <  1:228468/251‑1 (MQ=255)
                                                                                                                                                                             aGTCGCTGGTGCGGTGGCAACGGCGGGTTTTTTGCCGCGATTCGCGGGGAACGGG‑GCGACCGAAG‑CAGTCTTCTTTTTT‑TTTATCGGGG‑GAATCATTTTCATTTTTTTTGTTTGCTAATGCAATAGTTACTGGAATGATTCGGTGCGTTAATGTTGAACAAATCTCATGGTGAGTGGGGGTTGCTTTTTACAAC‑GATGCGGTTATGCCAGAATGGTTTGTTGGATTTTTTTTTAATATGGGtttggagggt                                                                  >  2:735326/1‑243 (MQ=255)
                                                                                                                                                                             aGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGG‑AATGGGTGCGACC‑ATGACAGTCT‑CCTTTTTCTGTATCGTGG‑AAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTT‑ACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCgtgt                                                                  <  2:645674/251‑1 (MQ=255)
                                                                                                                                                                                   tGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGATTTGCGGG‑GATGGGGGCGGCC‑CTGACAGTCT‑CCTTTTTTTGTATCGTGG‑AAATCATTTTTATTTTTTTTTTTAGCTTATGCAATAGGTACTTAACTGATCCCATTAGTTAATTTTTAAAAAATATCATGTTTCGTGGGGGGCGCTTTT‑ACCACAGATTCGTTTATGCCAGTATGGTTTTTTGAATTTTTTTTAAATCTTGGTTGAGAGTGTTggggg                                                            >  1:147290/1‑249 (MQ=255)
                                                                                                                                                                                               aaCGGCGGGTTTTCTGCCGCGTTTCGCGGT‑AATGGGTGCGACC‑ATGACAGTCT‑CCTTTTTCTGTATCGTGG‑AAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATTTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTT‑ACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGGCGGGAGCAAGTGTGagca                                                <  1:217297/251‑1 (MQ=255)
                                                                                                                                                                                                aCGGCGGGTTTTCTGCCGCGTTTCGCGGG‑AATGGGTGCGACC‑ATGACAGTCT‑CCTTTTTCTGTATCGTGG‑AAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTT‑ACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGagcag                                               <  2:411194/251‑1 (MQ=255)
                                                                                                                                                                                                      ggTTTTCTGCCGCGTTTCGCGGG‑AATGGGTGCGACC‑ATGACAGTCT‑CCTTTTTCTGTATCGTGG‑AAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTT‑ACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGt                                         <  2:402425/251‑1 (MQ=255)
                                                                                                                                                                                                      ggTTTTCTGCCGCGTTTCGCGGG‑AATGGGTGCGACC‑ATGACAGTCT‑CCTTTTTCTGTATCGTGG‑AAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTT‑ACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGt                                         <  1:337176/251‑1 (MQ=255)
                                                                                                                                                                                                          ttCTGCCGCGTTTCGCGGG‑AATGGGTGCGACC‑ATGACAGTCT‑CCTTTTTCTGTATCGTGG‑AAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTT‑ACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGa                                      <  1:503427/250‑1 (MQ=255)
                                                                                                                                                                                                           tCTGCCGCGTTTCGCGGG‑AATGGGTGCGACC‑ATGACAGTCT‑CCTTTTTCTGTATCGTGG‑AAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTT‑ACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGaac                                    >  1:622864/1‑251 (MQ=255)
                                                                                                                                                                                                           tCTGCCGCGTTTCGCGGG‑AATGGGTGCGACC‑ATGACAGTCT‑CCTTTGTCTGTATCGTGG‑AAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTT‑ACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGATCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGaa                                     <  2:140508/250‑1 (MQ=255)
                                                                                                                                                                                                               ccGCGTTTCGCGGG‑AATGGGTGCGACC‑ATGACATTCT‑CCTGTTTCTGTATCGTGG‑AAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTGGCGTGGTGGTCGCTTTT‑ACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACt                                <  2:89049/251‑1 (MQ=255)
                                                                                                                                                                                                                          ggg‑AATGGGTGCGACC‑ATGACAGTCT‑CCTTTTTCTGTATCGTGG‑GAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTT‑ACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGt                     >  2:587662/1‑251 (MQ=255)
                                                                                                                                                                                                                                  ggTGCGACC‑ATGACAGTCT‑CCTTTTTCTGTATCGTGG‑AAATCATTTTCATTTTTATTGTTAGCTACTGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTT‑ACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTaattg               >  2:169342/1‑247 (MQ=255)
                                                                                                                                                                                                                                          c‑ATGACAGTCT‑CCTTTTTCTGTATCGTGG‑AAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTT‑ACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTc      <  1:206528/251‑1 (MQ=255)
                                                                                                                                                                                                                                              gACAGTCT‑CCTTTTTCTGTATCGTGG‑AAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTT‑ACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTgg   <  2:579557/251‑1 (MQ=255)
                                                                                                                                                                                                                                               aCAGTCT‑CCTTTTTCTGTATCGTGG‑AAATCATGTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTT‑ACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGt  <  2:452283/251‑1 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    
AGTTCGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGG‑AATGGGTGCGACC‑ATGACAGTCT‑CTTTTTTCTGTATCGTGG‑AAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTT‑ACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGT  >  NC_000913/4223393‑4223880

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