Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 682,733 C→T 11.6% A221T (GCC→ACC)  hscC ← 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_000913682,7330CT11.6% 78.0 / 4.4 35A221T (GCC→ACC) hscCHsp70 family chaperone Hsc62; RpoD‑binding transcription inhibitor
Reads supporting (aligned to +/- strand):  ref base C (18/12);  new base T (2/2);  total (21/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.34e-02

TCTTCCAGTTCGTTCTCGTAAAATTCGCATTCCCGCGTTTCTTCCTGATACTGCCAGCGAATGTGCAATGGCGATTGATT‑GCTACATTTTGCCGCTTCCACACAGGCGTACAAGGCTGCCAGTTCACTCTCGTTAAGCGTGGTCCTGGCGACATCCGCGCGTTTTAAAACCTCATCGACCAGCATATGGGTAAAATCTTCGCCACCAAGAAAGTTGTCGCCAGCGGAGGCGTGAACTTCAATCACCGGCGTGGCGTACTCAAGCACCGTAACGTCAAACGTGCCGCCACCGAGATCA  >  NC_000913/682585‑682881
                                                                                                                                                     |                                                                                                                                                    
tcctCCAGTTCGTTCTCGTAAAATTCGCATTCCCGCGTTTCTTCCTGATACTGCCAGCGAATGTGCAATGGCGATTGATT‑GCTACATTTTGCCGCTTCCACACAGGCGTACAAGGCTGCCAGTTCACTCTCGTTAAGCGTGGTCCTGGc                                                                                                                                                      <  1:277034/146‑1 (MQ=255)
tCTTCCCGTCCGTTCTCGTAAAATTCGCATTCCCGCGTTTCTTCCTGATACTGCCAGCGAATGTGCAATGGCGATTGATT‑GCTACATTTTGCCGCTTCCACACAGGCGTACAAGGCTGCCAGTTCACTCTCGTTAAGCGTGGTCCTGGc                                                                                                                                                      <  1:499986/149‑1 (MQ=255)
ccTTCCAGGTCGTTCTCGTGAGACTCGCATTCCCGCGTTTCTTCCTGATACTGCCAGCGAATGTGCAATGGGGATTGCTTGGCGGCAGTCT‑TCTCTTGCCCACAGCGATACAAGGCTGCCAGTTCCCTCTCGTTAAGCGTGGTCCTGGc                                                                                                                                                      <  2:768845/148‑1 (MQ=255)
           gttcTCGTAAAATTCGTATTCCCGCGTTTCTTCCTGATACTGCCAGCGAATGTGCAATGGCGATTGATT‑GCTACATTTTGCCGCTTCCACACAGGCGTACAAGGCTGCCAGTTCACTCTCGTTAAGCGTGGTCCTGGCGACATCcgcgc                                                                                                                                           >  2:573609/1‑149 (MQ=255)
                cGTAAAATTCGCATTCCCGCGTTTCTTCCTGATACTGCCAGCGAATGTGCAATGGCGATTGATT‑GCTACATTTTGCCGCTTCCACACAGGCGTACAAGGCTGCCAGTTCACTCTCGTTAAGCGTGGTCCTGGCGACATCCGCGCGtttt                                                                                                                                      <  1:627388/149‑1 (MQ=255)
                    aaaTTCGCATTCCCGCGTTTCTTCCTGATACTGCCAGCGAATGTGCAATGGCGATTGATT‑GCTACATTTTGCCGCTTCCACACAGGCGTACAAGGCTGCCAGTTCACTCTCGTTAAGCGTGGTCCTGGCGACATCCGCGCGTTTTaaaa                                                                                                                                  <  1:555698/149‑1 (MQ=255)
                           cATTCCCGCGTTTCTTCCTGATACTGCCAGCGAATGTGCAATGGCGATTGATT‑GCTACATTTTGCCGCTTCCACACAGGCGTACAAGGCTGCCAGTTCACTCTCGTTAAGCGTGGTCCTGGCGACATCCGCGCGTTTTAAAACCTCATc                                                                                                                           <  1:590743/149‑1 (MQ=255)
                                     tttCTTCCTGATACTGCCAGCGAATGTGCAATGGCGATTGATT‑GCTACACTTTGCCGCTTCCACACAGGCGTACAAGGCTGCCAGTTCACTCTCGTTAAGCGTGGTCCTGGCGACATCCGCGCGTTTTAAAACCTCATCGACCAGCata                                                                                                                 >  1:332202/1‑149 (MQ=255)
                                        cttcCTGATACTGCCAGCGAATGTGCAATGGCGATTGATT‑GCTACATTTTGCCGCTTCCACACAGGCGTACAAGGCTGCCAGTTCACTCTCGTTAAGCGTGGTCCTGGCGACATCCGCGCGTTTTAAAACCTCATCGACCAGCATATgg                                                                                                              <  1:525744/149‑1 (MQ=255)
                                           ccTGATACTGCCAGCGAATGTGCAATGGCGATTGATT‑GCTACATTTTGCCGCTTCCACACAGGCGTACAAGGCTGCCAGTCCACTCTCGTTAAGCGTGGTCCTGGCGACATCCGCGCGTTTTAAAACCTCATCGACCAGCATATGGGTa                                                                                                           >  2:571577/1‑149 (MQ=255)
                                             tGATACTGCCAGCGAATGTGCAATGGCGATTGATT‑GCTACATTTTGCCGCTTCCACACAGGCGTACAAGGCTGCCAGTTCACTCTCGTTAAGCGTGGTCCTGGCGACATCCGCGCGTTTTAAAACCTCATCGACCAGCATATGGGTaaa                                                                                                         <  1:528072/149‑1 (MQ=255)
                                              gATACTGCCAGCGAATGTGCAATGGCGATTGATT‑GCTACATTTTGCCGCTTCCACACAGGCGTACAAGGCTGCCAGTTCACTCTCGTTAAGCGTGGTCCTGGCGACATCCGCGCGTTTTAAAACCTCATCGACCAGCATATGGGTaaaa                                                                                                        >  2:423261/1‑149 (MQ=255)
                                                       aGCGAATGTGCAATGGCGATTGATT‑GCTACATTTTGCCGCTTCCACACAGGCGTACAAGGCTGCCAGTTCACTCTCGTTAAGCGTGGTCCTGGCGACATCCGCGCGTTTTAAAACCTCATCGACCCGCATCTGGGTCAAATCTTCGccc                                                                                               >  1:24739/1‑148 (MQ=255)
                                                               tgCAATGGCGATTGATT‑GCTACATTTTGCCGCTTCCACACAGGCGTACAAGGCTGCCAGTTCACTCTCGTTAAGCGTGGTCCTGGCGACATCCGCGCGTTTTAAAACCTCATCGACCAGCATATGGGTAAAATCTTCGCCACCAAGaaa                                                                                       <  2:209329/149‑1 (MQ=255)
                                                                  aaTGGCGATTGATT‑GCTACATTTTGCCGCTTCCACACAGGCGTACAAGGCTGCCAGTTCACTCTCGTTAAGCGTGGTCCTGGCGACATCCGCGCGTTTTAAAACCTCATCGACCAGCATATGGGTAAAATCTTCGCCACCAAGAAAGtt                                                                                    >  1:154892/1‑149 (MQ=255)
                                                                     ggCGATTGATT‑GCTACATTTTGCCGCTTCCACACAGGCGTACAAGGCTGCCAGTTCACTCTCGTTAAGCGTGGTCCTGGCGACATCCGCGCGTTTTAAAACCTCATCGACCAGCATATGGGTAAAATCTTCGCCACCAAGAAAGTTGTc                                                                                 <  1:275960/149‑1 (MQ=255)
                                                                        gattgatt‑GCTACATTTTGCCGCTTCCACATAGGCGTACAAGGCTGCCAGTTCACTCTCGTTAAGCGTGGTCCTGGTGACATCCGCGCGTTTTAAAACCTCATCGACCAGCATATGGGTAAAATCTTCGCCACCAAGAAAGTTGTCGcc                                                                              >  1:415645/1‑149 (MQ=255)
                                                                        gattgatt‑GCTACATTTTGCCGCTTCCACACAGGCGTACAAGGCTGCCAGTTCACTCTCGTTAAGCGTGGTCCTGGTGACATCCGCGCGTTTTAAAACCTCATCGACCAGCATATGGGTAAAATCTTCGCCACCAAGAAAGTTGTCGcc                                                                              >  1:18366/1‑149 (MQ=255)
                                                                              tt‑gCTACATTTTGCCGCTTCCACACAGGCGTACAAGGCTGCCAGTTAACTCTCGTTTAGCGTGGTGCTGGCGACATCCGCGCGTTTTAAAACCTCATCGACCAGCATATGGGTAAAATCTTCTCCACCAAGCAAGTTGTGGCCAGCGGa                                                                        >  2:133625/1‑149 (MQ=255)
                                                                                            ccGCTTCCGCACAGGCGTACAAGGATGAAGGGGCACTCTCGTTAAGCGTGGTCCTGGCGACATCCGCGCGTTTTAAAACCTCATCGACCAGCATATGGGTAAAATCTTCGCCACCAAGAAAGTTGTCGGCAGCGGAGGCGTGAACTGCa                                                           >  2:370955/1‑149 (MQ=255)
                                                                                               cTTCCACACAGGCGTACAAGGCTGCCAGTTCACTCTCGTTAAGCGTGGTCCTGGCGACATCCGCGCGTTTTAAAACCTCATCGACCAGCATATGGGTAAAATCTTCGCCACCAAGTAAGTTGTCGCCAGCGGAGGCGTGAACTTCAATc                                                        >  1:800525/1‑149 (MQ=255)
                                                                                                  ccACACAGGCGTACAAGGCTGCCAGTTCACTCTCGTTAAGCGTGGTCCTGGCGACATCCGCGCGTTTTAAAACCTCATCGACCAGCATATGGGTAAAATCTTCGCCACCAAGAAAGTTGTCGCCAGCGGAGGCGTGAACTTCAATCAcc                                                     <  1:553554/149‑1 (MQ=255)
                                                                                                          gCGTACAAGGCTGCCAGTTCACTCTCGTTAAGCGTGGTCCTGGCGACATCCGCGCGTTTTAAAACCTCATCGACCAGCATATGGGTAAAATCTTCGCCACCAAGAAAGTTGTCGCCAGCGGAGGCGTGAACTTCAATCACCggcgtggc                                             >  2:691572/1‑149 (MQ=255)
                                                                                                          gCGTACAAGGCTGCCAGTTCACTCTCGTTAAGCGTGGTCCTGGAGACATCCGCGCGTGTTAAAACCTCATCGACCAGCATATGGGTAAAATCTTCGCCACCAAGAAAGTGGTCGCCAGGGGAGGCGTGAACTTCAATCACCggcgtggc                                             >  1:321992/1‑149 (MQ=255)
                                                                                                           cGTACAAGGCTGCCAGTTCACTCTCGTTAAGCGTGGTCCTGGTGACATCCGCGCGTTTTAAAACCTCATCGACCAGCATATGGGTAAAATCTTCGCCACCAAGAAAGTTGTCGCCAGCGGAGGCGTGAACTTCAATCACCggcgtggcg                                            <  2:415645/149‑1 (MQ=255)
                                                                                                           cGTACAAGGCTGCCAGTTCACTCTCGTTAAGCGTGGTCCTGGTGACATCCGCGCGTTTTAAAACCTCATCGACCAGCATATGGGTAAAATCTTCGCCACCAAGAAAGTTGTCGCCAGCGGAGGCGTGAACTTCAATCACCggcgtggcg                                            <  2:18366/149‑1 (MQ=255)
                                                                                                             tACAAGGCTGCCAGTTCACTCTCGTTAAGCGTGGTCCTGGCGACATCCGCGCGTTTTAAAACCTCATCGACCAGCACATGGGTAAAATCTTCGCCACCAAGAAAGTTGTCGCCAGCGGAGGCGTGAACTCCAATCACCggcggggcgtg                                          >  2:643447/1‑148 (MQ=255)
                                                                                                                    cTGCCAGTTCACTCTCGTTAAGCGTGACACGGGCGAGATCCGCGCGTTTTAAAACCTCATCGACCAGCATATGGGTAAAATCTTCGCCACCAAGAAAGTTGTCGCCAGCGGAGGCGTGAACTTCAATCACCGGCGTGGCGTACTCAAGc                                   >  2:570240/1‑149 (MQ=255)
                                                                                                                        cAGTTCACTCTCGTTAAGCGTGGTCCTGGCTACATCCGCGCGTTTTAAAACCTCATCGACCAGCATATGGGTAAAATCTTCGCCACCAAGAAAGTTGTCGCCAGCGGAGGCGTGAACTTCAATCACCGGCGTGGCGTACTCAAGCACCg                               >  1:411878/1‑149 (MQ=255)
                                                                                                                              aCTCTCGTTAAGCGTGGTCCTGGCGACATCCGCGCGTTTTAAAACCTCATCGACCAGCATATGGGTAAAATCTTCGCCACCAAGAAAGTTGTCGCCAGCGGAGGCGTGAACTTCAATCACCGGCGTGGCGTACTCAAGCACCGTAACGt                         >  1:454779/1‑149 (MQ=255)
                                                                                                                                      tAAGCGTGGTCCTGGCGACATCCGCGCGTTTTAAAACCTCATCGACCAGCATATGGGTAAAATCTTCGCCACCAAGCAAGTTGTCGCGAGCGGAGGCGTAAACCACAATCACGGGCGTGGCTTCCCCCAGCGCCGTAACGCCAAACGTg                 >  2:726779/1‑149 (MQ=255)
                                                                                                                                        cgCGTGGTCCTGGCGACATCCGCGCGTTTTAAAACCTCCTCGACCAGCATATGGGTAAAATCTTCGCCACCAAGAAAGTTGTCGCCAGCGGAGGCGTGAACTTCAATCACCGGCGTGGCGTACTCAAGCACCGTAACGTCAAACGTgcc               <  2:463875/148‑1 (MQ=255)
                                                                                                                                        aGCGTGGTCCTGGCGACATCCGCGCGTTTTAAAACCTCATCGACCAGCATATGGGTAAAATCTTCGCCACCAAGAAAGTTGTCGCCAGCGGAGGCGTGAACTTCAATCACCGGCGTGGCGTACTCAAGCACCGTAACGTCAAACGTgcc               <  2:678502/149‑1 (MQ=255)
                                                                                                                                         gCGTGGTCCTGGCGACATCCGCGCGTTTTAAAACCTCATCGACCAGCATATGGGTAAAATCTTCGCCACCAAGAAAGTTGTCGCCAGCGGAGGCGTGAACTTCAATCACCGGCGTGGCGTACTCAAGCACCGTAACGTCAAACGTgccg              <  2:800525/149‑1 (MQ=255)
                                                                                                                                            gggTCCTGGCGACATCCGCGCGTTTTAAAACCTCATCGACCAGCATATGGGTAAAATCTTCGCCACCAAGAAAGTTGTCGCCAGCGGAGGCGTGAACTTCAATCACCGGCGTGGCGGCGTCGAGCACCGTAACGTCAAACGGgcccgca           >  2:794697/2‑145 (MQ=255)
                                                                                                                                               tCCTGGCGACATCCGCGCGTTTTAAAACCTCATCGACCAGCATATGGGTAAAATCTTCGCCACCAAGAAAGTTGTCGCCAGCGGAGTCGTGAACTTCAATCACCGGCGTGGCGTACTCAAGCACCGTCACGTCCAACGCGCCGCccccc        >  2:550983/1‑148 (MQ=255)
                                                                                                                                               tCCTGGCGACATCCGCGCGTTTTAAAACCTCATCGACCAGCATATGGGGAAAATCTTCGCCACCAAGAAAGTTGTCGCCAGCGGAGGCGTGAACTTCAATCCCCGGCGTGGCGGACTCAAGCAGCGTAACGTCAAACGTGCCGACACCg        >  2:548598/1‑149 (MQ=255)
                                                                                                                                                   ggCGACATCCGCGCGTTTTAAAACCTCATCGACCAGCATATGGGTAAAATCTTCGCCACCAAGAAAGTTGTCGCCAGCGGAGGCGTGACCTTCAATCACCGGCGTGGCGTACTCAAGCACCGTAACGTCAAACGTGCCGCCACCGAGAt    <  1:423261/149‑1 (MQ=255)
                                                                                                                                                     cGACATCCGCGCGTTTTAAAACCTCATCGACCAGCATATGGGTAAAATCTTCGCCACCAAGAAAGTTGTCGCCAGCGGAGGCGTGAACTTCAATCACCGGCGTGGCGTACTCAAGCACCGTAACGTCAAACGTGCCGCCACCGAGATCa  <  1:691572/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TCTTCCAGTTCGTTCTCGTAAAATTCGCATTCCCGCGTTTCTTCCTGATACTGCCAGCGAATGTGCAATGGCGATTGATT‑GCTACATTTTGCCGCTTCCACACAGGCGTACAAGGCTGCCAGTTCACTCTCGTTAAGCGTGGTCCTGGCGACATCCGCGCGTTTTAAAACCTCATCGACCAGCATATGGGTAAAATCTTCGCCACCAAGAAAGTTGTCGCCAGCGGAGGCGTGAACTTCAATCACCGGCGTGGCGTACTCAAGCACCGTAACGTCAAACGTGCCGCCACCGAGATCA  >  NC_000913/682585‑682881

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑