Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,653,104 G→C 10.4% W204S (TGG→TCG) ‡ sseA → 3‑mercaptopyruvate sulfurtransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,653,1040GC10.4% 118.8 / 6.6 46W204S (TGG→TCG) ‡sseA3‑mercaptopyruvate sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base G (27/14);  new base C (3/2);  total (30/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.24e-01

TATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCAC  >  NC_000913/2652975‑2653244
                                                                                                                                 |                                                                                                                                            
tattGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTCCACGGAAcattca                                                                                                                                 <  1:631070/143‑6 (MQ=255)
tattCGCAACCCATGAAACTACAGCGCACATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGCTCATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTCCACGGAAcattca                                                                                                                                 <  1:14297/143‑6 (MQ=255)
  ttGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGGTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCGGGAGCACTGAATGTTACGTGGACGGAACTgggacg                                                                                                                               >  2:268407/1‑139 (MQ=255)
          ccATGAAAATACTGCGCAAAGTAGTGATGCGCGCCCGGCTGCACGTCTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCg                                                                                                                       <  2:126595/143‑1 (MQ=255)
             tGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAAc                                                                                                                    >  2:624600/1‑143 (MQ=255)
              gAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACt                                                                                                                   >  2:799032/1‑143 (MQ=255)
              gAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACt                                                                                                                   >  1:612989/1‑143 (MQ=255)
                   tACGGCGCAAATTATGGATGCCCGCCCGGCTGCACGTTTTAGCGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTaaaaa                                                                                                              <  2:507443/143‑1 (MQ=255)
                                 gtgctcgccgcccgGCTGGACGTTTTGACGCAGGAGTTCATGAACCTCGCCCAGGCGTACGTCGCGGACATATTCCCGGAGCATTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAACCGACCGATGAACTg                                                                                                <  2:707775/136‑1 (MQ=255)
                                    aTGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGAt                                                                                             >  2:616589/1‑143 (MQ=255)
                                    aTGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGAGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGAt                                                                                             >  1:238631/1‑143 (MQ=255)
                                     tGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATg                                                                                            <  1:439307/143‑1 (MQ=255)
                                     tGCCCGCCCGGCTGCACGTTTCAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATg                                                                                            <  1:443436/143‑1 (MQ=255)
                                               gCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATAttttt                                                                                  <  2:738738/143‑1 (MQ=255)
                                               gCTGCACGTTGTATCGCAGAAGTGGCGGGACCTCGCCCAGGCTTACGTCGCGGACATATTCCCGGAGCACTGAATGTCCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATAttttt                                                                                  <  2:801876/143‑1 (MQ=255)
                                                    aCGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGCAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAACACGACCGCTGAACTGGCTGCGATATTTTTTGGTc                                                                             >  2:510227/1‑143 (MQ=255)
                                                         ttAACGCAGAAGCTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCgcggc                                                                        <  2:209868/143‑1 (MQ=255)
                                                           aaCGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGGACTGGATGCGATATTTTTTGGTCGCGGCGt                                                                      >  1:686841/1‑143 (MQ=255)
                                                           aaCGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTTGTCGCGGCGt                                                                      >  1:178619/1‑143 (MQ=255)
                                                                aGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGGCGGAACTGGTGCGCGAAGGCGAACCAAACACGCTCGAGGAACTGGATGCGATAGTGTTTGCTCGCGGCGTCAGCt                                                                 >  1:670184/1‑143 (MQ=255)
                                                                aGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTCGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGCTGCGATATTTTTTGGTCGCGGCGTCAGCt                                                                 >  1:240024/1‑143 (MQ=255)
                                                                  aaGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGCACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATGTTTTTTGGTCGCGGCGTCGGCTAc                                                               >  2:75907/1‑143 (MQ=255)
                                                                    gTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTAAGGGGTCGCGGCGTCAGCTACGa                                                             >  2:522646/1‑143 (MQ=255)
                                                                        atgaACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACaaa                                                         >  1:417598/1‑143 (MQ=255)
                                                                         tgaACCTCGCCCAGCTTTACGTCGCGGACATATTCCCGGAGCACTGCATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACAGATGAACTGGATGCGACATTTTTTGGTCGCGGCGTCCGCTACGACAAAc                                                        >  2:349005/1‑143 (MQ=255)
                                                                            aCCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGGCAGCTACGACAAGCCag                                                     >  1:424994/1‑142 (MQ=255)
                                                                                 gcgcAGGTTTACGTCGCCGATATATTCCCGGAGCACTGACGGTTTCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATc                                                <  2:122600/140‑1 (MQ=255)
                                                                                 gCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTTAATGTTCCGTGGACGGAACTGGTGCGAGACCGGGAACTAAAAACGACCGATGAACTGGATGCAATATTGTGTCGTCGCGCCGTCTGCTACGACAAACCAATTTTc                                                >  2:674613/1‑143 (MQ=255)
                                                                                     nggTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGAGAAGGCGAACTAAAAACGACCGATGAACTGGACGCGCTATCCTTTGGTCGCGGCGTCAGGTACGCAAAACCAGTGATCggcg                                            >  1:67809/2‑140 (MQ=255)
                                                                                       gTTTACGTCGCGGACATAGTCCCGGAGCACAGAATGTTCCGTGGACGGAACTGGTGAGCGAAGGCGAACTAAAAACGACAGATGAACTGGACGCGATATTATTTGGTCGCGGCGGCAGCTACGACCAACCCATTAGCGGCacc                                          >  2:42080/1‑141 (MQ=255)
                                                                                        tttACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCt                                         >  1:691990/1‑143 (MQ=255)
                                                                                            cGTCGCGGACATATGCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCgg                                     <  2:69423/143‑1 (MQ=255)
                                                                                                gcgGACATATTCCCGGAGCACTGAATGGTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGctct                                 <  1:560907/143‑1 (MQ=255)
                                                                                                         ttCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCGCTGGTGTAGCg                        >  1:252863/1‑143 (MQ=255)
                                                                                                           cccGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGc                      <  2:417598/143‑1 (MQ=255)
                                                                                                            ccGGAGCACTGAATGTTCCGTGGACGGAACTCGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCa                     <  1:624600/143‑1 (MQ=255)
                                                                                                                  cACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGACCGCCGta               >  2:381085/1‑142 (MQ=255)
                                                                                                                    cTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGt             <  1:396326/143‑1 (MQ=255)
                                                                                                                    cTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGt             >  2:83166/1‑143 (MQ=255)
                                                                                                                    cTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGAGGAACTGGACGCGACAGTTTTTGGTCGCGGCGTCAGCTACGAGCAACCAATTATCGTCAGCTGCGCCGCGGGTGAAACGCCCGCCGCGGt             >  2:343238/1‑143 (MQ=255)
                                                                                                                     tGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTGACCGCAGCCGTGGtt            >  1:59905/1‑143 (MQ=255)
                                                                                                                       ccagTTCCGTCCACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGtttt          >  2:415556/4‑143 (MQ=255)
                                                                                                                       ccagTTCCGTCCACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGtttt          >  2:122317/4‑143 (MQ=255)
                                                                                                                           ttCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGCCTCTGGTGTCACGGCAGCCGTGGTTTTGATa      >  2:343134/1‑143 (MQ=255)
                                                                                                                              cGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCa   >  1:599744/1‑143 (MQ=255)
                                                                                                                               gTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTGTGTTAGCAc  >  1:317842/1‑143 (MQ=255)
                                                                                                                                 |                                                                                                                                            
TATTGGCAAGCCATGAAAATACGGCGCAAATTATTGATGCCCGCCCGGCTGCACGTTTTAACGCAGAAGTTGATGAACCTCGCCCAGGTTTACGTCGCGGACATATTCCCGGAGCACTGAATGTTCCGTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCAC  >  NC_000913/2652975‑2653244

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