Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 642,943 A→G S286G (AGC→GGC)  ybdR → putative Zn‑dependent NAD(P)‑binding oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913642,9430AG97.4% 133.2 / ‑6.6 39S286G (AGC→GGC) ybdRputative Zn‑dependent NAD(P)‑binding oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (17/21);  minor base C (1/0);  total (18/21)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGATGAAGACAGCGATCCGGCACAGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCAGCAGCGGTAAAGCGTTGCGTCAGTGTATTGCGGCGGTCAGGCGTGGCGGCATTGTTAGCGTACCGGGCGTCTACGCTGGATTTATTCACGGTTTCCTGTTTGGCGACGCCTTTGATAAAGGGTTGTCGTTTAAAATGGGAC  >  NC_000913/642801‑643084
                                                                                                                                              |                                                                                                                                             
tgatgaAGACAGAGATCCGGCACAGTCAATTATTGAACAAACGGCAGGTCAGCGGGGGGTGGGTGCAGTAATAGACGCCGTGGGTTTTGAAGCGAAAGCGAGCACCACGGAAACGGTGCTTACGAACCTGAAACTGGATGGCCGCAGTGGt                                                                                                                                       >  1:55305/1‑151 (MQ=255)
   tgaAGACAGCGATCCGGCACAGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCGGCAGCGGTaa                                                                                                                                     >  1:132937/1‑150 (MQ=255)
      aGACAGCGATCCGGCACAGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCGGCAGCGGTAAAGCg                                                                                                                                 >  1:733721/1‑151 (MQ=255)
         cAGCGATCCGGCACAGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCGGCAGCGGTAAAGCGTTg                                                                                                                              <  1:3417/151‑1 (MQ=255)
          aGCGATCCGGCACAGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCGGCAGCGGTAAAGCGTTGc                                                                                                                             <  2:6420/151‑1 (MQ=255)
            cGATCCGGCACAGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCGGCAGCGGTAAAGCGTTGCGt                                                                                                                           <  1:720812/151‑1 (MQ=255)
             gATCCGGCACAGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCGGCAGCGGTAAAGCGTTGCGTc                                                                                                                          >  2:516549/1‑151 (MQ=255)
                ccGGCACAGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCGGCAGCGGTAAAGCGTTGCGTCAgt                                                                                                                       <  2:286191/151‑1 (MQ=255)
                  ggCACAGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCGGCAGCGGTAAAGCGTTGCGTCAgtg                                                                                                                      >  2:354071/1‑150 (MQ=255)
                      caGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCGGCAGCGGTAAAGCGTTGCGTCAGTGTATTg                                                                                                                 <  1:657107/151‑1 (MQ=255)
                        gTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCGGCAGCGGTAAAGCGTTGCGTCAGTGTATTgcg                                                                                                               <  2:81765/151‑1 (MQ=255)
                        gTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCGGCAGCGGTAAAGCGTTGCGTCAGTGTATTgc                                                                                                                <  1:692608/150‑1 (MQ=255)
                                   aaCAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCGGCAGCGGTAAAGCGTTGCGTCAGTGTATTGCGGCGGTCAGGCg                                                                                                    <  1:626594/151‑1 (MQ=255)
                                    aCAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCGGCAGCGGTAAAGCGTTGCGTCAGTGTATTGCGGCGGTCAGGCg                                                                                                    >  2:235428/1‑150 (MQ=255)
                                        aCGGCAGGACACTGGGGCGTGGGTGAAGTAATAGACGCCGTTGGGTTTGAAGCGAAAGGCAGCACCAGGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCGGCAGCGGTAAAGCGGTGCGTCAGTGGATTGcggcggctacgggtgggg                                                                                               >  1:57656/1‑139 (MQ=255)
                                          ggCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCGGCAGCGGTAAAGCGTTGCGTCAGTGTATTGCGGCGGTCAGGCGTggcgg                                                                                              >  2:745666/1‑150 (MQ=255)
                                                 cACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCGGCAGCGGTAAAGCGTTGCGTCAGTGTATTGCGGCGGTCAGGCGTGGCGGCATTGTTa                                                                                      <  2:258515/151‑1 (MQ=255)
                                                   ccGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCGGCAGCGGTAAAGCGTTGCGTCAGTGTATTGCGGCGGTCAGGCGTGGCGGGATTGtt                                                                                       >  1:14751/1‑148 (MQ=255)
                                                   ccGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCGGCAGCGGTAAAGCGTTGCGTCAGTGTATTGCGGCGGTCAGGCGTGGCGGCATTGTTa                                                                                      <  1:528874/149‑1 (MQ=255)
                                                            ggATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCGGCAGCGGTAAAGCGTTGCGTCAGTGTATTGCGGCGGTCAGGCGTGGCGGCATTGTTAGCGTACCggg                                                                            >  2:574359/1‑150 (MQ=255)
                                                                gCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCGGCAGCGGTAAAGCGTTGCGTCAGTGTATTGCGGCGGTCAGGCGTGGCGGCATTGTTAGCGTACCGGGCGt                                                                         >  1:798922/1‑149 (MQ=255)
                                                                gCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCGGCAGCGGTAAAGCGTTGCGTCAGTGTATTGCGGCGGTCAGGCGTGGCGGCATTGTTAGCGTACCGGGCGTCt                                                                       <  1:607879/151‑1 (MQ=255)
                                                                  aGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCGGCAGCGGTAAAGCGTTGCGTCAGTGTATTGCGGCGGTCAGGCGTGGCGGCATTGTTAGCGTACCGGGCGTCTa                                                                      >  1:760679/1‑150 (MQ=255)
                                                                  aGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCGGCAGCGGTAAAGCGTTGCGTCAGTGTATTGCGGCGGTCAGGCGTGGCGGCATTGTTAGCGTACCGGGCGTCTAc                                                                     <  2:846930/151‑1 (MQ=255)
                                                                     aaTAGACGCCGTAGGTTTTGAAGCGAAAGGCAGCACCACCGAAACGGTGGTGACTAACCTGAAACTGGCGGGCGGCAGCGGTAAAGCGATGCGTCAGTGTATTGCGGCGGGCCGGCGCGGCGGGATTGTTAGGGTAGCGGGCggtttagc                                                                   >  1:125879/1‑143 (MQ=255)
                                                                        aGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCGGCAGCGGTAAAGCGTTGCGTCAGTGTATTGCGGCGGTCAGGCGGGGCGGCATTGTTAGCGTACCGGGCGTCTACGCTgg                                                                >  1:519315/1‑150 (MQ=255)
                                                                         gACTCCGTCGGTTTTGAAGCGGAAGGCGGCGCCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCGGCAGCGGTAAAGCGTTGCGTCAGTGTATTGCGGCGGTCAGGCGTGGCGGCATTGTTAGCGTACCGGGCGTCTACGCTGGAt                                                              <  1:827815/151‑1 (MQ=255)
                                                                         gACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCGGCAGCGGTAAAGCGTTGCGTCAGTGTATTGCGGCGGTCAGGCGTGGCGGCATTGTTAGCGTACCGGGCGTCTACGCTGGAt                                                              <  2:747318/151‑1 (MQ=255)
                                                                           cGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCGGCAGCGGTAAAGCGTTGCGTCAGTGTATTGCGGCGGTCAGGCGTGGCGGCATTGTTAGCGTACCGGGCGTCTACGCTGGAttt                                                            <  2:14772/151‑1 (MQ=255)
                                                                                             gAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCGGCAGCGGTAAAGCGTTGCGTCAGTGTATTGCGGCGGTCAGGCGTGGCGGCATTGTTAGCGTACCGGGCGTCTACGCTGGATTTATTCACGGTTTCCTGttt                                          >  1:162504/1‑151 (MQ=255)
                                                                                               aaGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCGGCAGCGGTAAAGCGTTGCGTCAGTGTATTGCGGCGGTCAGGCGTGGCGGCATTGTTAGCGTACCGGGCGTCTACGCTGGATTTATTCACGGTTTCCTGttttg                                        >  2:787723/1‑149 (MQ=255)
                                                                                                 ggCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCGGCAGCGGTAAAGCGTTGCGTCAGTGTATTGCGGCGGTCAGGCGTGGCGGCATTGTTAGCGTACCGGGCGTCTACGCTGGATTTATTCACGGTTTCCTGTTTGGCg                                      <  1:751910/151‑1 (MQ=255)
                                                                                                   cagcaCCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCGGCAGCGGTAAAGCGTTGCGTCAGTGTATTGCGGCGGTCAGGCGTGGCGGCATTGTTAGCGTACCGGGCGTCTACGCTGGATTTATTCACGGTTTCCTGTTTGGCGAc                                    <  1:35301/151‑1 (MQ=255)
                                                                                                       accacGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCGGCAGCGGTAAAGCGTTGCGTCAGTGTATTGCGGCGGTCAGGCGTGGCGGCATTGTTAGCGTACCGGGCGTCTACGCTGGATTTATTCACGGTTTCCTGTTTGGCGACGCCt                                <  1:713384/151‑1 (MQ=255)
                                                                                                                 cGGTGCTGACTAACCTGAAACTGGAGGGCGGCAGCGGTAAAGCGTTGCGTCAGTGTATTGCGGCGGTCAGGCGTGGCGGCATTGTTAGCGTACCGGGCGTCTACGCTGGATTTATTCACGGTTTCCTGTTTGGCGACGCCTTTGATAAAgg                      <  1:497674/151‑1 (MQ=255)
                                                                                                                   gTGCTGACTAACCTGAAACTGGAGGGCGGCAGCGGTAAAGCGTTGCGTCAGTGTATTGCGGCGGTCAGGCGTGGCGGCATTGTTAGCGTACCGGGCGTCTACGCTGGATTTATTCACGGTTTCCTGTTTGGCGACGCCTTTGATAAAggg                     >  2:610123/1‑150 (MQ=255)
                                                                                                                            aaCCTGAAACTGGAGGGCGGCAGCGGTAAAGCGTTGCGTCAGTGTATTGCGGCGGTCAGGCGTGGCGGCATTGTTAGCGTACCGGGCGTCTACGCTGGATTTATTCACGGTTTCCTGTTTGGCGACGCCTTTGATAAAGGGTTGTCGttt            <  1:708673/150‑1 (MQ=255)
                                                                                                                                 gAAACTGGAGGGCGGCAGCGGTAAAGCGTTGCGTCAGTGTATTGCGGCGGTCAGGCGTGGCGGCATTGTTAGCGTACCGGGCGTCTACGCTGGATTTATTCACGGTTTCCTGTTTGGCGACGCCTTTGATAAAGGGTTGTCGTTTAAAAt       <  2:297453/150‑1 (MQ=255)
                                                                                                                                      tGGAGGGCGGCAGCGGTAAAGCGTTGCGTCAGTGTATTGCGGCGGTCAGGCGTGGCGGCATTGTTAGCGTACCGGGCGTCTACGCCGGGTTTATTCACGGTTTCCTGTTTGGCGACGCCTTTGATAAAGGGGTGTCGTTTAAAATGGGAc  >  2:829985/1‑150 (MQ=255)
                                                                                                                                              |                                                                                                                                             
TGATGAAGACAGCGATCCGGCACAGTCAATTATTGAACAAACGGCAGGTCACCGGGGCGTGGATGCAGTAATAGACGCCGTCGGTTTTGAAGCGAAAGGCAGCACCACGGAAACGGTGCTGACTAACCTGAAACTGGAGGGCAGCAGCGGTAAAGCGTTGCGTCAGTGTATTGCGGCGGTCAGGCGTGGCGGCATTGTTAGCGTACCGGGCGTCTACGCTGGATTTATTCACGGTTTCCTGTTTGGCGACGCCTTTGATAAAGGGTTGTCGTTTAAAATGGGAC  >  NC_000913/642801‑643084

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