Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 764,243 A→G T22A (ACT→GCT)  sucD → succinyl‑CoA synthetase subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913764,2430AG96.4% 85.2 / ‑4.6 28T22A (ACT→GCT) sucDsuccinyl‑CoA synthetase subunit alpha
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (10/17);  total (10/18)
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.89e-01

TGCAGCAAAAGGTCTGACGGATGCAGCTCAGCAGGTTGTTGCCGCAGTGGAGGGGAAATAATGTCCATTTTAATCGATAAAAACACCAAGGTTATCTGCCAGGGCTTTACCGGTAGCCAGGGGACTTTCCACTCAGAACAGGCCATTGCATACGGCACTAAAATGGTTGGCGGCGTAACCCCAGGTAAAGGCGGCACCACCCACCTCGGCCTGCCGGTGTTCAACACCGTGC  >  NC_000913/764120‑764351
                                                                                                                           |                                                                                                            
tGCAGCAAAAGGTCTGACGGATGCAGCTCAGCAGGTTGTTGCCGCAGTGGAGGGGAAATAATGTCCATTTTAATCGATAAAAACACCAAGGTTATCTGCCAGGGCTTTACCGGTAGCCAGGGGGGTTTCCACt                                                                                                     <  2:66695/133‑1 (MQ=255)
tGCACCAAAAGGTCTGACGGATGCAGCTCAGCAGGTTCTTGCCGCAGTGGAGGGGAAATAATGTCCATTTTAATCGATAAAAACACCAAGGTTATCTGCCAGGGCTTTACCGGTAGCCAGGGGGCTTTCCACt                                                                                                     <  1:84807/133‑1 (MQ=255)
  cagcaAAAGGTCTGACGGATGCAGCTCAGCAGGTTGTTGCCGCAGTGGAGGGGAAATAATGTCCATTTTAATCGATAAAAACACCAAGGTTATCTGCCAGGGCTTTACCGGTAGCCAGGGGGCTTTCCACt                                                                                                     <  1:1269/131‑1 (MQ=255)
  cagcaAAAGGTCTGACGGATGCAGCTCAGCAGGTTGTTGCCGCAGTGGAGGGGAAATAATGTCCATTTTAATCGATAAAAACACCAAGGTTATCTGCCAGGGCTTTACCGGTAGCCAGGGGGCTTTCCACt                                                                                                     >  2:1269/1‑131 (MQ=255)
  cagcaAAAGGTCTGACGGATGCAGCTCAGCAGGTTGTTGCCGCAGTGGAGGGGAAATAATGTCCATTTTAATCGATAAAAACACCAAGGTTATCTGCCAGGGCTTTACCGGTAGCCAGGGGGCTTTCCACTCa                                                                                                   >  1:420223/1‑133 (MQ=255)
  cagcaAAAGGTCTGACGGATGCAGCTCAGCAGGTTGTTGCCGCAGTGGAGGGGAAATAATGTCCATTTTAATCGATAAAAACACCAAGGTTATCTGCCAGGGCTTTACCGGTAGCCAGGGGGCTTTCCACTCa                                                                                                   >  1:496256/1‑133 (MQ=255)
       aaaGGTCTGATGGATCCAGCTCAGCAGGGTGTGGCCGCAGTGGAGGGGAAATAATGTCCATTTTAATCGATAAAAACACCAAGGTTATCTGCCAGGGCTTTACCGGTAGCCAGGGGGCTTTCCACTCAGAACa                                                                                              <  2:149090/133‑1 (MQ=255)
          ggTCTGACGGATGCAGCTCAGCAGGTTGTTGCCGCAGTGGAGGGGAAATAATGTCCATTTTAATCGATAAAAACACCAAGGTTATCTGCCAGGGCTTTACCGGTAGCCAGGGGGATTTCCACTCAGAACAGGc                                                                                           >  1:221107/1‑133 (MQ=255)
           gTCTGACGGATGCAGCTCAGCAGGTTGTTGCCGCAGTGGAGGGGAAATAATGTCCATTTTAATCGATAAAAACACCAAGGTTATCTGCAAGGGCTTTACCGGTAGCCAGGGGACTTTCCACTCAGAACAGGcc                                                                                          <  2:86477/133‑1 (MQ=255)
                   gATGCAGCTCCGCAGGTTGTTGCCTCAGTGGAAGGGAAATACTGTCCATTTTAATCGATAAAAACACCAAGGTTATCTGCCAGGGCTTTACCGGTAGCCAGGGGGCTTTCCACTCAGAACAGGCCATTGCATa                                                                                  <  2:221107/133‑1 (MQ=255)
                     tGCAGCTCAGCAGGTTGTTGCCGCAGTGGAGGGGAAATAATGTCCATTTTAATCGATAAAAACACCAAGGTTATCTGCCAGGGCTTTACCGGTAGCCAGGGGGCTTTCCACTCAGAACAGGCCATTGCATACg                                                                                >  2:79312/1‑133 (MQ=255)
                                  gttgttGCCGCAGTGGAGGGGAAATAATGTCCATTTTAATCGATAAAAACACCAAGGTTATCTGCCAGGGCTTTACCGGTAGCCAGGGGGCTTTCCACTCAGAACAGGCCATTGCATACGGCACTAAAATGGt                                                                   >  2:361259/1‑133 (MQ=255)
                                           gCAGTGGAGGGGAAATAATGTCCATTTTAATCGATAAAAACACCAAGGTTATCTGCCAGGGCTTTACCGGTAGCCAGGGGGCTTTCCACTCAGAACAGGCCATTGCATACGGCACTAAAATGGTTGGCGGCGt                                                          <  2:66199/133‑1 (MQ=255)
                                                                   ttttAATCGATAAAAACACCAAGGTTATCTGCCAGGGCTTTACCGGTAGCCAGGGGGCTTTCCACTCAGAACAGGCCATTGCATACGGCACTAAAATGGTTGGCGGCGTAACCCCAGGTAAAGGCGGcaccac                                  >  1:514366/1‑133 (MQ=255)
                                                                   ttttAATCGATAAAAACACCAAGGTTATCTGCCAGGGCTTTACCGGTAGCCAGGGGGCTTTCCACTCAGAACAGGCCATTGCATACGGCACTAAAATGGTTGGCGGCGTAACCCCAGGTAAAGGCGGcaccac                                  >  2:406130/1‑133 (MQ=255)
                                                                       aaTCGATAAAAACACCAAGGTTATCTGCCAGGGCTTTACCGGTAGCCAGGGGGCTTTCCACTCAGAACAGGCCATTGCATACGGCACTAAAATGGTTGGCGGCGTAACCCCAGGTAAAGGCGGCAccacccac                              <  1:406130/133‑1 (MQ=255)
                                                                         tCGATAAAAACACCAAGGTTATCTGCCAGGGCTTTACCGGTAGCCAGGGGGCTTTCCACTCAGAACAGGCCATTGCATACGGCACTAAAATGGTTGGCGGCGTAACCCCAGGTAAAGGCGGCACCACCCACCt                            <  2:514366/133‑1 (MQ=255)
                                                                           gATAAAAACACCAAGGTTATCTGCCAGGGCTTTACCGGTAGCCAGGGGGCTTTCCACTCAGAACAGGCCATTGCATACGGCACTAAAATGGTTGGCGGCGTAACCCCAGGTAAAGGCGGCACCACCCACCTCg                          >  2:504864/1‑133 (MQ=255)
                                                                                   cacCAAGGTTGTCTGCCAGGGCTTCACCTGTATCCAGGGGGCTTTCCACTCAGAACAGGCCATTGCATACGGCACTAAAATGGTTGGCGGCGTAACCCCAGGTAAAGGCGGCACCACCCACCTCGGCCTGCCg                  <  2:219632/133‑1 (MQ=255)
                                                                                   cacCAAGGTTATCTGCCAGGGCTTTACCTGTACCCAGGGGGGTTTCCACTCAGAACAGGCCATTGCATACGGCACTAAAATGGTTGGCGGCGTAACCCCAGGTAAAGGCGGCACCACCCACCTCGGCCTGCCg                  <  2:176083/133‑1 (MQ=255)
                                                                                   cacCAAGGTTATCTGCCAGGGCTTTACCGGTAGCCAGGGGGCTTTCCACTCAGAACAGGCCATTGCATACGGCACTAAAATGGTTGGCGGCGTACCCCCAGGTAAAGGCGGCACCACCCACCTCGGCCTGCCg                  <  2:335067/133‑1 (MQ=255)
                                                                                   cacCAAGGTTATCTGCCAGGGCTTTACCGGTAGCCAGGGGGCTTTCCACTCAGAACAGGCCATTGCATACGGCACTAAAATGGTTGGCGGCGTAACCCCAGGTAAAGGCGGCACCACCCACCTCGGCCTGCCg                  <  2:339912/133‑1 (MQ=255)
                                                                                   cacCAAGGTTATCTGCCAGGGCTTTACCGGTAGCCAGGGGGCTTTCCACTCAGAACAGGCCATTGCATACGGCACTAAAATGGTTGGCGGCGTAACCCCAGGTAAAGGCGGCACCACCCACCTCGGCCTGCCg                  <  2:420223/133‑1 (MQ=255)
                                                                                   cacCAAGGTTATCTGCCAGGGCTTTACCGGTAGCCAGGGGGCTTTCCACTCAGAACAGGCCATTGCATACGGCACTAAAATGGTTGGCGGCGTAACCCCAGGTAAAGGCGGCACCACCCACCTCGGCCTGCCg                  <  2:496256/133‑1 (MQ=255)
                                                                                   cacCAAGGTTATCTGCCAGGGCTTTACCGGTAGCCAGGGGGCTTTCCACTCAGAACAGGCCATTGCATACGGCACTAAAATGGTTGGCGGCGTAACCCCAGGTAAAGGCGGCACCACCCACCTCGGCCTGCCg                  <  1:73541/133‑1 (MQ=255)
                                                                                   cacCAAGGTTATCTGCCAGGGCTTTACCGGTAGCCAGGGGGCTTTCCACTCAGAACAGGCCATTGCATACGGCACTAAAATGGTTGGCGGCGTAACCCCAGGTAAAGGCGGCACCACCCACCTCGGCCTGCCg                  <  1:361259/133‑1 (MQ=255)
                                                                                      cAAGGTTATCTGCCAGGGCTTTACCGGTAGCCAGGGGGCTTTCCACTCAGAACAGGCCATTGCATACGGCACTAAAATGGTTGGCGGCGTAACCCCAGGTAACGGCGGCAACACCCAGCTCGGCCTGCCGgtg               >  1:202642/1‑133 (MQ=255)
                                                                                              tCTGCCAGGGTTTTACCGGTAGCCAGGGCGCTGTCCGCTCAAAACAGGCCATTGCATACGGCACTAAAATGGTTGGCGGCGTAACCCCAGGTAAAGGCGGCCCCACCCACCTCGGCCTGCCGGTGTTCAacac       <  2:202642/133‑1 (MQ=255)
                                                                                              tCTGCCAGGGCTTTACCGGTAGCCAGGGGGCTTTCCACTCAGAACAGGCCATTGCATACGGCACTAAAATGGTTGGCGGCGTAACCCCAGGTAAAGGCGGCACCACCCACCTCGGCCTGCCGGTGTTCAacac       <  1:79312/133‑1 (MQ=255)
                                                                                                   cAGGGCTTTACCGGTAGCCAGGGGGCTTTCCACTCAGAACAGGCCATTGCATACGGCACTAAAATGGTTGGCGGCGTAACCCCAGGTAAAGGCGGCACCACCCACCTCGGCCTGCCGGTGTTCAACACcgtgc  >  2:314410/1‑133 (MQ=255)
                                                                                                                           |                                                                                                            
TGCAGCAAAAGGTCTGACGGATGCAGCTCAGCAGGTTGTTGCCGCAGTGGAGGGGAAATAATGTCCATTTTAATCGATAAAAACACCAAGGTTATCTGCCAGGGCTTTACCGGTAGCCAGGGGACTTTCCACTCAGAACAGGCCATTGCATACGGCACTAAAATGGTTGGCGGCGTAACCCCAGGTAAAGGCGGCACCACCCACCTCGGCCTGCCGGTGTTCAACACCGTGC  >  NC_000913/764120‑764351

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