Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 319,152 C→T 27.4% A284T (GCA→ACA)  rclA ← putative pyridine nucleotide‑disulfide oxidoreductase RclA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913319,1520CT27.4% 19.1 / 5.8 18A284T (GCA→ACA) rclAputative pyridine nucleotide‑disulfide oxidoreductase RclA
Reads supporting (aligned to +/- strand):  ref base C (6/7);  new base T (2/3);  total (8/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 = 5.40e-01

TTTTTCCGATCATCAGTACTACGTTTGCCTTCACCCAGTAACTCATCACGTACAATGCGGTAATCATCCAGTGATATGTAAGTAAATTGCAGCCCGCCGGTAACATCTCCCATCGCCCAAATATTGTCTGCGGTGGTATGTAATCGCTTGTCAACGACAATTGCCCCGCGCTCGTTTACTGCGATACCGGCATTTTCTGGATGTAACGAAGCGGTAGCCGGTTGACGACCGGAAGCTATTAACAGTGCATCCACCGCCAGTTGGGCGTG  >  NC_000913/319022‑319290
                                                                                                                                  |                                                                                                                                          
tttttCCGATCATCAGTACTACGTTTACCTTCACCCAGTAACTCATCACGTACAATGCGGTAATCATCGAGTGATATGTAAGTAAATTGCAGCCCGCCGGTAATATCTCCCATCGCCCAAATCTTGTCTGCGGTGGTATGTaa                                                                                                                                >  2:347010/1‑143 (MQ=255)
 ttttCCGACCATCAGTACTACGATAGCGTTCACCCGGCAATTAATCACATACAATGCGGTAATCATCCCGGGATATGTAAGTAAATTGCAGCCCGCCGGTAACATCTCCCATCGCCCAAATATTGTCTGCGGTGGTATGTAAt                                                                                                                               <  1:80173/143‑1 (MQ=255)
      ctatcatcAGTACCCCCTTCTCCATCACCTACTAATTCACCACTCACATTGCGCTAATCATCCAGTGATATTCACGTAAATTGCAGTACGCCTCCAACATCTCCCATCTCCCAAATATTGTCTGCGGTGGTATGTAATCGCtt                                                                                                                          <  2:8803/141‑1 (MQ=255)
              aGTACTACGTTTGCCTTCACCCAGTAACTCATCACGTATAATGCGGTAATCATCCAGTGATATGTAAGTAAATTGCAGCCCGCCGGTAACATCTCCCATCGCCCAAATATTGTCTGCGGTGGTATGTAATCGCTTGTCAACGa                                                                                                                  <  1:177064/143‑1 (MQ=255)
                                               aCGTACAATGCGGTAATCATCCAGTGATATGTAAGTAAATTGCAGCCCGCCGGTAACAGCTCCCATCGCCCAAATATTGTCTGCGGTGGTATGTAATCGCTTGTCAACGACAATTGCCCCGCGCTCGTTTACTGCGATACCgg                                                                                 <  1:147153/143‑1 (MQ=255)
                                                     aaTGCGGTAATCATCCAGCGATATGTAAGTAAATTGCAGCCCGCCGGTAACATCTCCCATCGCCCAACTATTGTCTGCGGTGGTATGTAATCACTTGTGAACGACAATAGCCACGCGATCAGATACTGCGATAACGACagatt                                                                           >  1:65683/1‑139 (MQ=255)
                                                                       tGATATGTAAGTAAATTGCAGCCCGCCGGTAACATCTCCCATCGCCCAAAAATTGTCTGCGGTGGTAGGTAATCGCATGGCAACGACAATTACCACGCGCTCGTGTGCTACGATACCGGCATTTTCTGGATAGAACGGAGCgg                                                         >  2:109679/1‑143 (MQ=255)
                                                                       tGATATGTAAGTAAATTGCAGCCCGCCGGTAACATCTCCCATCGCACAAATATTGTCTGCGGTGGTATGGAGTCGCTTGTCAAAGACAATTGCCCCGTGCTCGTGTACTCCGATACAGGCAATAACTGGATATAAAGaccagg                                                         >  2:80980/1‑138 (MQ=255)
                                                                             gtcggtacATTGCGGCCCGCCGGTAACATCTCCCATCGCCCGAATATTGTCTGCGGTGGTATGTAATCGCTTGTCAACGACAATTGCCCCGCGCTCGTGTACTGCGATACCGGCATTTTCTGGATGTAACGAAGCGGTAGCCg                                                   <  2:68252/139‑1 (MQ=255)
                                                                                       ttcAGCCCGCCGGTAACATCCCCCACCGCCCAAATATTGTCTGCGGTGGTATGTACTCGCTTGTCAACGACAATTGCCCCGCGCTCGTTTACTGCGATACCGGCATTTTCTGGATGTAACGAAGCGGTAGCCGGTTGACGAcc                                         <  1:68892/141‑1 (MQ=255)
                                                                                                       gagctcCCATCGCCCAAAAATTCTCTGTGGTGGTAGGAAAGCGCTTGTCAACGACAATTGCCCCGCGCTCGTTTAGTGCGATACCGGAATTTTCTGGATGTAACGAAGCGGTAGCCGGTTGACGACCGGAAGCTATTAACAGt                         <  2:254669/140‑1 (MQ=255)
                                                                                                       cATCTCCCATCGCCCAAATATTGTCTGTGGTGGTATGTAATCGCTTGTCAACGACAATTGCCCCGCGCTCGTTTACTGCGATACCGGCATTTTCTGGATGTAACGAAGCGGTAGc                                                     >  1:173595/1‑115 (MQ=255)
                                                                                                       cATCTCCCATCGCCCAAATATTGTCTGTGGTGGTATGTAATCGCTTGTCAACGACAATTGCCCCGCGCTCGTTTACTGCGATACCGGCATTTTCTGGATGTAACGAAGCGGTAGc                                                     <  2:171788/115‑1 (MQ=255)
                                                                                                       cATCTCCCATCGCCCAAATATTGTCTGTGGTGGTATGTAATCGCTTGTCAACGACAATTGCCCCGCGCTCGTTTACTGCGATACCGGCATTTTCTGGATGTAACGAAGCGGTAGc                                                     <  2:173595/115‑1 (MQ=255)
                                                                                                       cATCTCCCATCGCCCAAATATTGTCTGTGGTGGTATGTAATCGCTTGTCAACGACAATTGCCCCGCGCTCGTTTACTGCGATACCGGCATTTTCTGGATGTAACGAAGCGGTAGc                                                     >  1:171788/1‑115 (MQ=255)
                                                                                                          ctcCCATCGCCCAAATATTGTCTGCGGTGGTATGTAATCGCTTGTCAACGACAATTGCCCCGCGCTCGTTTACTGCGATACCGGCCTTTTCTGGATGTACCGTAGGGCCACCCGCTTGGTGACCGGGAGCTAGTAACAGTGCa                      >  2:255105/1‑143 (MQ=255)
                                                                                                                   cccAAATATTGTCTGCGGTGGTATGTAATCGCTTGTCAACGACAATTGCCCCGCGCTCGTTTACTGCGATACCGGCATTTTCTGGATGTAACGAAGCGGTAGCCGGTTGACGACCGGAAGCTATTAACAGTGCAACCACCGcc             >  2:199935/1‑143 (MQ=255)
                                                                                                                              gccgCGGTGGTATGTAATCGCTTGCCCACGACAATTGCCCCGCGCTCGTTTACTGCGATACCGGCATTTTCTGGATGTAACGAAGCGGTAGCCGGTTGACGACCGGAAGCTATTAACAGTGCATCCACCGCCAGTTGGGCGTg  <  1:128316/140‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                          
TTTTTCCGATCATCAGTACTACGTTTGCCTTCACCCAGTAACTCATCACGTACAATGCGGTAATCATCCAGTGATATGTAAGTAAATTGCAGCCCGCCGGTAACATCTCCCATCGCCCAAATATTGTCTGCGGTGGTATGTAATCGCTTGTCAACGACAATTGCCCCGCGCTCGTTTACTGCGATACCGGCATTTTCTGGATGTAACGAAGCGGTAGCCGGTTGACGACCGGAAGCTATTAACAGTGCATCCACCGCCAGTTGGGCGTG  >  NC_000913/319022‑319290

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