Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,657,874 A→G R354R (CGT→CGC)  hscA ← DnaK‑like molecular chaperone specific for IscU

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,657,8740AG93.5% 100.0 / ‑1.6 31R354R (CGT→CGC) hscADnaK‑like molecular chaperone specific for IscU
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (13/16);  total (14/17)
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.61e-01

CGTTTCGAGGCCCAGCGACAGTGGGATCACATCAAGCAACAGCATTTCGCTGTCTGGCTTGTTACCCACCAGAATATCCGCCTGAATCGCCGCGCCAATAGCGACGACTTTATCCGGGTCGATGGAAGTCAGCGGTGGACGACCGAAAAATTCGCCTACCCGTTCACGCACCAGCGGCACGCGAGTAGAACCGCCCACCATCACCACTTCCAGCACTTCATCAGCTTCTACACCCGCGTCTTTCAGCGCGCGACGACA  >  NC_000913/2657736‑2657993
                                                                                                                                          |                                                                                                                       
cGTTTCGAGGCCCAGCGACAGTGGGATCACATCAAGCAACAGCATTTCGCTGTCTGGCTTGTTACCCACCAGAATATCCGCCTGAATCGCCGCGCCAATAGCGACGACTTTATCCGGGTCGATGGAAGTCAGCGGTGga                                                                                                                         <  1:19031/139‑1 (MQ=255)
 gTTTCGAGGCCCAGCGACAGTGGGATCACATCAAGCAACAGCATTTCGCTGTCTGGCTTGTTACCCACCAGAATATCCGCCTGAATCGCCGCGCCAATAGCGACGACTTTATCCGGGTCGATGGAAGTCAGCGGTGgac                                                                                                                        >  2:260951/1‑139 (MQ=255)
             aGCGACAGTGGGATCACATCAAGCAACAGCATTTCGCTGTCTGGCTTGTTACCCACCAGAATATCCGCCTGAATCGCCGCGCCAATAGCGACGACTTTATCCGGGTCGATGGAAGTCAGCGGTGGGCGACCGaaaaa                                                                                                              >  2:133336/1‑137 (MQ=255)
             aGCGACAGTGGGATCAAATCAAGCAACAGCATTTCGCTGTCTGGCTTGTTACCCACCAGAATATCCGCCTGAATCGCCGCGCCAATAGCGACGACTTTATCCGGGTCGATGGAAGTCAGCGGTGGGCGACCGaaaaa                                                                                                              <  1:133336/137‑1 (MQ=255)
               cGACAGTGGGATCACATCAAGCAACAGCATTTCGCTGTCTGGCTTGTTACCCACCAGAATATCCGCCTGAATCGCCGCGCCAATAGCGACGACTTTATCCGGGTCGATGGAAGTCAGCGGTGGGCGACCGAAAAATTCg                                                                                                          >  2:206409/1‑139 (MQ=255)
                    gTGGGATCACATCAAGCAACAGCATTTCGCTGTCTGGCTTGTTACCCACCAGAATATCCGCCTGAATCGCCGCGCCAATAGCGACGACTTTATCCGGGTCGATGGAAGTCAGCGGTGGGCGACCGAAAAATTCGCCTAc                                                                                                     <  1:206409/139‑1 (MQ=255)
                             caTCAAGCAACAGCATTTCGCTGTCTGGCTTGTTACCCACCAGAATATCCGCCTGAATCGCCGCGCCAATAGCGACGACTTTATCCGGGTCGATGGAAGTCAGCGGTGGGCGACCGAAAAATTCGCCTACCCGTTCa                                                                                              >  1:10880/1‑137 (MQ=255)
                             caTCAAGCAACAGCATTTCGCTGTCTGGCTTGTTACCCACCAGAATATCCGCCTGAATCGCCGCGCCAATAGCGACGACTTTATCCGGGTCGATGGAAGTCAGCGGTGGGCGACCGAAAAATTCGCCTACCCGTTCa                                                                                              <  2:10880/137‑1 (MQ=255)
                                 aaGCAACAGCATTTCGCTGTCTGGCTTGTTACCCACCAGAATATCCGCCTGAATCGCCGCGCCAATAGCGACGACTTTATCCGGGTCGATGGAAGTCAGCGGTGGGCGACCGAAAAATTCGCCTACCCGTTCa                                                                                              >  2:658107/1‑133 (MQ=255)
                                 aaGCAACAGCATTTCGCTGTCTGGCTTGTTACCCACCAGAATATCCGCCTGAATCGCCGCGCCAATAGCGACGACTTTATCCGGGTCGATGGAAGTCAGCGGTGGGCGACCGAAAAATTCGCCTACCCGTTCa                                                                                              <  1:658107/133‑1 (MQ=255)
                                  aGCAACAGCATTTCGCTGTCTGGCTTGTTACCCACCAGAATATCCGCCTGAATCGCCGCGCCAATAGCGACGACTTTATCCGGGTCGATGGAAGTCAGCGGTGGGCGACCGAAAAATTCGCCTACCCGTTCACGCACCa                                                                                       >  1:450487/1‑139 (MQ=255)
                                      aCAGCATTTCGCTGTCTGGCTTGTTACCCACCAGAATATCCGCCTGAATCGCCGCGCCAATAGCGACGACTTTATCCGGGTCGATGGAAGTCAGCGGTGGGCGACCGAAAAATTCGCCTACCCGTTCACGCACCAGCgg                                                                                   >  1:34323/1‑139 (MQ=255)
                                         gCATTTCGCTGTCTGGCTTGTTACCCACCAGAATATCCGCCTGAATCGCCGCGCCAATAGCGACGACTTTATCCGGGTCGATGGAAGTCAGCGGTGGGCGACCGAAAAATTCGCCTACCCGTTCACGCACCAGCGGCAc                                                                                <  1:258259/139‑1 (MQ=255)
                                          cATTTCGCTGTCTGGCTTGTTACCCACCAGAATATCCGCCTGAATCGCCGCGCCAATAGCGACGACTTTATCCGGGTCGATGGAAGTCAGCGGTGGGCGACCGAAAAATTCGCCTACCCGTTCACGCACCAGCGGCAcg                                                                               >  2:658206/1‑139 (MQ=255)
                                                 cTGTCTGGCTTGTTACCCACCAGAATATCCGCCTGAATCGCCGCGCCAATAGCGACGACTTTATCCGGGTCGATGGAAGTCAGCGGTGGGCGACCGAAAAATTCGCCTACCCGTTCACGCACCAGCGGCACGCGAGTAg                                                                        <  1:220059/139‑1 (MQ=255)
                                                       ggCTTGTTACCCACCAGAATATCCGCCTGAATCGCCGCGCCAATAGCGACGACTTTATCCGGGTCGATGGAAGTCAGCGGTGGGCGACCGAAAAATTCGCCTACCCGTTCACGCACCAGCGGCACGCGAGTAGAACCGc                                                                  <  2:49415/139‑1 (MQ=255)
                                                           tGTTACCCACCAGAATATCCGCCTGAATCGCCGCGCCAATAGCGACGACTTTATCCGGGTCGATGGAAGTCAGCGGTGGGCGACCGAAAAATTCGCCTACCCGTTCACGCACCAGCGGCACGCGAGTAGAACCGCccac                                                              <  2:41981/139‑1 (MQ=255)
                                                            gTTACCCACCAGAATATCCGCCTGAATCGCCGCGCCAATAGCGACGACTTTATCCGGGTCGATGGAAGTCAGCGGTGGGCGACCGAAAAATTCGCCTACCCGTTCACGCACCAGCGGCACGCGAGTAGAACCGCccacc                                                             <  2:501726/139‑1 (MQ=255)
                                                                  caccaGAATATCCGCCTGAATCGCCGCGCCAATAGCGACGACTTTATCCGGGTCGATGGAAGTCAGCGGTGGGCGACCGAAAAATTCGCCTACCCGTTCACGCACCAGCGGCACGCGAGTAGAACCGCCCACCATcacc                                                       <  1:228836/139‑1 (MQ=255)
                                                                     caGAATATCCGCCTGAATCGCCGCGCCAATAGCGACGACTTTATCCGGGTCGATGGAAGTCAGCGGTGGGCGACCGAAAAATTCGCCTACCCGTTCACGCACCAGCGGCACGCGAGTAGAACCGCCCACCATCACCACt                                                    >  1:37968/1‑139 (MQ=255)
                                                                       gAATATCCGCCTGAATCGCCGCGCCAATAGCGACGACTTTATCCGGGTCGATGGAAGTCAGCGGTGGGCGACCGAAAAATTCGCCTACCCGTTCACGCACCAGCGGCACGCGAGTAGAACCGCCCACCATCACCACTTc                                                  >  1:495523/1‑139 (MQ=255)
                                                                          tatCCGCCTGAATCGCCGCGCCAATAGCGACGACTTTATCCGGGTCGATGGAAGTCAGCGGTGGGCGACCGAAAAATTCGCCTACCCGTTCACGCACCAGCGGCACGCGAGTAGAACCGCCCACCATCACCACTTCCAg                                               <  2:564687/139‑1 (MQ=255)
                                                                                 cTGAATCGCCGCGCCAATAGCGACGACTTTATCCGGGTCGATGGAAGTCAGCGGTGGACGACCGAAAAATTCGCCTACCCGTTCACGCACCAGCGGCACGCGAGTAGAACCGCCCACCATCACCACTTCCAGCACTtca                                        <  1:260951/139‑1 (MQ=255)
                                                                                           gcgcCAATAGCGACGACTTTATCCGGGTCGATGGAAGTCAGCGGTGGGCGACCGAAAAATTCGCCTACCCGTTCACGCACCAGCGGCACGCGAGTAGAACCGCCCACCATCACCACTTCCAGCACTTCATCAGCTTCTa                              <  2:495523/139‑1 (MQ=255)
                                                                                                                     gTCGATGGAAGTCAGCGGTGGGCGACCGAAAAATTCGCCTACCCGTTCACGCACCAGCGGCACGCGAGTAGAACCGCCCACCATCACCACTTCCAGCACTtcatca                                     >  2:470049/1‑106 (MQ=255)
                                                                                                                     gTCGATGGAAGTCAGCGGTGGGCGACCGAAAAATTCGCCTACCCGTTCACGCACCAGCGGCACGCGAGTAGAACCGCCCACCATCACCACTTCCAGCACTtcatca                                     <  1:470049/106‑1 (MQ=255)
                                                                                                                     gTCGATGGAAGTCAGCGGTGGGCGACCGAAAAATTCGCCTACCCGTTCACGCACCAGCGGCACGCGAGTAGAACCGCCCACCATCACCACTTCCAGCACTTCATCAGCTTCTACACCCGCGTCTTTCAGCGCGcgacga    >  1:177563/1‑139 (MQ=255)
                                                                                                                     gTCGATGGAAGTCAGCGGTGGACGACCGAAAAATTCGCCTACCCGTTCACGCACCAGCGGCACGCGAGTAGAACCGCCCACCATCACCACTTCCAGCACTTCATCAGCTTCTACACCCGCGTCTTTCAGCGCGcgacga    >  2:364475/1‑139 (MQ=255)
                                                                                                                      tCGATGGAAGTCAGCGGTGGGCGACCGAAAAATTCGCCTACCCGTTCACGCACCAGCGGCACGCGAGTa                                                                         >  2:457869/1‑69 (MQ=255)
                                                                                                                      tCGATGGAAGTCAGCGGTGGGCGACCGAAAAATTCGCCTACCCGTTCACGCACCAGCGGCACGCGAGTa                                                                         <  1:457869/69‑1 (MQ=255)
                                                                                                                       cGATGGAAGTCAGCGGTGGGCGACCGAAAAATTCGCCTACCCGTTCACGCACCAGCGGCACGCGAGTAGAACCGCCCACCATCACCACTTCCAGCACTTCATCAGCTTCTACACCCGCGTCTTTCAGCGCGCGACGACa  <  2:34323/139‑1 (MQ=255)
                                                                                                                                   gCGGTGGGCGACCGAAAAATTCGCCTACCCGTTCACGCACCAGCGGCACGCGAGTAGAACCGCCCACCATCACCACTTCCAGCACTTCATCAGCTTCTACACCCGCGTCTTTCAGCGCGcgacg     <  2:73161/124‑1 (MQ=255)
                                                                                                                                   gCGGTGGGCGACCGAAAAATTCGCCTACCCGTTCACGCACCAGCGGCACGCGAGTAGAACCGCCCACCATCACCACTTCCAGCACTTCATCAGCTTCTACACCCGCGTCTTTCAGCGCGcgacg     >  1:73161/1‑124 (MQ=255)
                                                                                                                                          |                                                                                                                       
CGTTTCGAGGCCCAGCGACAGTGGGATCACATCAAGCAACAGCATTTCGCTGTCTGGCTTGTTACCCACCAGAATATCCGCCTGAATCGCCGCGCCAATAGCGACGACTTTATCCGGGTCGATGGAAGTCAGCGGTGGACGACCGAAAAATTCGCCTACCCGTTCACGCACCAGCGGCACGCGAGTAGAACCGCCCACCATCACCACTTCCAGCACTTCATCAGCTTCTACACCCGCGTCTTTCAGCGCGCGACGACA  >  NC_000913/2657736‑2657993

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