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

AGTTCGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGT  >  NC_000913/4223393‑4223880
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    
agtgcGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCcttttt                                                                                                                                                                                                                                                 <  2:85369/247‑1 (MQ=255)
aGTTCGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGGTTCGCGGGGATGGGGGGGGCCATGACAGTCTCcttttt                                                                                                                                                                                                                                                 >  1:383012/1‑251 (MQ=255)
    cGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTg                                                                                                                                                                                                                                              <  1:422901/250‑1 (MQ=255)
    cGCCAAACTGACGCACGAAACCCTGCTGTTGCATCGCGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCAGTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGCTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGt                                                                                                                                                                                                                                             <  2:168135/251‑1 (MQ=255)
      ccAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACCCTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGCTTTGCGGGGAAGGGGGGGACCAAGAGAGTCTCCTTTTTTTGTAt                                                                                                                                                                                                                                           >  2:678593/1‑251 (MQ=255)
      ccAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGTTTCGCCGGAATGGGGGCGACCATGACAGTCACCTTTTTTTGTAt                                                                                                                                                                                                                                           >  1:76513/1‑251 (MQ=255)
          ccggACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGGGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACTCTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGTCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGGATCGTg                                                                                                                                                                                                                                       <  1:482704/248‑1 (MQ=255)
          ccTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTg                                                                                                                                                                                                                                       <  2:257949/251‑1 (MQ=255)
           cTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTCTGCCGCGTTTCGCGGGGATGGGTGCGGCCATGACAGTCTCCTTTTTCTGTATCGTgg                                                                                                                                                                                                                                      >  1:411194/1‑251 (MQ=255)
             gACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATCCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTgg                                                                                                                                                                                                                                      <  1:453508/249‑1 (MQ=255)
             gACGCACGAAACCCTGCTGGTGCAGCGTGGTGAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGGTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCTCTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAGCCTGTCCAGTCGCTGTTGCGGTGGCATCGGCGTGTTTTCTGCCGCGGGTCGCGGTAATGGGTGCGACCATTACATTCTCCTTTCTCTGGAGCGTGGaa                                                                                                                                                                                                                                   <  2:483339/251‑1 (MQ=255)
                        cccTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGGGGAAAtcagtttcat                                                                                                                                                                                                                        <  2:67032/251‑1 (MQ=255)
                            gcttttcaatagTGGTTAACAGGCGGTGGGTCGCGGAATTAGGTATCCCGGCTTGTTGCTCCAGTCCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCGGGCCACTCGTGAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCATCGGCGGTTTTTCTGCCGCGTTTGGCGGGAATTGGTGCGATCGTGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTGTCAttttt                                                                                                                                                                                                                    <  2:631754/240‑1 (MQ=255)
                            gctgTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAAtcattttcattt                                                                                                                                                                                                                      <  2:207351/249‑1 (MQ=255)
                               gTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTAtt                                                                                                                                                                                                                 <  1:276263/251‑1 (MQ=255)
                                      cGTGGTTAGCAGGCGGTGGGTCGTGTAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCt                                                                                                                                                                                                          <  1:442270/251‑1 (MQ=255)
                                         ggTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGTTTCGCGGGGATGGGGGCGGCCATGACAGGCTCCTTTTTTTGTATCGTGGAAATCATTTTTATTTTTTTTTTTTGCTAAt                                                                                                                                                                                                       >  1:402425/1‑251 (MQ=255)
                                         ggTTAGCAGGCGGTGGGTCGTGGAATGGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAAt                                                                                                                                                                                                       <  2:383633/251‑1 (MQ=255)
                                            tAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGATTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGc                                                                                                                                                                                                     <  2:413511/250‑1 (MQ=255)
                                             aGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGTCGCGTTTTGCGGGAGTGGGGGCGGCCAAGACAGTCTCCTTTTTTTGTATCGTGGGAAtcatttttctttttttttttttcttattcct                                                                                                                                                                                                   >  2:468457/1‑227 (MQ=255)
                                                       gggTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTTTGCCGCGTTTCGCGGGGAAGGGGGCGGCCATGACCGGCTCCTTTTTTTTTTTTGTGGGAAACCTTTTTCTTTTTTTTTTTAGCTAATGCTATAGtttttt                                                                                                                                                                                          >  2:337176/1‑246 (MQ=255)
                                                        ggTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTactgaa                                                                                                                                                                                        <  1:270008/251‑1 (MQ=255)
                                                              ggAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCAAGTACTTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTTTGCCCCGTTTCGCGGGAAAGGGGGCGACCCTGCCAGTCTCCTTTTTCTGTATCGTTGAAATCCTTTTCATTTTTTTTGTTAGCTACAGCAATAGTTACTGCACTGATc                                                                                                                                                                                  >  1:725039/1‑251 (MQ=255)
                                                                aaTTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTTTGCCGCGGTTCGCGGGGATGGGGGCGACCACGACAGACTCCTTTTTTTTTTTCGTGGGAATCATTTTTCTTTTTTTTTTTAGCTAATGCAATAGTTTCTGATATGATccc                                                                                                                                                                                >  2:609771/1‑250 (MQ=255)
                                                                     ggTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCCGTTCGCGGGGATGGGGGCGGCCCTGGCAGTCTCCCTTTTCTGTTTAGTGGGAATCATTTTTTCTTTTATTTTTTTCTTATTCAATAGTTacttaacttgtccgcatgg                                                                                                                                                                           >  1:452283/1‑240 (MQ=255)
                                                                       tAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGCCTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGATTTCTGCCCCGTTTCGCGGGGATGGGGGCGACCCCGACAGGCTCCTTTTTCTGTATTGTGGAGATCATTTTCCTTTTTATTGTTTTCTAATGCCATAATTTCTGTAATGATCCGATGGGtt                                                                                                                                                                         >  2:330363/1‑251 (MQ=255)
                                                                       tAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTCTTCCGCGTTTCGCGGGAAAGGGTGCGGCCCTGACAGGCTCCTTTTTTTTTATTGTGGAAATCATTTTCATTTTTTTTTTTAGCTTATTCAATAGTTAATGAACTGATCCGATGGGtt                                                                                                                                                                         >  1:579557/1‑251 (MQ=255)
                                                                         aCCCGTCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTaa                                                                                                                                                                       <  2:860/251‑1 (MQ=255)
                                                                           ccGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATg                                                                                                                                                                     <  1:408481/251‑1 (MQ=255)
                                                                           ccGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCAGTTTCATTTTTATTGTTAGCTAATGCAATATTTACTGAACTGATCCGAGGAGTTAATg                                                                                                                                                                     <  2:477488/251‑1 (MQ=255)
                                                                               cTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGc                                                                                                                                                                                                     <  1:58778/214‑1 (MQ=255)
                                                                                 tgttgCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTCGaa                                                                                                                                                                <  1:141872/250‑1 (MQ=255)
                                                                                      cgcCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTGGAACAAAtc                                                                                                                                                          <  2:116867/251‑1 (MQ=255)
                                                                                         cAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGAt                                                                                                                                                                              <  1:25273/228‑1 (MQ=255)
                                                                                         cAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTCTGCCGCGTTTTGCGGGGATGGGGGGGGCCAAGACAGTCTCCTTTTTTTGTATCGTGGGAATCATTTTTCTTTTTTTTGTTTGCTAATGCAATAGTTTCTGAACTTATCCggtagatcggaagagcgtttggtagg                                                                                                                                                       >  2:25273/1‑226 (MQ=255)
                                                                                              ccGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTTGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGtt                                                                                                                                                   <  1:678593/250‑1 (MQ=255)
                                                                                              ccGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTCTGGCGCGGTTTGCGGGGATGGGGGGGAGCATGACAGTCTCCTTTTTCTGTATCGGGGGAATCATTTTCATTTTTTTTTTTAGATAATGCAATAGTTACTGAACTGATTCGATTAGGTAATGTTTGACAAATCTCATGTTg                                                                                                                                                  >  2:206528/1‑251 (MQ=255)
                                                                                              ccGTGAGTGCCACACTGACATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTGTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTg                                                                                                                                                  <  2:119353/251‑1 (MQ=255)
                                                                                                 tGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTCTGGCGCGTTTCGCGGGGATGGGGGGGGCCAAGACAGTCTCCATTTTCTGTTTCGTGGGAATAATTTTCATTTTTTTTGTTTGCTTATTCGATGATTTCTGAAATCATCCGCTGAGTTTATGTTTAAAAAATCTCATTTTgtgg                                                                                                                                               >  2:125078/1‑248 (MQ=255)
                                                                                                  gAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGGCCTGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGTAATGGGTGCGACCATGTCAGTCTCCTTTTGCTGTATCGTGGAAATCATTGTCATTTTTATTGTTAGCTAATGCAATAGTGACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCgtg                                                                                                                                              <  1:635706/251‑1 (MQ=255)
                                                                                                      gCCAGACTGCCAGTGGATTCGGCAAGCCACTCCAGTAATTTCAGGCCACGCGTGAAAGACTGAACCTGTCCAGTCGCTGGTGTGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTCCTGGATCGTGGAGATGATTTTCAGTTTTATTGTTATCTAATGCAATATGTACGGAACTGATCCGAGGAGTTAATTTTGAACAAATCGCATTTTGGTtggtgg                                                                                                                                          <  2:294941/251‑1 (MQ=255)
                                                                                                        cacaCTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTTCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTc                                                                                                                                        <  1:705480/251‑1 (MQ=255)
                                                                                                              gCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTTTGCCGCGTTTCGCGGGGATGGGTGCGGCCATGACAGTCTCCTTTTTTTGTATCGTGGGAATCATTTTTATTTTTTTTGTTAGCTAATGCCATAGTTACTGAACTGGTCCGATGAGTTAATGTTGAAAAAATCTCATGTTGCGTGGTGGGCGCcttt                                                                                                                                  >  1:359558/1‑251 (MQ=255)
                                                                                                                                     ccAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCa                                                                                                            <  1:468457/249‑1 (MQ=255)
                                                                                                                                     ccAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGt                                                                                                          <  1:441358/251‑1 (MQ=255)
                                                                                                                                      cAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCTGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGTGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGGTAGCTAAGGCAATAGTTACTGAACTGAGCCGATGAGTTAATGTTGAACAAAGCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTa                                                                                                         <  1:98471/251‑1 (MQ=255)
                                                                                                                                       aGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCGGGTGCGGGGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTGAGCTAATGCAATAGTTACTGAACTGATCTGATGAGTGAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTAt                                                                                                        <  2:238310/251‑1 (MQ=255)
                                                                                                                                            tttCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTTCCGCGTTTCGCGGGGATGGGTGCGACCATGACAGGCTCCTTTTTTTGTTTCGTGGGAATCATTTTCATTTTTTTTGTTAGCTAATGCAATAGTTTCTGAACTGATCCGCTGAGTTAATTTTGGACAAATCTCATGTTTCGTGGGGGTTGCTTTTTCCACAGATGCGTTTTTTCCCGTATGGttt                                                                                                   >  2:362030/1‑251 (MQ=255)
                                                                                                                                               cAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCAGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTttgt                                                                                                 <  1:647630/250‑1 (MQ=255)
                                                                                                                                               cAGGCCACGCGTTAAAGACTGAACCTGTCCAGGCGCTGGTGCGGGGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGTAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAATAAATCTCATGTTGCGTGGTGGTCGTTTTTACCACAGATGCGTTTATGCCAGTATGgtgtgt                                                                                                 <  2:616836/250‑1 (MQ=255)
                                                                                                                                                  gCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGACCCGGTTTGCGGGGATCGGTGGGGCCATGGCCATCTCCTTTTTCTTTTTCGTGGAAATTCTTTTCATTTTTTTTGTTTGCTAATGCAATTGTTTCTTAAATAATCCGATGAGTTAATGTTGTACAAATCTCATGTTTCGTGGTGGTCGCTTTTACCCCAGATGAGTTTATGCCAGTATGGTttgttta                                                                                              >  2:314127/1‑248 (MQ=255)
                                                                                                                                                            aaaGACTTAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTTGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGAGGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTaa                                                                                   <  1:604736/251‑1 (MQ=255)
                                                                                                                                                               gACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATc                                                                                <  1:314127/251‑1 (MQ=255)
                                                                                                                                                                      ccTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTg                                                                         <  1:228468/251‑1 (MQ=255)
                                                                                                                                                                             aGTCGCTGGTGCGGTGGCAACGGCGGGTTTTTTGCCGCGATTCGCGGGGAACGGGGCGACCGAAGCAGTCTTCTTTTTTTTTATCGGGGGAATCATTTTCATTTTTTTTGTTTGCTAATGCAATAGTTACTGGAATGATTCGGTGCGTTAATGTTGAACAAATCTCATGGTGAGTGGGGGTTGCTTTTTACAACGATGCGGTTATGCCAGAATGGTTTGTTGGATTTTTTTTTAATATGGGtttggagggt                                                                  >  2:735326/1‑243 (MQ=255)
                                                                                                                                                                             aGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCgtgt                                                                  <  2:645674/251‑1 (MQ=255)
                                                                                                                                                                                   tGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGATTTGCGGGGATGGGGGCGGCCCTGACAGTCTCCTTTTTTTGTATCGTGGAAATCATTTTTATTTTTTTTTTTAGCTTATGCAATAGGTACTTAACTGATCCCATTAGTTAATTTTTAAAAAATATCATGTTTCGTGGGGGGCGCTTTTACCACAGATTCGTTTATGCCAGTATGGTTTTTTGAATTTTTTTTAAATCTTGGTTGAGAGTGTTggggg                                                            >  1:147290/1‑249 (MQ=255)
                                                                                                                                                                                               aaCGGCGGGTTTTCTGCCGCGTTTCGCGGTAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATTTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGGCGGGAGCAAGTGTGagca                                                <  1:217297/251‑1 (MQ=255)
                                                                                                                                                                                                aCGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGagcag                                               <  2:411194/251‑1 (MQ=255)
                                                                                                                                                                                                      ggTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGt                                         <  2:402425/251‑1 (MQ=255)
                                                                                                                                                                                                      ggTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGt                                         <  1:337176/251‑1 (MQ=255)
                                                                                                                                                                                                          ttCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGa                                      <  1:503427/250‑1 (MQ=255)
                                                                                                                                                                                                           tCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGaac                                    >  1:622864/1‑251 (MQ=255)
                                                                                                                                                                                                           tCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTGTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGATCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGaa                                     <  2:140508/250‑1 (MQ=255)
                                                                                                                                                                                                               ccGCGTTTCGCGGGAATGGGTGCGACCATGACATTCTCCTGTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTGGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACt                                <  2:89049/251‑1 (MQ=255)
                                                                                                                                                                                                                          gggAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGGAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGt                     >  2:587662/1‑251 (MQ=255)
                                                                                                                                                                                                                                  ggTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTACTGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTaattg               >  2:169342/1‑247 (MQ=255)
                                                                                                                                                                                                                                          cATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTc      <  1:206528/251‑1 (MQ=255)
                                                                                                                                                                                                                                              gACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTgg   <  2:579557/251‑1 (MQ=255)
                                                                                                                                                                                                                                               aCAGTCTCCTTTTTCTGTATCGTGGAAATCATGTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGt  <  2:452283/251‑1 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    
AGTTCGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGT  >  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: