Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,493,127 T→C D43G (GAC→GGC)  frc ← formyl‑CoA transferase, NAD(P)‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,493,1270TC88.2% 97.9 / 3.8 34D43G (GAC→GGC) frcformyl‑CoA transferase, NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base T (2/2);  new base C (15/15);  total (17/17)
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.32e-01

CGGATCAGCTTTTCCATTACCTCTTTGCCTTCCGCTGTTTTGGTATTTAACTCAATAGAACGTTTGTTACTGTTAAGCATGGTGAAGTAAAGCGCATCGATATCAGGAATATCTCGCAGCTGGTGACGCGTTACGTCACCAACGCCGGGACGTTCAATTTTAATAACGTCAGCGCCAAACCAGGCCAGCATTTGAGTACAAGATGGGCCAGATTGCACACCGGTGAAATCGAGAACTTTAATTCCTTGAAGTGGAGTTGACATAGCAAT  >  NC_000913/2492992‑2493260
                                                                                                                                       |                                                                                                                                     
cgGATCAGCTTTTCCATTACCTCTTTGCCTTCCGCTGTTTTGGTATTTAACTCAATAGAACGTTTGTTACTGTTAAGCATGGTGAAGTAAAGCGCATCGATATCAGGAATATCTCGCAGCTGGTGACGCGTTACGCCAc                                                                                                                                    <  1:125440/139‑1 (MQ=255)
                                    gTTGTGGTATTTAACTCAATAGAACTTTTGTTACTGTTAAGCATGGTGAAGTAAAGCGCATCGATATCAGGAATATCTCGCAGCTGGTGACGCGTTACGCCACCAACGCCGAGACGTTCAATTTTAATAACGTCAgcgc                                                                                                <  2:398364/139‑1 (MQ=255)
                                     ttttGGTATTTAACTCAATAGAACGTTTGTTACTGTTAAGCATGGTGAAGTAAAGCGCATCGATATCAGGAATATCTCGCAGCTGGTGACGCGTTACGCCACCAACGCCGGGACGt                                                                                                                      >  2:434157/1‑116 (MQ=255)
                                     ttttGGTATTTAACTCAATAGAACGTTTGTTACTGTTAAGCATGGTGAAGTAAAGCGCATCGATATCAGGAATATCTCGCAGCTGGTGACGCGTTACGCCACCAACGCCGGGACGt                                                                                                                      <  1:434157/116‑1 (MQ=255)
                                       ttGGTATTTAACTCAATAGAACGTTTGTTACTGTTAAGCATGGTGAAGTAAAGCGCATCGATATCAGGAATATCTCGCAGCTGGTGACGCGTTACGCCACCAACGCCGGGACGTTCAATTTTAATAACGTCAGCGCCaa                                                                                             <  1:447732/139‑1 (MQ=255)
                                                      aTAGAACGTTTGTTACTGTTAAGCATGGTGAAGTAAAGCGCATCGATATCAGGAATATCTCGCAGCTGGTGACGCGTTACGCCACCAACGCCGGGACGTTCAATTTTAATAACGTCAGCGCCAAACCAGGCCAGCAttt                                                                              <  2:403295/139‑1 (MQ=255)
                                                               ttGTTACTGTTAAGCATGGTGAAGTAAAGCGCATCGATATCAGGAATATCTCGCAGCTGGTGACGCGTTACGCCACCAACGCCGGGACGTTCAATTTTAATAACGTCAGCGCCAAACCAg                                                                                        <  2:480694/120‑1 (MQ=255)
                                                               ttGTTACTGTTAAGCATGGTGAAGTAAAGCGCATCGATATCAGGAATATCTCGCAGCTGGTGACGCGTTACGCCACCAACGCCGGGACGTTCAATTTTAATAACGTCAGCGCCAAACCAg                                                                                        >  1:480694/1‑120 (MQ=255)
                                                                 gTTACTGTTAAGCATGGTGAAGTAAAGCGCATCGATATCAGGAATATCTCGCAGCTGGTGACGCGTTACGCCACCAACGCCGGGACGTTCAATTTTAATAACGTCAGCGCCAAACCAGGCCAGCATTTGAGTACAAGAt                                                                   >  1:52801/1‑139 (MQ=255)
                                                                          aaGCATGGTGAAGTAAAGCGCATCGATATCAGGAATATCTCGCAGCTGGTGACGCGTTACGTCACCAACGCCGGGACGTTCAATTTTAATAACGTCAGCGCCAAACCAGGCCAGCATTTGAGTa                                                                         >  2:608828/1‑124 (MQ=255)
                                                                          aaGCATGGTGAAGTAAAGCGCATCGATATCAGGAATATCTCGCAGCTGGTGACGCGTTACGTCACCAACGCCGGGACGTTCAATTTTAATAACGTCAGCGCCAAACCAGGCCAGCATTTGAGTa                                                                         <  1:608828/124‑1 (MQ=255)
                                                                           aGCATGGTGAAGTAAAGCGCATCGATATCAGGAATATCTCGCAGCTGGTGACGCGTTACGCCACCAACGCCGGGACGTTCAATTTTAATAACGTCAGCGCCAAACCAGGCCAGCATTTGAGTACAAGATGGGCCAGAtt                                                         >  2:332238/1‑139 (MQ=255)
                                                                            gCATGGTGAAGTAAAGCGCATCGATATCAGGAATATCTCGCAGCTGGTGACGCGTTACGTCACCAACGCCGGGACGTTCAATTTTAATAACGTCAGCGCCAAACCAGGCCAGCAttt                                                                              >  2:455428/1‑117 (MQ=255)
                                                                            gCATGGTGAAGTAAAGCGCATCGATATCAGGAATATCTCGCAGCTGGTGACGCGTTACGTCACCAACGCCGGGACGTTCAATTTTAATAACGTCAGCGCCAAACCAGGCCAGCAttt                                                                              <  1:455428/117‑1 (MQ=255)
                                                                                 gtataagagacagcATCGATATCAGGAATATCTCGCAGCTGGTGACGCGTTACGCCACCAACGCCGGGACGTTCAATTTTACTAACGTCAGCGCCAAACCAGGCCAGCATTTGAGTACAAGATGGGCCAGATTGCacac                                                   >  2:72130/13‑139 (MQ=255)
                                                                                  tGAAGTAAAGCGCATCGATATCAGGAATATCTCGCAGCTGGTGACGCGTTACGCCACCAACGCCGGGACGTTCAATTTTAATAACGTCAGCGCCAAACCAGGCCAGCAttt                                                                              >  2:63743/1‑111 (MQ=255)
                                                                                  tGAAGTAAAGCGCATCGATATCAGGAATATCTCGCAGCTGGTGACGCGTTACGCCACCAACGCCGGGACGTTCAATTTTAATAACGTCAGCGCCAAACCAGGCCAGCAttt                                                                              <  1:63743/111‑1 (MQ=255)
                                                                                  tGAAGTAAAGCGCATCGATATCAGGAATATCTCGCAGCTGGTGACGCGTTACGCCACCAACGCCGGGACGTTCAATTTTAATAACGTCAGCGCCAAACCAGGCCAGCATTTGAGTACAAGATGGGCCAGATTGCACAcc                                                  >  1:84832/1‑139 (MQ=255)
                                                                                     aGTAAAGCGCATCGATATCAGGAATATCTCGCAGCTGGTGACGCGTTACGCCACCAACGCCGGGACGTTCAATTTTAATAACGTCAGCGCCAAACCAGGCCAGCATTTGAGTACAAGATGGGCCAGATTGCACACCGGt                                               >  2:148186/1‑139 (MQ=255)
                                                                                       tAAAGCGCATCGATATCAGGAATATCTCGCAGCTGGTGACGCGTTACGCCACCAACGCCGGGACGTTCAATTTTAATAACGTCAGCGCCAAACCAGGCCAGCATTTGAGTACAAGATGGGCCAGATTGCACACCGGTGa                                             >  1:453044/1‑139 (MQ=255)
                                                                                          aGCGCATCGATATCAGGAATATCTCGCAGCTGGTGACGCGTTACGCCACCAACGCCGGGACGTTCAATTTTAATAACGTCAGCGCCAAACCAGGCCAGCATTTGAGTACAAGATGGGCCAGATTGCACACCGGTGAAAt                                          >  2:327808/1‑139 (MQ=255)
                                                                                            cgcATCGATATCAGGAATATCTCGCAGCTGGTGACGCGTTACGCCACCAACGCCGGGACGTTCAATTTTAATAACGTCAGCGCCAAACCAGGCCAGCAttt                                                                              >  2:252012/1‑101 (MQ=255)
                                                                                            cgcATCGATATCAGGAATATCTCGCAGCTGGTGACGCGTTACGCCACCAACGCCGGGACGTTCAATTTTAATAACGTCAGCGCCAAACCAGGCCAGCAttt                                                                              <  1:252012/101‑1 (MQ=255)
                                                                                            cgcATCGATATCAGGAATATCTCGCAGCTGGTGACGCGTTACGCCACCAACGCCGGGACGTTCAATTTTAATAACGTCAGCGCCAAACCAGGCCAGCATTTGAGTACAAGATGGGcc                                                              <  1:497720/117‑1 (MQ=255)
                                                                                            cgcATCGATATCAGGAATATCTCGCAGCTGGTGACGCGTTACGCCACCAACGCCGGGACGTTCAATTTTAATAACGTCAGCGCCAAACCAGGCCAGCATTTGAGTACAAGATGGGcc                                                              >  2:497720/1‑117 (MQ=255)
                                                                                            cgcATCGATATCAGGAATATCTCGCAGCTGGTGACGCGTTACGCCACCAACGCCGGGACGTTCAATTTTAATAACGTCAGCGCCAAACCAGGCCAGCATTTGAGTACAAGATGGGCCAGATTGCACACCGGTGAAATCg                                        >  1:335859/1‑139 (MQ=255)
                                                                                                 cGATATCAGGAATATCTCGCAGCTGGTGACGCGTTACGCCACCAACGCCGGGACGTTCAATTTTAATAACGTCAGCGCCAAACCAGGCCAGCATTTGAGTACAAGATGGGCCAGATTGCACACCGGTGAAATCGAGAAc                                   >  2:121012/1‑139 (MQ=255)
                                                                                                             tatCTCGCAGCTGGTGACGCGTTACGCCACCAACGCCGGGACGTTCAATTTTAATAACGTCAGCGCCAAACCAGGCCAGCATTTGAGTACAAGATGGGCCAGATTGCACACCGGTGAAATCGAGAACTTTAATTCCt                         >  1:545628/1‑137 (MQ=255)
                                                                                                             tatCTCGCAGCTGGTGACGCGTTACGCCACCAACGCCGGGACGTTCAATTTTAATAACGTCAGCGCCAAACCAGGCCAGCATTTGAGTACAAGATGGGCCAGATTGCACACCGGTGAAATCGAGAACTTTAATTCCt                         <  2:545628/137‑1 (MQ=255)
                                                                                                                      gCTGGTGACGCGTTACGCCACCAACGCCGGGACGTTCAATTTTAATAACGTCAGCGCCAAACCAGGCCAGCATTTGAGTACAAGATGGGCCAGATTGCACACCGGTGAAATCGAGAACTTTAATTCCTTGAAGTGGAGt              <  2:335859/139‑1 (MQ=255)
                                                                                                                          gTGACGCGTTACGCCACCAACGCCGGGACGTTCAATTTTAATAACGTCAGCGCCAAACCAGGCCAGCATTTGAGTACAAGATGGGCCAGATTGCACACCGGTGAAATCGAGAACTTTAATTCCTTGAAGTGGAGTTGAc          <  1:332238/139‑1 (MQ=255)
                                                                                                                           tGACGCGTTACGCCACCAACGCCGGGACGTTCAATTTTAATAACGTCAGCGCCAAACCAGGCCAGCATTTGAGTACAAGATGGGCCAGATTGCACACCGGTGAAATCGAGAACTTTAATTCCTTGAAGTGGAGTTGACa         <  2:453044/139‑1 (MQ=255)
                                                                                                                              cgcgTTACGCCACCAACGCCGGGACGTTCAATTTTAATAACGTCAGCGCCAAACCAGGCCAGCATTTGAGTACAAGATGGGCCAGATTGCACACCGGTGAAATCGAGAACTTTAATTCCTTGAAGTGGAGTTGACATAg      <  1:121012/139‑1 (MQ=255)
                                                                                                                                  ttACGCCACCAACGCCGGGACGTTCAATTTTAATAACGTCAGCGCCAAACCAGGCCAGCATTTGAGTACAAGATGGGCCAGATTGCACACCGGTGAAATCGAGAACTTTAATTCCTTGAAGTGGAGTTGACATAGCAAt  <  1:252873/139‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                     
CGGATCAGCTTTTCCATTACCTCTTTGCCTTCCGCTGTTTTGGTATTTAACTCAATAGAACGTTTGTTACTGTTAAGCATGGTGAAGTAAAGCGCATCGATATCAGGAATATCTCGCAGCTGGTGACGCGTTACGTCACCAACGCCGGGACGTTCAATTTTAATAACGTCAGCGCCAAACCAGGCCAGCATTTGAGTACAAGATGGGCCAGATTGCACACCGGTGAAATCGAGAACTTTAATTCCTTGAAGTGGAGTTGACATAGCAAT  >  NC_000913/2492992‑2493260

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