Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 634,746 2 bp→AT [ybdL][ybdM] [ybdL], [ybdM]

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913634,7460GA100.0% 105.0 / NA 34*387* (TAG→TAAybdLmethionine aminotransferase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (15/19);  total (15/19)
*NC_000913634,7470CT97.1% 104.2 / ‑6.3 34*210* (TAG→TAAybdMSpo0J family protein, ParB‑like nuclease domain
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (15/18);  minor base G (0/1);  total (15/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAGCGGCGATTCCGCTGTCGGTGTTTTGCGCCGATCCCTTCCCACATAAACTGATTCGTCTCTGTTTTGCCAAGAAGGAATCGACGTTGCTGGCAGCAGCTGAACGCCTGCGCCAGCTTTAGCTATTTAACCGTCCAGGCGCGGGAATATTGACGGTCGGCAAACAACTCTTGCAGGCCGTTAATTTGCTTCAGGCGCAGCACTTCGTCGCTGTCCATTCCCAGCTCTTTGC  >  NC_000913/634625‑634856
                                                                                                                         ||                                                                                                             
tAGCGGCGATTCCGCTGTCGGTGTTTTGCGCCGATCCCTTCCCACATAAACTGATTCGTCTCTGTTTTGCCAAGAAGGAATCGACGTTGCTGGCAGCATCTGAACGCCTGCGCCAGCTTTAATTATTTaa                                                                                                        >  2:1641845/1‑130 (MQ=255)
tAGCGGCGATTCCGCTGTCGGTGTTTTGCGCCGATCCCTTCCCACATAAACTGATTCGTCTCTGTTTTGCCAAGAAGGAATCGACGTTGCTGGCAGCATCTGAACGCCTGCGCCAGCTTTAATTATTTaa                                                                                                        <  1:1641845/130‑1 (MQ=255)
tAGCGGCGATTCCGCTGTCGGTGTTTTGCGCCGATCCCTTCCCACATAAACTGATTCGTCTCTGTTTTGCCAAGAAGGAATCGACGTTGCTGGCAGCAGCTGAACGCCTGCGCCAGCTTTAATTATTTAACCGTCCAgg                                                                                               <  1:608248/139‑1 (MQ=255)
tAGCGGCGATTCCGCTGTCGGTGTTTTGCGCCGATCCCTTCCCACATAAACTGATTCGTCTCTGTTTTGCCAAGAAGGAATCGACGTTGCTGGCAGCAGCTGAACGCCTGCGCCAGCTTTAATTATTTAACCGTCCAgg                                                                                               <  2:408734/139‑1 (MQ=255)
 aGCGGCGATTCCGCTGTCGGTGTTTTGCGCCGATCCCTTCCCACATAAACTGATTCGTCTCTGTTTTGCCAAGAAGGAATCGACGTTGCTGGCAGCAGCTGAACGCCTGCGCCAGCTTTAATTATTTAACCGTCCAGgc                                                                                              >  1:1494202/1‑139 (MQ=255)
    ggcgATTCCGCTGTCGGTGTTTTGCGCCGATCCCTTCCCACATAAACTGATTCGTCTCTGTTTTGCCAAGAAGGAATCGACGTTGCTGGCAGCAGCTGAACGCCTGCGCCAGCTTTAATTATTTAACCGTCCAGGcgcg                                                                                           <  1:1775727/139‑1 (MQ=255)
        aTTCCGCTGTCGGTGTTTTGCGCCGATCCCTTCCCACATAAACTGATTCGTCTCTGTTTTGCCAAGAAGGAATCGACGTTGCTGGCAGCAGCTGAACGCCTGCGCCAGCTTTAATTATTTAACCGTCCAGGCGCGGGaa                                                                                       >  1:722046/1‑139 (MQ=255)
                        tttGCGCCGATCCCTTCCCACATAAACTGATTCGTCTCTGTTTTGCCAAGAAGGAATCGACGTTGCTGGCAGCAGCTGAACGCCTGCGCCAGCTTTAATTATTTAACCGTCCAGGCGCGGGAATATTGACGGTCGGCaa                                                                       <  1:504166/139‑1 (MQ=255)
                           gcgcCGATCCCTTCTCACATAAACTGATTCGTCTCTGTTTTGCCAAGAAGGAATCGACGTTGCTGGCAGCAGCTGAACGCCTGCGCCAGCTTTAATTATTTAACCGTCCAGGCGCGGGAATATTGACGGTCg                                                                           <  2:1056063/132‑1 (MQ=255)
                           gcgcCGATCCCTTCCCACATCAACTGATTCGTCTCTGTTTTGCCAAGAAGGAATCGACGTTGCTGGCAGCAGCTGAACGCCTGCGCCAGCTTTAATTATTTAACCGTCCAGGCGCGGGAATATTGACGGTCg                                                                           >  1:1056063/1‑132 (MQ=255)
                                 aTCCCTTCCCACATAAACTGATTCGTCTCTGTTTTGCCAAGAAGGAATCGACGTTGCTGGCAGCAGCTGAACGCCTGCGCCAGCTTTAATTATTTAACCGTCCAGGCGCGGGAATATTGACGGTCGGCAAACAACTCtt                                                              >  2:54532/1‑139 (MQ=255)
                                   cccTTCCCACATAAACTGATTCGTCTCTGTTTTGCCAAGAAGGAATCGACGTTGCTGGCAGCAGCTGAACGCCTGCGCCAGCTTTAATTATTTAACCGTCCAGGCGCGGGAATATTGACGGTCGGCAAACAACTCTTGc                                                            <  1:403614/139‑1 (MQ=255)
                                        cccACATAAACTGATTCGTCTCTGTTTTGCCAAGAAGGAATCGACGTTGCTGGCAGCAGCTGAACGCCTGCGCCAGCTTTAATTATTTAACCGTCCAGGCGCGGGAATATTGACGGTCGGCAAACAACTCTTGCAGGcc                                                       <  2:1023070/139‑1 (MQ=255)
                                             aTAAACTGATTCGTCTCTGTTTTGCCAAGAAGGAATCGACGTTGCTGGCAGCAGCTGAACGCCTGCGCCAGCTTTAATTATTTAACCGTCCAGGcgcg                                                                                           >  1:1003551/1‑98 (MQ=255)
                                             aTAAACTGATTCGTCTCTGTTTTGCCAAGAAGGAATCGACGTTGCTGGCAGCAGCTGAACGCCTGCGCCAGCTTTAATTATTTAACCGTCCAGGcgcg                                                                                           <  2:1003551/98‑1 (MQ=255)
                                              tAAACTGATTCGTCTCTGTTTTGCCAAGAAGGAATCGACGTTGCTGGCAGCAGCTGAACGCCTGCGCCAGCTTTAATTATTTAACCGTCCAGGCGCGGGAATATTGACGGTCGGCAAACAACTCTTGCAGGCCGTTAAt                                                 >  2:337670/1‑139 (MQ=255)
                                                     aTTCGTCTCTGTTTTGCCAAGAAGGAATCGACGTTGCTGGCAGCAGCTGAACGCCTGCGCCAGCTTTAATTATTTAACCGTCCAGGCGCGGGAATATTGACGGTCGGCAAACAACTCTTGCAGGCCGTTAATTTGCTTc                                          <  1:337670/139‑1 (MQ=255)
                                                       tCGTCTCTGTTTTGCCAAGAAGGAATCGACGTTGCTGGCAGCAGCTGAACGCCTGCGCCAGCTTTAATTATTTAACCGTCCAGGCGCGGGAATATTGACGGTCGGCAAACAACTCTTGCAGGCCGTTAATTTGCTTCAg                                        >  1:444140/1‑139 (MQ=255)
                                                          tctcTGTTTTGCCAAGAAGGAATCGACGTTGCTGGCAGCAGCTGAACGCCTGCGCCAGCTTTAATTATTTAACCGTCCAGGCGCGGGAATATTGACGGt                                                                             >  1:1353901/1‑99 (MQ=255)
                                                          tctcTGTTTTGCCAAGAAGGAATCGACGTTGCTGGCAGCAGCTGAACGCCTGCGCCAGCTTTAATTATTTAACCGTCCAGGCGCGGGAATATTGACGGt                                                                             <  2:1353901/99‑1 (MQ=255)
                                                                     ccAAGAAGGAATCGACGTTGCTGGCAGCAGCTGAACGCCTGCGCCAGCTTTAATTATTTAACCGTCCAGGCGCGGGAATATTGACGGTCGGCAAACAACTCTTGCAGGCCGTTAATTTGCTTCAGGCGCAGCACTtcgt                          >  1:1694654/1‑139 (MQ=255)
                                                                     ccAAGAAGGAATCGACGTTGCTGGCAGCAGCTGAACGCCTGCGCCAGCTTTAATTATTTAACCGTCCAGGCGCGGGAATATTGACGGTCGGCAAACAACTCTTGCAGGCCGTTAATTTGCTTCAGGCGCAGCACTtcgt                          >  1:1392550/1‑139 (MQ=255)
                                                                       aagaagGAATCGACGTTGCTGGCAGCAGCTGAAAGCCTGCGCCAGCTTTAATTATTTAACCGTCCAGGCGCGGGAATATTGACGGTCGGCAAACAACTCTTGCAGGCCGTTAATTTGCTTCAGGCGCAGCACTtcgtcg                        <  2:229780/139‑1 (MQ=255)
                                                                             gAATCGACGTGGCTGGCAGCAGCTGAACGCCTGCGCCAGCTTTAAGTATTTAACCGTCCAGGCGCGGGAATATTGACGGTCGGCAAACAACTCTTGCAGGCCGTTAATTTGCTTCAGGCGCAGCACTTCGTCGCTGTcc                  <  2:1694654/139‑1 (MQ=255)
                                                                                    cGTTGCTGGCAGCAGCTGAACGCCTGCGCCAGCTTTAATTATTTAACCGTCCAGGCGCGGGAATATTGACGGTCGGCAAACAGCTCTTGCAGGCCGTTAATTTGCTTCAGGCGCAGCACTTCGTCGCTGTCCATTccc            <  1:1134786/138‑1 (MQ=255)
                                                                                    cGTTGCTGGCAGCAGCTGAACGCCTGCGCCAGCTTTAATTATTTAACCGTCCAGGCGCGGGAATATTGACGGTCGGCAAACAGCTCTTGCAGGCCGTTAATTTGCTTCAGGCGCAGCACTTCGTCGCTGTCCATTccc            >  2:1134786/1‑138 (MQ=255)
                                                                                     gTTGCTGGCAGCAGCTGAACGCCTGCGCCAGCTTTAATTATTTAACCGTCCAGGCGCGGGAATATTGACGGTCGGCAAACAACTCTTGCAGGCCGTTAATTTGCTTCAGGCGCAGCACTTCGTCGCTGTCCATTCCCAg          <  2:722046/139‑1 (MQ=255)
                                                                                        gCTGGCAGCAGCTGAACGCCTGCGCCAGCTTTAATTATTTAACCGTCCAGGCGCGGGAATATTGACGGTCGGCAAACAACTCTTGCAGGCCGTTAATTTGCTTCAGGCGCAGCACTTCGTCGCTGTCCATTCCCAGctc       >  2:644548/1‑139 (MQ=255)
                                                                                             cagcagCTGAACGCCTGCGCCAGCTTTAATTATTTAACCGTCCAGGCGCGGGAATATTGACGGTCGGCAAACAACTCTTGCAGGCCGTTAATTTGCTTCAGGCGCAGCACTTCGTCGCTGTCCATTCCCAGCTCTTTGc  <  2:1392550/139‑1 (MQ=255)
                                                                                             cagcagCTGAACGCCTGCGCCAGCTTTAATTATTTAACCGTCCAGGCGCGGGAATATTGACGGTCGGCAAACAACTCTTGCAGGCCGTTAATTTGCTTCAGGCGCAGCACTTCGTCGCTGTCCATTCCCAGCTCTTTGc  <  1:1556989/139‑1 (MQ=255)
                                                                                                 agcTGAACGCCTGCGCCAGCTTTAATTATTTAACCGTCCa                                                                                                 <  2:968323/40‑1 (MQ=37)
                                                                                                 agcTGAACGCCTGCGCCAGCTTTAATTATTTAACCGTCCa                                                                                                 >  1:968323/1‑40 (MQ=255)
                                                                                                   cTGAACGCCTGCGCCAGCTTTAATTATTTAACCGTCCAGGCGCGGGAATATTGACGGTCGGCa                                                                        >  2:1023665/1‑63 (MQ=255)
                                                                                                   cTGAACGCCTGCGCCAGCTTTAATTATTTAACCGTCCAGGCGCGGGAATATTGACGGTCGGCa                                                                        <  1:1023665/63‑1 (MQ=255)
                                                                                                                         ||                                                                                                             
TAGCGGCGATTCCGCTGTCGGTGTTTTGCGCCGATCCCTTCCCACATAAACTGATTCGTCTCTGTTTTGCCAAGAAGGAATCGACGTTGCTGGCAGCAGCTGAACGCCTGCGCCAGCTTTAGCTATTTAACCGTCCAGGCGCGGGAATATTGACGGTCGGCAAACAACTCTTGCAGGCCGTTAATTTGCTTCAGGCGCAGCACTTCGTCGCTGTCCATTCCCAGCTCTTTGC  >  NC_000913/634625‑634856

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: