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.9% 74.9 / ‑2.0 29A79A (GCC→GCT) creAPF05981 family protein CreA
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base T (16/10);  total (18/11)
Fisher's exact test for biased strand distribution p-value = 5.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.21e-01

GAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGATGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACCGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCCATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTATTCAAAAAACGCACGTCCCTGGTCTTTAAGTCGTTACAGGTCGTGCGCTTTTATGATGCCAAACGCAACGCGCTCGCTTATCTGGCTTACTCCGACAAAGTTGTAGAAGGTTCGCCGAAAAACGCGATTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAG  >  NC_000913/4635523‑4635999
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                  
gAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGAGTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTATGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCTCCTGTGATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAg                                                                                                                                                                                                                                    <  1:198143/251‑1 (MQ=255)
               gATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAAc                                                                                                                                                                                                                     >  1:78063/1‑251 (MQ=255)
               gATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAAc                                                                                                                                                                                                                     >  1:299502/1‑251 (MQ=255)
                agcctttcgttAAGACTGATAATGCGGGGTCCATTGGCTCATGCTGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGATGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGTGTCTGGCGGATGTTACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACt                                                                                                                                                                                                                    <  2:63421/242‑1 (MQ=255)
                        tttAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATc                                                                                                                                                                                                            >  1:331754/1‑251 (MQ=255)
                                 gATAATTCTGGGGTCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGGCACCTGTTATTTGAGCCGGGCGTAAACAGGTGGTAGTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAATTCGGGCCGATTGAACTGTCGGATCGTATTaa                                                                                                                                                                                                     <  1:46258/249‑1 (MQ=255)
                                 gATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTaaa                                                                                                                                                                                                    >  1:213279/1‑250 (MQ=255)
                                    aaTGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCTGATCACAAAATTGTTGGTGAAGCCTTTTACGATCCCGATGTGAAAAATTTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGTGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACg                                                                                                                                                                                                <  1:29135/251‑1 (MQ=255)
                                                aTTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAgg                                                                                                                                                                                    >  1:181463/1‑251 (MQ=255)
                                                          gCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGgtagt                                                                                                                                                                          >  1:227202/1‑251 (MQ=255)
                                                                       gTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCCCCTGTTATGTGAGCAGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGCGGTAGTATTCaaaaaaagc                                                                                                                                                             >  1:42663/1‑248 (MQ=255)
                                                                               gccacCGTATTTAAAATGATCGGACCGGAGCACAAAATTGTTGTGGAAGCCTTTGACAATCTCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGTTGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGCGGTAGAATTCAAAAAACTCACGTcccc                                                                                                                                                     >  2:62396/3‑250 (MQ=255)
                                                                                      tATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTATTCAAAAAAAGCACGTCCCTGGTCttt                                                                                                                                              >  1:274928/1‑251 (MQ=255)
                                                                                       aTTTAAAATGAACTGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGGCACCTGTTATGTGAGCCGGGCGAAAAAAGGTGGGATTAAAGTGGGATTGGGTCTGGCGGCAGATACCTCAGATGAGGATATTACTTGTCAACAAGTCGAGGCGAATGAACTGTCGGGTCGTATTAAAAAAGGCTAAGCTCAGGGCGAGGGAGTATTCACAAAACTCAAGTACCTGGGCtta                                                                                                                                              >  2:285272/1‑249 (MQ=255)
                                                                                                       ccGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGAGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTATTCAAAAAACGCACGTCCCTGGTCTTTAAGTCGTTACAGGGCGt                                                                                                                             >  1:102842/1‑251 (MQ=255)
                                                                                                            tCACAAAATTGTTGTGGAAGCCTTTGATGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACCGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCCATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTAGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGAAGTATTCAAAAAACGAACGTCCCTGGTCTTTAAGTCGTTACAGGTCGTGCGCt                                                                                                                        >  1:51162/1‑251 (MQ=255)
                                                                                                                  aaTTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGAATTCAAAAAACGCACGTCCCTGGTCTTTAAGTCGTTACAGGTCGTGCGATTTTatg                                                                                                                  >  1:232030/1‑251 (MQ=255)
                                                                                                                   aTTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTATTCAAAAAACGCACGTCCCTGGTCTTTAAGTCGTTACAGGTCGTGCGCTTTTatga                                                                                                                 <  1:308941/251‑1 (MQ=255)
                                                                                                                                   ttGACGATCCCGATGTGAATAATGTCAACAGTTAAGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTTTGGCGGAATATACCTCCGATGCTGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTATTCAAAAAACGCACGTCCCTGGTCTTTAAGTCGTTACAGGGCGTGCGCTTTTATGACGCCAaacgcaacgc                                                                                                  <  2:41117/250‑1 (MQ=255)
                                                                                                                                    tgacgatCCCGATGTGAAAAATGTCCCCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGTGCCGATTGAACTGTCGGATCGTATTAAAAACGTCAAAGCTCAGGTCGAGGTAGTATTCAAAAAACGCACGTCCCTGGTCTTTAAGTCGTTACAGTTCGTGCGCTTTTATGACGCCAAACGCAAcgcg                                                                                                 <  1:155492/246‑1 (MQ=255)
                                                                                                                                                            cACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGgtagt                                                                                                                                                                          <  2:263372/153‑1 (MQ=255)
                                                                                                                                                            cACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGgtagt                                                                                                                                                                          >  1:263372/1‑153 (MQ=255)
                                                                                                                                                                         aGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTATTCAAAAAACGCACGTCCCTGGTCTTTAAGTCGTTACAGGTCGTGCGCTTTTATGACGCCAAACGCAACGCGCTCGCTTATCTGGCTTACTCCGACAAAGTTGTAGAAgg                                                           >  1:33071/1‑251 (MQ=255)
                                                                                                                                                                           ccGGGCGAAACCAGGTGGTATTAAAGGTGGATTGGGTCTGGCGGAAGATACCTCAGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGCGGTAGTATTCAAAAAACGCACGTCCCTGGTCTTTAAGTCGTTACAGGTCGTGCGCTTTTATGACGCCAAACGCAACGCGCTCGCTAATCTGGCTTAATCCGACAAAGTTGtgtaagggt                                                         >  2:148652/1‑243 (MQ=255)
                                                                                                                                                                               gCGAAAACAGGTGGTATTAAAGGGGGATTTGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTATTCAAAAAACGCACGTCCCTGGTCTTTAAGTCGTTACAGGTCGTGCGCTTTTATGACGCCAAACGCAACGCGCTCGCTTATCTGGCTTACTCCGACAAAGTTGTAGAAGGTTCGcc                                                     <  1:204809/251‑1 (MQ=255)
                                                                                                                                                                                     acagagggtATTAAGTGGGTATTGGGTCGGGCGGAAGATTCCTACGATGAGTCGATATCATGTCAGGAAGTCGGGCCTATGTAATTGTCGGAGCGTATTAAAAACGGCAAAGCTCAGGACGAGGGAGTAGTCAAAAAACGCACGTCACTGGTCTTTAAGTCGTTACAGGTCGTGCGCTTGTATGACGTCAAACGCAGCGCGCTCGATAATCTGGCTTACTCCGACAAAGTTGTATAAGTTTCGCCGaaaaa                                               <  2:117753/245‑1 (MQ=255)
                                                                                                                                                                                                                 tGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTATTCAAAAAACGCACGTCCCTGGTCTTTAAGTCGTTACAGGTCGTGCGCTTTTATGACGCCAAACGCAACGCGCTCGCTTATCTGGCTTACTCCGACAAAGTTGTAGAAGGTTCGCCGAAAAACGCGATTATCGCGGTTCCTGTCATGCCg                   >  2:274046/1‑251 (MQ=255)
                                                                                                                                                                                                                             ccTCCGATGCGGCCGTTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAGACGGCAAAGCTCAGGGCGAGGTAGTATTCAAAAAACGAACGTCCCTGGTCTTTAAGTCGTTACAGGTCGTGCGCTTTTATGATGCCAAACGCAACGCGCTCGCTTATCTGGCTTACTCCGACAAAGTTGTAGTAGGTTCACCGAAAAACGCGATTAGCGCGGTTCCTGTCATGCCGTGGCGGCaaaa        >  2:143261/1‑248 (MQ=255)
                                                                                                                                                                                                                                  gATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTATTCAAAAAACGCACGTCCCTGGTCTTTAAGTCGTTACAGGTCGTGCGCTTTTATGACGCCAAACGCAACGCGCTCGCTTATCTGGCTTACTCCGACAAAGTTGTAGAAGGTTCGCCGAAAAACGCGATTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACagag  <  1:285272/251‑1 (MQ=255)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                  
GAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGATGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACCGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCCATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTATTCAAAAAACGCACGTCCCTGGTCTTTAAGTCGTTACAGGTCGTGCGCTTTTATGATGCCAAACGCAACGCGCTCGCTTATCTGGCTTACTCCGACAAAGTTGTAGAAGGTTCGCCGAAAAACGCGATTAGCGCGGTTCCTGTCATGCCGTGGCGGCAATAACAGAG  >  NC_000913/4635523‑4635999

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