Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,635,757 C→T A79A (GCC→GCT)  creA → PF05981 family protein CreA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,635,7570CT92.0% 64.1 / ‑3.3 25A79A (GCC→GCT) creAPF05981 family protein CreA
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (12/11);  total (13/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 = 7.45e-01

ATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATT‑GTTGTGGAAGCCTTTGATGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACCGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCCATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTATTCAAAAAACGCACGTCCCTGGTCTTTAAGTCGTTACAGGTCGTGCGCTTTTATGATGCCAAACGCAACGCGCTCGCTTATCTGGCTTACTCCGACAAAGTTGTAGAAGGTTCGCCGAAAAACGCGATTAGCGCGGTTCCTGTCATGCCGTGGC  >  NC_000913/4635510‑4635986
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                     
aTGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATT‑GTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGCTGCGGCTAtt                                                                                                                                                                                                                                    >  1:321813/1‑251 (MQ=255)
       cAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATT‑GTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTCAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTc                                                                                                                                                                                                                             >  1:443460/1‑251 (MQ=255)
             gAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCGCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATT‑GTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGTTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAg                                                                                                                                                                                                                       <  2:141413/251‑1 (MQ=255)
                        aTTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATT‑GTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTc                                                                                                                                                                                                                             <  2:238206/234‑1 (MQ=255)
                        aTTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATT‑GTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTc                                                                                                                                                                                                                             >  1:238206/1‑234 (MQ=255)
                                                                   tCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATT‑GTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAgg                                                                                                                                                                 >  1:99261/1‑251 (MQ=255)
                                                                                                   tATTTAAAATGATCGGCCCGGATCACAAATTT‑GTTGTGGCAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCACCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCAGATTGAACTGGAGGAACGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAATATTCAAATAACGCACGGCCCTGGTCttt                                                                                                                                 >  1:478109/1‑251 (MQ=255)
                                                                                                    aTTTAAAATGATCGGCCCGGCTCACAAAATT‑GTTGTGGAAGCCTTTAACGCTCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAGCAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGCTACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGGCGGATCGTATTATAAACGGCAAAGCACAGGGCGCGGAAGTATTCAAAAACCGCACGACCCTGGCCttt                                                                                                                                 >  2:92319/1‑250 (MQ=255)
                                                                                                         aaaTGATCGGACCGGAACACAAAATTTGTTG‑GGAAGCATTTGACGATCCCGATGTGAAAAATGGCCCCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTATTCAAAAAACGCACGTCCCTGGTCTTTAAGTCg                                                                                                                           <  2:310169/251‑1 (MQ=255)
                                                                                                                          cacaAAATT‑GTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTATTCAAAAAACGCACGTCCCTGGTCTTTAAGTCGTTACAGGTCGTGCGCtt                                                                                                          <  1:402848/251‑1 (MQ=255)
                                                                                                                                    gttgtGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCCCCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTATTCAAAAAACGCACGTCCCTGGTCTTTAAGTCGTTACAGGTCGTGCGCTTTTatgacg                                                                                                  <  1:221698/250‑3 (MQ=255)
                                                                                                                                            aGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTATTCAAAAAACGCACGTCCCTGGTCTTTAAGTCGTTACAGGTCGTGCGCTTTTATGACGCCAaacgca                                                                                         <  1:430948/251‑1 (MQ=255)
                                                                                                                                              ccTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTGCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTATTCAAAAAACGCACGTCCCTGGTCTTTAAGTCGTTACAGGTCGTGCGCTTTTATGACGCCAaacgcaac                                                                                       <  2:466636/251‑1 (MQ=255)
                                                                                                                                                tttGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTATTCAAAAAACGCACGTCCCTGGTCTTTAAGTCGTTCCAGGTCGTGCGCTTTTAAGACGCAAaacgcaccgc                                                                                     >  2:466181/1‑251 (MQ=255)
                                                                                                                                                   gatgatCCCCATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACCGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCCATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTATTCAAAAAACGCACGTCCCTGGTCTTTAAGTCGTTACAGGTCGTGCGCTTTTATGATGCCAAACGCAACGcgct                                                                                  >  1:258327/1‑251 (MQ=255)
                                                                                                                                                             aTGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTACGATGCGGCTATTTCTTGTCAGCACGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCCAAGCTCATGGCGAGGTAGTATTCAAAAAACGAACGTACCTGGTCTTTAAGTCGTTACAGGTAGTGCGCTTTTATGACGCCAAACGAAACGCGCTCGCTTATCTg                                                                        >  1:203630/1‑251 (MQ=255)
                                                                                                                                                             aTGTGAAAAAGGTCCCCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGTGGATTGGGTCTGGCGGAATATACCTCCGAGGCGGCTATTTCTGGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTAGTCAAAAAACGCACGTCCCTGGTCGTTAAGTCGTTACAGGTCGGGCGCTTTTATGACGCCAAACGCAACGCGCTCGCTTATCTg                                                                        <  2:468747/251‑1 (MQ=255)
                                                                                                                                                                               gTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTATTCAAAAAACGCACGTCCCTGGTCTTTAAGTCGTTACAGGTCGTGCGCTTTTATGACGCCAAACGCAACGCGCTCGCTTATCTGGCTTACTCCGACAAAg                                                        <  1:317383/249‑1 (MQ=255)
                                                                                                                                                                                       aGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGCTTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTATTCAAAAAACGCACGTCCCTGGTCTTTAAGTCGTTACAGGTCGTGCGCTTTTATGACGCCAAACGCAACGCGCTCGCTTATCTGGCTTACTACGACAAAGTTGTAGAAgg                                              >  1:467829/1‑251 (MQ=255)
                                                                                                                                                                                              cGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGTCAAAGCTCAGGGCGAGGTAGTATTCAAAAAACGCACGGCCCTGGTCTTTAAGTCGTTACAGGTCGTGCGCTTTTATGACGCCAAACGCAACGCGCTCGCTTATCTGGCTTACTCCGACAAAGTTGTAGAAGGTTCGCCg                                       <  2:439403/251‑1 (MQ=255)
                                                                                                                                                                                                          gtATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAATCGGCAAAGCTCAGGGCGAGGTAGTATTCAAAAAACGCACGTCCCTGGTCTTTAAGTCGTTACAGGTCGTGCGCTTTTATGACGCCAAAAGCAACGCGCTCGCTTATCTGGCTTACTCTGCCAAAGGTGTAGAAGGTTCGCCGAAAAACGCGAtt                           >  2:372913/1‑251 (MQ=255)
                                                                                                                                                                                                                        ttGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTATTCAAAAAACGCACGTCCCTGGTCTTTAAGTCGTTACAGGTCGTGCGCTTTTATGACGCCAAACGCAACGCGCTCGCTTATCTGGCTTACTCCGACAAAGTTGTAGAAGGTTCGCCGAAAAACGCGATTAGCGCGGTTCCTGt             >  1:162140/1‑251 (MQ=255)
                                                                                                                                                                                                                          gggTCTGGCGGAAGATACCTCCGATGCGGCCATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTATTCAAAAAACGCACGTCCCTGGTCTTTAAGTCGTTACAGGTCGTGCGCTTTTATGATGCCAAACGCAACGCGCTCGCTTATCTGGCTTACTCCGACAAAGTTGTAGAAGGTTCGCCGAAAAACGCGATTAGCGCGGTTCCTGTCa           <  2:258327/251‑1 (MQ=255)
                                                                                                                                                                                                                             tCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTATTCAAAAAACGCACGTCCCTGGTCTTTAAGTCGTTACAGGTCGTGCGCTTTTATGACGCCAAACGCAACGCGCTCGCTTATCTGGCTTCCTCCGACAACTTTGTAGAAGGTTCGCCGAAAAACGCGATTAGCTCGGTTCCTGTCATGc        >  2:374316/1‑251 (MQ=255)
                                                                                                                                                                                                                                   ggAAGAGACCTCCGATGCTGCTATTTCTTGTCAGCAAGTCTGGCCGATTGACCGGTCGGGTCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTATTATTCAAAAAACGCACGTCCCTGGTCTTTAAGTCTTTACAGGTCGTGCTCTTTTAGTACGCCAAACGCAACGCGCTCGCTTATCTGGCTTACTCCGACAAAGTTGTAGAAGGATCGCCGAAAAACTCTATTAGCGCGGTTCCTGTCATGCCGTggc  <  2:377725/251‑1 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                     
ATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATT‑GTTGTGGAAGCCTTTGATGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACCGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCCATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTATTCAAAAAACGCACGTCCCTGGTCTTTAAGTCGTTACAGGTCGTGCGCTTTTATGATGCCAAACGCAACGCGCTCGCTTATCTGGCTTACTCCGACAAAGTTGTAGAAGGTTCGCCGAAAAACGCGATTAGCGCGGTTCCTGTCATGCCGTGGC  >  NC_000913/4635510‑4635986

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