Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 605,124 G→A 13.2% R145H (CGC→CAC)  thiL → thiamin monophosphate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947605,1240GA13.2% 85.8 / 6.8 38R145H (CGC→CAC) thiLthiamin monophosphate kinase
Reads supporting (aligned to +/- strand):  ref base G (21/12);  new base A (2/3);  total (23/15)
Fisher's exact test for biased strand distribution p-value = 3.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.00e-01

TGGCGGCCTATGCCGACGGCCTTGGTCGCATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGAGCGCGAGGCGCCTGCCGCGCAGGCCGAGCCGCTGCTGGCTCATTACTGGTCACCGATG  >  NC_002947/604989‑605269
                                                                                                                                       |                                                                                                                                                 
tggcGGCCTATGCCGACGGCCTTGGTCGCATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCaa                                                                                                                                     >  1:232043/1‑150 (MQ=255)
tggcGGCCTATGCCGACGGCCTTGGTCGCATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCaa                                                                                                                                     <  1:529085/150‑1 (MQ=255)
tggcGGCCTATGCCGACGGCCTTGGTCGCATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGACTTGCGCCGCAGCGGTGCaa                                                                                                                                     >  1:135204/1‑150 (MQ=255)
 ggcggcCTATGCCGACGGCCTTGGTCGCATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCCGGCCTTGCACCGCTGCGGCGCAAg                                                                                                                                    <  2:140844/150‑1 (MQ=255)
       cTATGCCGACGGCCTTGGTCGCATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCCGGCCTTGCACCGCTGCGGCGCAAGGCCtgg                                                                                                                              <  1:262532/150‑1 (MQ=255)
                   ccTTGGTCGCATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTtgtg                                                                                                                  <  1:397422/150‑1 (MQ=255)
                   ccTTGGTCGCATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTtgtg                                                                                                                  >  1:423750/1‑150 (MQ=255)
                           gCATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCggtg                                                                                                          >  1:351401/1‑150 (MQ=255)
                            cATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCACCGCTGCGGCGCAAg                                                                                                                                    >  2:539901/1‑123 (MQ=255)
                            cATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCACCGCTGCGGCGCAAg                                                                                                                                    <  1:539901/123‑1 (MQ=255)
                                  ccACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCt                                                                                                   >  2:168849/1‑150 (MQ=255)
                                   cACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTg                                                                                                  >  1:132398/1‑150 (MQ=255)
                                   cACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTg                                                                                                  >  1:43323/1‑150 (MQ=255)
                                                    ccTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTAGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGc                                                                                 <  1:784866/150‑1 (MQ=255)
                                                              ggtgACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTg                                                                       <  2:43323/150‑1 (MQ=255)
                                                                  acacCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGc                                                                   >  2:791660/1‑150 (MQ=255)
                                                                      cacGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTtgtgt                                                                                                                 <  1:376872/100‑1 (MQ=255)
                                                                      cacGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTtgtgt                                                                                                                 >  2:376872/1‑100 (MQ=255)
                                                                      cacGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTggg                                                               >  2:381168/1‑150 (MQ=255)
                                                                      cacGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTggg                                                               >  1:355170/1‑150 (MQ=255)
                                                                      cacGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTggg                                                               >  1:217639/1‑150 (MQ=255)
                                                                            cgGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGAgcg                                                         <  2:423750/150‑1 (MQ=255)
                                                                              gCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGGGAGcgcg                                                       >  1:265856/1‑150 (MQ=255)
                                                                                   ttGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGc                                                                                                                                            <  2:359624/60‑1 (MQ=255)
                                                                                   ttGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGc                                                                                                                                            >  1:359624/1‑60 (MQ=255)
                                                                                    tGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGAGCGCGAGgcgc                                                 <  2:135204/150‑1 (MQ=255)
                                                                                     gAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGAGCGCGAGGCGcc                                                <  1:381168/150‑1 (MQ=255)
                                                                                                   ccctgtTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGAGCGCGAGGCGCCTGCCGCGCAGGCCg                                  >  1:797503/4‑150 (MQ=255)
                                                                                                    cgtgTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGAGCGCGAGGCGCCTGCCGCGCAGGCCGa                                 >  1:793894/1‑150 (MQ=255)
                                                                                                     gtgtTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGAGCGCGAGGCGCCTGCCGCGCAGGCCGAg                                >  1:734751/1‑150 (MQ=255)
                                                                                                        atCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGa                                                            >  2:304442/2‑119 (MQ=255)
                                                                                                        atCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGa                                                            <  1:304442/118‑1 (MQ=255)
                                                                                                                 gTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGAGCGCGAGGCGCCTGCCGCGCAGGCCGAGCCGCTGCTGGCt                    >  2:397764/1‑150 (MQ=255)
                                                                                                                 gTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGAGCGCGAGGCGCCTGCCGCGCAGGCCGAGCCGCTGCTGGCt                    >  1:652220/1‑150 (MQ=255)
                                                                                                                        ttGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGAGCGCGAGGTGCCTGCGGCGCAGGCCGAGCCGCTGCTGGCTCATTACt             <  2:232043/149‑1 (MQ=255)
                                                                                                                             catgCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGAGCGCGAGGCGCCTGCCGCGCAGGCCGAGCCGCTGCTGGCTCATTAc              <  1:631912/141‑1 (MQ=255)
                                                                                                                             cAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGAGCGCGAGGCGCCTGCCGCGCAGGCCGAGCCGCTGCTGGCTCATTAc              >  2:631912/1‑144 (MQ=255)
                                                                                                                                   ttGCACCGCTGCGGCGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGAGCGCGAGGCGCCTGCCGCGCAGGCCGAGCCGCTGCTGGCTCATTACTGGTCACCGATg  >  2:285561/1‑150 (MQ=255)
                                                                                                                                       |                                                                                                                                                 
TGGCGGCCTATGCCGACGGCCTTGGTCGCATGGCCCACCGTTGCCGCATGAGCCTGATCGGTGGTGACACCACGCGCGGCCCGTTGAGCATCACCGTGACCGTGTTCGGCCGCGTGCCGGCTGGCCAGGCCTTGCGCCGCAGCGGTGCAAGGCCGGGTGATCTGTTGTGTGTCGGTGGTTCGCTGGGCAAGGCCGCCGGTGCCTTGCCGCTGGTGCTGGGTGAGCGCGAGGCGCCTGCCGCGCAGGCCGAGCCGCTGCTGGCTCATTACTGGTCACCGATG  >  NC_002947/604989‑605269

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