Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 677,829 T→C W139R (TGG→CGG)  djlB → J domain‑containing protein DjlB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913677,8290TC92.0% 63.9 / ‑1.9 25W139R (TGG→CGG) djlBJ domain‑containing protein DjlB
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (4/19);  total (5/20)
Fisher's exact test for biased strand distribution p-value = 3.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

AGAATATGAGGTAGCAGAACATGAAATTCTGCTCGCCTTTCGTGCGTTACTTGCCTCTGATAGTGAACGTTTTCTGCCCTCCGCCTGGCAGCGATTCATTCAGCAATTAAATTATTGCTCGATGGAGGAGATTGATGAATTACGCTGGTCGCTGTGCACAATCGCCATGAACACTGCCCATTTATCCTTCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCGGTGGTT‑ACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCGAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGATTTTTACCAAATGTTCGCTCGGATTTGGTCATCGCATCCACAATGGCTGACATTGTATTTAGCGCAACATCGCGCAGTGATTATCCCCGATGATGCAAAACTGCACA  >  NC_000913/677606‑678072
                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                    
agaaTATGAGGTAGCAGAACATGAAATTCTGCTCGCCTTTCGTGCGTTACTTGCCTCTGATAGTGAACGTTTTCTGCCCTCCGCCTGGCAGCGATTCATTCAGCAATTAAATTATTGCTCGATGGAGGAGATTGATGAATTACGCTGGTCGCTGTGCACAATCGCCATGAACACTGCCCATTTATCCTTCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCggcggtagatatggggaaataaacgggga                                                                                                                                                                                                                            >  1:131731/1‑227 (MQ=255)
       gAGGTAGCAGAACATGAAATTCTGCTCGCCTTTCGTGCGTTACTTGCCTCTGATAGTGAACGTTTTCTGCCCTCCGCCTGGCAGCGATTCATTCAGCAATTAAATTATTGCTCGATGGAGGAGATTGATGAATTACGCTGGTCGCTGTGCACAATCGCCATGAACACTGCCCATTTATCCTTCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCGGCGGTT‑TTAGGAGGGAAAAAACGGGGGAATAGa                                                                                                                                                                                                                      >  1:370987/1‑248 (MQ=255)
          gTAGCAGAACATGAAATTCTGCTCGCCTTTCGTGCGTTACTTGCCTCTGATAGTGAACGTTTTCTGCCCTCCGCCTGGCAGCGATTCATTCAGCAATTAAATTATTGCTCGATGGAGGAGATTGATGAATTACGCTGGTCGCTGTGCACAATCGCCATGAACACTGCCCATTTATCCTTCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCGGCGGTT‑TTAGGAGGAGAATACCGGGGGAAAAGACga                                                                                                                                                                                                                   >  1:290648/1‑248 (MQ=255)
              cAGAACATGAAATTCTGCTCGCCTTTCGTGCGTTACTTGCCTCTGATAGTGAACGTTTTCTGCCCTCCGCCTGGCAGCGATTCATTCAGCAATTAAATTATTGCTCGATGGAGGAGATTGATGAATTACGCTGGTCGCTGTGCACAATCGCCATGAACACTGCCCATTTATCCTTCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCGGCGGTTTACAGGGGGaaaaaccaggggaataaaaggggaaga                                                                                                                                                                                                              >  1:234690/1‑227 (MQ=255)
              cAGAACATGAAATTCTGCTCGCCTTTCGTGCGTTACTTGCCTCTGATAGTGAACGTTTTCTGCCCTCCGCCTGGCAGCGATTCATTCAGCAATTAAATTATTGCTCGATGGAGGAGATTGATGAATTACGCTGGTCGCTGTGCACAATCGCCATGAACACTGCCCATTTATCCTTCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCGGCGGTT‑ACAGGAGGAAAATACCGGGGAAATAGACgaagaagaa                                                                                                                                                                                                            <  2:290648/251‑1 (MQ=255)
                     tGAAATTCTGCTCGCCTTTCGTGCGTTACTTGCCTCTGATAGTGAACGTTTTCTGCCCTCCGCCTGGCAGCGATTCATTCAGCAATTAAATTATTGCTCGATGGAGGAGATTGATGAATTACGCTGGTCGCTGTGCACAATCGCCATGAACACTGCCCATTTATCCTTCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCGGCGGTT‑ACAGGGGGGAAAGAACGGGGGAATAAACggagaagaaaaggg                                                                                                                                                                                                       >  1:240249/1‑244 (MQ=255)
                      taaaTTCTGCTCGCCTTTCGTGCGTTACTTGCCTCTGATAGTGAACGTTTTCTGCCCTCCGCCTGGCAGCGATTCATTCAGCAATTAAATTATTGCTCGATGGAGGAGATTGATGAATTACGCTGGTCGCTGTGCACAATCGCCATGAACACTGCCCATTTATCCTTCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCGGCGGTT‑ACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAAt                                                                                                                                                                                                     <  1:115223/249‑1 (MQ=255)
                      gAAATTCTGCTCGCCTTTCGTGCGTTACTTGCCGCTGATAGTGAACGTTTTCTGCCCTCCGCCTGGCAGCGATTCATTCAGCAATTAAATTATTGCTCGATGGAGGAGATTGATGAATTACGCTGGTCGCTGTGCACAATCGCCATGAACACTGCCCATTTATCCTTCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCGGTGGTT‑ACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAAt                                                                                                                                                                                                     <  2:48636/250‑1 (MQ=255)
                       aaaTTCTGCTCGCCTTTCGTGCGTTACTTGCCTCTGATAGTGAACGTTTTCTGCCCTCCGCCTGGCAGCGATTCATTCAGCAATTAAATTATTGCTCGATGGAGGAGATTGATGAATTACGCTGGTCGCTGTGCACAATCGCCATGAACACTGCCCATTTATCCTTCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCGGCGGTT‑ACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATcc                                                                                                                                                                                                   <  1:445554/251‑1 (MQ=255)
                                 cGCCTTTCGTGCGTTACTTGCCTCTGATAGTGAACGTTTTCTGCCCTCCGCCTGGCAGCGATTCATTCAGCAATTAAATTATTGCTCGATGGAGGAGATTGATGAATTACGCTGGTCGCTGTGCACAATCGCCATGAACACTGCCCATTTATCCTTCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCGGCGGTT‑ACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTtata                                                                                                                                                                                           <  1:297403/249‑1 (MQ=255)
                                        cGTGCTTTACTTGCCTCTGATAGTGAACGTTTTCTGCCCTCCGCCTGGCAGCGATTCATTCAGCAATTAAATTATTGCTCGATGGAGGAGATTGATGAATTACGCTGGTCGCTGTGCACAATCGCCATGAACACTGCCCATTTATCCTTCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCGGCGGTT‑ACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGc                                                                                                                                                                                  <  2:232912/251‑1 (MQ=255)
                                           gCGTTACTTGCCTCTGATAGTGAACGTTTTCTGCCCTCCGCCTGGCAGCGATTCATTCAGCAATTAAATTATTGCTCGATGGAGGAGATTGATGAATTACGCTGGTCGCTGTGCACAATCGCCATGAACACTGCCCATTTATCCTTCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCGGCGGTT‑ACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCGa                                                                                                                                                                                <  1:145457/250‑1 (MQ=255)
                                                   tGCCTCTGATAGTGAACGTTTTCTGCCCTCCGCCTGGCAGCGATTCATGCAGCAATTAAATTATTGCTCGATGGAGGAGATTGATGAATTACGCTGGTCGCTGTGCACAATCGCCATGAACACTGCCCATTTATCCTTCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCGGCGGTT‑ACATGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCGAAGGGGAAt                                                                                                                                                                        <  2:240249/250‑1 (MQ=255)
                                                           aTAGTGAACGTTTTCTGCCCTCCGCCTGGCAGCGATTCATTCAGCAATTAAATTATTGCTCGATGGAGGAGATTGATGAATTACGCTGGTCGCTGTGCACAATCGCCATGAACACTGCCCATTTATCCTTCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCGGCGGTT‑ACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCGAAGGGGAATGTTTTTAAc                                                                                                                                                               <  2:26836/251‑1 (MQ=255)
                                                                  aCGTTTTCTGCCCTCCGCCTGGCAGCGATTCATTCAGCAATTAGATTATTGCTCGATGGAGGAGATTGATGAATTACGATGGTCGCTGTGCACAATCGCCATGAACACTGCCCATTTATCCTTCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCGGTGGTT‑ACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCATTTTTATATGCCATTGCGaa                                                                                                                                                                               >  1:433301/1‑228 (MQ=255)
                                                                                                       cAATTAAATTATTGCTCGATGGAGGAGATTGATGAATTACGCTGGTCGCTGTGCACAATCGCCATGAACACTGCCCATTTATCCTTCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCGGCGGTT‑ACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCGAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGacac                                                                                                                   <  1:91747/251‑1 (MQ=255)
                                                                                                                                      atgaagtgcGCTGGTCGCTGTGCACAATCGCCATGAACACTGACCAGTGAACCTTCGAGTGCGTGGTGTTATTAACAGAAAGATGGCTGCGGGT‑ACAGGAGGAAAAGACCGGGGTAATAGACGAAGAAGAACTGGAATACGTTTTATATGCCATTGCGAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCGGTTGACGTGCAAAATGACACCATCGACTTTTACCCACTGTTCGCTCGGAtt                                                                                    <  2:300436/243‑1 (MQ=255)
                                                                                                                                                                                  cATTTATCCTTCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCGGCGGTT‑ACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCGAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGATTTTTACCAAATGTTCGCTCGGATTTGGTCATCGCATCCACAATGGCTGACATTGTATTTAGCGCAACa                                        <  1:350680/251‑1 (MQ=255)
                                                                                                                                                                                            tCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCGGCGGTT‑ACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCGAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGATTTTTACCAAATGTTCGCTCGGATTTGGTCATCGCATCCACAATGGCTGACATTGTATTTAGCGCAACATCGCGCAGt                               <  1:17645/250‑1 (MQ=255)
                                                                                                                                                                                             cGAGTTCGTGGTGTTATTAGCAGAAAGATTGCGGCGGTT‑ACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCGAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGATTTTTACCAAATGTTCGCTCGGATTTGGTCATCGCATCCACAATGGCTGACATTGTATTTAGCGCAACATCGCGCAGTGa                             <  2:485567/251‑1 (MQ=255)
                                                                                                                                                                                             cGAGTGCGTGGTGTTATTAGCAGAAAGATTGCGGCGGTT‑ACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCGAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGATTTTTACCAAATGTTCGCTCGGATTTGGTCATCGCATCCACAATGGCTGACATTGTATTTAGCGCAACATCGCGCAGTGa                             <  2:234690/251‑1 (MQ=255)
                                                                                                                                                                                                gTGCGTGGTGTTATTAGCAGAAAGATTGCGGCGGTT‑ACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCGAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGATTTTTACCAAATGTTCGCTCGGATTTGGTCATc                                                                            <  2:154722/201‑1 (MQ=255)
                                                                                                                                                                                                      ggtgTTATTAGCAGAAAGATTGCGGCGGTT‑ACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCGAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGATTTTTACCAAATGTTCGCTCGGATTTGGTCATCGCATCCACAATGGCTGACATTGTATTTAGCGCAACATCGCGCAGTGATTATCCCCg                    <  2:59655/251‑1 (MQ=255)
                                                                                                                                                                                                               aGCAGAAAGATTGCGGCGGTT‑ACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCGAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGATTTTTACCAAATGTTCGCTCGGATTTGGTCATCGCATCCACAATGGCTGACATTGTATTTAGCGCAACATCGCGCAGTGATTATCCCCGATGATGCa            <  1:228507/250‑1 (MQ=255)
                                                                                                                                                                                                                       gATTGCGGCGGTT‑ACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCGAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGATTTTTACCAAATGTTCGCTCGGATTTGGTCATCGCATCCACAATGGCTGACATTGTATTTAGCGCAACATCGCGCAGTGATTATCCCCGATGATGCAAAACTGcac   <  1:412984/251‑1 (MQ=255)
                                                                                                                                                                                                                        agtGCGGCGGTT‑ACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCGAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCAGCGATTTTTACCAAATGTTCGCTCGGATTTGGTCATCGCATCCACAATGGCTGACATTGTATTTAGCGCAACATCGCGCAGTGATTATCCCCGATGATGCAAAACTGcaca  <  2:348186/249‑1 (MQ=255)
                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                    
AGAATATGAGGTAGCAGAACATGAAATTCTGCTCGCCTTTCGTGCGTTACTTGCCTCTGATAGTGAACGTTTTCTGCCCTCCGCCTGGCAGCGATTCATTCAGCAATTAAATTATTGCTCGATGGAGGAGATTGATGAATTACGCTGGTCGCTGTGCACAATCGCCATGAACACTGCCCATTTATCCTTCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCGGTGGTT‑ACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCGAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGATTTTTACCAAATGTTCGCTCGGATTTGGTCATCGCATCCACAATGGCTGACATTGTATTTAGCGCAACATCGCGCAGTGATTATCCCCGATGATGCAAAACTGCACA  >  NC_000913/677606‑678072

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