Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 681,694 A→G 37.5% intergenic (+31/+29) djlC → / ← hscC J domain‑containing HscC co‑chaperone; Hsc56/Hsp70 family chaperone Hsc62; RpoD‑binding transcription inhibitor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913681,6940AG37.5% 8.7 / 13.7 16intergenic (+31/+29)djlC/hscCJ domain‑containing HscC co‑chaperone; Hsc56/Hsp70 family chaperone Hsc62; RpoD‑binding transcription inhibitor
Reads supporting (aligned to +/- strand):  ref base A (4/6);  new base G (2/4);  total (6/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.28e-01

AAAAATGATGAGGGTGAAAAACCACAGAGCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGG  >  NC_000913/681553‑681817
                                                                                                                                             |                                                                                                                           
aaaaaTGATGAGGGTGAAAAACCACAGAGCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGttt                                                                                                                       <  2:391393/148‑1 (MQ=255)
    atgatgAGGGTGAAAAACCACAGAGCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTAcgc                                                                                                                   <  1:10449/148‑1 (MQ=255)
          aGGGTGAAAAACCACAGAGCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCGGCGTAAACGCTGCATc                                                                                                             <  4:919158/148‑1 (MQ=255)
                               ttACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCGGCGTAAAcgctg                                                                                                                 <  3:1093973/123‑3 (MQ=255)
                               ttACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCGGCGTAAAcgctg                                                                                                                 >  4:1093973/1‑121 (MQ=255)
                                                 tCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAg                                                                                                                            >  3:418616/1‑94 (MQ=255)
                                                 tCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAg                                                                                                                            <  4:418616/94‑1 (MQ=255)
                                                                 ggTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTAcg                                                                                                                    >  1:383724/1‑86 (MQ=255)
                                                                 ggTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTAcg                                                                                                                    <  2:383724/86‑1 (MQ=255)
                                                                        ggTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAg                                                                                                                            <  3:103349/71‑1 (MQ=255)
                                                                        ggTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAg                                                                                                                            >  4:103349/1‑71 (MQ=255)
                                                                          ttATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATg                                                                                            >  1:638350/1‑101 (MQ=255)
                                                                          ttATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATg                                                                                            <  2:638350/101‑1 (MQ=255)
                                                                           tATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCGGCGTAAAcgc                                                                                                                   <  2:604903/77‑1 (MQ=38)
                                                                                                                                         aTGCGGCGTAAACGCTGCATCCGGCATCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCgg  >  1:441391/1‑128 (MQ=255)
                                                                                                                                         aTGCGGCGTAAACGCTGCATCCGGCATCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCgg  <  2:441391/128‑1 (MQ=255)
                                                                                                                                             |                                                                                                                           
AAAAATGATGAGGGTGAAAAACCACAGAGCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGG  >  NC_000913/681553‑681817

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