Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 754,446 A→G 100% L8P (CTC→CCC)  gltA ← citrate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913754,4460AG94.0% 38.6 / ‑5.5 15L8P (CTC→CCC) gltAcitrate synthase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (8/6);  total (9/6)
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.91e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CAGGATGCGGTTGAAGTGAAGCCTGGGTCAAAGGTGAACACACCTTTTGAACCGAGAGTACGGATATCAATAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGAGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTAAATTTGCCTTCGCACATCAACCTGGCTTTACCCGTTTTTTATTTGGCTCGCCGCTCTGTGAAAGAG  >  NC_000913/754311‑754579
                                                                                                                                       |                                                                                                                                     
cAGGATGCGGTTGAAGCGAAGCCTCGGTCAAAGCTGACCACACCTTTCGAACCGAGAGTACGGATATCAATAACATCTTGACCCAGCGTGCCTTTCCGCACATCCAGTCCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGc                                                                                                                                <  1:82504/143‑1 (MQ=255)
      gCGGTTGAAGTGAAGCCTGGGTCAAAGGTGAACACACCTTTTGAACCGAGAGTACGGATATCAATAACATCTTGACCCAGCGTGCCTTTCAGCACAGCCAGTTCAACATCTGTATCCCCGTTCAGGGAGAGTCTTGCCTTTGt                                                                                                                          >  1:120029/1‑143 (MQ=255)
                    ggatGGGTCAAAGGTGAACACACCTTTCGAACCGAGCGTACGGATATCAATAACATCTTGGCCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAg                                                                                                            <  2:221662/140‑1 (MQ=255)
                                         acCTTGCGAACCGAGATCACGGATCTCAGTAACAGCTTGGCCCAGCGTGGCTTTCAGCACATCCAGTTCAATAGCTGTATCCACGTTGATGTTGGGTTTTGCTTTTGTATCAGCTATTTAAGGTCTCCTTAGCGCCTTATTGc                                                                                       <  2:52827/143‑1 (MQ=255)
                                                 aaCCGAGAGTACGGATATCAATAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACt                                                                               >  2:257630/1‑143 (MQ=255)
                                                                  tCAATAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGAATGCCGGAACTTAAATTTg                                                              >  1:277617/1‑143 (MQ=255)
                                                                        tcATCTTGACCCATCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTCTCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTAAATTTGCCTTCg                                                        <  1:180564/142‑1 (MQ=255)
                                                                        aCATCTTGACCCGGTGTGCCTTGCAGCACATGCAGTGCAACGGCTGTTTCCCCACAAAGGGTGGGTGTTGCTTTGGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTAAATTTGCCTTCg                                                        <  1:125649/143‑1 (MQ=255)
                                                                                    aGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACCTAAATTTACCTTCTCACAtaaccatg                                            >  2:129021/1‑136 (MQ=255)
                                                                                        tGCCTTTCAGCACAGCCAGTGCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTAGCAGCCACCTAAGGTCCCCTTAGCGCCCAATTGCGTAAGACTGCCGGAACATAAATTTGCCTTCGCACATCAACCTCGCtt                                        >  1:182105/1‑143 (MQ=255)
                                                                                                  cacaTCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTAAATTTGCCTTCGCACATCAACCTGGCTTTACCCGtttt                              >  2:67807/1‑143 (MQ=255)
                                                                                                      tCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTAAATTTGCCTTCGCACATCAACATGGCTTTACCCGTTTTTTAt                          >  1:56348/1‑143 (MQ=255)
                                                                                                                 aGCTGTATCCCCGTTGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTCAGACTGCCGGAACTTAAATTTGCCTTCGCACATCAACCTGGCTTTACCCGTTTTTTATTTGGCtaggcg               >  2:144505/1‑138 (MQ=255)
                                                                                                                     gTATCTCCGTGTAGGGTCGTTTTTGCTTTTGTATCAGCCCTGTCAGGTCTCCTTAGCGCCTTATTGCGTAAGACTTCTGGACCTTAAATTTGCCTTCCCACATCAACCTGGCTTTTCCCGTTTTTTATTTGGCTCGCCGctct           <  2:120029/143‑1 (MQ=255)
                                                                                                                              ttGAGGGTGGGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGCACTTAAATTTGCCTTCGCACATCAACCTGGCTTTACCAGTTTTTGATTTGGCTCGCCGCTCTATGaacgcg  >  2:81894/1‑139 (MQ=255)
                                                                                                                                       |                                                                                                                                     
CAGGATGCGGTTGAAGTGAAGCCTGGGTCAAAGGTGAACACACCTTTTGAACCGAGAGTACGGATATCAATAACATCTTGACCCAGCGTGCCTTTCAGCACATCCAGTTCAACAGCTGTATCCCCGTTGAGGGTGAGTTTTGCTTTTGTATCAGCCATTTAAGGTCTCCTTAGCGCCTTATTGCGTAAGACTGCCGGAACTTAAATTTGCCTTCGCACATCAACCTGGCTTTACCCGTTTTTTATTTGGCTCGCCGCTCTGTGAAAGAG  >  NC_000913/754311‑754579

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