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,8290TC80.8% 39.6 / 6.0 26W139R (TGG→CGG) djlBJ domain‑containing protein DjlB
Reads supporting (aligned to +/- strand):  ref base T (2/3);  new base C (0/21);  total (2/24)
Fisher's exact test for biased strand distribution p-value = 3.08e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.15e-02

GCTAAATCTGTTTGGCAACCAGAAGAATATGAGGTAGCAGAACATGAAATTCTGCTCGCCTTTCGTGCGTTACTTGCCTCTGATAGTGAACGTTTTCTGCCCTCCGCCTGGCAGCGATTCATTCAGCAATTAAATTATTGCTCGATGGAGGAGATTGATGAATTACGCTGGTCGCTGTGCACAATCGCCATGAACACTGCCCATTTATCCTTCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCGGTGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCGAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGATTTTTACCAAATGTTCGCTCGGATTTGGTCATCGCATCCACAATGGCTGACATTGTATTTAGCGCAACATCGCGCAGTGATTATCCCCGATGATGCAAAACTGCACAGA  >  NC_000913/677583‑678074
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     
gCTAAATCTGTTTGGCAACCAGAAGAATATGAGGTAGCAGAACATGAAATTCTGCTCGCCTTTCGTGCGTTACTTGCCTCTGATAGTGAACGTTTTCTGCCCTCCGCCTGGCAGCGATTCATTCAGCAATTAAATTATTGCTCGATGGAGGAGATTGATGAATTACGCTGGTCGCTGTGCACAATCGCCATGAACACTGCCCATTTATCCTTCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCggcggt                                                                                                                                                                                                                                                    <  1:596068/250‑1 (MQ=255)
gCTAAATCTGTTTGGCAACCAGAAGAATATGAGGTAGCAGAACATGAAATTCTGATCGCCTTTCGTGCGTTACTTGCCTCTGATAGTGAACGTTTTCTGCCCTCCGCCTGGCAGCGATTCATTCAGCAATTAAATTATTGCTCGATGGAGGAGATTGATGAATTACGCTGGTCGCTGTGCACAATCGCCATGAACACTGCCCATTTATCCTTCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCggcggtg                                                                                                                                                                                                                                                   >  1:322457/1‑250 (MQ=255)
    aaTCTGTTTGGCAACCAGAAGAATATGAGGTAGCAGAACATGAAATTCTGCTCGCCTTTCGTGCGTTACTTACCTCTGATAGTGAAAGTTTTCTGCCCTCCGCCTGGCAGCGATTCATTCAGCAATTAAATTATTGCTCGATGGAGGAGATTGATGAATTACGCTGGTCGCTGTGCACAATCGCCATGAACACTGCCCATTTATCCTTCGAGTGCGTGGTGTTCTTAGCAGAAAGATTGCggcggtaat                                                                                                                                                                                                                                                 >  1:611050/1‑246 (MQ=255)
            cggcattcagtagaaTATGAGGTTGCAGATCATGAAATTCTGCTCGCCTTTCGGGCGTTACTTTCCTCTGATAGTTAACGTTTTCTGCCCTCGGCCTGGCATCGATTCATTCAGTAATTAAATTATTGCTCGATGGAGGAGATTGATGAATTACGCTGGTAGCGGTGCGTAATCGACATGAACACTGCTCATTTATCCTACGAGTGCGAGGTGTTATTATCAGAAAGATTACGGCGGTTTCAGGAGGaaaa                                                                                                                                                                                                                                       <  2:238218/240‑1 (MQ=255)
                  ccAGAAGAATATGAGGTAGCAGAACATGAAATTCTGCTCGCCTTTCGTGCGTTACTTGCCTCTGATAGTGAACGTTTTCTGCCCTCCGCCTGGCAGCGATTCATTCAGCAATTAAATTATTGCTCGATGGAGGAGATTGATGAATTACGCTGGTCGCTGTGCACAATCGCCATGAACACTGCCCATTTATCCTTCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCggcggtagcatggggggaaaaat                                                                                                                                                                                                                                   >  1:371924/1‑232 (MQ=255)
                     gaagaaTATGAGGTAGCAGAACATGAAATTCTGCTCGCCTTTCGTTCGTTACTTGCCTCTGATAGTGAACTTTTTCTGCCCTTCGCCTGGCAGCTATTCATTCAGCAACTAAAGTATTGCTCGATGGAGGAGATTGATGAATTACGCTGGTCGCTGTGCACAATCGCCATGAACACTGCCCATTTATCCTTCGAGTGCGTGGTGTTATTAGCATAAAGATTGCGGCGGTTACAGGAGGAACATACCGGGGa                                                                                                                                                                                                                              <  2:382695/251‑1 (MQ=255)
                     gaagaaTATGAGGTAGCAGAACATGAAATTCTGCTCGCCTTTCGTGCGTTACTTGCCTCTGATAGTGAACGTTTTCTGCCCTCCGCCTGGCAGCGATTCATTCAGCAATTAAATTATTGCTCGATGGAGGAGATTGATGAATTACGCTGGTCGCTGTGCACAATCGCCATGAACACTGCCCATTTATCCTTCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCGGCGGTTACAGGAGGAAAATACCGGGGa                                                                                                                                                                                                                              <  1:568426/251‑1 (MQ=255)
                       agaaTATGAGGAAGCAGAACATGAAATTCTGCTCGCCTTTCGTGCGTTACTTGCCTCTGATAGTGAACGTTTTCTGCCCTCCGCCTGGCAGCGATTCATTCAGCAATTAAATTATTGCTCGATGGAGGAGATTGATGAATTACGCTGGTCGCTGTGCACAATCGCCATGAACACTGCCCATTTATCCTTCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCggcggtataagggaggaaataaaggggga                                                                                                                                                                                                                             >  1:366913/1‑227 (MQ=255)
                                      aGAACATGAAATTCTGCTCGCCTTTCGTGCGTTACTTGCCTCTGATAGTGAACGTTTTCTGCCCTCCGCCTGGCAGCGATTCATTCAGCAATTAAATTATTGCTCGATGGAGGAGATTGATGAATTACGCTGGTCGCTGTGCACAATCGCCATGAACACTGCCCATTTATCCTTCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCGGCGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAAc                                                                                                                                                                                                             <  1:356815/251‑1 (MQ=255)
                                                     gCTCGCCTTTCGTGCGTTACTTGCCTCTGATAGTGAACGTTTTCTGCCCTCCGCCTGGCAGCGATTCATTCAGCAATTAAATTATTGCTCGATGGAGGAGATTGATGAATTACGCTGGTCGCTGTGCACAATCGCCATGAACACTGCCCATTTATCCTTCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCGGCGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTtat                                                                                                                                                                                              <  1:181213/251‑1 (MQ=255)
                                                        cGCCTTTCTTGCGTTACTTTCCTTTGATAGTGAACGTTTTCTGCCCTCCGCCTGGCAGCGATGCATTCAGCAATTAAATTATTGCTCGATGGAGGAGATTGATGAATTACGCTGGTCGCTGTGCACAATCGCCATGAACACTACCCATTTATCCTTCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCGGCGGTTACAGGAGGAAAATACCGGGGAAATATACGAAGAAGAACTGGAATCCTTTTTATATg                                                                                                                                                                                           <  2:616894/251‑1 (MQ=255)
                                                                           gcctctgagaagtaaCGTTTTCTGCCTGCCGTCTGGCAGCGATTCATTCAGCTATTAAATTATTGTGCGATGGAGGAGATTGATGAATTACGCTGGTCTCTGTTCTCAATCGCCATGAACATTGACCATTTATCCTTCGAGTGCGTTGTGTTATTAGTAGAAAGATTGCGGCGGGTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAAGCCTTTTTATATTCCATTGCGAAGGGGAATg                                                                                                                                                                         <  2:414351/237‑1 (MQ=255)
                                                                                 gATAGTGAACGTTTTCTGCCCTCCGCCTGGCAGCGATTCATTCAGCAATTAAATTATTGCTCGATGGAGGAGATTGATGAATTACGCTGGTCGCTGTGCACAATCGCCATGAACACTGCCCATTTATCCTTCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCGGCGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCGAAGGGGAATGttttt                                                                                                                                                                    <  1:202152/249‑1 (MQ=255)
                                                                                 gATAGTGAACGTTTTCTGCCCTCCGCCTGGCAGCGATTCATTCAGCAATTAAATTATTGCTCGATGGAGGAGATTGATGAATTACGCTGGTCGCTGTGCACAATCGCCATGAACACTGCCCATTTATCCTTCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCGGCGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCGAAGGGGAATGTTTTTa                                                                                                                                                                   <  1:673736/250‑1 (MQ=255)
                                                                                                          ccTGGCAGCGATTCATTCAGCAATTAAATTATTGCTCGATGGAGGAGATTGATGAATTACGCTGGTCGCTGGGCACAAACGCCATGAACACTGCCCATTTATCCCTCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCGGCGGTTACAGTAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCGAAGGGGAATGTTTTTAACTTCCAGATCATTCTGCATCTGccc                                                                                                                                         <  2:655982/251‑1 (MQ=255)
                                                                                                             gtacgagATTCATTCAGCAATTAAATTATTGCCCGATGTAGGAGAGTTATGAATGACGCTGGTAGCTGTGAACAATTGGCATGAACACTGCCCAGTTATCCTTCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCGGCGGTTACAGGAGGAAAATACCGGGGAAATAGACGAATAAGAACTGGAATCCTTTTTATATGGCATTGCGAAGGGGATTTTTTTTAAGTTCTAGTCCAGTCTTCATCTGCCCg                                                                                                                                        <  2:611050/243‑1 (MQ=255)
                                                                                                                                   aaattaTTGCTCGATGGAGGAGATTGATGAATTTCGCGGGTCGCTGTGCACAATCGCCATGAACACTGCCCATTTATCCTTCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCGGTGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCGAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCg                                                                                                                <  2:190360/251‑1 (MQ=255)
                                                                                                                                   aaagtattGCGCGATGGAGGAGATTGATGAATTTCGCTGGCCGCTGTGCACAATCGCCATTAACACTGCCCATTTATCCGTCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCGGCGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCGAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATc                                                                                                                 <  1:25698/246‑1 (MQ=255)
                                                                                                                                           gCTCTATGGAGGAGATTGATGAATTACGCTGGTCGCTGTGCACAATCGCCATGAACACTGCCCATTTATCCTTCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCGGCGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCGAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGATTTTTAc                                                                                                        <  1:648175/251‑1 (MQ=255)
                                                                                                                                              cGATGGAGGAGATTGATGAATTACGCTGGTCGCTGTGCACAATCGCCATGAACACTGCCCATTTATCCTTCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCGGCGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCGAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGATTTTTACCaa                                                                                                     <  1:380614/251‑1 (MQ=255)
                                                                                                                                                                      gCTGGTCGCTGTGCACAATCGCCATGAACACTGCCCATTTATCCTTCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCGGCGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCGAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGATTTTTACCAAATGTTCGCTCGGATTTGGTCATCg                                                                             <  2:182937/251‑1 (MQ=255)
                                                                                                                                                                                    acaATCGCCATGAACACTGCCCATTTATCCTTCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCGGCGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCGAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGATTTTTACCAAATGTTCGCTCGGATTTGGTCATCGCATCCACAATGGCt                                                               <  1:614953/251‑1 (MQ=255)
                                                                                                                                                                                       aTCGCCATGAACACTGCCCATTTATCCTTCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCGGTGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCGAAGGGGAATGTTTTTAAc                                                                                                                                                                 >  1:275856/1‑150 (MQ=255)
                                                                                                                                                                                       aTCGCCATGAACACTGCCCATTTATCCTTCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCGGTGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCGAAGGGGAATGTTTTTAAc                                                                                                                                                                 <  2:275856/150‑1 (MQ=255)
                                                                                                                                                                                            cATGAACACTGCCCATTTATCCTTCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCGGCGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTTGAATCCTTTTTATATGCCATTGCGAAGGGGaa                                                                                                                                                                           <  1:423834/135‑1 (MQ=255)
                                                                                                                                                                                                 acacTGCCCATTTATCCTTCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCGGTGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCGAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTACCCGTTGCCGTGCAAAATGACACCATCGATTTTTACCAAATGTTCGCTCGGATTTGGTCATCGCATCCACAATGGCTGACATTGTAttt                                                   >  1:339591/1‑250 (MQ=255)
                                                                                                                                                                                                      gCCCATTTATCCTTCGAGGTCGGGGTGTTATTAGCAGAAAGATTGCGGCGGTTACAGTAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCGAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGATTTTTACCAAATGTTCGCTCGGATTTGGTCATCGCATCCACAATGGCTGACATTGTATTTAgcgc                                              <  2:583306/250‑1 (MQ=255)
                                                                                                                                                                                                              aTCCTTCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCGGCGGTGACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAAGCCTTTTTATATGCCATTGCGAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGATTTTTTCCAAATGTTCTCTCGGATTTGTTCAGCGCATCCACAATGGCTGACATTGTATTTAGCGCAACATcgc                                      <  2:566460/250‑1 (MQ=255)
                                                                                                                                                                                                                      aGTGCGTGGGGTGAATAGCAGAAAGATTGCGGCGGTTTCAGGAGAAAAATACCTGGGAAATAGCCTAAGAAGAACTGGAATCCTGTTTATATGCCATTGCGAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTTCAAAATGACACCATCGATTTTTACCAAATGTTCGCTCGGATTTGGTCATCGGATCCACAATGGCGGACAGTGTATTTAGCGCAACATCTCGCAGTGAtt                             <  2:585181/251‑1 (MQ=255)
                                                                                                                                                                                                                      aGTGCGGGGGGTTATTAGCAGAAAGATTGCGGCGGTTACAGGAGGAAAATACCGGGGAAAGAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCGAAGGGGAATGTTTATAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGATTTTTACCAAATGTTCGCTCGGATTTGGTCATCGCATCCACAATGGCTGACATTGTATTTAGCGCATCATCGCGCTGTGAtt                             <  2:85547/251‑1 (MQ=255)
                                                                                                                                                                                                                                               aTTGCGGTGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCGAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTACCCGTTGCCGTGCAAAATGACACCATCGATTTTTACCAAATGTTCGCTCGGATTTGGTCATCGCATCCACAATGGCTGACATTGTATTTAGCGCAACATCGCGCAGTGATTATCCCCGATGATGCAAAACTGcaca    <  2:339591/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 ggcggcgGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTGTTTATATGCCATTGCGAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGATTTTTACCAAATGTTCGCTCGGATTTGGTCATCGCATCCACAATGGCTGACATTGTATTTATCGCAACATCGCGCAGTGATTATCCCCGATGATGCAAAACTGCACAGa  <  2:660229/250‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     
GCTAAATCTGTTTGGCAACCAGAAGAATATGAGGTAGCAGAACATGAAATTCTGCTCGCCTTTCGTGCGTTACTTGCCTCTGATAGTGAACGTTTTCTGCCCTCCGCCTGGCAGCGATTCATTCAGCAATTAAATTATTGCTCGATGGAGGAGATTGATGAATTACGCTGGTCGCTGTGCACAATCGCCATGAACACTGCCCATTTATCCTTCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCGGTGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCGAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGATTTTTACCAAATGTTCGCTCGGATTTGGTCATCGCATCCACAATGGCTGACATTGTATTTAGCGCAACATCGCGCAGTGATTATCCCCGATGATGCAAAACTGCACAGA  >  NC_000913/677583‑678074

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: