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,8290TC87.0% 48.5 / 0.7 23W139R (TGG→CGG) djlBJ domain‑containing protein DjlB
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base C (1/19);  total (3/20)
Fisher's exact test for biased strand distribution p-value = 3.44e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

GGCAACCAGAAGAATATGAGGTAGCAGAACATGAAATTCTGCTCGCCTTTCGTGCGTTACTTGCCTCTGATAGTGAACGTTTTCTGCCCTCCGCCTGGCAGCGATTCATTCAGCAATTAAATTATTGCTCGATGGAGGAGATTGATGAATTACGCTGGTCGCTGTGCACAATCGCCATGAACACTGCCCATTTATCCTTCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCGGTGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCGAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGATTTTTACCAAATGTTCGCTCGGATTTGGTCATCGCATCCACAATGGCTGACATTGTATTTAGCGCAACATCGCGCAGTGATTATCCCCGATGAT  >  NC_000913/677596‑678059
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                      
ggCAACCAGAAGAATATGAGGTAGCAGAACATGAAATTCTGCTCGCCTTTCGTGCGTTACTTGCCTCTGATAGTGAACGTTTTCTGCCCTCCGCCTGGCAGCGATTCATTCAGCAATTAACTTATTGCTCGATGGAGGAGATTGATGACTTACGCTGGTCGCTGTGCACAATCGCCATGAACACTGCCCATTTATCCTTCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCGGCGGTTatacgtggaaga                                                                                                                                                                                                                        >  1:486536/1‑239 (MQ=255)
       agaagaATATGAGGTAGCAGAACATGAAATTCTGCTCGCCTTTCGTGCGTTACTTGCCTCTGATAGTGAACGTTTTCTGCCCTCCGCCTGGCAGCGATTCATTCAGCAATTAAATTATTGCTCGATGGAGGAGATTGATGAATTACGCTGGTCGCTGTGCACAATCGCCATGAACACTGCCCATTTATCCTTCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCGGCGGTTACAGGAGGAAAATACCgggg                                                                                                                                                                                                                <  1:499273/251‑1 (MQ=255)
                    gTAGCAGAACATTAAATTCTGCTCGCCTTTCGTGCGTTACTTGCCTCTGATAGTGAACGTTTTCTGCCCTCCGCCTGGCAGCGATTCATTCAGCAATTAAATTATTGCTCGAGGGAGGAGATTGATGAATTACGCTGGTCGCTGTGCACAATCGCCATGAACACTTCCCATTTATCCTTCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCGGCGGTTACAGGAGGAAAATACCGGGGAAATAGACgaaga                                                                                                                                                                                                   <  2:486536/251‑1 (MQ=255)
                                         cTCGCCTTTCGTGCGTTACTTGCCTCTGATAGTGAACGTTTTCTGCCCTCCGCCTGGCAGCGATTAATTCAGCAATTAAATTATTGCTCGATGGAGGAGATTGATGAATTACGCTGGTCGCTGTGCACAATCGCCATGAACACTGCCCATTTATCCTTCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCGGTGGTTa                                                                                                                                                                                                                                   <  1:216187/198‑1 (MQ=255)
                                         cTCGCCTTTCGTGCGTTACTTGCCTCTGATAGTGAACGTTTTCTGCCCTCCGCCTGGCAGCGATTAATTCAGCAATTAAATTATTGCTCGATGGAGGAGATTGATGAATTACGCTGGTCGCTGTGCACAATCGCCATGAACACTGCCCATTTATCCTTCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCGGTGGTTa                                                                                                                                                                                                                                   >  2:216187/1‑198 (MQ=255)
                                             ccTTTCGTGCGTTACCTGCCTCTGATAGTGAACGTTTTCTGTCCTCCGCCTGGCAGCGATTCATTCAGCAATTAAAGTATTGCTCGATGGAGGAGATTGATGAATTTCGCTGGTCGCTGTGCACAATCGCCATGAACACTGCCCATTTATCCTTCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCGGCGGTTACAGGAGGAAAATCCCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGc                                                                                                                                                                           <  2:174694/250‑1 (MQ=255)
                                                         tACTTGCCTCTGATAGTGAACGTTTTCTGCCCTCCGCCTGGCAGCGATTCATTCAGCAATTAAATTATTGCTCGATGGAGGAGATTGATGAATTACGCTGGTCGCTGTGCACAATCGCCATGAACACTGCCCATTTATCCTTCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCGGCGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCGAAggg                                                                                                                                                               <  1:494197/250‑1 (MQ=255)
                                                             tGCCTCTGATAGTTAACGTTTTCTGCCCTCCGCCTGGCAGCGATTCATTCAGCAATTAAATTATTGCTCGATGGAGGAGATTGATGAATTACGCTGGTCGCTGTGCACAATCGCCATGAACACTGCCCATTTATCCTTCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCGGCGGTTACAGGAGGAAAATACCGGGGAAAGAGACGAAGAAGAACGGGAATCCTGTTTATATGCCATTGCGAAGGGGAAt                                                                                                                                                           <  1:139036/250‑1 (MQ=255)
                                                                          gAACGTTTTCTGCCCTCCGCCTGGCAGCGATTCATTAAGCAATTAAATTATTGCTCGATGGAGGAGATTGATGAATTACGCTGGTCGCTGTGCACAATCGCCATGAACACTGCCCATTTATCCGTCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCGGCGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCGAAGGGGAATGTTTTTAACTTcc                                                                                                                                              <  1:111878/250‑1 (MQ=255)
                                                                                      cccTCCGCCTGGCAGCGATTCATTCAGCAATTAAATTATTGCTCGATGGAGGAGATTGATGAATTACGCTGGTCGCTGTGCACAATCGCCATGAACACTGCCCATTTATCCTTCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCGGCGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCGAAGGGGAATGTTTTTAACTTCCAGACCTTTCTGCa                                                                                                                                 <  1:306443/251‑1 (MQ=255)
                                                                                                 gCAGGGATTCATTCTGCAATTAAATTATTGCTCGATGGAGGTGAGTGATTAATTACTCTGGTCGCTGTGCATAATCGCCATGATCGCAGCCGATTGATCCTTCGAGTGCGTGGTGTTATTAGCAGAAAGATGGCGGCTGTGACAGGAGGAAAAGACCGGGGAAATGGACGATGAAGATCTTGAATCTGTTTTGTAGGCCGTCGCGAAGGGGAATGTGTTGAACTTCCAGAGCATTATGCATCTGCCCGtgg                                                                                                                      <  2:328671/251‑3 (MQ=255)
                                                                                                         tcattcaGCAATTAAATTATTGCTCGATGGAGGAGATTGATGAATTACGCTGGTCGCTGTGCACAATCGCCATGAACACTGCCCATTTATCCTTCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCGGCGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCGAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCTTTGCCGTGCa                                                                                                               <  2:437410/250‑1 (MQ=255)
                                                                                                                                               gattaaTTACGCTGGTCGCTGTGCACAATCGCCATGAACACTGCCCATTTATCCTTCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCGGCGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCGAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGTCTCCATCGATTTTTACCAAATGTTCGCTCGGAt                                                                        <  1:198658/247‑1 (MQ=255)
                                                                                                                                                      tactcTGGTCGCTGTGCACAATCGCCATGAACACTGCCCATTTATCCTTCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCGGCGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCGAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCTTGCAAAATGACACCATCGATTTTTACCAAATGTTCGCTCGGATTTGGTc                                                                  <  1:367855/246‑1 (MQ=255)
                                                                                                                                                            ggTCGCTGTGCACAATCGCCATGAACACTGCCCATTTATCCTTCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCGGCGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCGAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGATTTTTACCAAATGTTCGCTCGGATTTGGTCATCGCAt                                                           <  1:199788/251‑1 (MQ=255)
                                                                                                                                                             gTCGCTGTGCACAATCGCCATGAACACTGCCCATTTATCCTTCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCGGCGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCGAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGAtt                                                                                              <  2:438709/215‑1 (MQ=255)
                                                                                                                                                                     gCACAATCGCCATGAACACTGCCCATTTATCCTTCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCGGCGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCGAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGATTTTTACCAAATGTTCGCTCGGATTTGGTCATCGCATCCACAATg                                                   <  1:12775/250‑1 (MQ=255)
                                                                                                                                                                         aaTCGCCATGAACACTGCCCATTTATCCTTCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCGGTGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGACCTGGAATCCTTTTTATATGCCATCGCGAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGATTTTTACCAAATGTTCGCTCGGCTTTGGTCATCGCATCCAAAATGGCTgc                                              >  2:27579/1‑250 (MQ=255)
                                                                                                                                                                                aTGAACACTGCCCATTTATCCTTCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCGGCGGTTACAGGAGGAAAATACCGGGGAAATAGACg                                                                                                                                                                                                       <  2:17061/91‑1 (MQ=255)
                                                                                                                                                                                     cacTGCCCATTTATCCTTCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCGGCGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCGAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGATTTTTACCAAATGTTCGCTCGGATTTGGTCATCGCATCCACAATGGCTGACa                                            <  2:3287/241‑1 (MQ=255)
                                                                                                                                                                                                         aGTGCGTGGTGTTATTAGCAGAAAGATTGCGGCGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCGAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGATTTTTACCAAATGTTCGCTCGGATTTGGTCATCGCATCCACAAGTGCTGACATTGTATTTAGCGCAGCATCGCGCAGTGAtt              <  1:136431/251‑1 (MQ=255)
                                                                                                                                                                                                                   gtTATTAGCAGAAAGATTGCGGCGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCGAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGATTTTTACCAAATGTTCGCTCGGATTTGGTCATCGCATCCACAATGGCTTACATTGGATTTAGCGCAACATCGCGCAGTGATTATCCCCgatg    <  1:446962/251‑1 (MQ=255)
                                                                                                                                                                                                                      attaGCAGAAAGATTGCGGCGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTAGAATCCTTTTTATATGCCATTGCGAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGATTTTTACCAAATGTTCGCTCGGATTTGGTCATCGCATCCACAATGGCTGACATTGTATTTAGCGCAACATCGCGCAGTGATTATCCCCgatgat  <  1:199437/250‑1 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                      
GGCAACCAGAAGAATATGAGGTAGCAGAACATGAAATTCTGCTCGCCTTTCGTGCGTTACTTGCCTCTGATAGTGAACGTTTTCTGCCCTCCGCCTGGCAGCGATTCATTCAGCAATTAAATTATTGCTCGATGGAGGAGATTGATGAATTACGCTGGTCGCTGTGCACAATCGCCATGAACACTGCCCATTTATCCTTCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCGGTGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCGAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGATTTTTACCAAATGTTCGCTCGGATTTGGTCATCGCATCCACAATGGCTGACATTGTATTTAGCGCAACATCGCGCAGTGATTATCCCCGATGAT  >  NC_000913/677596‑678059

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