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,6380TC96.4% 167.7 / ‑4.0 56intergenic (‑10/‑190)iclR/metHDNA‑binding transcriptional repressor IclR/cobalamin‑dependent methionine synthase
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base C (15/38);  total (18/38)
Fisher's exact test for biased strand distribution p-value = 9.16e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCT‑CTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGC‑TTTT‑ACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTG  >  NC_000913/4223395‑4223881
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                     
ttCGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGAAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCAACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTCTGCCCCGTTTCCCGGGGATGGGGGCGAGCATGCCCGTCT‑CCCTTTTCt                                                                                                                                                                                                                                                >  1:251594/1‑251 (MQ=255)
 tcggcAACCTGACGCACGAAACCCTGCTGGTGCATCGTGGTTAGCAGGCGGTGGGTCGTTGAATTGGGTAACCCGGCTTGGTGCGCCAGTTCCGTGAGTGCCACTCTGCCAGTGGATGCGGCAATCCTCTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGTAATGGGTGCGACCATGACAGTCT‑CCTTTTTCt                                                                                                                                                                                                                                                <  2:590300/246‑1 (MQ=255)
 tCGCCAACCTTACGCAAGAAGCCCTGGTGTTGCTTCGTGGTGAGCAGGCGGTTGGGCGTGTAATGGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGAGTCGGGAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGATCCGGTCCAGTCGCTGGTGCGGTTGCGATGGCGGTTTTGCTGCCGCGTTTCGCGGGAATGGGTGCGTTCATGACATTCT‑TCTTTTTCt                                                                                                                                                                                                                                                <  2:31199/250‑1 (MQ=255)
  cGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCT‑CCTTTTTCTGt                                                                                                                                                                                                                                              <  1:5670/251‑1 (MQ=255)
   gCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTTCCGCGTTTCGCGGGGAGGGGGGCGACCATGTCAGTCT‑CCTTTTTTTGtt                                                                                                                                                                                                                                             >  1:20271/1‑250 (MQ=255)
    ccAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCT‑CCTTTTTCTGTAt                                                                                                                                                                                                                                            <  2:506243/251‑1 (MQ=255)
        ccTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCT‑CCTTTTTCTGTATCGTg                                                                                                                                                                                                                                        <  1:127016/251‑1 (MQ=255)
           gACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCT‑CCTTTTTCTGTATCGTGGaa                                                                                                                                                                                                                                     <  1:352779/251‑1 (MQ=255)
             cGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCT‑CCTTTTTCTGTATCGTGGAAAt                                                                                                                                                                                                                                   <  1:258972/251‑1 (MQ=255)
               cACGAAACCGTGCTGTGGCAGCGTGGTTAGCGGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAAGTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCT‑CCTTTTTCTGTATCGTGGAAAtca                                                                                                                                                                                                                                 <  2:83811/251‑1 (MQ=255)
               cACGAAACCCTTATGTTGCACCGGGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCT‑CCTTTTTCTGTATCGTGGAAAtca                                                                                                                                                                                                                                 <  2:95680/251‑1 (MQ=255)
                     aCCCTGCTGGTGCATCGTGGTTAGCGGGAGGTGGGTCGTGGAAGTGGGTAACCCGGCTTGTTGCGCCAGGTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACGCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGGCT‑CCTTTTTCTGTATCGTGGAAAtcattttca                                                                                                                                                                                                                           <  2:201774/251‑1 (MQ=255)
                         ggctgTTGCAGCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATGCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCT‑CCTTTTTCTGTATCGTGGAAATcattttcatttt                                                                                                                                                                                                                       <  2:198253/250‑1 (MQ=255)
                               tGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACTGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCT‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTg                                                                                                                                                                                                                  <  2:623590/250‑1 (MQ=255)
                                  aTCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCCCGTTTCGCGGGAGAGGGGGCGGCCATGACAGACTTCTTTTTT‑TGTTTCCTGGAAATTCTTtttatttttattttttg                                                                                                                                                                                                              >  1:154994/1‑246 (MQ=255)
                                    cGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCATTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCT‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGc                                                                                                                                                                                                             <  2:352421/250‑1 (MQ=255)
                                          tAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCCGTAATTTCAGGCCACGCGTTAAAGCCTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGTTTCGCGGGGAAGGGGGCGCCCATGACAGTCT‑CCCTTTTTTTTTTCGTGGAGATCATTTTTTTTTTTATTTTTAGCTAATGc                                                                                                                                                                                                       >  1:73772/1‑250 (MQ=255)
                                           aGCAGGCGGTGGGGCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCT‑CCTTGTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCaa                                                                                                                                                                                                     <  2:234807/251‑1 (MQ=255)
                                                gCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACCCTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCT‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGtt                                                                                                                                                                                                <  2:449877/251‑1 (MQ=255)
                                                                              ttgttgCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTCTGCCGCGGTTCGCGGGGATGGGGGCGGCCCTGAAAGTCC‑CCTTTTTCTGTATCGTGGAAAACATTTTTTTTTTTATTTTTTGCTTATGCAATAGTTTCTGAAATTGTCCGAGGAGTTTATGTTgt                                                                                                                                                                   >  2:589892/1‑249 (MQ=255)
                                                                               tgttgCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCT‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGa                                                                                                                                                                   <  1:408768/249‑1 (MQ=255)
                                                                                tgtgCGCCAGTTCCGTGAGTGCCACACTGCCAGTGGATTCGTCAATCCAGTCCAGTAATTGCAGGCCACTCGTTAAAGACTGAACCTGGCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGTGTTTCGCGGGAATGGGTGCTACCATGACAGTCG‑CCTTTTTCTGTATCGGGGAAAACATTTTCATATTTATTTTTAGCTAATGCAATAGTTAGTGAACTGATCCGATGGGTGAATGGTGatc                                                                                                                                                                 <  2:587718/248‑3 (MQ=255)
                                                                                              gTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCT‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCg                                                                                                                                                  <  1:437512/251‑1 (MQ=255)
                                                                                                gAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCT‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCgt                                                                                                                                                 <  1:185157/250‑1 (MQ=255)
                                                                                                 atggCCACACTGCCATGTGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGTGAATGGGTGCGACCATGACAGTCT‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATTTTGAACAAATCTCATGTGGCgtgg                                                                                                                                               <  2:280485/248‑1 (MQ=255)
                                                                                                                           aaTCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCT‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGC‑TTTT‑ACCACAGATGCGttt                                                                                                                   <  2:471554/251‑1 (MQ=255)
                                                                                                                             tCCACTCCAGTCATTTCAGGCCACGCGTTAAAGCCTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTTTGCCGCGCTTCGCGGGGATGGGTGCGACCATGAAAGTGT‑CCTTTTTTTGTATCGTGGAAAACATTTTTATTTTTCTTGTTAGCTAATGCAATAGTTACTGAAATGATCCGATGGGGTAATTTTTAAAAAATATCATGTTGCGAGGTGGTCGCGTTTT‑Acaacgaagcggttat                                                                                                                  >  2:607643/1‑237 (MQ=255)
                                                                                                                              ccACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGTTTCGCGGGAATGGGTGCGACCATGAAAGTCT‑CCTTTTTCTGTATCGTGGAAATCATTTTTATTTTTTTTGTTTGCTAATGCAATAGTTTCTGAACTGATCCGGTGAGTTTATGTTTGACAAAACTCATGTTGCGTGGTGGTCGC‑TTTT‑TCCACAGAttcggttttg                                                                                                                >  1:496278/1‑242 (MQ=255)
                                                                                                                                   ccAGTAATTTCAGGCCACGCGTTAAAGACTGATCCTGTCCAGTCGCTGGTTGGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCT‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTGATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGC‑TTTT‑ACCACAGATGCGTTTATGcct                                                                                                             <  2:553499/249‑2 (MQ=255)
                                                                                                                                   ccAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTTTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCT‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATTCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGC‑TTTT‑ACCACAGAGGCGTTTATGCCAg                                                                                                            <  2:20774/250‑1 (MQ=255)
                                                                                                                                   ccAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCT‑CCTTTTTCTGt                                                                                                                                                                                                                                              <  1:250294/122‑1 (MQ=255)
                                                                                                                                   ccAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCT‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGC‑TTTT‑ACCACAGATGCGTTTATGCCAGt                                                                                                           <  1:240660/251‑1 (MQ=255)
                                                                                                                                            tCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCT‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGC‑TTTT‑ACCACAGATGCGTTTATGCCAGTATGGTttg                                                                                                   <  1:589892/250‑1 (MQ=255)
                                                                                                                                                     gcgTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTCTGCCCCGGTTCGCGGGGAAGGGGGGGGCCATTGCAGTCT‑CCTTTTTTTGTATCGTGGGAATCATTTTTATTTTTATTGTTAGCTAATGCCATAGTTACTGAACTGATCCGATGAGTTAATGTTTAACAAATTTCCTGTTGCGTTGGGGTCGC‑TTTT‑TCCACCGATGCGTTTATGCCAGTATGGTTTTTTGAAttttt                                                                                         >  1:208927/1‑251 (MQ=255)
                                                                                                                                                       gTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCAGGACAGTCT‑CCTTGTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGC‑TTTT‑ACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTttat                                                                                       <  1:184992/251‑1 (MQ=255)
                                                                                                                                                       gTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGCACTTTCTTCCGCGTTTCGCGGGAAAGGGGGCGACCATGACACGCT‑CCCTTTTCTGTATCGTGGGAATCAATTTCATTTTTTTTTTTAACTAATCCAATAGTTACTTAAATGCTCCCATGAGTTTATGTTTAACAAATCTCATGTTTCGTGGTGGTGGC‑TTTT‑ACCTCACCAGCGTTTATGCCAGTATGGTTTGTTTAATTTttat                                                                                       >  1:270629/1‑251 (MQ=255)
                                                                                                                                                            aGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGTTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCT‑CCTTTATCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGC‑TTTT‑ACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAAt                                                                                  <  2:51680/251‑1 (MQ=255)
                                                                                                                                                             gACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCT‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGC‑TTTT‑ACCACAGATGCGTTTAAGCCAGTATGGTTTGTTGAATTTTTATTAAATc                                                                                 <  2:118532/251‑1 (MQ=255)
                                                                                                                                                                     cTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTTTTCCGGGTTTTGCGGGAATGGGGGCGAACACGTCAGTCT‑CCCTTTTTTGTATTGTGGGAATCATTTTTATTTTTATTGTTTGCTACTGCAATAGTTACTGAAACGATCCGCTGCGTTAAAGTTGAAAAAATCTCATGTTTCATGGGGGTCGC‑TTTTTCCCACA‑ATGCGGTTATTTCAGTATTGTTTGTTGAATTTTTTTTAAATTTGGGTTgg                                                                         >  1:30800/1‑250 (MQ=255)
                                                                                                                                                                         ccAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCT‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGC‑TTTT‑ACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCg                                                                      <  1:226838/250‑1 (MQ=255)
                                                                                                                                                                         ccAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTCTGCCGCGTTTCGCGGGGAAGGGGGCGACCATGGCAGTCT‑CCTTTTTCTGGCTCGTGGGAATCATTTTTATTTTTATTTTTTGCTTATGCAATAGTTTCTTAAATGATCCGATGAGTTAAAGTTTAACAAATTTCCTGTTGCGTGGGGGTGGC‑CTTT‑TCCACAGATGCGTTTATGCCAGTCTGGGTTGTTGAATTTTTTTTAAATTTGGGGTGAGCgg                                                                     >  1:405490/1‑250 (MQ=255)
                                                                                                                                                                                 gggTGCGGGGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCT‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGC‑TTTT‑ACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAg                                                             <  2:508211/250‑1 (MQ=255)
                                                                                                                                                                                                    ggTTTACTCCCGCGTTCCGCGGGAATGGGGGCGACCATGACCGCCC‑CCTCTTTCTGTATCGTGGCAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGCTGAGTTCATGCTGAACAACTCTCATGTTGCGTGGTGGTCGC‑TTTT‑ACCACAGATGCGATTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGt                                          >  2:314059/1‑251 (MQ=255)
                                                                                                                                                                                                       tttCTGCCGCGTTTCAAGGGAATGGGTGCGACCATGACAGTCT‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGACCTGATCCGATGGGTTAATGTTGAACAAATGTCATGTTGCGTGGTGGTCGA‑TTTT‑ACCACAGATGCGTTTATGCAAGTATGGTTTGTTGAATTTTTATAAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAgagg                                        >  1:218270/1‑247 (MQ=255)
                                                                                                                                                                                                        ttCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCT‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGC‑TTTT‑ACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGa                                       <  2:595812/250‑1 (MQ=255)
                                                                                                                                                                                                                cgTTTCGCGGGAATGGGTGCGACCATGACAGTCT‑CCTTTTTCTGTATCGTAGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGC‑TTTT‑ACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAg                                                                        >  2:279281/1‑209 (MQ=255)
                                                                                                                                                                                                                cgTTTCGCGGGAATGGGTGCGACCATGACAGTCT‑CCTTTTTCTGTATCGTAGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGC‑TTTT‑ACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAg                                                                        <  1:279281/209‑1 (MQ=255)
                                                                                                                                                                                                                   ttCGCGGGAATGGGTGCGACCATGACAGTCT‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGC‑TTTT‑ACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGTTTGCTTTTTTCTGGAGTAAGTGTGAGCAGCAAAGTGGAACAACtgcgtt                            <  1:355152/250‑2 (MQ=255)
                                                                                                                                                                                                                           aaTGGGTGCGACCATGACAGTCT‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCTTGGTGGTCGC‑TTTT‑ACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTaa                   <  1:228791/251‑1 (MQ=255)
                                                                                                                                                                                                                                    gACCATGACAGTCT‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGC‑TTTT‑ACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTa          <  2:496278/251‑1 (MQ=255)
                                                                                                                                                                                                                                     aCCATGACAGTCT‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGC‑TTTT‑ACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTAt         <  1:470364/251‑1 (MQ=255)
                                                                                                                                                                                                                                      ccATGACAGTCT‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGC‑TTTT‑ACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTAtt        >  2:131820/1‑251 (MQ=255)
                                                                                                                                                                                                                                          gACAGTCT‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGC‑TTTT‑ACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTgg    >  2:233485/1‑251 (MQ=255)
                                                                                                                                                                                                                                           aCAGTCT‑CTTTTTTCTGTATCGTGAAAATCATTTTCATATTTATGGAAAGAGAAGGCAATAGTTAAGGAAGGGAACAGACGAGTTAAAGATGAACAAATCTCAAGTTGAGTGGAGGTCGC‑AAAA‑AACACAGATGCGAATATGCCAGTATGGAGGGAGGAATGATAAAAAAAGCTGGGACGAGCGTGACGGGAGCAAGTGCGAGCAGCAAAGTGGAACAACTGCGTGCGCAGATAAAGGAACGTATTCTGGt   >  1:424980/1‑251 (MQ=255)
                                                                                                                                                                                                                                           aCAGTCT‑CATTTTTCTGTATAGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGC‑TTTT‑ACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGt   >  2:219976/1‑251 (MQ=255)
                                                                                                                                                                                                                                            cAGTCT‑CCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGC‑TTTT‑ACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTg  <  2:594166/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                     
TTCGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCT‑CTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGC‑TTTT‑ACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTG  >  NC_000913/4223395‑4223881

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