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

TTCGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTG  >  NC_000913/4223395‑4223881
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                     
ttCGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGAAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCAACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTCTGCCCCGTTTCCCGGGGATGGGGGCGAGCATGCCCGTCTCCCTTTTCt                                                                                                                                                                                                                                                >  1:251594/1‑251 (MQ=255)
 tcggcAACCTGACGCACGAAACCCTGCTGGTGCATCGTGGTTAGCAGGCGGTGGGTCGTTGAATTGGGTAACCCGGCTTGGTGCGCCAGTTCCGTGAGTGCCACTCTGCCAGTGGATGCGGCAATCCTCTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGTAATGGGTGCGACCATGACAGTCTCCTTTTTCt                                                                                                                                                                                                                                                <  2:590300/246‑1 (MQ=255)
 tCGCCAACCTTACGCAAGAAGCCCTGGTGTTGCTTCGTGGTGAGCAGGCGGTTGGGCGTGTAATGGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGAGTCGGGAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGATCCGGTCCAGTCGCTGGTGCGGTTGCGATGGCGGTTTTGCTGCCGCGTTTCGCGGGAATGGGTGCGTTCATGACATTCTTCTTTTTCt                                                                                                                                                                                                                                                <  2:31199/250‑1 (MQ=255)
  cGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGt                                                                                                                                                                                                                                              <  1:5670/251‑1 (MQ=255)
   gCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTTCCGCGTTTCGCGGGGAGGGGGGCGACCATGTCAGTCTCCTTTTTTTGtt                                                                                                                                                                                                                                             >  1:20271/1‑250 (MQ=255)
    ccAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTAt                                                                                                                                                                                                                                            <  2:506243/251‑1 (MQ=255)
        ccTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTg                                                                                                                                                                                                                                        <  1:127016/251‑1 (MQ=255)
           gACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGaa                                                                                                                                                                                                                                     <  1:352779/251‑1 (MQ=255)
             cGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAAt                                                                                                                                                                                                                                   <  1:258972/251‑1 (MQ=255)
               cACGAAACCGTGCTGTGGCAGCGTGGTTAGCGGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAAGTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAAtca                                                                                                                                                                                                                                 <  2:83811/251‑1 (MQ=255)
               cACGAAACCCTTATGTTGCACCGGGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAAtca                                                                                                                                                                                                                                 <  2:95680/251‑1 (MQ=255)
                     aCCCTGCTGGTGCATCGTGGTTAGCGGGAGGTGGGTCGTGGAAGTGGGTAACCCGGCTTGTTGCGCCAGGTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACGCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGGCTCCTTTTTCTGTATCGTGGAAAtcattttca                                                                                                                                                                                                                           <  2:201774/251‑1 (MQ=255)
                         ggctgTTGCAGCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATGCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATcattttcatttt                                                                                                                                                                                                                       <  2:198253/250‑1 (MQ=255)
                               tGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACTGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTg                                                                                                                                                                                                                  <  2:623590/250‑1 (MQ=255)
                                  aTCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCCCGTTTCGCGGGAGAGGGGGCGGCCATGACAGACTTCTTTTTTTGTTTCCTGGAAATTCTTtttatttttattttttg                                                                                                                                                                                                              >  1:154994/1‑246 (MQ=255)
                                    cGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCATTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGc                                                                                                                                                                                                             <  2:352421/250‑1 (MQ=255)
                                          tAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCCGTAATTTCAGGCCACGCGTTAAAGCCTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGTTTCGCGGGGAAGGGGGCGCCCATGACAGTCTCCCTTTTTTTTTTCGTGGAGATCATTTTTTTTTTTATTTTTAGCTAATGc                                                                                                                                                                                                       >  1:73772/1‑250 (MQ=255)
                                           aGCAGGCGGTGGGGCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTGTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCaa                                                                                                                                                                                                     <  2:234807/251‑1 (MQ=255)
                                                gCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACCCTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGtt                                                                                                                                                                                                <  2:449877/251‑1 (MQ=255)
                                                                              ttgttgCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTCTGCCGCGGTTCGCGGGGATGGGGGCGGCCCTGAAAGTCCCCTTTTTCTGTATCGTGGAAAACATTTTTTTTTTTATTTTTTGCTTATGCAATAGTTTCTGAAATTGTCCGAGGAGTTTATGTTgt                                                                                                                                                                   >  2:589892/1‑249 (MQ=255)
                                                                               tgttgCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGa                                                                                                                                                                   <  1:408768/249‑1 (MQ=255)
                                                                                tgtgCGCCAGTTCCGTGAGTGCCACACTGCCAGTGGATTCGTCAATCCAGTCCAGTAATTGCAGGCCACTCGTTAAAGACTGAACCTGGCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGTGTTTCGCGGGAATGGGTGCTACCATGACAGTCGCCTTTTTCTGTATCGGGGAAAACATTTTCATATTTATTTTTAGCTAATGCAATAGTTAGTGAACTGATCCGATGGGTGAATGGTGatc                                                                                                                                                                 <  2:587718/248‑3 (MQ=255)
                                                                                              gTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCg                                                                                                                                                  <  1:437512/251‑1 (MQ=255)
                                                                                                gAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCgt                                                                                                                                                 <  1:185157/250‑1 (MQ=255)
                                                                                                 atggCCACACTGCCATGTGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGTGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATTTTGAACAAATCTCATGTGGCgtgg                                                                                                                                               <  2:280485/248‑1 (MQ=255)
                                                                                                                           aaTCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGttt                                                                                                                   <  2:471554/251‑1 (MQ=255)
                                                                                                                             tCCACTCCAGTCATTTCAGGCCACGCGTTAAAGCCTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTTTGCCGCGCTTCGCGGGGATGGGTGCGACCATGAAAGTGTCCTTTTTTTGTATCGTGGAAAACATTTTTATTTTTCTTGTTAGCTAATGCAATAGTTACTGAAATGATCCGATGGGGTAATTTTTAAAAAATATCATGTTGCGAGGTGGTCGCGTTTTAcaacgaagcggttat                                                                                                                  >  2:607643/1‑237 (MQ=255)
                                                                                                                              ccACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGTTTCGCGGGAATGGGTGCGACCATGAAAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTTATTTTTTTTGTTTGCTAATGCAATAGTTTCTGAACTGATCCGGTGAGTTTATGTTTGACAAAACTCATGTTGCGTGGTGGTCGCTTTTTCCACAGAttcggttttg                                                                                                                >  1:496278/1‑242 (MQ=255)
                                                                                                                                   ccAGTAATTTCAGGCCACGCGTTAAAGACTGATCCTGTCCAGTCGCTGGTTGGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTGATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGcct                                                                                                             <  2:553499/249‑2 (MQ=255)
                                                                                                                                   ccAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTTTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATTCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGAGGCGTTTATGCCAg                                                                                                            <  2:20774/250‑1 (MQ=255)
                                                                                                                                   ccAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGt                                                                                                                                                                                                                                              <  1:250294/122‑1 (MQ=255)
                                                                                                                                   ccAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGt                                                                                                           <  1:240660/251‑1 (MQ=255)
                                                                                                                                            tCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTttg                                                                                                   <  1:589892/250‑1 (MQ=255)
                                                                                                                                                     gcgTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTCTGCCCCGGTTCGCGGGGAAGGGGGGGGCCATTGCAGTCTCCTTTTTTTGTATCGTGGGAATCATTTTTATTTTTATTGTTAGCTAATGCCATAGTTACTGAACTGATCCGATGAGTTAATGTTTAACAAATTTCCTGTTGCGTTGGGGTCGCTTTTTCCACCGATGCGTTTATGCCAGTATGGTTTTTTGAAttttt                                                                                         >  1:208927/1‑251 (MQ=255)
                                                                                                                                                       gTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCAGGACAGTCTCCTTGTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTttat                                                                                       <  1:184992/251‑1 (MQ=255)
                                                                                                                                                       gTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGCACTTTCTTCCGCGTTTCGCGGGAAAGGGGGCGACCATGACACGCTCCCTTTTCTGTATCGTGGGAATCAATTTCATTTTTTTTTTTAACTAATCCAATAGTTACTTAAATGCTCCCATGAGTTTATGTTTAACAAATCTCATGTTTCGTGGTGGTGGCTTTTACCTCACCAGCGTTTATGCCAGTATGGTTTGTTTAATTTttat                                                                                       >  1:270629/1‑251 (MQ=255)
                                                                                                                                                            aGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGTTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTATCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAAt                                                                                  <  2:51680/251‑1 (MQ=255)
                                                                                                                                                             gACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTAAGCCAGTATGGTTTGTTGAATTTTTATTAAATc                                                                                 <  2:118532/251‑1 (MQ=255)
                                                                                                                                                                     cTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTTTTCCGGGTTTTGCGGGAATGGGGGCGAACACGTCAGTCTCCCTTTTTTGTATTGTGGGAATCATTTTTATTTTTATTGTTTGCTACTGCAATAGTTACTGAAACGATCCGCTGCGTTAAAGTTGAAAAAATCTCATGTTTCATGGGGGTCGCTTTTTCCCACAATGCGGTTATTTCAGTATTGTTTGTTGAATTTTTTTTAAATTTGGGTTgg                                                                         >  1:30800/1‑250 (MQ=255)
                                                                                                                                                                         ccAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCg                                                                      <  1:226838/250‑1 (MQ=255)
                                                                                                                                                                         ccAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTCTGCCGCGTTTCGCGGGGAAGGGGGCGACCATGGCAGTCTCCTTTTTCTGGCTCGTGGGAATCATTTTTATTTTTATTTTTTGCTTATGCAATAGTTTCTTAAATGATCCGATGAGTTAAAGTTTAACAAATTTCCTGTTGCGTGGGGGTGGCCTTTTCCACAGATGCGTTTATGCCAGTCTGGGTTGTTGAATTTTTTTTAAATTTGGGGTGAGCgg                                                                     >  1:405490/1‑250 (MQ=255)
                                                                                                                                                                                 gggTGCGGGGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAg                                                             <  2:508211/250‑1 (MQ=255)
                                                                                                                                                                                                    ggTTTACTCCCGCGTTCCGCGGGAATGGGGGCGACCATGACCGCCCCCTCTTTCTGTATCGTGGCAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGCTGAGTTCATGCTGAACAACTCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGATTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGt                                          >  2:314059/1‑251 (MQ=255)
                                                                                                                                                                                                       tttCTGCCGCGTTTCAAGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGACCTGATCCGATGGGTTAATGTTGAACAAATGTCATGTTGCGTGGTGGTCGATTTTACCACAGATGCGTTTATGCAAGTATGGTTTGTTGAATTTTTATAAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAgagg                                        >  1:218270/1‑247 (MQ=255)
                                                                                                                                                                                                        ttCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGa                                       <  2:595812/250‑1 (MQ=255)
                                                                                                                                                                                                                cgTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTAGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAg                                                                        >  2:279281/1‑209 (MQ=255)
                                                                                                                                                                                                                cgTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTAGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAg                                                                        <  1:279281/209‑1 (MQ=255)
                                                                                                                                                                                                                   ttCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGTTTGCTTTTTTCTGGAGTAAGTGTGAGCAGCAAAGTGGAACAACtgcgtt                            <  1:355152/250‑2 (MQ=255)
                                                                                                                                                                                                                           aaTGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCTTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTaa                   <  1:228791/251‑1 (MQ=255)
                                                                                                                                                                                                                                    gACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTa          <  2:496278/251‑1 (MQ=255)
                                                                                                                                                                                                                                     aCCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTAt         <  1:470364/251‑1 (MQ=255)
                                                                                                                                                                                                                                      ccATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTAtt        >  2:131820/1‑251 (MQ=255)
                                                                                                                                                                                                                                          gACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTgg    >  2:233485/1‑251 (MQ=255)
                                                                                                                                                                                                                                           aCAGTCTCTTTTTTCTGTATCGTGAAAATCATTTTCATATTTATGGAAAGAGAAGGCAATAGTTAAGGAAGGGAACAGACGAGTTAAAGATGAACAAATCTCAAGTTGAGTGGAGGTCGCAAAAAACACAGATGCGAATATGCCAGTATGGAGGGAGGAATGATAAAAAAAGCTGGGACGAGCGTGACGGGAGCAAGTGCGAGCAGCAAAGTGGAACAACTGCGTGCGCAGATAAAGGAACGTATTCTGGt   >  1:424980/1‑251 (MQ=255)
                                                                                                                                                                                                                                           aCAGTCTCATTTTTCTGTATAGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGt   >  2:219976/1‑251 (MQ=255)
                                                                                                                                                                                                                                            cAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTg  <  2:594166/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                     
TTCGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTG  >  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: