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,2170CA100.0% 88.7 / NA 27W14M (TGG→CTG) bioAadenosylmethionine‑8‑amino‑7‑oxononanoate aminotransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (13/14);  total (13/14)
*NC_000913809,2180AG96.3% 91.4 / NA 27W14M (TGG→CTG) bioAadenosylmethionine‑8‑amino‑7‑oxononanoate aminotransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (13/13);  minor base T (0/1);  total (13/14)
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.97e-01

TGCGGCACAGCTCAATGGCTGGCGCATGGGTGATACCGCCAAACATCACATGCGACATGGCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCCAGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTCGCAAGTCACAGAATTATTTGAAAAACCGTTGCTGGATCTGCTGTTTGAAGCGCAGCAGGTGCATCGCCAGCATTTCGAT  >  NC_000913/808978‑809451
                                                                                                                                                                                                                                               ||                                                                                                                                                                                                                                          
tGCGGCACAGCTCAATGGCTGGCGCATGGGTGATACCGCCAAACATCACATGCGACATGGCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCg                                                                                                                                                                                                                                   <  2:15098/250‑1 (MQ=255)
                        cATGGGTGATACCGCCAAACATCACATGCGACATGGCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGtt                                                                                                                                                                                                         >  2:125892/1‑251 (MQ=255)
                         aTGGGTGATACCGCCAAACATCCCATGCGACATGGCATCAATTGGCGACTTCATCGCCGCATTAAGCTGCGGGTGAGTGTAGCCGTGGATCGCCGCCCACCAGGACGATATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCGCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCATGATCTGGCGGTTGGCAACGGCAATATCGTCCGTTg                                                                                                                                                                                                        <  2:114011/251‑1 (MQ=255)
                                        aaaCATCACATGCGACATGGCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCCCCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGt                                                                                                                                                                                         <  1:54833/251‑1 (MQ=255)
                                            aTCACATGCGACATGGCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAAc                                                                                                                                                                                     <  2:88379/251‑1 (MQ=255)
                                              cacaTGCGACATGGCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAAcc                                                                                                                                                                                    <  2:15811/250‑1 (MQ=255)
                                                       gaTGGCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCGGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGaaaa                                                                                                                                                                          <  2:4639/250‑1 (MQ=255)
                                                          ggCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGaaa                                                                                                                                                                           >  1:217539/1‑247 (MQ=255)
                                                                      cGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACaa                                                                                                                                                           <  1:193343/251‑1 (MQ=255)
                                                                        aCTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCCGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCCTAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGt                                                                                                                                                         >  2:145594/1‑251 (MQ=255)
                                                                             aTCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAACTTGAAAAGATTTAGGTTTACAAGTCTac                                                                                                                                                     >  1:117729/1‑250 (MQ=255)
                                                                               cgccgcATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACAcc                                                                                                                                                  >  1:181805/1‑251 (MQ=255)
                                                                                  cgcATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGaa                                                                                                                                               <  1:93749/251‑1 (MQ=255)
                                                                                             gCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAAcaacaa                                                                                                                                     <  1:120625/250‑1 (MQ=255)
                                                                                                   gATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAacac                                                                                                                              <  1:116995/251‑1 (MQ=255)
                                                                                                        tAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGtttt                                                                                                                         >  2:138511/1‑251 (MQ=255)
                                                                                                                                                     cGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTCGCa                                                                              >  2:8637/1‑249 (MQ=255)
                                                                                                                                                           ccGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTCGCAAGTCACAg                                                                      <  1:8637/251‑1 (MQ=255)
                                                                                                                                                               cAGACAAAATCAGCTCGTAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATCTGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAACGATTCAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTCGCAAGCCACAGAatt                                                                  >  2:73260/1‑251 (MQ=255)
                                                                                                                                                                    aaaaTCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTCGCAAGTcaca                                                                       <  2:221635/241‑1 (MQ=255)
                                                                                                                                                                    aaaaTCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGAGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTCGCAAGTcaca                                                                       >  1:221635/1‑241 (MQ=255)
                                                                                                                                                                      aaTCAGCTCGCAACCTTAGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCAACGCTGGACATTGTCGCAAGTCACAGAATTATTTGa                                                            >  1:271339/1‑250 (MQ=255)
                                                                                                                                                                                   ccTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTCGCAAGTCACAGAATTATTTGAAAAACCGTTGCTgg                                              >  1:6051/1‑251 (MQ=255)
                                                                                                                                                                                                 caccGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTCGCAAGTCACAGAATTATTTGAAAAACCGTTGCTGGATCTGCTGTTTGa                                 >  2:228362/1‑250 (MQ=255)
                                                                                                                                                                                                                             aTGGATGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTCGCAAGTCACAGAATTATTTGAAAAACCGTTGCTGGATCTGCTGTTTGAAGCGCAGCAGGTGCATCGCCAGCATTTCg    <  2:14380/251‑1 (MQ=255)
                                                                                                                                                                                                                             aTGGATGAGTATGGGTGCAGGATATGGCGTTGGTAAAAGGCAAGATCGTCCGTTGTCATAATCGAATTGTAAACCAAATTGCAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTCGCAAGTCCCAGAATTATTTGCAAAACCGATGCTGGATCTGCTGTTTGAAGCGAAGCAGGTGACTCGCCAGCATTTCg    >  2:172620/1‑251 (MQ=255)
                                                                                                                                                                                                                               gtaTGTGTATGGGTGCAGGATATGGCGTTGGTCAAAGGCAGTATCTTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTATGTTTACAAGTCTACACCGAATTAACACCAAAAAACTCGTGTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTCGCAAGTCCCAGAATTATTTGAAAAACCGTTGCTGGAGCGGCTGTTTGAAGCGCAGCAGGTGCATCGCCAGCATTTCGAt  <  1:126848/249‑1 (MQ=255)
                                                                                                                                                                                                                                               ||                                                                                                                                                                                                                                          
TGCGGCACAGCTCAATGGCTGGCGCATGGGTGATACCGCCAAACATCACATGCGACATGGCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGATTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCCAGATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTCGCAAGTCACAGAATTATTTGAAAAACCGTTGCTGGATCTGCTGTTTGAAGCGCAGCAGGTGCATCGCCAGCATTTCGAT  >  NC_000913/808978‑809451

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