Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 75,070 C→T P137P (CCG→CCA thiB ← thiamine ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091375,0700CT95.5% 66.4 / NA 22P137P (CCG→CCAthiBthiamine ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (11/10);  minor base G (1/0);  total (12/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGGTGGTGTAACTCAGTACCAGATCGCTTTCACCTTTTAAAAACAGGCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGGCGTCTAACAGGTTGTTATCCAGCCCCAGCACCACATCGGCTTTACTGTTTTTGCCTTCCATCCGTAGACGGTTG  >  NC_000913/74828‑75289
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                           
gagGTGGTGTAACTCAGTACCAGACCGCTTTCCCCTTTTAAAAACAGGCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTGTGTGGTGGGttttt                                                                                                                                                                                                                     <  1:63746/251‑1 (MQ=255)
                                aCCTTTAAATAACAGGCCTTACGCTTCGCTCCAGCCTTTGGTGACCGGGTCCGTTTTCTTCGCCAGTTTCTGCCAGGATTGTGGGGCGTCATCGCCATCGTCCTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTTCTGGTGCGCGGATCCTGATAAATCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGTGGGTTTTTCAGTTTGTTCGTGTCATAAACGAAGGCGAAGt                                                                                                                                                                                     <  1:51556/251‑1 (MQ=255)
                                  cTTTTAAACACAGGCCGTCGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGTGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGt                                                                                                                                                                                     >  1:65096/1‑249 (MQ=255)
                                                 cGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGTGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGt                                                                                                                                                                    >  2:202870/1‑251 (MQ=255)
                                                            cTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGTGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAACCGGTCCGAAAGTGTc                                                                                                                                                         >  2:190674/1‑251 (MQ=255)
                                                                  ccTTTGGTGACCGTGACTGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATAGCCTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGTTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGTGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCattatt                                                                                                                                                   <  1:159340/251‑1 (MQ=255)
                                                                     ttGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGTGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGGGGGTTTTTCAGTTTGTTCTTGTCACAAACGAAGGCGAGGTAGCCATAATCAAACGGTACGCAAGTGTCATTATTcc                                                                                                                                                 >  2:1054/1‑250 (MQ=255)
                                                                                                        ccAGGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGTGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCa                                                                                                                                                        >  1:6346/1‑208 (MQ=255)
                                                                                                        ccAGGCTTGTGGGGCGTCAGCGCCATAGACTTTTTGCATCCATAGCAACAGCCCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCTCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGTGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCa                                                                                                                                                        <  2:6346/208‑1 (MQ=255)
                                                                                                                               cATAGCCTTTTTGCACCCATAGCAACAGACCCAGCCCCGGTGTACGGGGGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGTGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGtt                                                                                      <  1:194504/251‑1 (MQ=255)
                                                                                                                                       ttttGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGTGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTggc                                                                              >  2:7071/1‑251 (MQ=255)
                                                                                                                                          ggcagccATAGCAACAGTCCCAGCCCCGGGGTACTGGTGCGCGGATCCTGAAAAGTCCCCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGTGGGTTTTTCAGTTTGTTCTTGTCATAAACTAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTggcggc                                                                           <  1:151252/246‑1 (MQ=255)
                                                                                                                                            cacccATAGCAACAGCCCCAGCCCCGGTGTACTGGTGCGCGGACCCTGATAAATCCCCCGCCATTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGTGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATaa                                                                                                                                                                             <  1:127723/148‑1 (MQ=255)
                                                                                                                                            cATCCCTAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGTGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATaa                                                                                                                                                                             >  2:127723/1‑151 (MQ=255)
                                                                                                                                                     cAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGTGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGGCGTCTAACAGgt                                                                <  1:141052/251‑1 (MQ=255)
                                                                                                                                                      atggaaCCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGCTCTTTCAGGCTTTGTGGTGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCCCACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGGCGTCTAACAGgtt                                                               >  2:44877/6‑251 (MQ=255)
                                                                                                                                                             ccAGCCCCGGTGTGCTGGTGCGCGGATCCTGATAATTCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGTGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGGCGTCTAACAGGTTGTTATcc                                                        <  1:210859/251‑1 (MQ=255)
                                                                                                                                                                       tgtACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGTGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGGCGTCTAACAGGTTGTTATCCAGCCCCAGc                                               >  1:145225/1‑250 (MQ=255)
                                                                                                                                                                             ggTGCGCGGATCCTGATAAATCAACCGCCAGTTTTGCTCGATCTCAACCAGTTCTTTCAGGCTTTGTGGTGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAAAGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTACCACACCGCTTTTTGCAAACAGTCCGGTTTTACTGGAGGCGGCTAACCGGTTGTTATCCAGCCCCAGCACCACAt                                        >  2:162718/1‑251 (MQ=255)
                                                                                                                                                                                  gcgGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGTGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGGCGTCTAACAGGTTGTTATCCAGCCCCAGCACCACATCGGCt                                   >  1:221623/1‑251 (MQ=255)
                                                                                                                                                                                      aTCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGTGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGGCGTCTAACAGGTTGTTATCCAGCCCCAGCACCACATCGGCTTTAc                               >  1:225037/1‑251 (MQ=255)
                                                                                                                                                                                                                   cGCTCTCAACCAGTTCTTTCAGGCTTTGTGGTGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGGCGTCTAACAGGTTGTTATCCAGCCCCAGCACCACATCGGCTTTACTGTTTTTGCCTTCCATCCGTAGACGGTTg  <  1:202870/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                           
GAGGTGGTGTAACTCAGTACCAGATCGCTTTCACCTTTTAAAAACAGGCCGTAGGCTTCGCTCCAGCCTTTGGTGACCGTGACCGTTTTCTTCGCCAGTTTCTGCCAGGCTTGTGGGGCGTCATCGCCATAGACTTTTTGCATCCATAGCAACAGACCCAGCCCCGGTGTACTGGTGCGCGGATCCTGATAAATCACCCGCCAGTTTTGATCGCTCTCAACCAGTTCTTTCAGGCTTTGTGGCGGGTTTTTCAGTTTGTTCTTGTCATAAACGAAGGCGAAGTAGCCATAATCAAACGGTACGAAAGTGTCATTATTCCAGCCGCCGGGAACGTTAACGGCATCCGCTGCCACACCGCTTTTGGCAAACAGTCCGGTTTTACTGGCGGCGTCTAACAGGTTGTTATCCAGCCCCAGCACCACATCGGCTTTACTGTTTTTGCCTTCCATCCGTAGACGGTTG  >  NC_000913/74828‑75289

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