Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,925,529 G→A R30R (AGG→AGA)  yobB → putative carbon‑nitrogen hydrolase family protein YobB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,925,5290GA90.5% 50.2 / ‑3.0 21R30R (AGG→AGA) yobBputative carbon‑nitrogen hydrolase family protein YobB
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (8/11);  total (9/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

GCTTCGGTCAATCTGTCACGTTTACCTTACGAGCCCAAACTTAAATAAAACTTATACAGAGTTACACTTTCTTACATAACGCCTGCTAAATTATGAGTATTTTCTAAACCGCACTCATAATTTGCAGTCATTTTGAAAAGGAAGTCATTATGTCGTTCTGGAAAGTTGCAGCAGCGCAATATGAGCCTCGTAAAACCTCACTCACCGAACAGGTAGCTCATCATCTGGAGTTTGTGAGGGCCGCCGCCAGGCAACAATGTCAGCTGTTAGTTTTTCCATCGCTTTCTTTACTGGGGTGTGATTATTCGCGACGTGCCCTTCCTGCCCCACCCGATCTCTCACTGTTAGACCCGCTTTGTTATGCCGCAACGACCTGGCGAATGACCATCATTGCCGGCCTTCCCGTTGAATATAACGATCGCTTTATTCGTGGAATTGCGGTGTTCGCACCGTGGCGAAAAACGCCT  >  NC_000913/1925291‑1925757
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                    
gCTTCGGTCAATCTGTCACGTTTACCTTACGAGCCCAAACTTAAATAAAACTTATACAGAGTTACACGGTCTTACAAAGCGCCTGCTAAATTATGAGTATTTTCTAAAGCGCTCTCATAATTTGGAGGCATAGTGAAAAGGAATTCATAATGTCGGGCTGGAAAGTTGCATCAGCACAATATGAGCCGCGTAAAACCTCACTCACCGAGCAGGTAGCTCATCATCTGGAGTTTGTGAGAGCCGCTGCCAgg                                                                                                                                                                                                                          <  1:44836/251‑1 (MQ=255)
                   gTTTACCTTACGAGCCCAAACTTAAAGAAAACTTATCCAGAGTTACACTTTCGTACATAACGCCTGCTAAATGATGAGTATTTTCTAAACCGCACTCATAATTTGCAGTCATTTTGAAAAGGAAGTCATTATGTCGTTCTGGAAAGTTGCAGCAGCGGAATATGAGCCTGGTAAAACCTCACTCACCGAACAGGGAGCTCGTCATCTGGAGTTTGTGAGAGCCGCCGCCAGGCAACAATGTCAGCTGTTAg                                                                                                                                                                                                       <  1:168177/251‑1 (MQ=255)
                     ttACCTTACGAGCCCAAACTTAAATAAAACTTATACAGAGTTACACTTTCTTACATAACGCCTGCTAAATTATGAGTATTTTCTAAACCGCACTCATAATTTGCAGTCATTTTGAAAAGGTAGTCATTATGTCGTTCTGGAAAGTTGCAGCAGCGCAATATGAGCCTCGTAAAACCTCACTCACCGAACAGGTAGCTCATCATCTGGAGTTTGTGAGAGCCGCCGCCAGGCAACAATGTCAGCTGTTAGt                                                                                                                                                                                                      <  1:294827/250‑1 (MQ=255)
                                   cAAACTTAAATAAAACTTATACAGAGTTACACTTTCTTACATAACGCCTGCTACATTATGAGTATTTTCTAAACCGCACTCATAATTTGCAGTCATTTTGAAAAGGAAGTCATTATGTCGTTCTGGAAAGTTGCAGCAGCGCAATATGAGCCTCGTAAAACCTCACTCACCGAACAGGTAGCTCATCATCTGGAGTTTGTGAGAGCCGCCGCCAGGCAACAATGTCAGCTGTTAGTTTTTCCATCGctttc                                                                                                                                                                                       >  1:175063/1‑251 (MQ=255)
                                   cAAACTTAAATAAAACTTATACAGAGTTACACTTTCTTACATAACGCCTGCTAAATTATGAGTATTTTCTAAACCGCACTCATAATTTGCAGTCATTTTGAAAAGGAAGTCATTATGTCGTTCTGGAAAGCTGCAGCAGCGCAATATGAGCCTCGTAAAACCTCACTCACCGAACAGGGAGCGCATCATCTGGAGTTTGTGAGAGCCGCCGCCAGGCAACAATGTCAGCTGTTAGTTTTTCCATCGctttc                                                                                                                                                                                       <  1:190275/251‑1 (MQ=255)
                                                              tACACTTTCTTACATAACGCCTGCTAAATTATGAGTATTTTCTAAACCGTACTCATAATTTGCAGTCATTTTGAAAAGGAAGTCATTATGTCGTTCTGGAAAGTTGCAGCAGCGCAATATGAGCCTCGTAAAACCTCACTCACCGAACAGGTAGCTCATCATCTGGAGTTTGTGAGAGCCGCCGCCAGGCAACAATGTCAGCTGTTAGTTTTTCCATCGCTTTCTTTACTGGGGAGTGATTATTCGCGACg                                                                                                                                                            >  1:25444/1‑251 (MQ=255)
                                                                          naTAACGCCTGCTAAATTCTGAGTATTTTCTAAACCGCACTCATAATTTGCAGTCATTTTGCAAAGGAAGTAATTATGTCGTTCTGGAAAGTTGCAGCAGCGCAATATGAGCCTCGTAAAACCTCACTCACCGAACAGGTAGCTCATCATCTGGTGTTTGTGAGAGCCGCCGCCAGGCAACAATGTCAGCTGTTAGTTTTTCCATCGCTTTCTTTACTGGGGTGAGATTATTCGCGACGTGCCCTTcccgc                                                                                                                                                >  1:199159/2‑248 (MQ=255)
                                                                                              gAGTATTTTCTAAACCGCACTCATAATTTGCAGTCATTTTGAAAAGGAAGTCATTATGTCGTTCTGGAAAGTTGCAGCAGCGCAATATGAGCCTCGTAAAACCTCACTCACCGAACAGGTAGCTCATCATCTGGAGTTTGTGAGAGCCGCCGCCAGGCAACAATGTCAGCTGTTAGTTTGTCCATCGCTTTCTTTACTGGGGTGTGATTATTCGCGACGTGCCCTTCCTGCCCCACCCGATCTCTCACTGt                                                                                                                            <  1:283459/251‑1 (MQ=255)
                                                                                                                            caggcATTTTGAAAAGGAAGTCATTATGTCGTTCGGGAAAGTTGCATCAGCGCAATAGGAGTCTCGTCAAGCCTCACTCACCGAACAGCTAGCTCATCATCTGTAGTTTTTGAGAGCCGCCGCCAGGCAACAATGTCAGCTGtt                                                                                                                                                                                                         <  1:166128/140‑1 (MQ=255)
                                                                                                                            cAGTCATTTTGAAAAGGAAGTCATTATGTCGTTCTGGAAAGTTGCAGCAGCGCAATATGAGCCTCGTAAAACCTCACTCACCGAACAGGTAGCTCCTCATCTGGAGTTTGTGAGAgccgccgc                                                                                                                                                                                                                              >  1:98804/1‑123 (MQ=255)
                                                                                                                              gTCATTTTGAAAGGGAAGTCATTATGTCGTTCTGGAAAGTTGCAGCAGCGCAATATGAGCCTCGTAAAGCCGCACTCACCGAACAGGTATCTCATCATCTGGAGTGTGTTAGGGCCGCCGCCAGGCAACAATATCAGCTGTTAGTTTTTCCATCTCTTTCTTTACTGGGGTGTTATTATTCGCGACTTTCCCTTCCTGCCGCACCCGATCTCTCACTGTTAGACCCGCTTTGTTATGCCGCAACGACCTgg                                                                                            <  1:299939/251‑1 (MQ=255)
                                                                                                                                  ttttGAAAAGGAAGTCATTATGTCGTTCTGGAAAGTTGCAGCAGCGCAATATGAGCCTCGTAAAACCTCACTCACCGAACAGGTAGCTCATCATCTGGAGTTTGTGAGAgccg                                                                                                                                                                                                                                  <  1:209296/113‑1 (MQ=255)
                                                                                                                                                     aTGTCGTTCTGGAAAGTTGCAGCAGCGCAATATGAGCCTCGTAAAACCGCACTCACCGAACAGGTAGCTCATCATCTGGAGTTTGTGAGAgccgc                                                                                                                                                                                                                                 <  1:278441/95‑1 (MQ=255)
                                                                                                                                                                    gTTGCAGCAGCGCAATATGAGCCTCGTAAAACCTCACTCACCGAACAGGTAGCGCATCATCTGGACTTTGTGAGAGCCGCCGCCAGGCAACAATGTCAGCTGt                                                                                                                                                                                                          <  1:299652/103‑1 (MQ=255)
                                                                                                                                                                       gcagcaGCGCAATATGCGCCTCGTAAAACCTCACTCACCGAACAGGTAGCTCATCATCTGGAGTTTGTGAGAGCCGACGCCAGGCAACAATGTCAGCTGTTAGTTTTTCCATCGCTTTCTTTACTGGGGtgtg                                                                                                                                                                         >  1:169391/1‑133 (MQ=255)
                                                                                                                                                                                      ggnnnnngTANAACCTCATTCACCGAACAGGTAGCTCATCATCTGAAGTTTGTGAGAGCCGCCGCCAGGCAACAATGTCAGCTGTTAGTTTTTCCATCGCTTTCTTTACTGGGGTGTGATGATTCGCGACGTGCTCTTCCTGcc                                                                                                                                               <  2:199159/137‑1 (MQ=255)
                                                                                                                                                                                         ncTCGTAAAACCTCACTCACCGAACAGGTAGCTCATCATCTGGAGTTTGTAAGAGCCGCCGCCAGGCAACAATGTCAGCTGTTAGTTTTTCCATCGCTTTCTTTACTGGGGTGTGATTATTCGCGACGTGCCCTTCCTGCCCCACCCGATCTCTCACTGTTAGACCCGCTTTGTTATGCCGCAACGACCTGGCGAATGACcat                                                                                 >  1:63101/2‑203 (MQ=255)
                                                                                                                                                                                                     antcncacCGNACAGNTAGCTCATCATCTGGNGTGTGTGATAGCCGCCTCCNNGNAACGANGTCAg                                                                                                                                                                                                              <  2:120697/63‑1 (MQ=255)
                                                                                                                                                                                                                   ngTAGCTCATCATCTGGCGTTTGTGAGAGCCGCCGCCAGGCAACAATGTCAGCTGTTAGTTTTTCCATCGCTTTCTTTACTGGGGTGTGATTATACGCGACGTGCCCTTCCTGCCCCACCCGATCTCTCACTGTTAGACCCGCTTTGTTATGCCGCAACGACCTGGCGAATGAc                                                                                    >  1:195558/2‑174 (MQ=255)
                                                                                                                                                                                                                        ctcataatatgggggttgagAGGGCCGGCGCAAGGCAACAATGTCAGCTGCTAGATTTTCCATCGCTTCCTTTACTGGGGTGTGATTATTAGCGAAGTGCCCTTCCTGAACCACCCGATCTCTCACTGTTAGACCGGCTTTGTTATGCCGCAACGACCTGGAGAAGGACCATCATTGCCGGCCTTCCCGTTGAATATAACGAACGCTTTATTCGGGGAATTGCGGTGTTCGCACCGTGGCGAAAAACGCCt  >  1:119323/20‑251 (MQ=255)
                                                                                                                                                                                                                                         gtgAGAGCCGCCGCCAGGCAACAATGTCAGCTGTTAGTTTTTCCATCGCTTTCTTTAATGGGGTGTGATTATTCGCGACGTGCCCTTCCTGCCCCACCCGATCTCTCACTGTTAGACCCCCTTTGTTATGCCGCACCGACCTGGAGAATGACCATCATTGCCGGCCTTCCCGTTGAATATAACGATCGCTTTTTTCGTGGCATTGCGGTGTTc                       >  1:172018/1‑213 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                    
GCTTCGGTCAATCTGTCACGTTTACCTTACGAGCCCAAACTTAAATAAAACTTATACAGAGTTACACTTTCTTACATAACGCCTGCTAAATTATGAGTATTTTCTAAACCGCACTCATAATTTGCAGTCATTTTGAAAAGGAAGTCATTATGTCGTTCTGGAAAGTTGCAGCAGCGCAATATGAGCCTCGTAAAACCTCACTCACCGAACAGGTAGCTCATCATCTGGAGTTTGTGAGGGCCGCCGCCAGGCAACAATGTCAGCTGTTAGTTTTTCCATCGCTTTCTTTACTGGGGTGTGATTATTCGCGACGTGCCCTTCCTGCCCCACCCGATCTCTCACTGTTAGACCCGCTTTGTTATGCCGCAACGACCTGGCGAATGACCATCATTGCCGGCCTTCCCGTTGAATATAACGATCGCTTTATTCGTGGAATTGCGGTGTTCGCACCGTGGCGAAAAACGCCT  >  NC_000913/1925291‑1925757

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑