Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 649,194 A→G G126G (GGT→GGC)  citF ← citrate lyase alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913649,1940AG97.1% 115.9 / NA 34G126G (GGT→GGC) citFcitrate lyase alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (17/16);  minor base T (1/0);  total (18/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCACAGGACGGGACGCCGAGGAAAGCCACGTCGATATTCAGTTCGCCGCTCTGTACCAGATGCACACGACCGCCGTGAGAGTGGATCTGCACCGGTTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGACCACGCAGGCCGGAGGTATAAATGCGGGTAACCACGCCCTGGCGAATGTGTTCTACCAGCGGCGCATGGCAATCACTCAGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATC  >  NC_000913/649059‑649329
                                                                                                                                       |                                                                                                                                       
atcaCAGGACGGGACGCCGAGGAAAGCCACGTCGATATTCAGTTCGCCGCTCTGTACCAGATGCACACGACCGCCGTGAGAGTGGATCTGCACCGGTTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGGCCAc                                                                                                                                     >  1:863416/1‑140 (MQ=255)
      ggacgggacgCCGAGGAAAGCCACGTCGATATTCAGTTCGCCGCTCTGTACCAGATGCACACGACCGCCGTGAGAGTGGATCTGCACCGGTTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGGCCACGCAGGc                                                                                                                               <  2:138772/140‑1 (MQ=255)
         cgggacgCCGAGGAAAGCCACGTCGATATTCAGTTCGCCGCTCTGTCCCAGATGCACACGACCGCCGTGAGAGTGGATCTGCACCGGTTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGGCCACGCAGGCCgg                                                                                                                            <  2:33467/140‑1 (MQ=255)
                      aaaGCCACGTCGATATTCAGTTTGCCGCTGTGTACCACATGCACACGACCGCCGTGAGAGTGGATCTGCACCGGTTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGGCCACGCAGGCCGGAGGTATAAATGCg                                                                                                               <  1:358874/140‑1 (MQ=255)
                      aaaGCCACGTCGATATTCAGTTCGCCGCTCTGTACCAGATGCACACGACCGCCGTGAGAGTGGATCTGCACCGGTTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGGCCACGCAGGCCGGAGGTATAAATGCg                                                                                                               >  2:525748/1‑140 (MQ=255)
                                 gATATTCAGTTCCCCGCTCTGTCCAAGATGCACACGACCGCCGTGAGAGTGGATCTGCACCGGTTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGGCCACGCAGGCCGGAGGTATAAATGCGGGTAACCACGc                                                                                                    <  2:863416/140‑1 (MQ=255)
                                 gATATTCAGTTCACCGCTCTGTACCAGATGCACACGACCGCCGTGAGAGTGGATCTGCACCGGTTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGGCCACGCAGGCCGGAGGTATAAATGCGGGTAACCACGc                                                                                                    <  1:685645/140‑1 (MQ=255)
                                   tatTCAGTTCGCCGCCCTGTCCGAGTTGCCCACGACCGCCGTGAGAGTGGATCTGCACCGGTTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGGCCACGCAGGCCGGAGGTATAAATGCGGGTAACCACGccc                                                                                                  <  1:340118/140‑1 (MQ=255)
                                          ttCGCCGCTCTGTACCAGATGCACACGACCGCCGTGAGAGTGGATCTGCACCGGTTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGGCCACGCAGGCCGGAGGTATAAATGCGGGTAACCACGCCCTGGCGaa                                                                                           >  2:519631/1‑140 (MQ=255)
                                           tCGCCGCTCTGTACCAGATGCACACGACCGCCGTGAGAGTGGATCTGCACCGGTTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGGCCACGCAGGCCGGAGGTATAAATGCGGGTAACCACGCCCTGGCGAAt                                                                                          >  1:765753/1‑140 (MQ=255)
                                            cgccgcTCTGTACCAGATGCACACGACCGCCGTGAGAGTGGATCTGCACCGGTTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGTCCACGCAGGCCGGAGGTATAAATGCAGGTAGCCACGCCCTGGCGAAtg                                                                                         >  2:340118/1‑140 (MQ=255)
                                               cgcTCTGTACCAGATGCACACGACCGCCGTGAGAGTGGATCTGCACCGGTTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGGCCACGCAGGCCGGAGGTATAAATGCGGGTAACCACGCCCTGGCGAATgtgt                                                                                      >  2:177844/1‑140 (MQ=255)
                                                       aCCAGATGCACACGACCGCCGTGAGAGTGGATCTGCACCGGTTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGGCCACGCAGGCCGGAGGTATAAATGCGGGTAACCACGCCCTGGCGAATGTGTTCTACCAg                                                                              <  2:765753/140‑1 (MQ=255)
                                                           gATGCACACGACCGCCGTGAGAGTGGATCTGCACCGGTTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGGCCACGCAGGCCGGAGGTATAAATGCGGGTAACCACGCCCTGGCGAATGTGTTCTACCAgcggc                                                                          >  1:837218/1‑140 (MQ=255)
                                                            aTGCACACGACCGCCGTGAGAGTGGATCTGCACCGGTTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGGCCACGCAGGCCGGAGGTATAAATGCGGGTAACCACGCCCTGGCGAATGTGTTCTACCAgcggcg                                                                         >  1:426468/1‑140 (MQ=255)
                                                            aTGCACACGACCGCCGTGAGAGTGGATCTGCACCGGTTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGGCCACGCAGGCCGGAGGTATAAATGCGGGTAACCACGCCCTGGCGAATGTGTTCTACCAgcggcg                                                                         >  1:116545/1‑140 (MQ=255)
                                                              gCACACGACCGCCGTGAGAGTGGATCTGCACCGGTTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGGCCACGCAGGCCGGAGGTATAAATGCGGGTAACCACGCCCTGGCGAATGTGTTCTACCAGCGGCGCa                                                                       >  1:772068/1‑140 (MQ=255)
                                                              gCACACCACCGCCGTGAGAGTGGATCTGCACCGGTTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGGCCACGCAGGCCGGAGGTATAAATGCGGGTAACCACGCCCTGGCGAATGTGTTCTACCAGCGGCGCa                                                                       <  2:747281/140‑1 (MQ=255)
                                                                acacGACCGCCGTGAGAGTGGATCTGCACCGGTTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGGCCACGCAGGCCGGAGGTATAAATGCGGGTAACCACGCCCTGGCGAATGTGTTCTACCAGCGGCGCATg                                                                     >  2:28046/1‑140 (MQ=255)
                                                                acacGACCGCCGTGAGAGTGGATCTGCACCGGTTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGGCCACGCAGGCCGGAGGTATAAATGCGGGTAACCACGCCCTGGCGAATGTGTTCTACCAGCGGCGCATg                                                                     >  1:675775/1‑140 (MQ=255)
                                                                   cGACCGCCGTGAGAGTGGATCTGCACCGGTTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGGCCACGCAGGCCGGAGGTATAAATGCGGGTAACCACGCCCTGGCGAATGTGTTCTACCAGCGGCGCATGGCa                                                                  <  1:177844/140‑1 (MQ=255)
                                                                   acaCCGCCGGGAGAGTGGATCTGCTCCGGTTCTGCAAGCAGACCACGGGAGATCTCTTCCGCCAGTGGGCCACGCAGGCCGGAGGTATAAATGCGGGTAACCACGCCCTGGCGAATGTGTTCTACCAGCGGCGCATGGCa                                                                  <  2:426468/138‑1 (MQ=255)
                                                                    gACCGCCGTGAGAGTGGATCTGCACCGGTTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGGCCACGCAGGCCGGAGGTATAAATGCGGGTAACCACGCCCTGGCGAATGTGTTCTACCAGCGGCGCATGGCaa                                                                 >  1:693816/1‑140 (MQ=255)
                                                                     aCCGCCGTGAGAGTGGATCTGCACCGGTTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGGCCACGCAGGCCGGAGGTATAAATGCGGGTAACCACGCCCTGGCGAATGTGTTCTACCAGCGGCGCATGGCAAt                                                                <  1:525748/140‑1 (MQ=255)
                                                                                         gCACCGGTTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGGCCACGCAGGCCGGAGGTATAAATGCGGGTAACCACGCCCTGGCGAATGTGTTCTACCAGCGGCGCATGGCAATCACTCAGGGAGCTGGACGcc                                            <  2:772068/140‑1 (MQ=255)
                                                                                              ggTTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGGCCACGCAGGCCGGAGGTATAAATGCGGGTAACCACGCCCTGGCGAATGTGTTCTACCAGCGGCGCATGGCAATCACTCAGGGAGCTg                                                  <  1:598759/129‑1 (MQ=255)
                                                                                              ggTTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGGCCACGCAGGCCGGAGGTATAAATGCGGGTAACCACGCCCTGGCGAATGTGTTCTACCAGCGGCGCATGGCAATCACTCAGGGAGCTg                                                  >  2:598759/1‑129 (MQ=255)
                                                                                                                cGGGAGATCTCTTCCGCCAGTGGGCCACGCAGGCCGGAGGTATAAATGCGGGTAACCACGCCCTGGCGAATGTGTTCTACCAGCGGCGCATGGCAATCACTCAGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCAt                     >  1:701015/1‑140 (MQ=255)
                                                                                                                      aTCTCTTCCGCCAGTGGGCCACGCAGGCCGGAGGTATAAATGCGGGTAACCACGCCCTGGCGAATGTGTTCTACCAGCGGCGCATGGCAATCACTCAGGGAGGTGGACGCAAGGTTCTGGTTTTTAAAGCCCATCTTcgc               >  1:318278/1‑140 (MQ=255)
                                                                                                                               gCCAGTGGGCCACGCTGGCCGGCGGCATACATGCGGGTAACCACGCCCTGGCGAATGTGTTCTACCAGCGGCGCATGGCAATCACTCAGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTc      <  1:267727/140‑1 (MQ=255)
                                                                                                                               gCCAGTGGGCCACGCAGGCCGGAGGTATAAATGCGGGTAACCACGCCCTGGCGAATGTGTTCTACCAGCGGCGCATGGCAATCACTCAGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTc      <  2:693816/140‑1 (MQ=255)
                                                                                                                               gCCAGTGGGCCACGCAGGCCGGAGGTATAAATGCGGGTAACCACGCCCTGGCGAATGTGTTCTACCAGCGGCGCATGGCAATCACTCAGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTc      <  1:449220/140‑1 (MQ=255)
                                                                                                                               gCCAGTGGGCCACGCAGGCCGGAGGTATAAATGCGGGTAACCACGCCCTGGCGAATGTGTTCTACCAGCGGCGCATGGCAATCACTCAGGGAGCTGGACGACAGGGTCGGGTTTTTAGAGCCCATCTTCGCGATGACGTc      >  2:258852/1‑140 (MQ=255)
                                                                                                                                ccAGTGGGCCACGCAGGCCGGAGGTATAAATGCGGGTAACCACGCCCTGGCGAATGTGTTCTACCAGCGGCGCATGGCAATCACTCAGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTcc     >  2:865486/1‑140 (MQ=255)
                                                                                                                                   gTGGGCCACGCAGGCCGGAGGTATAAATGCGGGTAACCACGCCCTGGCGAATGTGTTCTACCAGCGGCGCATGGCAATCACTCAGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATc  <  2:544609/140‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                       
ATCACAGGACGGGACGCCGAGGAAAGCCACGTCGATATTCAGTTCGCCGCTCTGTACCAGATGCACACGACCGCCGTGAGAGTGGATCTGCACCGGTTCTGCCAGCAGACCACGGGAGATCTCTTCCGCCAGTGGACCACGCAGGCCGGAGGTATAAATGCGGGTAACCACGCCCTGGCGAATGTGTTCTACCAGCGGCGCATGGCAATCACTCAGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATC  >  NC_000913/649059‑649329

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑