Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB143 5,390,023 Δ1 bp 24.3% coding (117/1666 nt) ligH → 10‑formyltetrahydrofolate synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1435,390,0230G.24.3% 102.9 / 29.9 37coding (117/1666 nt)ligH10‑formyltetrahydrofolate synthetase
Reads supporting (aligned to +/- strand):  ref base G (18/10);  new base . (4/5);  total (22/15)
Fisher's exact test for biased strand distribution p-value = 4.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAAGGACGAGGAGGACACATGGCACAAGGCTCCGATATCGAGATCGCCCGCGAAGCGAAGATGCAGAATATCGCGGAGGTCGGCGCGAAAGTCGGCATCCCGCAGGACGCGCTGCTGAACTACGGCCCCTACAAGGCCAAGCTGAGCTGGGACTTCATCAACAGCGTTCAGGGCAACCAGGACGGCAAGCTGATCCTCGTCACCGCGATCAACCCGACGCCGGCCGGCGAAGGCAAGACCACGACCACCGTGGGCCTGGCCGACGGCC  >  ACB143/5389888‑5390156
                                                                                                                                       |                                                                                                                                     
gaAAGGTCGAGGAGGACACATGGCACAAGGCTCCGATATCGAGATCGCCCGCGAAGCGAAGATGCAGAATATCGCGGAGGTCGGCGCGAAAGTCGGCATCCCGCAGGACGCGCTGCTGAACTACGGCCCCTACAAg                                                                                                                                       <  2:145416/136‑1 (MQ=255)
    ggACGAGGAGGACACATGGCACAAGGCTCCGATATCGAGATCGCCCGCGAAGCGAAGATGCAGAATATCGCGGAGGTCGGCGCGAAAGTCGGCATCCCGCAGGACGCGCTGCTGAACTACGGCCCCTACAAgccaa                                                                                                                                   <  3:282758/136‑5 (MQ=255)
      aCGAGGAGGACACATGGCACAAGGCTCCGATATCGAGATCGCCCGCGAAGCGAAGATGCAGAATATCGCGGAGGTCGGCGCGAAAGTCGGCATCCCGCAGGACGCGCTGCTGAACTACGGCCCCTAcaa‑gccaagc                                                                                                                                >  3:91116/1‑136 (MQ=255)
                   aTGGCACAAGGCTCCGATATCGAGATCGCCCGCGAAGCGAAGATGCAGAATATCGCGGAGGTCGGCGCGAAAGTCGGCATCCCGCAGGACGCGCTGCTGAACTACGGCCCCTACAAGGCCAAGCTGAGCTGGGACt                                                                                                                    >  1:231218/1‑136 (MQ=255)
                        acaAGGCTCCGATATCGAGATCGCCCGCGAAGCGAAGATGCAGAATATCGCGGAGGTCGGCGCGAAAGTCGGCATCCCGCAGGACGCGCTGCTGAACTACGGCCCCTACAAGGCCAAGCTGAGCTGGGACTtcat                                                                                                                >  5:344781/1‑135 (MQ=255)
                           aGGCTCCGATATCGAGATCGCCCGCGAAGCGAAGATGCAGAATATCGCGGAGGTCGGCGCGAAAGTCGGCATCCCGCAGGACGCGCTGCTGAACTACGGCCCCTACAA‑GCCAAGCTGAGCTGGGACTTCATCaa                                                                                                             >  2:131261/1‑134 (MQ=255)
                             gCTCCGATATCGAGATCGCCCGCGAAGCGAAGATGCAGAATATCGCGGAGGTCGGCGCGAAAGTCGGCATCCCGCAGGACGCGCTGCTGAACTACGGCCCCTACAAGGCCAAGCTGAGCTGGGACTTCATCAACAg                                                                                                          <  2:303944/136‑1 (MQ=255)
                             gCTCCGATATCGAGATCGCCCGCGAAGCGAAGATGCAGAATATCGCGGAGGTCGGCGCGAAAGTCGGCATCCCGCAGGACGCGCTGCTGAACTACGGCCCCTACAA‑GCCAAGCTGAGCTGGGACTTCATCAACAGc                                                                                                         <  3:386314/136‑1 (MQ=255)
                              cTCCGATATCGAGATCGCCCGCGAAGCGAAGATGCAGAATATCGCGGAGGTCGGCGCGAAAGTCGGCATCCCGCAGGACGCGCTGCTGAACTACGGCCCCTACAAGGCCAAGCTGAGCTGGGACTTCATCAACAGc                                                                                                         <  6:231861/136‑1 (MQ=255)
                              cTCCGATATCGAGATCGCCCGCGAAGCGAAGATGCAGAATATCGCGGAGGTCGGCGCGAAAGTCGGCATCCCGCAGGACGCGCTGCTGAACTACGGCCCCTACAAGGCCAAGCTGAGCTGGGACTTCATCAACAGc                                                                                                         <  1:323637/136‑1 (MQ=255)
                                 cGATATCGAGATCGCCCGCGAAGCGAAGATGCAGAATATCGCGGAGGTCGGCGCGAAAGTCGGCATCCCGCAGGACGCGCTGCTGAACTACGGCCCCTACAA‑GCCAAGCTGAGCTGGGACTTCATCAACAGCGTTc                                                                                                     <  2:33734/136‑1 (MQ=255)
                                  gATATCGAGATCGCCCGCGAAGCGAAGATGCAGAATATCGCGGAGGTCGGCGCGAAAGTCGGCATCCCGCAGGACGCGCTGCTGAACTACGGCCCCTACAAGGCCAAGCTGAGCTGGGACTTCATCAACAGCGtt                                                                                                      >  3:93597/1‑135 (MQ=255)
                                  gATATCGAGATCGCCCGCGAAGCGAAGATGCAGAATATCGCGGAGGTCGGCGCGAAAGTCGGCATCCCGCAGGACGCGCTGCTGAACTACGGCCCCTACAAGGCCAAGCTGAGCTGGGACTTCATCAACAGCGTTc                                                                                                     <  3:341691/136‑1 (MQ=255)
                                   atatCGAGATCGCCCGCGAAGCGAAGATGCAGAATATCGCGGAGGTCGGCGCGAAAGTCGGCATCCCGCAGGACGCGCTGCTGAACTACGGCCCCTACAA‑GCCAAGCTGAGCTGGGACTTCATCAACAGCGTTCAg                                                                                                   >  7:246916/1‑136 (MQ=255)
                                    tatCGAGATCGCCCGCGAAGCGAAGATGCAGAATATCGCGGAGGTCGGCGCGAAAGTCGGCATCCCGCAGGACGCGCTGCTGAACTACGGCCCCTACAAGGCCAAGCTGAGCTGGGACTtcatca                                                                                                              <  3:129337/125‑1 (MQ=255)
                                    tatCGAGATCGCCCGCGAAGCGAAGATGCAGAATATCGCGGAGGTCGGCGCGAAAGTCGGCATCCCGCAGGACGCGCTGCTGAACTACGGCCCCTACAAGGCCAAGCTGAGCTGGGACTtcatca                                                                                                              >  4:129337/1‑125 (MQ=255)
                                         agaTCGCCCGCGAAGCGAAGATGCAGAATATCGCGGAGGTCGGCGCGAAAGTCGGCATCCCGCAGGACGCGCTGCTGAACTACGGCCCCTACAA‑GCCAAGCTGAGCTGGGACTTCATCAACAGCGTTCAGGGCAAc                                                                                             >  6:367115/1‑136 (MQ=255)
                                           aTCGCCCGCGAAGCGAAGATGCAGAATATCGCGGAGGTCGGCGCGAAAGTCGGCATCCCGCAGGACGCGCTGCTGAACTACGGCCCCTACAAGGCCAAGCTGAGCTGGGACTTCATCAACAGCGTTCAGGGCacc                                                                                             >  3:14479/1‑133 (MQ=255)
                                              gCCCGCGAAGCGAAGATGCAGAATATCGCGGAGGTCGGCGCGAAAGTCGGCATCCCGCAGGACGCGCTGCTGAACTACGG‑CCCTACAAGGCCAAGCTGAGCTGGGACTTCATCAACAGCGTTCAGGGCAACCAGGa                                                                                        >  5:367104/1‑136 (MQ=255)
                                                  gcgaagcgaaGATGCAGAATATCGCGGAGGTCGGCGCGAAAGTCGGCATCCCGCAGGACGCGCTGCTGAACTACGGCCCCTACAAGGCCAAGCTGAGCTGGGACTTCATCAACAGCGTTCAGGGCAACCAGGACgg                                                                                     >  8:124970/1‑136 (MQ=255)
                                                  gcgaagcgaaGATGCAGAATATCGCGGAGGTCGGCGCGAAAGTCGGCATCCCGCAGGACGCGCTGCTGAACTACGGCCCCTACAAGGCCAAGCTGAGCTGGGACTTCATCAACAGCGTTCAGGGCAACCAGGACgg                                                                                     >  8:190373/1‑136 (MQ=255)
                                                    gaagcgaagATGCAGAATATCGCGGAGGTCGGCGCGAAAGTCGGCATCCCGCAGGACGCGCTGCTGAACTACGGCCCCTACAAGGCCAAGCTGAGCTGGGACTTCATCAACAGCGTTCAGGGCAACCAGGACGGCa                                                                                   <  3:56523/136‑1 (MQ=255)
                                                                 aGAATATCGCGGAGGTCGGCGCGAAAGTCGGCATCCCGCAGGACGCGCTGCTGAACTACGGCCCCTACAAGGCCAAGCTGAGCTGGGACTTCATCAACAGCGTTCAGGGCAACCAGGACGGCAAGCTGATCCTCGt                                                                      <  8:103782/136‑1 (MQ=255)
                                                                       tCGCGGAGGTCGGCGCGAAAGTCGGCATCCCGCAGGACGCGCTGCTGAACTACGGCCCCTACAA‑GCCAAGCTGAGCTGGGACTTCATCAACAGCGTTCAGGGCAACCAGGACGGCAAGCTGATCCTCGTCACcgc                                                                <  1:342710/135‑1 (MQ=255)
                                                                       tCGCGGAGGGCGGCGCGAAAGTCGGCATCCCGCAGGACGCGCTGCTGAACTACGGCCCCTACAA‑GCCAAGCTGAGCTGGGACTTCATCAACAGCGTTCAGGGCAACCAGGACGGCAAGCTGATCCTCGTCACcgc                                                                <  6:194346/135‑1 (MQ=255)
                                                                           ggAGGTCGGCGCGAAAGTCGGCATCCCGCAGGACGCGCTGCTGAACTACGGCCCCTACAAGGCCAAGCTGAGCTGGGACTTCATCAACAGCGTTCAGGGCAACCAGGACGGCAAGCTGATCCTCGTCACCGCGAt                                                             >  7:92936/1‑135 (MQ=255)
                                                                                tCGGCGCGAAAGTCGGCATCCCGCAGGACGCGCTGCTGAACTACGGCCCCTCCAAGGCCAAGCTGAGCTGGGACTTCATCAACAGCGTTCAGGGCAACCAGGACGGCAAGCTGATCCTCGTCACCGCGATCATccc                                                       >  7:80313/1‑136 (MQ=255)
                                                                                  ggCGCGAAAGTCGGCATCCCGCAGGACGCGCTGCTGAACTACGG‑CCCTACAAGGCCAAGCTGAGCTGGGACTTCATCAACAGCGTTCAGGGCAACCAGGACGGCAAGCTGATCCTCGTCACCGCGATCAACCCgcc                                                    >  7:292874/1‑134 (MQ=255)
                                                                                                cATCCCGCAGGACGCGCTGCTTAACTACGGCCCCTACAA‑GCCAAGCTGTGCTGGTACTTCATCAACAGCGTTCAGGGAAACCAGGACGGCATGCTTATCCTCGTCACCGCTATCAACCCGACGCCGGCCGGCGAAg                                      <  6:41722/136‑1 (MQ=255)
                                                                                                    ccGCAGGACGCGCTGCTGAACTACGGCCCCTACAAGGCCAAGCTGAGCTGGGACTTCATCAACAGCGTTCAGGGCAACCAGGACGGCAAGCTGATCCTCGTCACCGCGATCAACCCGACGCCGGCCGGCGAAGGc                                    <  1:282545/135‑1 (MQ=255)
                                                                                                     cGCAGGACGCGCTGCTGAACTACGGCCCCTACAAGGCCAAGCTGAGCTGGGACTTCATCAACAGCGTTCAGGGCAACCAGGACGGCAAGCTGATCCTCGTCACCGCGATCAACCCGACGCCGGCCGGCGAAGGCaa                                  <  7:152144/136‑1 (MQ=255)
                                                                                                           aCGCGCTGCTGAACTACGGCCCCTACAAGGCCAAGCTGAGCTGGGACTTCATCAACAGCGTTCAGGGCAACCAGGACGGCAAGCTGATCCTCGTCACCGCGATCAACCCGACGCCGGCCGGCGAAGGCAAgaccac                            <  1:405323/136‑1 (MQ=255)
                                                                                                             gcgcTGCTGAACTACGGCCCCTACAAGGCCAAGCTGAGCTGGGACTTCATCAACAGCGTTCAGGGCAACCAGGACGGCAAGCTGATCCTCGTCACCGCGATCAACCCGACGCCGGCCGGCGAAGGCAAgaccacga                          >  1:271169/1‑136 (MQ=255)
                                                                                                                              ccccTACAAGGCCAAGCTGAGCTGGGACTTCATCAACAGCGTTCAGGGCAACCAGGACGGCAAGCTGATCCTCGTCACCGCGATCAACCCGACGCCGGCCGGCGAAGGCAAGACCACGACCACCGTGGGCCTGGcc         >  1:101465/1‑136 (MQ=255)
                                                                                                                              ccccTACAAGGCCAAGCTGAGCTGGGACTTCATCAACAGCGTTCAGGGCAACCAGGACGGCAAGCTGATCCTCGTCACCGCGATCAACCCGACGCCGGCCGGCGAAGGCAAGACCACGACCACCGTGGGCCTGGcc         >  1:165693/1‑136 (MQ=255)
                                                                                                                              ccccTACAAGGCCAAGCTGAGCTGGGACTTCATCAACAGCGTTCAGGGCAACCAGGACGGCAAGCTGATCCTCGTCACCGCGATCAACCCGACGCCGGCCGGCGAAGGCAAGACCACGACCACCGTGGGCCTGGcc         >  7:229669/1‑136 (MQ=255)
                                                                                                                              ccccTACAA‑GCCAAGCTGAGCTGGGACTTCATCAACAGCGTTCAGGGCAACCAGGACGGCAAGCTGATCCTCGTCACCGCGATCAACCCGACGCCGGCCGGCGAAGGCAAGACCACGACCACCGTGGGCCTGGcc         >  1:263013/1‑135 (MQ=255)
                                                                                                                                   aCAAGGCCAAGCTGAGCTGGGACTTCATCAACAGCGTTCAGGGCAACCAGGACGGCAAGCTGATCCTCGTCACCGCGATCAACCCGACGCCGGCCGGCGAAGGCAAGACCACGACCACCGTGGGCCTGGCCGACgg    >  7:307918/1‑136 (MQ=255)
                                                                                                                                     aaGGCCAAGCTGAGCTGGGACTTCATCAACAGCGTTCAGGGCAACCAGGACGGCAAGCTGATCCTCGTCACCGCGATCAACCCGACGCCGGCCGGCGAAGGCACGACCACGCCCCCCGTGGGCCAGGCCGCCGccc  >  3:296278/1‑136 (MQ=255)
                                                                                                                                     aaGGCCAAGCTGAGCTGGGACTTCATCAACAGCGTTCAGGGCAACCAGGACGGCAAGCTGATCCTCGTCACCGCGATCAACCAGACGCCGGCCGGCGAAGGCAAGACCACGACCACCGTGGGCCTGGCCGACGGcc  >  5:115829/1‑136 (MQ=255)
                                                                                                                                       |                                                                                                                                     
GAAAGGACGAGGAGGACACATGGCACAAGGCTCCGATATCGAGATCGCCCGCGAAGCGAAGATGCAGAATATCGCGGAGGTCGGCGCGAAAGTCGGCATCCCGCAGGACGCGCTGCTGAACTACGGCCCCTACAAGGCCAAGCTGAGCTGGGACTTCATCAACAGCGTTCAGGGCAACCAGGACGGCAAGCTGATCCTCGTCACCGCGATCAACCCGACGCCGGCCGGCGAAGGCAAGACCACGACCACCGTGGGCCTGGCCGACGGCC  >  ACB143/5389888‑5390156

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