Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,369,793 G→A 100% E35E (GAG→GAA pspE → thiosulfate:cyanide sulfurtransferase (rhodanese)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,369,7930GA98.5% 214.6 / ‑4.3 67E35E (GAG→GAApspEthiosulfate:cyanide sulfurtransferase (rhodanese)
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (36/30);  total (36/31)
Fisher's exact test for biased strand distribution p-value = 4.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.37e-01

AAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTATTGC  >  NC_000913/1369696‑1369889
                                                                                                 |                                                                                                
aaaaaGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACa                                                                                                <  2:984183/100‑1 (MQ=255)
  aaaGGCTTACTTGCGCTGGCACTGGTGTTTTCACTGCTCGTTTTCGCAGCTGTACACTGGATCAATGTTCGTGTTCCAGAGCAGTATCAGCAAGAacac                                                                                               <  2:586813/99‑1 (MQ=255)
     ggCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTc                                                                                           <  2:2268807/100‑1 (MQ=255)
       cTTACTTGCTCTGGCACTGGTGTTTTCACTGCGCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAg                                                                                         >  1:1393913/1‑100 (MQ=255)
        ttACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAgg                                                                                        <  1:609550/100‑1 (MQ=255)
           cTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAggggg                                                                                     <  1:475471/100‑1 (MQ=255)
            ttGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGc                                                                                    <  1:2488976/100‑1 (MQ=255)
             tGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGCACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGcc                                                                                   >  2:1741425/1‑100 (MQ=255)
                tctGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATc                                                                                <  1:2493171/100‑1 (MQ=255)
                    gCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAata                                                                            >  1:414373/1‑100 (MQ=255)
                       cTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTc                                                                         >  1:1668744/1‑100 (MQ=255)
                         ggTGTTTTCACTGCCCGGTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTcc                                                                        <  2:2158664/99‑1 (MQ=255)
                          gtgtTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTcccc                                                                      <  2:338231/100‑1 (MQ=255)
                          gtgtTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTcccc                                                                      <  2:1257514/100‑1 (MQ=255)
                              tttCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTgaa                                                                  <  2:383077/100‑1 (MQ=255)
                                   cTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTgaaggaag                                                             >  2:1247561/1‑100 (MQ=255)
                                       ccGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGaa                                                         >  1:435902/1‑100 (MQ=255)
                                        cGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGaaa                                                        >  1:606713/1‑100 (MQ=255)
                                        cGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGAACCATCAATATTCCCCTGAAGGAAGTGaaa                                                        >  2:2297903/1‑100 (MQ=255)
                                         gTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAag                                                       <  1:154228/100‑1 (MQ=255)
                                             tCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAgcg                                                   <  2:1473436/100‑1 (MQ=255)
                                             tCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAgcg                                                   >  1:2241535/1‑100 (MQ=255)
                                              cgccgcTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAgcgc                                                  >  1:1106867/1‑100 (MQ=255)
                                               gccgcTGAACACTGGATCGATGTTCGTGTTCCAGATCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCa                                                 <  2:399854/100‑1 (MQ=255)
                                                ccgcTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCAt                                                >  1:1132787/1‑100 (MQ=255)
                                                    tGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGcc                                            >  1:1841872/1‑100 (MQ=255)
                                                     gAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCa                                           <  1:1832021/100‑1 (MQ=255)
                                                     gAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCa                                           <  1:1827790/100‑1 (MQ=255)
                                                     gAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCa                                           >  2:441250/1‑100 (MQ=255)
                                                      aaCACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTTCCCTGAAGGAAGTGAAATAGCGCATTGCCAc                                          <  1:1725207/100‑1 (MQ=255)
                                                         acTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCAccgc                                       >  2:2118076/1‑100 (MQ=255)
                                                            ggATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGt                                    >  1:2008671/1‑100 (MQ=255)
                                                            ggATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGt                                    >  2:2118571/1‑100 (MQ=255)
                                                             gATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGtt                                   >  2:1131229/1‑100 (MQ=255)
                                                                  aTGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGa                              <  1:1586762/100‑1 (MQ=255)
                                                                   tGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGAt                             >  2:2091894/1‑100 (MQ=255)
                                                                   tGTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGAt                             <  1:1821771/100‑1 (MQ=255)
                                                                    gTTCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATa                            >  1:573557/1‑100 (MQ=255)
                                                                      tCGTGTTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCCATATTCACCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATaa                           <  2:2154325/99‑1 (MQ=255)
                                                                         tgtTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAAc                       <  1:1685/100‑1 (MQ=255)
                                                                         tgtTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAAc                       <  2:1942115/100‑1 (MQ=255)
                                                                          gtTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGACGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACg                      <  2:1563693/100‑1 (MQ=255)
                                                                          gtTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACg                      <  1:167271/100‑1 (MQ=255)
                                                                          gtTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACg                      >  1:529453/1‑100 (MQ=255)
                                                                          gtTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACg                      <  2:2008219/100‑1 (MQ=255)
                                                                          gtTCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACg                      >  2:2499485/1‑100 (MQ=255)
                                                                            tCCAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGac                    >  2:1593664/1‑100 (MQ=255)
                                                                             ccAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGaca                   >  1:72564/1‑100 (MQ=255)
                                                                             ccAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGaca                   >  2:255928/1‑100 (MQ=255)
                                                                             ccAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGaca                   >  1:209848/1‑100 (MQ=255)
                                                                              cAGAGCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGacac                  >  1:1542587/1‑100 (MQ=255)
                                                                               agagCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACAcc                 >  2:479393/1‑100 (MQ=255)
                                                                                gagCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCg                <  1:1593675/100‑1 (MQ=255)
                                                                                  gCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTg              <  1:1035475/100‑1 (MQ=255)
                                                                                  gCAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTg              >  2:520679/1‑100 (MQ=255)
                                                                                   cAGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATCCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGa             <  1:11861/100‑1 (MQ=255)
                                                                                    aGTATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGaa            >  2:73018/1‑100 (MQ=255)
                                                                                      tATCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAg          <  1:1131237/100‑1 (MQ=255)
                                                                                       aTCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAg          >  2:2045109/1‑99 (MQ=255)
                                                                                        tCAGCAAGAACACGTTCAGGGGGCCATCAATATTCGCCTGAAGGAAGTGAAAGAGCGCATTGCAACCGCCGTTCCGGATAAAAACGACACAGTGAAAgt         >  2:2242626/1‑99 (MQ=255)
                                                                                        tCAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAgtg        <  2:99833/100‑1 (MQ=255)
                                                                                         cAGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCCGATAAACACGACACCGTGAAAatg        >  2:669800/1‑99 (MQ=255)
                                                                                         cAGCAAGAACACGTTCAGGGAGCCATCAATATTCCGATGAAGGAAGTGAAAGAGCGCATTGCCACTGCCCATCCGGATAAAACCGACACCATGAAAgtgt       >  2:874638/1‑100 (MQ=255)
                                                                                          aGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTa      <  1:1240752/100‑1 (MQ=255)
                                                                                          aGCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTa      >  1:482938/1‑100 (MQ=255)
                                                                                           gCAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTAt     >  2:1685341/1‑100 (MQ=255)
                                                                                            cAAGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTAtt    >  1:324661/1‑100 (MQ=255)
                                                                                             aaGAACACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTATTg   <  2:1706333/100‑1 (MQ=255)
                                                                                              agaaCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTATTGc  >  1:1465891/1‑100 (MQ=255)
                                                                                              agaaCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTATTGc  <  1:2118585/100‑1 (MQ=255)
                                                                                                 |                                                                                                
AAAAAGGCTTACTTGCTCTGGCACTGGTGTTTTCACTGCCCGTTTTCGCCGCTGAACACTGGATCGATGTTCGTGTTCCAGAGCAGTATCAGCAAGAGCACGTTCAGGGGGCCATCAATATTCCCCTGAAGGAAGTGAAAGAGCGCATTGCCACCGCCGTTCCGGATAAAAACGACACCGTGAAAGTGTATTGC  >  NC_000913/1369696‑1369889

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