Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 72,253 C→T S225S (TCG→TCA thiQ ← thiamine ABC transporter ATP binding subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091372,2530CT96.8% 101.8 / NA 31S225S (TCG→TCAthiQthiamine ABC transporter ATP binding subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (18/12);  minor base G (0/1);  total (18/13)
Fisher's exact test for biased strand distribution p-value = 4.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATCCGCGCCGCATCTTCCACGCTGTGCGACACCATCAATAGCGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAACATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTG  >  NC_000913/72019‑72486
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                         
gATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTa                                                                                                                                                                                                                            <  1:182661/250‑1 (MQ=255)
       gctcgcggctgggtggCGCTGGTGGTTTTAAGTCGCCACCCGTTGATGCCGGTGTATATCGATATTTGGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTGGGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTca                                                                                                                                                                                                                     <  2:198240/239‑1 (MQ=255)
         cATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTcaac                                                                                                                                                                                                                   >  1:133743/1‑250 (MQ=255)
             ggCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTc                                                                                                                                                                                                              >  2:211640/1‑251 (MQ=255)
                  ggtggCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTg                                                                                                                                                                                                          <  1:2690/250‑1 (MQ=255)
                     ggCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTct                                                                                                                                                                                                      >  1:57819/1‑250 (MQ=255)
                                        cGCCCCCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATgc                                                                                                                                                                                    <  1:626/250‑1 (MQ=255)
                                                   gATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCg                                                                                                                                                                         >  1:222680/1‑250 (MQ=255)
                                                                atCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAgcg                                                                                                                                                           >  1:95342/1‑251 (MQ=255)
                                                                      atTTTGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTCTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGc                                                                                                                                                     >  1:25656/1‑251 (MQ=255)
                                                                         ttGCGTAAAGTGGTTGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACCAAACGCACTGTCCGTAGGACGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCAGTCGGCGACCACCACCGCGCGCGTGGCgct                                                                                                                                                  >  2:198313/1‑249 (MQ=255)
                                                                                  gTGGTTGGGGTTTAGGGGATTGTTTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCgcgcggc                                                                                                                                         <  2:177614/251‑4 (MQ=255)
                                                                                       tGGGGTTTAGGGGATTGTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCtt                                                                                                                                    >  1:70482/1‑251 (MQ=255)
                                                                                                ggggattgTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACtgt                                                                                                                           >  1:108338/8‑251 (MQ=255)
                                                                                                   gattgTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTTAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACCCtgtgtg                                                                                                                        >  1:121472/5‑249 (MQ=255)
                                                                                                    attgTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACtgtgtg                                                                                                                        >  1:222099/4‑248 (MQ=255)
                                                                                                     ttttCTTGAGCGGGTTGCCTGGTGCGACGCTGACCGGTTCTTGTCGGGCCTACAAAACGCACTGTCCGTAGTTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCGTAAGCACTCGCCTTACCGCTCAACAACTCATGAGTCTTTCCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTAGCGATCCGCGCGGCATCTTCCACACtgtgtga                                                                                                                       <  2:1101/247‑4 (MQ=255)
                                                                                                                gggTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATa                                                                                                            >  1:221030/1‑250 (MQ=255)
                                                                                                                                                    gCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCACTAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCAcc                                                                         >  1:32660/1‑249 (MQ=255)
                                                                                                                                                        aCAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCg                                                                     <  1:58646/249‑1 (MQ=255)
                                                                                                                                                        aCAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTc                                                                   >  1:93806/1‑251 (MQ=255)
                                                                                                                                                              cgcactgtccGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGTATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCAtc                                                             <  2:133743/243‑1 (MQ=255)
                                                                                                                                                               gcactgtccGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGCTCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCAtct                                                            >  1:17188/8‑251 (MQ=255)
                                                                                                                                                                cactgtccGTAGGTCGGATAAGGTGTTCACGCCGCATCGGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTGACCCTGCCAGGCGATGCGCCCGGCGGCGACCACCACCGAGCTCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACTTCTGTTTTTGCTGCGGGCAGCTCGTGCTCACCAGCGTCAGCAtctc                                                           <  1:167862/245‑1 (MQ=255)
                                                                                                                                                                                   aaGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATc                                         <  1:39902/250‑1 (MQ=255)
                                                                                                                                                                                   aaGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATc                                         <  2:70463/250‑1 (MQ=255)
                                                                                                                                                                                                         cAGTGCATACTAACCCGTAATACCCAATAGTGCTGAAGCACTCGCCTTACCGCTCCACAACTCATTGTTCTTACCCTGCCAGGCGATGCTCCCGTCGGCGACCACCACCGAGCGCGTGACGATCCGCGCGGCCTCTTCCACACTGTGTGCCACCCTCAATACCGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCg                   >  2:42762/1‑250 (MQ=255)
                                                                                                                                                                                                                tACTAACCCGTAATACCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGTTCGAGCGCCGATAACGGTTCactgagaaata            >  2:215519/1‑240 (MQ=255)
                                                                                                                                                                                                                     aCCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCg      <  2:217954/251‑1 (MQ=255)
                                                                                                                                                                                                                        cGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCt   >  2:24194/1‑251 (MQ=255)
                                                                                                                                                                                                                         gTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGGCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACCCCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTg  <  2:166446/251‑1 (MQ=255)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                         
GATTTCTGCCATCGGCGTGGTGGCGCTGGTGGTGTTAATTCGCCACCCGTTGATGCCGGTGTATATCGATATTTTGCGTAAAGTGGTTGGGGTTTAGGAGATAGTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATCCGCGCCGCATCTTCCACGCTGTGCGACACCATCAATAGCGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAACATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTG  >  NC_000913/72019‑72486

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: