Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,555,324 A→G S116P (TCT→CCT)  maeA ← malate dehydrogenase, decarboxylating, NAD‑requiring; malic enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,555,3240AG97.7% 147.1 / ‑4.3 44S116P (TCT→CCT) maeAmalate dehydrogenase, decarboxylating, NAD‑requiring; malic enzyme
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (28/15);  total (28/16)
Fisher's exact test for biased strand distribution p-value = 3.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.59e-01

TATTATGGTTCGGCACGTTTTGCAGAATATCGTCCATATTGTGCCGGTTCTGGTAAGAGATAAACACGCCGCGTGAACGGCGGTAGATCTCAGAAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTATAAATAACAGGCATCATCTCATCAAGATGATTGTTTACCAGACGGTAGAAGAGGGTTTCGTT  >  NC_000913/1555231‑1555420
                                                                                             |                                                                                                
tattatGGTTCGGCACGTTTTGCAGAATATCGTCCATATTGTGCCGGTTCTGGTAAGAGATAAACACGCCGCGTGAACGGCGGTAGATCTCAGGAAAACg                                                                                            >  1:315410/1‑100 (MQ=255)
       gTTCGGCACGTTTTGCAGAATATCGTCCATATTGTGCCGGTTATGGTAAGAGATAAACACGCCGCGTGAAAGGCGGTAGATCTCAGGAAAACGCTCACAg                                                                                     <  2:921811/100‑1 (MQ=255)
        ttCGGCACGTTTTGCAGAATATCGTCCATATTGTGCCGGTTCTGGTAAGAGATAAACACGCCGCGTGAACGGCGGTAGATCTCAGGAAAACGCTCACAgg                                                                                    >  1:1508438/1‑100 (MQ=255)
         tCGGCACGTTTTGCAGAATATCGTCCATATTGTGCCGGTTCTGGTAAGAGATAAACACGCCGCGTGAACGGCGGTAGATCTCAGGAAAACGCTCACAGGc                                                                                   >  1:1549550/1‑100 (MQ=255)
              aCGTTTTGCAGAATATCGTCCATATTGTGCCGGTTCTGGTAAGAGATAAACACGCCGCGTGAACGGCGGTAGATCTCAGGAAAACGCTCACAGGCTgcgc                                                                              >  1:112620/1‑100 (MQ=255)
                gTTTTGCAGAATATCGTCCATATTGTGCCGGTTCTGGTAAGAGATAAACACGCCGCGTGAACGGCGGTAGATCTCAGGAAAACGCTCACAGGCTGCGCCg                                                                            >  2:378849/1‑100 (MQ=255)
                   ttGCAGAATATCGTCCATATTGTGCCGGTTCTGGTAAGAGATAAACACGCCGCGTGAACGGCGGTAGATCTCAGGAAAACGCTCACAGGCTGCGCCGACg                                                                         >  1:624258/1‑100 (MQ=255)
                     gCAGAATATCGTCCATATTGTGCCGGTTCTGGTAAGAGATAAACACGCCGCGTGAACGGAGGTAGATCTCAGGAAAACGCTCACAGGCTGCGCCGACGGt                                                                       >  1:438643/1‑100 (MQ=255)
                         aaTATCGTCCATATTGTGCCGGTTCTGGTAAGAGATAAACACGCCGCGTGAACGGCGGTAGATCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTggg                                                                   >  1:381984/1‑100 (MQ=255)
                          atatCGTCCATATTGTGCCGGTTCTGGTAAGAGATAAACACGCCGCGTGAACGGCGGTAGATCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTgggg                                                                  <  2:149947/100‑1 (MQ=255)
                            atCGTCCATATTGTGCCGGTTCTGGTAAGAGATAAACACGCCGCGTGAACGGCGGTAGATCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGta                                                                >  1:16077/1‑100 (MQ=255)
                              cGTCCATATTGTGCCGGTTCTGGTAAGAGATAAACACGCCGCGTGAACGGCGGTAGATCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGtata                                                              >  2:1460015/1‑100 (MQ=255)
                                tCCATATTGTGCCGGTTCTGGTAAGAGATAAACACGCCGCGTGAACGGCGGTAGATCTCAGAAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTataaa                                                            <  2:427664/100‑1 (MQ=255)
                                    tatTGTGCCGGTTCTGGTAAGAGATAAACACGCCGCGTGAACGGCGGTAGATCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTATAAATAAc                                                        >  1:598877/1‑100 (MQ=255)
                                        gtgCCGGTTCTGGTAAGAGATAAACACGCCGCGTGAACGGCGGTAGATCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTATAAATAACAGGc                                                    >  1:134387/1‑100 (MQ=255)
                                         tgCCGGTTCTGGTAAGAGATAAACACGCCGCGTGAACGGCGGTAGATCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTATAAATAACAGGca                                                   >  1:1225898/1‑100 (MQ=255)
                                          gCCGGTTCTGGTCAGAGATAAACACGCCGCGTGAACGGCGGTAGATCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTATAAATAACAGGcat                                                  >  2:87617/1‑100 (MQ=255)
                                          gCCGGTTCTGGTAAGAGATAAACACGCCGCGTGAACGGCGGTAGATCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTATAAATAACAGGcat                                                  >  2:1052079/1‑100 (MQ=255)
                                               ttCTGGTAAGAGATAAACACGCCGCGTGAACGGCGGTAGATCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTATAAATAACAGGCAtcatct                                             >  1:898035/1‑100 (MQ=255)
                                                  tGGTAAGAGATAAACACGCCGCGTGAACGGCGGTAGATCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTATAAATAACAGGCAtcatctcat                                          >  1:197221/1‑100 (MQ=255)
                                                  tGGTAAGAGATAAACACGCCGCGTGAACGGCGGTAGATCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTATAAATAACAGGCAtcatctcat                                          >  1:1168949/1‑100 (MQ=255)
                                                   ggTAAGAGATAAACACGCCGCGTGAACGGCGGTAGATCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTATAAATAACAGGCAtcatctcatc                                         <  1:119752/100‑1 (MQ=255)
                                                    gTAAGAGTTAAACACGCCGCGTGAACGGCGGTAGATCTCAGGATACCGCTCACTGGCTGCGCCGACGGTTGGGGTATAAATTACAGGCATCATCtcatca                                        <  1:929433/100‑1 (MQ=255)
                                                     tAAGAGATAAACACGCCGCGTGAACGGCGGTAGATCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTATAAATAACAGGCATCATCTCATCaa                                       >  2:1141390/1‑100 (MQ=255)
                                                      aaGAGATAAACACGCCGCGTGAACGGCGGTAGATCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTATAAATAACAGGCATCATCTCATCAAg                                      >  1:1288658/1‑100 (MQ=255)
                                                        gagaTAAACACGCCGCGTGAACGGCGGTAGATCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTATAAATAACAGGCATCATCTCATCACgat                                    >  2:271015/1‑100 (MQ=255)
                                                          tccgatctcgccgcGTGAACGGCGGTAGATCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTATAAATAACAGGCATCATCTCATCAAgatga                                  <  2:789480/92‑1 (MQ=255)
                                                             aaaCACGCCGCGTGAACGGCGGTAGATCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTATAAATAACAGGCATCATCTCATCAAGATGATTg                               >  2:1344468/1‑100 (MQ=255)
                                                             aaaCACGCCGCGTGAACGGCGGTAGATCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTATAAATAACAGGCATCATCTCATCAAGATGATTg                               <  2:1473341/100‑1 (MQ=255)
                                                               acacGCCGCGTGAACGGCGGTAGATCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTATAAATAACAGGCATCATCTCATCAAGATGATTGtt                             <  2:550588/100‑1 (MQ=255)
                                                                 acGCCGCGTGAACGGCGGTAGATCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTATAAATAACAGGCATCATCTCATCAAGATGATTGTTTa                           <  2:1225893/100‑1 (MQ=255)
                                                                  cgccgcGTGAACGGCGGTAGATCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTATAAATAACAGGCATCATCTCATCAAgatgaagatcgga                          >  1:789485/1‑92 (MQ=255)
                                                                      gcgTGAACGGCGGTAGATCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTATAAATAACAGGCATCATCTCATCAAGATGATTGTTTACCAGa                      <  1:24792/100‑1 (MQ=255)
                                                                          gAACGGCGGTAGATCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTATAAATAACAGGCATCATCTCATCAAGATGATTGTTTACCAGACGGt                  <  2:19189/100‑1 (MQ=255)
                                                                          gAACGGCGGTAGATCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTATAAATAACAGGCATCATCTCATCAAGATGATTGTTTACCAGACGGt                  <  1:1052743/100‑1 (MQ=255)
                                                                            aCGGCGGTAGATCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTATAAATAACAGGCATCATCTCATCAAGATGATTGTTTACCAGACGGTag                <  2:718957/100‑1 (MQ=255)
                                                                              ggcggTAGATCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTATAAATAACAGGCATCATCTCATCAAGATGATTGTTTACCAGACGGTagaa              <  1:1485648/100‑1 (MQ=255)
                                                                              ggcggTAGATCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTATAAATAACAGGCATCATCTCATCAAGATGATTGTTTACCAGACGGTagaa              <  2:1549540/100‑1 (MQ=255)
                                                                               gcggTAGATCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTATAAATAACAGGCATCATCTCATCAAGCTGATTGTTTACCAGACGGTagaag             >  1:692334/1‑100 (MQ=255)
                                                                                   tAGATCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTATAAATAACAGGCATCATCTCATCAAGATGATTGTTTACCAGACGGTAGAAGAggg         <  2:982435/100‑1 (MQ=255)
                                                                                      aTCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTATAAATAACAGGCATCATCTCATCAAGATGATTGTTTACCAGACGGTAGAAGAGGGttt      >  2:174995/1‑100 (MQ=255)
                                                                                      aTCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTATAAATAACAGGCATCATCTCATCAAGATGATTGTTTACCAGACGGTAGAAGAGGGttt      >  1:1440229/1‑100 (MQ=255)
                                                                                      aTCTCAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGACTAAATAACAGTCATCATCTCATCAATGTGATCGTTTAcacac                      >  2:1336447/1‑80 (MQ=255)
                                                                                          cAGGAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTATAAATAACAGGCATCATCTCATCAAGATGATTGTTTACCAGACGGTAGAAGAGGGTTTCGtt  >  2:765088/1‑100 (MQ=255)
                                                                                             |                                                                                                
TATTATGGTTCGGCACGTTTTGCAGAATATCGTCCATATTGTGCCGGTTCTGGTAAGAGATAAACACGCCGCGTGAACGGCGGTAGATCTCAGAAAAACGCTCACAGGCTGCGCCGACGGTTGGGGTATAAATAACAGGCATCATCTCATCAAGATGATTGTTTACCAGACGGTAGAAGAGGGTTTCGTT  >  NC_000913/1555231‑1555420

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