Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 809,217 2 bp→AG coding (40‑41/1290 nt) bioA ← adenosylmethionine‑8‑amino‑7‑oxononanoate aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913809,2170CA92.3% 58.4 / ‑2.6 26W14M (TGG→CTG) bioAadenosylmethionine‑8‑amino‑7‑oxononanoate aminotransferase
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base A (9/15);  total (9/17)
Fisher's exact test for biased strand distribution p-value = 5.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_000913809,2180AG96.0% 76.5 / ‑3.9 26W14M (TGG→CTG) bioAadenosylmethionine‑8‑amino‑7‑oxononanoate aminotransferase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (9/15);  total (9/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 = 2.30e-01

TGCGGCACAGCTCAATGGCTGGCGCATGGGTGATACCGCCAAACATCACATGCGACATGGCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCCAGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTCGCAAGTCACAGAATTATTTGAAAAACCGTTGCTGGATCTGCTGTTTGAAGCGCAGCAGGTGCATCGCCAGCATTTCGATCCTCGTCA  >  NC_000913/808978‑809459
                                                                                                                                                                                                                                               ||                                                                                                                                                                                                                                                 
tGCGGCACAGCTCAATGGCTGGCGCATGGGTGATACCGCCAAACATCACATGCGACATGGCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCg                                                                                                                                                                                                                                          <  1:5450/250‑1 (MQ=255)
      acaGCTCAATGGCTGGCGCATGGGTGATACCGCCAAACATCACATGCGACATGGCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGCGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTc                                                                                                                                                                                                                                    >  1:266088/1‑250 (MQ=255)
                     gcgcATGGGTGACACCGCCAAACATCACATGCGACATGGCATCAATTTGCGCCTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATCCCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCCCCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTcc                                                                                                                                                                                                                    <  2:69077/251‑1 (MQ=255)
                       gcATGGGTGATACCGCCAAACAGCACATGCGACATGGCATCAATTTGCGACTGCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGt                                                                                                                                                                                                                  <  1:144878/251‑1 (MQ=255)
                                       cAAACATCACATGTGACATGGCATCAATTTGCGACTTCATCGCTGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCTCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCGTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAAACGACTTg                                                                                                                                                                                                  <  2:104405/251‑1 (MQ=255)
                                         aacatctcagtcGACATGGCATCAAGTGGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATTGCGTTGGTCAAAGGCAAGATCGTCCGTGGTCATAATCGACTTGTa                                                                                                                                                                                                <  2:266088/240‑1 (MQ=255)
                                               caaggCGACATGGCCTCAATTTGCGCCTTCAGCGCCGCACTAAGCTGTGGGTGATTGTAGCCGTGGATCGCGGCCCACCAGGACGACATACCGTCACCCAGGCGTCTGCCGTCAGATAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAATGCAAGAGCGTCCGTTGTCAGAATCGACTTGTAAACCac                                                                                                                                                                                          <  2:112423/247‑2 (MQ=255)
                                                        aTGGCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAACATGaaaa                                                                                                                                                                                  >  2:154827/1‑250 (MQ=255)
                                                            cATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGGAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGAttt                                                                                                                                                                             <  1:271971/251‑1 (MQ=255)
                                                              ctaaTTGGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGGAGCCGTGGACCGCCGCCCACCAGGACGACATTCCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCGCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCTTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGAttt                                                                                                                                                                             <  2:43424/247‑1 (MQ=255)
                                                                                     aTTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGTTGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATAGACTTGTAAACCAAATTGACAAGATTTAGGTTTCCAAGTCTACACCGAAtt                                                                                                                                                     >  1:12925/1‑250 (MQ=255)
                                                                                        aaGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTaac                                                                                                                                                  >  1:97599/1‑250 (MQ=255)
                                                                                                                      cgccCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAACTCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCAGCAGAGGGGAGGGCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAACCCAAATTGAAAAGATTTAGGTTTACAAATCTACACCGAATTAACAACAAAAAACACGATTTg                                                                                                                                >  2:220074/1‑238 (MQ=255)
                                                                                                                      cgccCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTg                                                                                                                                <  1:220074/238‑1 (MQ=255)
                                                                                                                                    acATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACg                                                                                                     >  2:41982/1‑251 (MQ=255)
                                                                                                                                               aCCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGt                                                                                          <  1:240317/251‑1 (MQ=255)
                                                                                                                                                           ccGTCAGACAAAATCCGCTCACAAACGTCTGCGCTCACCACCGGATAAACCGGCAGCGGCGGGGTCATGGGTGTGTATGGAGGCAGGATATGGCGTTGGGCAAAGGCAAGATCGTCCGTTGTCATAATCGACATGTAAACAAAATGGCACAGATTTAGGTTTACAAGTCTCCAACGACTTAAAAACAACAAACACGTTTTGGAGACGCCTCATGGGTCACAGCCCACGCTAGACGTTGTCggacttcaca                                                                               >  2:177944/1‑241 (MQ=255)
                                                                                                                                                             gTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCCAGATATGGTGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTa                                                                                                                                                                      <  1:146187/161‑1 (MQ=255)
                                                                                                                                                                aGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTCGCAAGTCACAGAatta                                                                         <  1:154827/251‑1 (MQ=255)
                                                                                                                                                                     aaaTCAGCTCGCAACCTTCGGCGCACACCTCCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCGATGGCTCACCGCCCACGCTGGACATTGTCGCAAGTCACAGAATTATTTGa                                                                    <  2:161638/251‑1 (MQ=255)
                                                                                                                                                                               gCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACACGTCTACACCGAATTCACAACAAAAAACACGTTTTGGAGAAGCCCCACAGCTCACCGCCCACGCTGGGCATTGGCGCAAGTCACAGAATTATTTGaaaaaacgatg                                                          >  2:64588/1‑245 (MQ=255)
                                                                                                                                                                                cAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTCGCAAGTCACAGAATTATTTGAAAAACCGTTGc                                                         <  1:282939/251‑1 (MQ=255)
                                                                                                                                                                                           gcacgcctccGGATAAACTGGCAGAGGGGAGATCAAGGATGTGCATAGGTGCCCGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAAACGACTTGTAAACCAAATTGAAAAGATTTATGTTTACGAGTCTACAACAAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTTGACATTGTCGCAAGTCACAGAATTATTTGAAAAACCGTTGCTGGATctgctg                                              <  1:159764/243‑1 (MQ=255)
                                                                                                                                                                                           gcTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTCGCAAGTCACAGAATTATTTGAAAAACCGTTGCTGGATctgctg                                              <  1:140902/251‑1 (MQ=255)
                                                                                                                                                                                                                 agagGGGAGGTCATGGATGTGAATGGGTGCAGGATATGGCGTTGGTCAAAGGCACGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTCGCAAGTCACAGAATTATTTGAAAAACCGTTGCTGGATCTGCTGTTTGAAGCGCACCAGGAGCATc                        >  2:272736/1‑251 (MQ=255)
                                                                                                                                                                                                                                       aTGGGTGCAGGATAGGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGGCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTCGCAAGTCACAGAATTATTTGAAAAACCGTTGCTGGATCTGCTGTTTGAAGCGCAGCAGGTGCATCGCCAGCATTTCGATCCTCGTCa  <  1:17771/251‑1 (MQ=255)
                                                                                                                                                                                                                                               ||                                                                                                                                                                                                                                                 
TGCGGCACAGCTCAATGGCTGGCGCATGGGTGATACCGCCAAACATCACATGCGACATGGCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCCAGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTCGCAAGTCACAGAATTATTTGAAAAACCGTTGCTGGATCTGCTGTTTGAAGCGCAGCAGGTGCATCGCCAGCATTTCGATCCTCGTCA  >  NC_000913/808978‑809459

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