Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,679,735 T→A 14.9% L60H (CTC→CAC)  ydgI → putative arginine/ornithine antiporter transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,679,7350TA14.9% 55.0 / 4.2 29L60H (CTC→CAC) ydgIputative arginine/ornithine antiporter transporter
Reads supporting (aligned to +/- strand):  ref base T (13/11);  new base A (2/2);  total (16/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.37e-01

GCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGACGGCGGTATCTTTACCTATGCCCGCGAAGGGTTTGGCGAGCTAATCGGCTTTTGTTCCGCATGGGGATACTGGCTGTGCGCAGTCATCGCCAACGTCTCCTATCTGGTGATCGTTTTTTCCG  >  NC_000913/1679592‑1679881
                                                                                                                                               |                                                                                                                                                  
tcTGTTATTAAGCTCAATCCTGTCCCCGGCTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCAcgc                                                                                                                                               <  2:509411/148‑1 (MQ=255)
  tGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCCCACgcgc                                                                                                                                             >  2:641551/1‑149 (MQ=255)
  tGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCACACgcgc                                                                                                                                             >  1:605502/1‑149 (MQ=255)
  tGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCACACgcgc                                                                                                                                             >  1:677805/1‑149 (MQ=255)
                      ggCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGAcgg                                                                                                                         <  1:782703/149‑1 (MQ=255)
                                ttttCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTTAcggcgggcgcttt                                                                                                               >  2:81729/1‑142 (MQ=255)
                                ttttCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGACGGCGGTATCttt                                                                                                               >  1:142500/1‑149 (MQ=255)
                                             gCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGGTCCTCACGCGCATTCGTCCCGAACTTGACGGCGGTATCTTTACCTATGCCagcg                                                                                                  >  1:437644/1‑149 (MQ=255)
                                                         ggcAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGACGGCGGGATCTTTACCTATGCCCGCGAAGGGTTTGGCg                                                                                      >  1:515308/1‑149 (MQ=255)
                                                                 ccAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGACGGCGGTATCTTTACCTATGCCCGCGAAGGGTTTGGCGCGCTAATc                                                                              >  1:66672/1‑149 (MQ=255)
                                                                 ccAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGACGGCGGTATCTTTACCTATGCCCGCGAAGGGTTTGGCGAGCTAATc                                                                              <  1:458569/149‑1 (MQ=255)
                                                                  cAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGATGATGCTCACGCGAATTCGTCACGAACTTGACGGCGGGATCTTTACCTATGCCCGCGAAGGGGTTGGCGAGCTAATCg                                                                             >  2:671480/1‑149 (MQ=255)
                                                                      ccgtgagcaCTGCTCATCGGCTGGAGTATTCCTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGACGGCGGTATCTTTACCTATGCCCGCGAAGGGTTTGGCGAGCTAATCGGCtt                                                                         <  2:712387/144‑1 (MQ=255)
                                                                      ccGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGACGGCGGTATCTTTACCTATGCCCGCGAAGGGTTTGGCGAGCTAATCGGCtt                                                                         <  1:655968/149‑1 (MQ=255)
                                                                             cagtGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTGACGCGCATTCGTCCCGAACTTGACGGCGGTATCTTTACCTATGCCCGCGAGGGGGTTGGCGAGCTAATCGGCTTTTGtttc                                                                  >  1:158808/4‑147 (MQ=255)
                                                                             caCTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCTATGCTGATCCTCACGCGCATTCGTCCCGAACTTGTCGGCGGTATCTTTACCTATGCCCGCGAAGGGTTTGGCGAGCTAATCGGCTTTTGTTcc                                                                  <  2:458075/149‑1 (MQ=255)
                                                                                     aTCGGCTGGCGTAGTTCTGGCGCTGGCAGTTTATTGCTCGCCTTTGCCAGGCTGATTCTCACTCGCATTCGTCCCGAACTTGACGGCGGTATCTTTACCTATGCCCGCGAAGGGTTTGGCGAGCTAATCGGCTTTTGTTCCGCATgggg                                                          <  2:607070/149‑1 (MQ=255)
                                                                                           tGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGACGGCGGTATCTTTACCTATGCCCGCGAAGGGTTTGGCGAGCTAATCGGCTTTTGTTCCGCATGGGGATACTg                                                    >  1:280390/1‑149 (MQ=255)
                                                                                                tATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGACGGCGGTATCTTTACCTATGCCCGCGAAGGGTTTGGCGAGCTAATCGGCTTTTGTTCCGCATGGGGATACTGGCtgt                                               <  1:668584/149‑1 (MQ=255)
                                                                                                  ttACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCACACGCGCATTCGTCCCGAACTTGACGGCGGTATCTTTACCTATGCCCGCGAAGGGTTTGGCGAGCTAATCGGCTTTTGTTCCGCATGGGTTTACTGGCTGTgc                                             <  2:677805/149‑1 (MQ=255)
                                                                                                  ttACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCACACGCGCATTCGTCCCGAACTTGACGGCGGTATCTTTACCTATGCCCGCGAAGGGTTTGGCGAGCTAATCGGCTTTTGTTCCGCATGGGGATACTGGCTGTgc                                             <  2:605502/149‑1 (MQ=255)
                                                                                                     cTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGACGGCGGTATCTTTACCTATGCCCGCGAAGGGTTTGGCGAGCTAATCGGCTTTTGTTCCGCATGGGGATACTGGCTGTGCGCa                                          <  1:234140/149‑1 (MQ=255)
                                                                                                        gcggTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGACGGCGGTATCTTTACCTATGCCCGCGAAGGGTTTGGCGAGCTAATCGGCTTTTGTTCCGCATGGGGATACTGGCTGTGCGCAGTc                                       <  2:689440/149‑1 (MQ=255)
                                                                                                             ggCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGACGGCGGTATCTTTACCTATGCCCGCGAAGGGTTTGGCGAGCTAATCGGCTTTTGTTCCGCATGGGGATACTGGCTGTGCGCAGTCATCGc                                  >  2:773331/1‑149 (MQ=255)
                                                                                                               cATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGACGGCGGTATCTTTACCTATGCCCGCGAAGGGTTTGGCGAGCTAATCGGCTTTTGTTCCGCATGGGGATACTGGCTGTGCGCAGTCATCGCCa                                >  1:235897/1‑149 (MQ=255)
                                                                                                                  tttATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGACGGCGGTATCTTTCCCTATGCCCGCGAAGGGTTTGGCGAGCTAATCGGCTTTTGTTCCGCATGGGGATACTGGCTGTGCGCAGTCATCGCCAACg                             >  2:91654/1‑149 (MQ=255)
                                                                                                                               tttGCCATGCTGATCCTCACGCGCAGTCGTCCCGAACTTGACGGCGGTATCTTTACCTATGCCCGCGAAGGGTTTGGCGAGCTAGTCGGCTTTTGTGCCGCATGGGGATACTGGCTGTGCGCGGTCATCGCCAACGTCTCCTCTCTGGt                >  1:132523/1‑149 (MQ=255)
                                                                                                                                    cATGCTGATCCTCACGCGCATTCGTCCCGAACGTGACGGCGGTATCTTTACCTATGCCCGCGAAGGGTTTGGCGAGCTAATCGGCTTTTGTTCCGCATGGGGATACTGGCTGTGCGCAGTCATCGCCAACGTCTCCTATCTGGTGATCg           <  1:717458/149‑1 (MQ=255)
                                                                                                                                      tGCTGATCCTCACGCGCATTCGTCCCGAACTTGACGGCGGTATCTTTACCTATGCCCGCGAAGGGTTTGGCGAGCTAATCGGCTTTTGTTCCGCATGGGGATACTGGCTGTGCGCAGTCATCGCCAACGTCTCCTATCTGGTGATCGtt         >  1:488418/1‑149 (MQ=255)
                                                                                                                                             cctccctcGCATTCGTCCCGAACTTGACGGCGGTATCTTTACCTATGCCCGCGAAGGGTTTGGCGAGCTAATCGGCTTTTGTTCCGCATGGGGATACTGGCTGTGCGCAGTCATCGCCAACGTCTCCTATCTGGTGATCGTTTTTTCcg  <  2:757988/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
GCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGACGGCGGTATCTTTACCTATGCCCGCGAAGGGTTTGGCGAGCTAATCGGCTTTTGTTCCGCATGGGGATACTGGCTGTGCGCAGTCATCGCCAACGTCTCCTATCTGGTGATCGTTTTTTCCG  >  NC_000913/1679592‑1679881

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