Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 677,889 G→A A159T (GCC→ACC)  djlB → putative HscC co‑chaperone, uncharacterized J domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913677,8890GA97.8% 132.9 / ‑6.1 45A159T (GCC→ACC) djlBputative HscC co‑chaperone, uncharacterized J domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (21/23);  total (22/23)
Fisher's exact test for biased strand distribution p-value = 4.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCGCTGTGCACAATCGCCATGAACACTGCCCATTTATCCTTCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCGGTGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCGAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGATTTTTACCAAATGTTCGCTCGGATTTGGTCATCGCATCCACAATGGCTG  >  NC_000913/677753‑678014
                                                                                                                                        |                                                                                                                             
gTCGCTGTGCACAATCGCCATGAACACTGCCCATTTATCCTTCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCGGTGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATAcc                                                                                                                             >  2:422697/1‑139 (MQ=255)
             aTCGCCATGAACACTGCCCATTTATCCTTCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCGGTGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATACCATTGCGAAGGGGa                                                                                                                <  1:129077/139‑1 (MQ=255)
              tCGCCATGAACACTGCCCATTTATCCTTCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCGGTGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATACCATTGCGAAGGGGaa                                                                                                               >  2:1332885/1‑139 (MQ=255)
                     gAACACTGCCCATTTATCCTTCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCGGTGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATACCATTGCGAAGGGGAATGttttt                                                                                                        <  1:130833/139‑1 (MQ=255)
                     gAACACTGCCCATTTATCCTTCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCGGTGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATACCATTGCGAAGGGGAATGttttt                                                                                                        <  1:223241/139‑1 (MQ=255)
                      aaCACTGCCCATTTATCCTTCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCGGTGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATACCATTGCGAAGGGGAATGTTTTTa                                                                                                       <  1:562233/139‑1 (MQ=255)
                           tGCCCATTTATCCTTCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCGGTGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATACCATTGCGAAGGGGAATGTTTTTa                                                                                                       <  1:1054092/134‑1 (MQ=255)
                           tGCCCATTTATCCTTCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCGGTGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATACCATTGCGAAGGGGAATGTTTTTa                                                                                                       >  2:1054092/1‑134 (MQ=255)
                           tGCCCATTTATCCTTCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCGGTGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATACCATTGCGAAGGGGAATGTTTTTAACTTc                                                                                                  >  2:1060142/1‑139 (MQ=255)
                            gCCCATTTATCCTTCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCGGTGGTTACAGGAGGAAAATACCGGGGCAATAGACGAAGAAGAACTGGAATCCTTTTTATATACCATTGCGAAATGGAATGTTTTTAACTTcc                                                                                                 >  2:746258/1‑139 (MQ=255)
                            gCCCACTTATCCTTCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCGGTGTTTACAGGAGGAAAATACCGGGGATGTAGACGAAGAAGAACTGGAATCCTTTTTATATACCATTGCGGAGGGTAACGTGTTTACCTTcc                                                                                                 >  2:988298/1‑139 (MQ=255)
                               cATTTATCCTTCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCGGTGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATACCATTGCGAAGGGGAATGTTTTTAACTTCCAGa                                                                                              <  1:988298/139‑1 (MQ=255)
                                aTTTATCCTTCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCGGTGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCa                                                                                                                            >  1:861547/1‑108 (MQ=255)
                                aTTTATCCTTCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCGGTGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATACCa                                                                                                                            <  2:861547/108‑1 (MQ=255)
                                    aTCCTTCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCGGTGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATACCATTGCGAAGGGGAATGTTTTTAACTTCCAGACCAtt                                                                                         <  1:383169/139‑1 (MQ=255)
                                     tCCTTCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCGGTGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATACCATTGCGAAGGGGAATGTTTTTAACTTCCAGACCATTc                                                                                        >  1:127292/1‑139 (MQ=255)
                                         tCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCGGTGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATACCATTGCGAAGGGGAATGttttt                                                                                                        >  1:68171/1‑119 (MQ=255)
                                         tCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCGGTGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATACCATTGCGAAGGGGAATGttttt                                                                                                        <  2:68171/119‑1 (MQ=255)
                                            aGTGCGTGGTGTTATTAGCAGAAAGATTGCGGTGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATACCATTGCGAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCAt                                                                                   <  2:217411/137‑1 (MQ=255)
                                            aGTGCGTGGTGTTATTAGCAGAAAGATTGCGGTGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATACCATTGCGAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCAt                                                                                   >  1:217411/1‑137 (MQ=255)
                                                cGTGGTGTTATTAGCAGAAAGATTGCGGTGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATACCATTGCGAAGGGGAATGTTTTTAACtt                                                                                                   <  1:882902/117‑1 (MQ=255)
                                                cGTGGTGTTATTAGCAGAAAGATTGCGGTGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATACCATTGCGAAGGGGAATGTTTTTAACtt                                                                                                   >  2:882902/1‑117 (MQ=255)
                                                cGTGGTGTTATTAGCAGAAAGATTGCGGTGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATACCATTGCGAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGccc                                                                             >  2:926611/1‑139 (MQ=255)
                                                     tgtTATTAGCAGAAAGATTGCGGTGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATACCATTGCGAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGc                                                                        <  1:1332885/139‑1 (MQ=255)
                                                         attaGCAGAAAGATTGCGGTGGTTACAGGAGGAAACTACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATACCATTGCGAAGGGGAATGTTTTTAACTTCCAGa                                                                                              <  2:851133/113‑1 (MQ=255)
                                                         attaGCAGAAAGATTGCGGTGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATACCATTGCGAAGGGGAATGTTTTTAACTTCCAGa                                                                                              >  1:851133/1‑113 (MQ=255)
                                                               aGAAAGATTGCGGTGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATACCATTGCGAAGGGGAATGTTTTTAACTTCCAGACCAtt                                                                                         <  1:464056/112‑1 (MQ=255)
                                                               aGAAAGATTGCGGTGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATACCATTGCGAAGGGGAATGTTTTTAACTTCCAGACCAtt                                                                                         >  2:464056/1‑112 (MQ=255)
                                                                 aaaGATTGCGGTGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATACCATTGCGAAGGGGAATGttttt                                                                                                        <  2:1154301/95‑1 (MQ=255)
                                                                 aaaGATTGCGGTGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATACCATTGCGAAGGGGAATGttttt                                                                                                        >  1:1154301/1‑95 (MQ=255)
                                                                    gATTGCGGTGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATACCATTGCGAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGacac                                                         >  1:724105/1‑139 (MQ=255)
                                                                        gCGGTGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATACCATTGCGAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATc                                                     >  2:1020887/1‑139 (MQ=255)
                                                                        gCGGTGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATACCATTGCGAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATc                                                     >  1:131634/1‑139 (MQ=255)
                                                                            tggtTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATACCATTGCGAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGAtt                                                 <  1:1020887/139‑1 (MQ=255)
                                                                                       ggaAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATACCATTGGGAGGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTg                                                                    >  2:889525/1‑109 (MQ=255)
                                                                                       ggaAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATACCATTGCGAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTg                                                                    <  1:889525/109‑1 (MQ=255)
                                                                                        gaAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATACCATTGCGAAGGGGAATGTTTTTAACTTCCAGa                                                                                              <  2:356297/82‑1 (MQ=255)
                                                                                        gaAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATACCATTGCGAAGGGGAATGTTTTTAACTTCCAGa                                                                                              >  1:356297/1‑82 (MQ=255)
                                                                                        gaAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATACCATTGCGAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCAt                                                                                   >  2:1100346/1‑93 (MQ=255)
                                                                                        gaAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATACCATTGCGAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCAt                                                                                   <  1:1100346/93‑1 (MQ=255)
                                                                                        gaAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATACCATTGCGAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGATTTTTACCAAATGt                                     >  1:485691/1‑139 (MQ=255)
                                                                                             tACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATACCATTGCGAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGATTTTTACCAAATGTTCGCt                                <  1:422697/139‑1 (MQ=255)
                                                                                                        tAGACGAAGAAGAACTGGAATCCTTTTTATATACCATTGCGAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGATTTTTACCAAATGTTCGCTCGGATTTGGTc                     <  2:131634/139‑1 (MQ=255)
                                                                                                                      cTGGAATCCTTTTTATATACCATTGCGAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGATTTTTACCAAATGTTCGCTCGGATTTGGTCATCGCATCCACAAt       <  2:724105/139‑1 (MQ=255)
                                                                                                                           aTCCTTTTTATATACCATTGCGAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGATTTTTACCAAATGTTCGCTCGGATTTGGTCATCGCATCCACAATGGCTg  <  1:926611/139‑1 (MQ=255)
                                                                                                                                     tataCCATTGCGAAGGGGAATGTTTTTAACTTCCa                                                                                                <  1:956715/35‑1 (MQ=255)
                                                                                                                                     tataCCATTGCGAAGGGGAATGTTTTTAACTTCCa                                                                                                >  2:956715/1‑35 (MQ=255)
                                                                                                                                        |                                                                                                                             
GTCGCTGTGCACAATCGCCATGAACACTGCCCATTTATCCTTCGAGTGCGTGGTGTTATTAGCAGAAAGATTGCGGTGGTTACAGGAGGAAAATACCGGGGAAATAGACGAAGAAGAACTGGAATCCTTTTTATATGCCATTGCGAAGGGGAATGTTTTTAACTTCCAGACCATTCTGCATCTGCCCGTTGCCGTGCAAAATGACACCATCGATTTTTACCAAATGTTCGCTCGGATTTGGTCATCGCATCCACAATGGCTG  >  NC_000913/677753‑678014

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