Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 849,617 C→T A161T (GCT→ACT)  rhtA ← L‑threonine/L‑homoserine exporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913849,6170CT90.0% 54.2 / ‑3.2 20A161T (GCT→ACT) rhtAL‑threonine/L‑homoserine exporter
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base T (14/4);  total (14/6)
Fisher's exact test for biased strand distribution p-value = 7.89e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGGTGAGGGCAATCATTTCCAGCGAATAAGGCAGAGCGGTCGAGAGAATAGCGACAGCCAGACCCAATGGAATAACCGACCAGTGCCAGAGTGCTTCACCAGCCTGAAGCGCTCCAATTGGCACGAAAATTAACGCTGCAATCAACGAACCAATTGCCACCGTCGCAGGGCCATGTTCCGCTCCTGCGCGTTGCCCACTTAAAATGTAAATAGCCCAACAAGCCCCGGCCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAACGTCTTGCCCCAGCGGTAGCAGGAACCACAGACCAAGAACCGCCAGCACAACCCAGACGAAATCTACCGGGCGACGAGAAGAGAACAGCGCCACCGCCAGTGGTCCGGTGAACTCCAGCGCCACCGCAATACCCAGCGGTACTGTCTGAATAGAAAGATAAAAAAGATAATTCATCCCA  >  NC_000913/849403‑849852
                                                                                                                                                                                                                       |                                                                                                                                                                                                                                            
cGGGGGAGGGCAATCATTTCCAGCGAATAAGGCAGAGCGGTCGAGAGAATAGCGACAGCCAGACCCAATGGAATAACCGACCAGTGCCAGAGTGCTTCACCAGCCTGAAGCGCTCCAATTGGCACGAAAATTAACGCTGCAATCAACGAACCAATTGCCACCGTCGCAGGGCCATGTTCCGCTACTGCGCGTTGCCCACTTAAAATGTAAATAGTCCAACAAGCCCCGGCCACCAGTGTCAGAGCACAgtc                                                                                                                                                                                                         >  1:305960/1‑249 (MQ=255)
  ggTGAGGGCAATCATTTCCAGCGAATAAGGCAGAGCGGTCGAGAGAATAGCGACAGCCAGACCCAATGGAATAACCGACCAGTGCCAGAGTGCTTCACCAGCCTGAAGCGCTCCAATTGGCACGAAGATTAACGCTGCAATCAACGAACCAATTGCCACCGTCGCAGGGCCATGTTCCGCTCCTGCGCGTTGCCCACTTAAAATGTAAATAGTCCAACAAGCCCCGGCCCCCAGTGCCAGCGCACAGCCgg                                                                                                                                                                                                       <  2:289183/251‑1 (MQ=255)
                        gAATAAGGCAGCGCGGTCGAGAGAATAGCGACAGCCAGACCCAATGGAATAACCGACCAGTGCCAGAGTGCTTCACCAGCCTGAAGCGCTCCAATTGGCACGACAATTAACGCTGCAATCAACGAACCAATTGCCACCGTCGCAGGGCCATGTTCCGCTCCTGCGCGTTGCCCACTTAAAATGTAAATAGTCCAACAAGCCCCGGCCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAACGTc                                                                                                                                                                                 >  2:178427/1‑251 (MQ=255)
                                       gTCGAGAGAATAGCGACAGCCAGACCCAATGGAATAACCGACCAGTGCCAGAGTGCTTCACCAGCCTGAAGCGCTCCAATTGGCACGAAAATTAACGCTGCAATCAACGAACCAATTGCCACCGTCGCAGGGCCATGTTCCGCTCCTGCGCGTTGCCCACTTAAAATGTAAATAGTCCAACAAGCCCCGGCCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAACGTCTTGcccc                                                                                                                                                                          <  1:335805/243‑1 (MQ=255)
                                       gTCGAGAGAATAGCGACAGCCAGACCCAATGGAATAACCGACCAGTGCCAGAGTGCTTCACCAGCCTGAAGCGCTCCAATTGGCACGAAAATTAACGCTGCAATCAACGAACCAATTGCCACCGTCGCAGGGCCATGTTCCGCTCCTGCGCGTTGCCCACTTAAAATGTAAATAGTCCAACAAGCCCCGGCCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAACGTCTTGcccc                                                                                                                                                                          >  2:335805/1‑243 (MQ=255)
                                                  aGCGACAGCCAGACCCAATGGAATAACCGACCAGTGCCAGAGTGCTTCACCAGCCTGAAGCGCTCCAATTGGCACGAAAATTAACGCTGCAATCAACGAACCAATTGCCACCGTCGCAGGGCCATGTTCCGCTCCTGCGCGTTGCCCACTTAAAATGTAAATAGTCCAACAAGCCCCGGCCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAACGTCTTGCCCCAGCGGCAGCAGGACCcac                                                                                                                                                        >  2:251121/1‑250 (MQ=255)
                                                      aCAGCCAGACCCAATGGAATAACCGACCAGTGCCAGAGTGCTTCACCAGCCTGAAGCGCTCCAATTGGCACGAAAATTAACGCTGCAATCAACGAACCAATTGCCACCGTCGCAGGGCCATGTTCCGCTCCTGCGCGTTGCCCACTTAAAATGTAAATAGTCCAACAAGCCCCGGCCCCCAGTGCCAGCGcac                                                                                                                                                                                                             >  1:44525/1‑193 (MQ=255)
                                                      aCAGCCAGACCCAATGGAATAACCGACCAGTGCCAGAGTGCTTCACCAGCCTGAAGCGCTCCAATTGGCACGAAAATTAACGCTGCAATCAACGAACCAATTGCCACCGTCGCAGGGCCATGTTCCGCTCCTGCGCGTTGCCCACTTAAAATGTAAATAGTCCAACAAGCCCCGGCCCCCAGTGCCAGCGcac                                                                                                                                                                                                             <  2:44525/193‑1 (MQ=255)
                                                                      gAATAACCGACCAGTGCCAGAGTGCTTCACCAGCCTGAAGCGCTCCAATTGGCACGAAAATTAACGCTGCAATCAACGAACCAATTGCCACCGTCGCAGGGCCATGTTCCGCTCCTGCGCGTTGCCCACTTAAAATGTAAATAGTCCAACAAGCCCCGGCCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAACGTCTTGCCCCAGCGGTAGCAGGAACCACAGACCAAGAACCGCCAGcac                                                                                                                                    >  2:139090/1‑250 (MQ=255)
                                                                                       cAGAGTGCTTCACCAGCCTGAAGCGCTCCAATTGGCACGAAAATTAACGCTGCAATCAACGAACCAATTGCCACCGTCGCAGGGCCATGTTCCGCTCCTGCGCGTTGCCCACTTAAAATGTAAATAGTCCAACAAGCCCCGGCCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAACGTCTTg                                                                                                                                                                              <  2:88972/191‑1 (MQ=255)
                                                                                       cAGAGTGCTTCACCAGCCTGAAGCGCTCCAATTGGCACGAAAATTAACGCTGCAATCAACGAACCAATTGCCACCGTCGCAGGGCCATGTTCCGCTCCTGCGCGTTGCCCACTTAAAATGTAAATAGTCCAACAAGCCCCGGCCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAACGTCTTg                                                                                                                                                                              >  1:88972/1‑191 (MQ=255)
                                                                                                 ccccAGCCGGAACCGCTCCAATTGACACGAAATTACCCGCTGAACACACCGAACAAATTGCCACCGCCGCAGGCCCGTGTCCCGCCCCCGCGCGTTGCCCACTAAAAATGTAAATAGCCCAAAAGCCCCCGCCCCCCAGGCCCAGCCCCCAGCCGGTTAAACGGACAGGGAAAACGGCTTCCCCCAGGGGGAGCAGGAACCACAAACAAAGAACCGCCAGCACAACCCAGACGAAATCTACCGGGcgacga                                                                                                        <  2:338681/249‑1 (MQ=255)
                                                                                                   ccagccagAAGCGCTCCAATGGGCACAAAAATGAACGCGGCAACAAACGAACAAAGTCCCACCGTCGCGGGGCCTTTTCCCGCCCCGGCGCGTTGCCCACTAAAAAGGAAAATGGCCCAACAAGCCCCGCCCCCCAGGCCCAGCCCACGGCCGGTTAAACCGACAGGGGAAACGTCTCCCCCCGGCGGTAGCAGAACCCACAGACAAAAAACCGCCAGCACAACCCAGACGAAATCTACCGGGCGACgaga                                                                                                      <  2:13070/251‑1 (MQ=255)
                                                                                                                                   tAACGCTGCAATCAACGAACCAATTGCCACCGTCGCAGGGCCATGTTCCGCTCCTGCGCGTTGCCCACTTAAAATGTAAATAGTCCAACAAGCCCCGGCCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAACGTCTTGCCCCAGCGGTAGCAGGAACCACAGACCAAGAACCGCCAGCACAACCCAGACGAAATCTACCGGGCGACGAGAAGAGAACAGCGCCACCGCCAGTGGTCCGGTg                                                                       >  1:170715/1‑250 (MQ=255)
                                                                                                                                                       ccAATTGCCACCGTCGCAGGGCCATGTTCCGCTCCTGCGCGTTGCCCACTTAAAATGTAAATAGTCCAACAAGCCCCGGCCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAACGTCTTGCCCCAGCGGTAGCAGGAACCACAGACCAAGAACCGCCAGCACAACCCAGACGAAATCTACCGGGCGACGAGAAGAGAACAGCGCCACCGCCAGTGGTCCGGTGAACTCCAGCGCCACCGCAAt                                                   >  1:365763/1‑251 (MQ=255)
                                                                                                                                                             gCCACCGTCGCAGGGCCATGTTCCGCTCCTGCGCGTTGCCCACTTAAAATGTAAATAGTCCAACAAGCCCCGGCCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAACGTCTTGCCCCAGCGGTAGCAGGAACCACAGACCAAGAACCGCCAGCACAACCCAGACGAAATCTACCGGGCGACGAGAAGAGAACAGCGCCACCGCCAGTGGTCCGGTGAACTCCAGCGCCACCGCAATAccc                                               >  1:184589/1‑249 (MQ=255)
                                                                                                                                                                        aGGGCCATGTTCCGCTCCTGCGCGTTGCCCCCTTAAAATGTAACTAGTCCAACACGCCCCGGCCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAACGTCTTGCCCCAGCGGTAGCAGGAACCACAGACCAAGAACCGCCAGCACAACCCAGACGAAATCTACCGGGCGACGAGAAGAGAACAGCGCCACCGCCAGTGGTCCGGTGAACTCCAGCGCCACCGCAATACCCAGCGGTACTGTCt                                  >  1:45598/1‑251 (MQ=255)
                                                                                                                                                                                        ccTGCGCGTTGCCCACTTAACATGTAAATAGTCCACCAATCCCCTGCCCCCAGTGCCCGCGCACAGCCGGTTAAATCGACATGGGCAACGTCTTGCCCCAGCGGTAGCCGGAACCACAGACCAAGAACCGCCAGCACAACCCAGACGAAATCTCCCGGGCGACGCGAAGAGAACAGCGCCACCGCCAGTGGTCCGGTGAACTCCATCGCCACCGCACTACCCCGCGGTACTGTCTGAATAgcacgataca                   >  2:148448/1‑241 (MQ=255)
                                                                                                                                                                                                 tGCCCACTTAAAATGTAAATAGTCCAACAAGCCCCGGCCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAACGTCTTGCCCCAGCGGTAGCAGGAACCACAGACCAAGAACCGCCAGCACAACCCAGACGAAATCTACCGGGCGACGAGAAGAGAACAgcgc                                                                                          >  2:76937/1‑170 (MQ=255)
                                                                                                                                                                                                        ttAAAATGTAAATAGTCCAACAAGCCCCGGCCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAACGTCTTGCCCCAGCGGTAGCAGGAACCACAGACCAAGAACCGCCAGCACAACCCAGACGAAATCTACCGGGCGACGAGAAGAGAACAGCGCCACCGCCAGTGGTCCGGTGAACTCCAGCGCCACCGCAATACCCAGCGGTACTGTCTGAATAGAAAGATAAAAAAGATAATTCATcccc  >  2:322271/1‑250 (MQ=255)
                                                                                                                                                                                                                       |                                                                                                                                                                                                                                            
CGGGTGAGGGCAATCATTTCCAGCGAATAAGGCAGAGCGGTCGAGAGAATAGCGACAGCCAGACCCAATGGAATAACCGACCAGTGCCAGAGTGCTTCACCAGCCTGAAGCGCTCCAATTGGCACGAAAATTAACGCTGCAATCAACGAACCAATTGCCACCGTCGCAGGGCCATGTTCCGCTCCTGCGCGTTGCCCACTTAAAATGTAAATAGCCCAACAAGCCCCGGCCCCCAGTGCCAGCGCACAGCCGGTTAAATCGACATGGGAAACGTCTTGCCCCAGCGGTAGCAGGAACCACAGACCAAGAACCGCCAGCACAACCCAGACGAAATCTACCGGGCGACGAGAAGAGAACAGCGCCACCGCCAGTGGTCCGGTGAACTCCAGCGCCACCGCAATACCCAGCGGTACTGTCTGAATAGAAAGATAAAAAAGATAATTCATCCCA  >  NC_000913/849403‑849852

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