Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,193,973 A→C 8.5% W106G (TGG→GGG)  nudJ ← bifunctional thiamine pyrimidine pyrophosphate hydrolase and thiamine pyrophosphate hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,193,9730AC 8.5% 155.0 / 7.2 47W106G (TGG→GGG) nudJbifunctional thiamine pyrimidine pyrophosphate hydrolase and thiamine pyrophosphate hydrolase
Reads supporting (aligned to +/- strand):  ref base A (19/24);  new base C (2/2);  total (21/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.76e-01

TTAGATGACACCCTTTGTAAAAGGCCAGTTAAAATCGCCAATCATCTCCAGCGGATAACGTTGCCCGCTTTGATAACAACGAATACTTTCCGCCACCAGCGGCGAACGAAGATTTGACGCCTGTAAAATTTCTTCGGCGCTGACCCAACGGCAGCAGTCGATATCGCTGTCATGAGGCTGGGTCGGGCATATTTGCTCAAGCTCAATGGCAAAGAGGAAACGCAAAAACGGCGTTTTATCTGGCGCAATCCACTGATGCATACGAATAAAGTGTTGCGGCTGCGCGCTGATGCC  >  NC_000913/1193827‑1194120
                                                                                                                                                  |                                                                                                                                                   
ttAGATGACACCCTTTGTAAAAGGCCAGTTAAAATCGCCAATCATCTCCAGCGGATAACGTTGCCCGCTTTGATAACAACGAATACTTTCCGCCACCAGCGGCGAACGAAGATTTGACGCCTGTAAAATTTCTTCGGCGCTGACCCAAc                                                                                                                                                   <  2:413591/149‑1 (MQ=255)
    aTGACACCCTTTGTAAAAGGCCAGTTAAAATCGCCAATCATCTCCAGCGGATAACGTTGCCCGCTTTGATAACAACGAATACTTTCCGCCACCAGCGGCGAACGAAGATTTGACGCCTGTAAAATTTCTTCGGCGCTGACCCAACGgca                                                                                                                                               <  4:190260/149‑1 (MQ=255)
      gACACCCTTTGTAAAAGGCCAGTTAAAATCGCCAATCATCTCCAGCGGATAACGTTGCCCGCTTTGATAACAACGAATACTTTCCGCCACCAGCGGCGAACGAAGATTTGACGCCTGTAAAATTTCTTCGGCGCTGACCCAACGgcagc                                                                                                                                             >  4:38225/1‑149 (MQ=255)
       acacCCTTTGTAAAAGGCCAGTTAAAATCGCCAATCATCTCCAGCGGATAACGTTGCCCGCTTTGATAACAACGAATACTTTCCGCCACCAGCGGCGAACGAAGATTTGACGCCTGTAAAATTTCTTCGGCGCTGACCCAACGgcagca                                                                                                                                            >  3:500393/1‑149 (MQ=255)
         acCCTTTGTAAAAGGCCAGTTAAAATCGCCAATCATCTCCAGCGGATAACGTTGCCCGCTTTGATAACAACGAATACTTTCCGCCACCAGCGGCGAACGAAGATTTGACGCCTGTAAAATTTCTTCGGCGCTGACCCCACGGCAGCAGt                                                                                                                                          >  3:503259/1‑149 (MQ=255)
         acCCTTTGTAAAAGGCCAGTTAAAATCGCCAATCATCTCCAGCGGATAACGTTGCCCGCTTTGATAACAACGAATACTTTCCGCCACCAGCGGCGAACGAAGATTTGACGCCTGTAAAATTTCTTCGGCGCTGACCCCACGGCAGCAGt                                                                                                                                          >  3:412201/1‑149 (MQ=255)
          cccTTTGTAAAAGGCCAGTTAAAATCGCCAATCATCTCCAGCGGATAACGTTGCCCGCTTTGATAACAACGAATACTTTCCGCCACCAGCGGCGAACGAAGATTTGACGCCTGTAAAATTTCTTCGGCGCTGACCCAACGGCAGCAGTc                                                                                                                                         >  1:22521/1‑149 (MQ=255)
               tGTAAAAGGCCAGTTAAAATCGCCAATCATCTCCAGCGGATAACGTTGCCCGCTTTGATAACAACGAATACTTTCCGCCACCAGCGGCGAACGAAGATTTGACGCCTGTAAAATTTCTTCGGCGCTGACCCAACGGCAGCAGTCGatat                                                                                                                                    <  3:227238/149‑1 (MQ=255)
                   aaaGGCCAGTTAAAATCGCCAATCATCTCCAGCGGATAACGTTGCCCGCTTTGATAACAACGAATACTTTCCGCCACCAGCGGCGAACGAAGATTTGACGCCTGTAAAATTTCTTCGGCGCTGACCCAACGGCAGCAGTCGATATCGCt                                                                                                                                >  1:26269/1‑149 (MQ=255)
                       gCCAGTTAAAATCGCCAATCATCTCCAGCGGATAACGTTGCCCGCTTTGATAACAACGAATACTTTCCGCCACCAGCGGCGAACGAAGATTTGACGCCTGTAAAATTTCTTCGGCGCTGACCCAACGGCAGCAGTCGATATCGCTGTCa                                                                                                                            <  1:17306/149‑1 (MQ=255)
                        ccAGTTAAAATCGCCAATCATCTCCAGCGGATAACGTTGCCCGCTTTGATAACAACGAATACTTTCCGCCACCAGCGGCGAACGAAGATTTGACGCCTGTAAAATTTCTTCGGCGCTGACCCAACGGCAGCAGTCGATATCGCTGTCAt                                                                                                                           <  2:26269/149‑1 (MQ=255)
                           gTTAAAATCGCCAATCATCTCCAGCGGATAACGTTGCCCGCTTTGATAACAACGAATACTTTCCGCCACCAGCGGCGAACGAAGATTTGACGCCTGTAAAATTTCTTCGGCGCTGACCCAACGGCAGCAGTCGATATCGCTGTCATGAg                                                                                                                        >  2:321925/1‑149 (MQ=255)
                             tAAAATCGCCAATCATCTCCAGCGGATAACGTTGCCCGCTTTGATAACAACGAATACTTTCCGCCACCAGCGGCGAACGAAGATTTGACGCCTGTAAAATTTCTTCGGCGCTGACCCAACGGCAGCAGTCGATATCGCTGTCATGAGGc                                                                                                                      <  1:321925/149‑1 (MQ=255)
                                   cGCCAATCATCTCCAGCGGATAACGTTGCCCGCTTTGATAACAACGAATACTTTCCGCCACCAGCGGCGAACGAAGATTTGACGCCTGTAAAATTTCTTCGGCGCTGACCCAACGGCAGCAGTCGATATCGCTGTCATGAGGCTGGGTc                                                                                                                <  1:139762/149‑1 (MQ=255)
                                       aaTCATCTCCAGCGGATAACGTTGCCCGCTTTGATAACAACGAATACTTTCCGCCACCAGCGGCGAACGAAGATTTGACGCCTGTAAAATTTCTTCGGCGCTGACCCAACGGCAGCAGTCGATATCGCTGTCATGAGGCTGGGTCGGGc                                                                                                            >  2:357738/1‑149 (MQ=255)
                                          catcTCCAGCGGATAACGTTGCCCGCTTTGATAACAACGAATACTTTCCGCCACCAGCGGCGAACGAAGATTTGACGCCTGTAAAATTTCTTCGGCGCTGACCCAACGGCAGCAGTCGATATCGCTGTCATGAGGCTGGGTCGGGCata                                                                                                         <  3:360982/149‑1 (MQ=255)
                                           atcTCCAGCGGATAACGTTGCCCGCTTTGATAACAACGAATACTTTCCGCCACCAGCGGCGAACGAAGATTTGACGCCTGTAAAATTTCTTCGGCGCTGACCCAACGGCAGCAGTCGATATCGCTGTCATGAGGCTGGGTCGGGCatat                                                                                                        >  1:20552/1‑149 (MQ=255)
                                             ctcCAGCGGATAACGTTGCCCGCTTTGATAACAACGAATACTTTCCGCCACCAGCGGCGAACGAAGATTTGACGCCTGTAAAATTTCTTCGGCGCTGACCCAACGGCAGCAGTCGATATCGCTGTCATGAGGCTGGGTCGGGCATAttt                                                                                                      <  2:337071/149‑1 (MQ=255)
                                                   cGGATAACGTTGCCCGCTTTGATAACAACGAATACTTTCCGCCACCAGCGGCGAACGAAGATTTGACGCCTGTAAAATTTCTTCGGCGCTGACCCAACGGCAGCAGTCGATATCGCTGTCATGAGGCTGGGTCGGGCATATTTgctcaa                                                                                                <  1:621016/149‑1 (MQ=255)
                                                          cGTTGCCCGCTTTGATAACAACGAATACTTTCCGCCACCAGCGGCGAACGAAGATTTGACGCCTGTAAAATTTCTTCGGCGCTGACCCAACGGCAGCAGTCGATATCGCTGTCATGAGGCTGGGTCGGGCATATTTGCTCAAGCTCAAt                                                                                         <  4:500393/149‑1 (MQ=255)
                                                                   cTTTGATAACAACGAATACTTTCCGCCACCAGCGGCGAACGAAGATTTGACGCCTGTAAAATTTCTTCGGCGCTGACCCAACGGCAGCAGTCGATATCGCTGTCATGAGGCTGGGTCGGGCATATTTGCTCAAGCTCAATGGCAAagag                                                                                <  2:22521/149‑1 (MQ=255)
                                                                          aacaacGAATACTTTCCGCCACCAGCGGCGAACGAAGATTTGACGCCTGTAAAATTTCTTCGGCGCTGACCCAACGGCAGCAGTCGATATCGCTGTCATGAGGCTGGGTCGGGCATATTTGCTCAAGCTCAATGGCAAAGAGGAAACGc                                                                         >  3:34933/1‑149 (MQ=255)
                                                                             aacGAATACTTTCCGCCACCAGCGGCGAACGAAGATTTGACGCCTGTAAAATTTCTTCGGCGCTGACCCAACGGCAGCAGTCGATATCGCTGTCATGAGGCTGGGTCGGGCATATTTGCTCAAGCTCAATGGCAAAGAGGAAACGCaaa                                                                      >  1:29276/1‑149 (MQ=255)
                                                                                gAATACTTTCCGCCACCAGCGGCGAACGAAGATTTGACGCCTGTAAAATTTCTTCGGCGCTGACCCAACGGCAGCAGTCGATATCGCTGTCATGAGGCTGGGTCGGGCATATTTGCTCAAGCTCAATGGCAAAGAGGAAACGCAAAAAc                                                                   >  3:98792/1‑149 (MQ=255)
                                                                                     cTTTCCGCCACCAGCGGCGAACGAAGATTTGACGCCTGTAAAATTTCTTCGGCGCTGACCCAACGGCAGCAGTCGATATCGCTGTCATGAGGCTGGGTCGGGCATATTTGCTCAAGCTCAATGGCAAAGAGGAAACGCAAAAACGGCGt                                                              <  3:119303/149‑1 (MQ=255)
                                                                                      tttCCGCCACCAGCGGCGAACGAAGATTTGACGCCTGTAAAATTTCTTCGGCGCTGACCCAACGGCAGCAGTCGATATCGCTGt                                                                                                                              >  3:265645/1‑84 (MQ=255)
                                                                                      tttCCGCCACCAGCGGCGAACGAAGATTTGACGCCTGTAAAATTTCTTCGGCGCTGACCCAACGGCAGCAGTCGATATCGCTGt                                                                                                                              <  4:265645/84‑1 (MQ=255)
                                                                                      tttCCGCCACCAGCGGCGAACGAAGATTTGACGCCTGTAAAATTTCTTCGGCGCTGACCCAACGGCAGCAGTCGATATCGCTGTCATGAGGCTGGGTCGGGCATATTTGCTCAAGCTCAATGGCAAAGAGGAAACGCAAAAACGGCGtt                                                             >  1:206766/1‑149 (MQ=255)
                                                                                      tttCCGCCACCAGCGGCGAACGAAGATTTGACGCCTGTAAAATTTCTTCGGCGCTGACCCAACGGCAGCAGTCGATATCGCTGTCATGAGGCTGGGTCGGGCATATTTGCTCAAGCTCAATGGCAAAGAGGAAACGCAAAAACGGCGtt                                                             <  2:357500/149‑1 (MQ=255)
                                                                                                caGCGGCGAACGAAGATTTGACGCCTGTAAAATTTCTTCGGCGCTGACCCAACGGCAGCAGTCGATATCGCTGTCATGAGGCTGGGTCGGGCATATTTGCTCAAGCTCAATGGCAAAGAGGAAACGCAAAAACGGCGTTTTATCTGgcg                                                   <  3:38225/149‑1 (MQ=255)
                                                                                                        aacgaaGATTTGACGCCTGTAAAATTTCTTCGGCGCTGACCCAACGGCAGCAGTCGATATCGCTGTCATGAGGCTGGGTCGGGCATATTTGCTCAAGCTCAATGGCAAAGAGGAAACGCAAAAACGGCGTTTTATCTGGCGCAATCCAc                                           <  2:399950/149‑1 (MQ=255)
                                                                                                         acgaaGATTTGACGCCTGTAAAATTTCTTCGGCGCTGACCCAACGGCAGCAGTCGATATCGCTGTCATGAGGCTGGGTCGGGCATATTTGCTCAAGCTCAATGGCAAAGAGGAAACGCAAAAACGGCGTTTTATCTGGCGCAATCCACt                                          <  1:235853/149‑1 (MQ=255)
                                                                                                           gaaGATTTGACGCCTGTAAAATTTCTTCGGCGCTGACCCAACGGCAGCAGTCGATATCGCTGTCATGAGGCTGGGTCGGGCATATTTGCTCAAGCTCAATGGCAAAGAGGAAACGCAAAAACGGCGTTTTATCTGGCGCAATCCACTGa                                        >  1:538548/1‑149 (MQ=255)
                                                                                                                 ttGACGCCTGTAAAATTTCTTCGGCGCTGACCCAACGGCAGCAGTCGATATCGCTGTCATGAGGCTGGGTCGGGCATATTTGCTCAAGCTCAATGGCAAAGAGGAAACGCAAAAACGGCGTTTTATCTGGCGCAATCCACTGATGCATa                                  <  4:520535/149‑1 (MQ=255)
                                                                                                                 ttGACGCCTGTAAAATTTCTTCGGCGCTGACCCAACGGCAGCAGTCGATATCGCTGTCATGAGGCTGGGTCGGGCATATTTGCTCAAGCTCAATGGCAAAGAGGAAACGCAAAAACGGCGTTTTATCTGGCGCAATCCACTGATGCATa                                  <  2:29276/149‑1 (MQ=255)
                                                                                                                     cGCCTGTAAAATTTCTTCGGCGCTGACCCAACGGCAGCAGTCGATATCGCTGTCATGAGGCTGGGTCGGGCATATTTGCTCAAGCTCAATGGCAAAGAGGAAACGCAAAAACGGCGTTTTATCTGGCGCAATCCACTGATGCATACGaa                              <  2:20552/149‑1 (MQ=255)
                                                                                                                      gCCTGTAAAATTTCTTCGGCGCTGACCCAACGGCAGCAGTCGATATCGCTGTCATGAGGCTGGGTCGGGCATATTTGCTCAAGCTCAATGGCAAAGAGGAAACGCAAAAACGGCGTTTTATCTGGCGCAATCCACTGATGCATACGAAt                             >  2:442883/1‑149 (MQ=255)
                                                                                                                              aaTTTCTTCGGCGCTGACCCAACGGCAGCAGTCGATATCGCTGTCATGAGGCTGGGTCGGGCATATTTGCTCAAGCTCAATGGCAAAGAGGAAACGCAAAAACGGCGTTTTATCTGGCGCAATCCACTGATGCATACGAATAAAGTGtt                     >  2:256250/1‑149 (MQ=255)
                                                                                                                                    ttcGGCGCTGACCCAACGGCAGCAGTCGATATCGCTGTCATGAGGCTGGGTCGGGCATATTTGCTCAAGCTCAATGGCAAAGAGGaaa                                                                            >  3:222666/1‑88 (MQ=255)
                                                                                                                                    ttcGGCGCTGACCCAACGGCAGCAGTCGATATCGCTGTCATGAGGCTGGGTCGGGCATATTTGCTCAAGCTCAATGGCAAAGAGGaaa                                                                            <  4:222666/88‑1 (MQ=255)
                                                                                                                                    ttcGGCGCTGACCCAACGGCAGCAGTCGATATCGCTGTCATGAGGCTGGGTCGGGCATATTTGCTCAAGCTCAATGGCAAAGAGGAAACGCAAAAACGGCGTTTTATCTGGCGCAATCCACTGATGCATACGAATAAAGTGTTGCGGCt               >  4:576283/1‑149 (MQ=255)
                                                                                                                                      cGGCGCTGACCCAACGGCAGCAGTCGATATCGCTGTCATGAGGCTGGGTCGGGCATATTTGCTCAAGCTCAATGGCAAAGAGGAAACGCAAAAACGGCGTTTTATCTGGCGCAATCCACTGATGCATACGAATAAAGTGTTGCGGCTgc             <  3:75746/149‑1 (MQ=255)
                                                                                                                                           cTGACCCAACGGCAGCAGTCGATATCGCTGTCATGAGGCTGGGTCGGGCATATTTGCTCAAGCTCAATGGCAAAGAGGAAACGCAAAAACGGCGTTTTATCTGGCGCAATCCACTGATGCATACGAATAAAGTGTTGCGGCTGCGCGCt        <  4:34933/149‑1 (MQ=255)
                                                                                                                                             gACCCCACGGCAGCAGTCGATATCGCTGTCATGAGGCTGGGTCGGGCATATTTGCTCAAGCTCAATGGCAAAGAGGAAACGCAAAAACGGCGTTTTATCTGGCGCAATCCACTGATGCATACGAATAAAGTGTTGCGGCTGCGCGCTGa      <  4:412201/149‑1 (MQ=255)
                                                                                                                                             gACCCCACGGCAGCAGTCGATATCGCTGTCATGAGGCTGGGTCGGGCATATTTGCTCAAGCTCAATGGCAAAGAGGAAACGCAAAAACGGCGTTTTATCTGGCGCAATCCACTGATGCATACGAATAAAGTGTTGCGGCTGCGCGCTGa      <  4:503259/149‑1 (MQ=255)
                                                                                                                                                ccAACGGCAGCAGTCGATATCGCTGTCATGAGGCTGGGTCGGGCATATTTGCTCAAGCTCAATGGCAAAGAGGAAACGCAAAAACGGCGTTTTATCTGGCGCAATCCACTGATGCATACGAATAAAGTGTTGCGGCTGCGCGCTGATGc   >  1:153943/1‑149 (MQ=255)
                                                                                                                                                 cAACGGCAGCAGTCGATATCGCTGTCATGAGGCTGGGTCGGGCATATTGGCTCAAGCTCAATGGCAAAGAGGAAACGCAAAAACGGCGTTTTATCTGGCGCAATCCACTGATGCATACGAATAAAGTGTTGCGGCTGCGCGCTGATGcc  >  4:452342/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TTAGATGACACCCTTTGTAAAAGGCCAGTTAAAATCGCCAATCATCTCCAGCGGATAACGTTGCCCGCTTTGATAACAACGAATACTTTCCGCCACCAGCGGCGAACGAAGATTTGACGCCTGTAAAATTTCTTCGGCGCTGACCCAACGGCAGCAGTCGATATCGCTGTCATGAGGCTGGGTCGGGCATATTTGCTCAAGCTCAATGGCAAAGAGGAAACGCAAAAACGGCGTTTTATCTGGCGCAATCCACTGATGCATACGAATAAAGTGTTGCGGCTGCGCGCTGATGCC  >  NC_000913/1193827‑1194120

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